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.     

Does Joint Physical Custody “Cause” Children’s Better Outcomes?

ABSTRACT

Policymakers and researchers are concerned with whether joint physical custody (JPC) produces better outcomes for children than sole custody. Although several review articles summarizing up to 61 empirical articles demonstrate very positive answers, many of the research designs used compromise the ability to claim that it is JPC per se—and not selection effects—that causes the effect. We discuss several research design issues, such as propensity score analysis, that can more powerfully probe the question of causality. Some studies have already been conducted employing these strategies and more are recommended and likely to soon be forthcoming. On the basis of this comprehensive review we conclude that JPC probably does cause benefits to children on average, and that social scientists can now provisionally recommend rebuttably presumptive JPC to policymakers.     

Varieties of political support in emerging democracies: A cross-national analysis

Political support is a multidimensional construct encompassing evaluations of political leaders and institutions (specific support) as well as adherence to basic regime principles (diffuse support). Scholars have traditionally assumed that diffuse and specific support are driven by different forces and evolve largely independently. Prior empirical work, however, has struggled to untangle the two support dimensions and focused predominantly on cross-national differences, ignoring their dynamics. This paper develops an analytical and empirical approach to examine the levels and dynamic interplay of both support dimensions and estimate their determinants, applying it to South American democracies between 1996 and 2015. Contrary to received wisdom, we show that both dimensions are quite volatile and closely linked in this region. In particular, negative economic shocks not only undermine support for government actors, but also fuel democratic disenchantment. Nonetheless, while regime support is rather fickle in South America, it can be ultimately resistant to performance fluctuations.     

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.     

The t-information and its use in multivariate problems and time-series analysis

The second moment of recently introduced scalar inference function can be viewed as generalized Fisher information of the continuous probability distributions. In this paper we call it the t-information and give some possible applications of the new concept.     

R×C ecological inference: bounds, correlations, flexibility and transparency of assumptions

Summary.  Despite its potential pitfalls, ecological inference is an unavoidable part of some quantitative settings, including US voting rights litigation. In such applications, the analyst will typically encounter two-way tables with more than two rows and columns. Although several ecological inference methods are currently available for 2×2 tables, there are fewer options for analysing general R × C tables, and virtually none that model counts as opposed to fractions. We propose a count R × C method that respects the bounds deterministically, that allows for complex relationships between internal cell quantities, that is easily extensible and that results from transparent assumptions. We study the method via simulation, and then apply it to an example that is drawn from the state of Texas relevant to recent redistricting litigation there.     

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.     

从释义到推论——余英时《朱熹的历史世界》史料解读献疑

余英时先生《朱熹的历史世界》一书重构两宋（尤其是南宋）政坛面貌,旁征博引,新见迭出,深得学界好评。然而通读全书,作者对若干史料的运用,似不无令人生疑之处。书中史料处理的疑点。主要可分为：字词含义、语段脉络、史料性质、推论程序四类。     

Adaptive multiple importance sampling for Gaussian processes

In applications of Gaussian processes (GPs) where quantification of uncertainty is a strict requirement, it is necessary to accurately characterize the posterior distribution over Gaussian process covariance parameters. This is normally done by means of standard Markov chain Monte Carlo (MCMC) algorithms, which require repeated expensive calculations involving the marginal likelihood. Motivated by the desire to avoid the inefficiencies of MCMC algorithms rejecting a considerable amount of expensive proposals, this paper develops an alternative inference framework based on adaptive multiple importance sampling (AMIS). In particular, this paper studies the application of AMIS for GPs in the case of a Gaussian likelihood, and proposes a novel pseudo-marginal-based AMIS algorithm for non-Gaussian likelihoods, where the marginal likelihood is unbiasedly estimated. The results suggest that the proposed framework outperforms MCMC-based inference of covariance parameters in a wide range of scenarios.