Higher-order Bayesian Approximations for Pseudo-posterior Distributions

The theory of higher-order asymptotics provides accurate approximations to posterior distributions for a scalar parameter of interest, and to the corresponding tail area, for practical use in Bayesian analysis. The aim of this article is to extend these approximations to pseudo-posterior distributions, e.g., posterior distributions based on a pseudo-likelihood function and a suitable prior, which are proved to be particularly useful when the full likelihood is analytically or computationally infeasible. In particular, from a theoretical point of view, we derive the Laplace approximation for a pseudo-posterior distribution, and for the corresponding tail area, for a scalar parameter of interest, also in the presence of nuisance parameters. From a computational point of view, starting from these higher-order approximations, we discuss the higher-order tail area (HOTA) algorithm useful to approximate marginal posterior distributions, and related quantities. Compared to standard Markov chain Monte Carlo methods, the main advantage of the HOTA algorithm is that it gives independent samples at a negligible computational cost. The relevant computations are illustrated by two examples.     

Approximate Bayesian computation with composite score functions

Both approximate Bayesian computation (ABC) and composite likelihood methods are useful for Bayesian and frequentist inference, respectively, when the likelihood function is intractable. We propose to use composite likelihood score functions as summary statistics in ABC in order to obtain accurate approximations to the posterior distribution. This is motivated by the use of the score function of the full likelihood, and extended to general unbiased estimating functions in complex models. Moreover, we show that if the composite score is suitably standardised, the resulting ABC procedure is invariant to reparameterisations and automatically adjusts the curvature of the composite likelihood, and of the corresponding posterior distribution. The method is illustrated through examples with simulated data, and an application to modelling of spatial extreme rainfall data is discussed.     

持续改善你的服务

服务不能建立库存,以后再用;如果需求不稳定, 也没办法进行完美预测服务行业在运营方面有一些其他行业所没有的特征,正是这些特征,使管理服务比管理其他行业更为困难。那么,怎样才能提高服务的质量和价值?这需要提供服务的组织从战略的角度去看问题。     

Identity Agents: Parents as Active and Reflective Participants in Their Children's Identity Formation

The paper introduces the concept of identity agents. This concept refers to those individuals who actively interact with children and youth with the intention of participating in their identity formation, and who reflectively mediate larger social influences on identity formation. This contrasts with the focus of mainstream research in the identity field that tends to portray adolescents as the sole reflective agents involved in mature identity development. The paper presents a theoretical analysis presenting the importance of the concept for the formulation of a comprehensive contextual theory of identity formation. The particulars of this concept are illustrated through the presentation of a qualitative report of religious parents actively encouraging their children's processes of identification, co‐participating in their children's identity's formation, and reflectively deliberating their parental roles and goals in regards to this process.     

Quasi-likelihood from<Emphasis Type="Italic">M</Emphasis>-estimators: A numerical comparison with empirical likelihood

In this paper we compare two robust pseudo-likelihoods for a parameter of interest, also in the presence of nuisance parameters. These functions are obtained by computing quasi-likelihood and empirical likelihood from the estimating equations which define robustM-estimators. Application examples in the context of linear transformation models are considered. Monte Carlo studies are performed in order to assess the finite-sample performance of the inferential procedures based on quasi-and empirical likelihood, when the objective is the construction of robust confidence regions.     

Bootstrap adjustments of signed scoring rule root statistics

Scoring rules give rise to methods for statistical inference and are useful tools to achieve robustness or reduce computations. Scoring rule inference is generally performed through first-order approximations to the distribution of the scoring rule estimator or of the ratio-type statistic. In order to improve the accuracy of first-order methods even in simple models, we propose bootstrap adjustments of signed scoring rule root statistics for a scalar parameter of interest in presence of nuisance parameters. The method relies on the parametric bootstrap approach that avoids onerous calculations specific of analytical adjustments. Numerical examples illustrate the accuracy of the proposed method.     

Between stability and higher-order asymptotics

We discuss the effects of model misspecifications on higher-order asymptotic approximations of the distribution of estimators and test statistics. In particular we show that small deviations from the model can wipe out the nominal improvements of the accuracy obtained at the model by second-order approximations of the distribution of classical statistics. Although there is no guarantee that the first-order robustness properties of robust estimators and tests will carry over to second-order in a neighbourhood of the model, the behaviour of robust procedures in terms of second-order accuracy is generally more stable and reliable than that of their classical counterparts. Finally, we discuss some related work on robust adjustments of the profile likelihood and outline the role of computer algebra in this type of research.     

Modified quasi-profile likelihoods from estimating functions

We discuss higher-order adjustments for a quasi-profile likelihood for a scalar parameter of interest, in order to alleviate some of the problems inherent to the presence of nuisance parameters, such as bias and inconsistency. Indeed, quasi-profile score functions for the parameter of interest have bias of order O(1)$\mathrm{O}(1)$, and such bias can lead to poor inference on the parameter of interest. The higher-order adjustments are obtained so that the adjusted quasi-profile score estimating function is unbiased and its variance is the negative expected derivative matrix of the adjusted profile estimating equation. The modified quasi-profile likelihood is then obtained as the integral of the adjusted profile estimating function. We discuss two methods for the computation of the modified quasi-profile likelihoods: a bootstrap simulation method and a first-order asymptotic expression, which can be simplified under an orthogonality assumption. Examples in the context of generalized linear models and of robust inference are provided, showing that the use of a modified quasi-profile likelihood ratio statistic may lead to coverage probabilities more accurate than those pertaining to first-order Wald-type confidence intervals.