Approximate composite marginal likelihood inference in spatial generalized linear mixed models

Non-Gaussian spatial responses are usually modeled using spatial generalized linear mixed model with spatial random effects. The likelihood function of this model cannot usually be given in a closed form, thus the maximum likelihood approach is very challenging. There are numerical ways to maximize the likelihood function, such as Monte Carlo Expectation Maximization and Quadrature Pairwise Expectation Maximization algorithms. They can be applied but may in such cases be computationally very slow or even prohibitive. Gauss–Hermite quadrature approximation only suitable for low-dimensional latent variables and its accuracy depends on the number of quadrature points. Here, we propose a new approximate pairwise maximum likelihood method to the inference of the spatial generalized linear mixed model. This approximate method is fast and deterministic, using no sampling-based strategies. The performance of the proposed method is illustrated through two simulation examples and practical aspects are investigated through a case study on a rainfall data set.     

Bayesian inference for stable Lévy–driven stochastic differential equations with high‐frequency data

In this paper, we consider parametric Bayesian inference for stochastic differential equations driven by a pure‐jump stable Lévy process, which is observed at high frequency. In most cases of practical interest, the likelihood function is not available; hence, we use a quasi‐likelihood and place an associated prior on the unknown parameters. It is shown under regularity conditions that there is a Bernstein–von Mises theorem associated to the posterior. We then develop a Markov chain Monte Carlo algorithm for Bayesian inference, and assisted with theoretical results, we show how to scale Metropolis–Hastings proposals when the frequency of the data grows, in order to prevent the acceptance ratio from going to zero in the large data limit. Our algorithm is presented on numerical examples that help verify our theoretical findings.     

Integrated nested Laplace approximation for the analysis of count data via the combined model: A simulation study

The combined model accounts for different forms of extra-variability and has traditionally been applied in the likelihood framework, or in the Bayesian setting via Markov chain Monte Carlo. In this article, integrated nested Laplace approximation is investigated as an alternative estimation method for the combined model for count data, and compared with the former estimation techniques. Longitudinal, spatial, and multi-hierarchical data scenarios are investigated in three case studies as well as a simulation study. As a conclusion, integrated nested Laplace approximation provides fast and precise estimation, while avoiding convergence problems often seen when using Markov chain Monte Carlo.     

话语标记“要不”简

“要不”有多个义项,到目前为止前人对“要不”表示“建议义”的用法还很少涉及。表示“建议义”的“要不”已经虚化为一个话语标记,对话语的生成和理解起到组织和调控功能;在语句中主要起到语篇衔接、间接否定、开启话题和人际互动功能：其生成动因是语言表达的主观性和交互主观性;生成机制是语境吸收和语用推理。     

试从关联理论角度看语用翻译

文章根据关联理论的明示-推理交际模式,提出了语用翻译实际上是双明示-推理的交际过程,是涉及到原文作者、译者、译文读者三个交际主体的动态的三元翻译,并探讨了认知语用学的关联理论对语用翻译的启发和指导意义,进而尝试将关联理论的研究成果应用于解决翻译实践中的文化、语用类问题。     

多焊接电源控制技术

现代的膜式水冷壁制造工艺要求多达２０台焊接电源同时工作,当外网原因使车间供电电压波动的时候,焊接电源群同时进行调整。群焊接电源的冲击性载荷对电网进行反激,导致电网的更大波动。文章引入优先级的概念,将自适应神经模糊推理系统(ANFIS)应用到焊接电源群的优先级控制上,设计了自适应神经模糊推理网络,并利用MATLAB平台进行整个焊接系统的建模和仿真。结果表明应用该优先级群控策略,各焊机调整能力加强,调整时间缩短,动态性能得到改善。     

Averaging Lorenz curves

A large number of functional forms has been suggested in the literature for estimating Lorenz curves that describe the relationship between income and population shares. The traditional way of overcoming functional-form uncertainty when estimating a Lorenz curve is to choose the function that best fits the data in some sense. In this paper we describe an alternative approach for accommodating functional-form uncertainty, namely, how to use Bayesian model averaging to average the alternative functional forms. In this averaging process, the different Lorenz curves are weighted by their posterior probabilities of being correct. Unlike a strategy of picking the best-fitting function, Bayesian model averaging gives posterior standard deviations that reflect the functional-form uncertainty. Building on our earlier work (Chotikapanich and Griffiths, 2002), we construct likelihood functions using the Dirichlet distribution and estimate a number of Lorenz functions for Australian income units. Prior information is formulated in terms of the Gini coefficient and the income shares of the poorest 10% and poorest 90% of the population. Posterior density functions for these quantities are derived for each Lorenz function and are averaged over all the Lorenz functions.     

Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations

Summary.  Structured additive regression models are perhaps the most commonly used class of models in statistical applications. It includes, among others, (generalized) linear models, (generalized) additive models, smoothing spline models, state space models, semiparametric regression, spatial and spatiotemporal models, log-Gaussian Cox processes and geostatistical and geoadditive models. We consider approximate Bayesian inference in a popular subset of structured additive regression models, latent Gaussian models , where the latent field is Gaussian, controlled by a few hyperparameters and with non-Gaussian response variables. The posterior marginals are not available in closed form owing to the non-Gaussian response variables. For such models, Markov chain Monte Carlo methods can be implemented, but they are not without problems, in terms of both convergence and computational time. In some practical applications, the extent of these problems is such that Markov chain Monte Carlo sampling is simply not an appropriate tool for routine analysis. We show that, by using an integrated nested Laplace approximation and its simplified version, we can directly compute very accurate approximations to the posterior marginals. The main benefit of these approximations is computational: where Markov chain Monte Carlo algorithms need hours or days to run, our approximations provide more precise estimates in seconds or minutes. Another advantage with our approach is its generality, which makes it possible to perform Bayesian analysis in an automatic, streamlined way, and to compute model comparison criteria and various predictive measures so that models can be compared and the model under study can be challenged.     

SEMIPARAMETRIC REGRESSION AND GRAPHICAL MODELS

Semiparametric regression models that use spline basis functions with penalization have graphical model representations. This link is more powerful than previously established mixed model representations of semiparametric regression, as a larger class of models can be accommodated. Complications such as missingness and measurement error are more naturally handled within the graphical model architecture. Directed acyclic graphs, also known as Bayesian networks, play a prominent role. Graphical model-based Bayesian 'inference engines', such as bugs and vibes , facilitate fitting and inference. Underlying these are Markov chain Monte Carlo schemes and recent developments in variational approximation theory and methodology.     

On reliability estimation based on fuzzy lifetime data

Real lifetime data are never precise numbers but more or less non-precise, also called fuzzy. This kind of imprecision is connected with all measurement results of continuous variables, therefore also with time observations. Imprecision is different from errors and variability. Therefore estimation methods for reliability characteristics have to be adapted to the situation of fuzzy lifetimes in order to obtain realistic results.