Multilevel maximum likelihood estimation with application to covariance matrices

The asymptotic variance of the maximum likelihood estimate is proved to decrease when the maximization is restricted to a subspace that contains the true parameter value. Maximum likelihood estimation allows a systematic fitting of covariance models to the sample, which is important in data assimilation. The hierarchical maximum likelihood approach is applied to the spectral diagonal covariance model with different parameterizations of eigenvalue decay, and to the sparse inverse covariance model with specified parameter values on different sets of nonzero entries. It is shown computationally that using smaller sets of parameters can decrease the sampling noise in high dimension substantially.     

Providing a secure base for LGBTQ young people in foster care: The role of foster carers

The experiences and needs of lesbian, gay, bisexual, trans and queer/questioning (LGBTQ) young people in care have been overlooked in England, in both policy and research. This paper reports on findings from the first study of LGBTQ young people in care in England and focuses on the nature of foster carers' experiences and perspectives on caring for LGBTQ young people. Qualitative interviews regarding the fostering role in caring for LGBTQ young people were conducted with a sample of foster carers (n = 26) and analysed thematically. Foster carers described the importance of offering LGBTQ young people not only the nurturing relationships that all children in care need but also availability, sensitivity and acceptance to help young people manage stigma and other challenges associated with minority sexual orientation and gender identity. The Secure Base caregiving model provided a framework for analysing the different dimensions of these relationships. Understanding caregiving roles and relationships for LGBTQ young people in care has important implications for recruiting, training, matching and supporting foster carers to care for LGBTQ young people effectively.     

Conditional and marginal estimates in case-control family data – extensions and sensitivity analyses

This work considers two specific estimation techniques for the family-specific proportional hazards model and for the population-averaged proportional hazards model. So far, these two estimation procedures were presented and studied under the gamma frailty distribution mainly because of its simple interpretation and mathematical tractability. Modifications of both procedures for other frailty distributions, such as the inverse Gaussian, positive stable and a specific case of discrete distribution, are presented. By extensive simulations, it is shown that under the family-specific proportional hazards model, the gamma frailty model appears to be robust to frailty distribution mis-specification in both bias and efficiency loss in the marginal parameters. The population-averaged proportional hazards model, is found to be robust under the gamma frailty model mis-specification only under moderate or weak dependency within cluster members.     

Bootstrapping estimators for the seemingly unrelated regressions model

In this paper we consider the possibility of using the bootstrap to estimate the finite sample variability of feasible generalized least squares and improved estimators applied to the seemingly unrelated regressions model. The improved estimators we employ include members of the Stein-rule family and a hierarchical Bayes estimator proposed by Blattberg and George (1991). Simulation experiments are carried out using several SUR examples as well as a very large example based on the price-promotion model, and data, from marketing research.     

Analytic strategies for longitudinal networks with missing data

Missing data are often problematic when analyzing complete longitudinal social network data. We review approaches for accommodating missing data when analyzing longitudinal network data with stochastic actor-based models. One common practice is to restrict analyses to participants observed at most or all time points, to achieve model convergence. We propose and evaluate an alternative, more inclusive approach to sub-setting and analyzing longitudinal network data, using data from a school friendship network observed at four waves (N = 694). Compared to standard practices, our approach retained more information from partially observed participants, generated a more representative analytic sample, and led to less biased model estimates for this case study. The implications and potential applications for longitudinal network analysis are discussed.     

Marginal likelihood estimation from the Metropolis output: tips and tricks for efficient implementation in generalized linear latent variable models

The marginal likelihood can be notoriously difficult to compute, and particularly so in high-dimensional problems. Chib and Jeliazkov employed the local reversibility of the Metropolis–Hastings algorithm to construct an estimator in models where full conditional densities are not available analytically. The estimator is free of distributional assumptions and is directly linked to the simulation algorithm. However, it generally requires a sequence of reduced Markov chain Monte Carlo runs which makes the method computationally demanding especially in cases when the parameter space is large. In this article, we study the implementation of this estimator on latent variable models which embed independence of the responses to the observables given the latent variables (conditional or local independence). This property is employed in the construction of a multi-block Metropolis-within-Gibbs algorithm that allows to compute the estimator in a single run, regardless of the dimensionality of the parameter space. The counterpart one-block algorithm is also considered here, by pointing out the difference between the two approaches. The paper closes with the illustration of the estimator in simulated and real-life data sets.     

Accounting for Uncertainty in Heteroscedasticity in Nonlinear Regression

Toxicologists and pharmacologists often describe toxicity of a chemical using parameters of a nonlinear regression model. Thus estimation of parameters of a nonlinear regression model is an important problem. The estimates of the parameters and their uncertainty estimates depend upon the underlying error variance structure in the model. Typically, a priori the researcher would not know if the error variances are homoscedastic (i.e., constant across dose) or if they are heteroscedastic (i.e., the variance is a function of dose). Motivated by this concern, in this paper we introduce an estimation procedure based on preliminary test which selects an appropriate estimation procedure accounting for the underlying error variance structure. Since outliers and influential observations are common in toxicological data, the proposed methodology uses M-estimators. The asymptotic properties of the preliminary test estimator are investigated; in particular its asymptotic covariance matrix is derived. The performance of the proposed estimator is compared with several standard estimators using simulation studies. The proposed methodology is also illustrated using a data set obtained from the National Toxicology Program.     

A comparison of maximum likelihood and least squares for the estimation of a cumulative distribution

Monte Carlo methods are used to compare the methods of maximum likelihood and least squares to estimate a cumulative distribution function. When the probabilistic model used is correct or nearly correct, the two methods produce similar results with the MLE usually slightly superior When an incorrect model is used, or when the data is contaminated, the least squares technique often gives substantially superior results.     

An integrated model based on statistical process control and maintenance for two-stage dependent processes

The close relationship between quality and maintenance of manufacturing systems has contributed to the development of integrated models which use the concept of statistical process control (SPC) and maintenance. This article demonstrates the integration of the Shewhart individual-residual (Z_X ? Z_e) joint control chart and maintenance for two-stage dependent processes by jointly optimizing their policies to minimize the expected total costs associated with quality, maintenance and inspection. To evaluate the effectiveness of the proposed model, two stand-alone models—a maintenance model and an SPC model—are proposed. Then a numerical example is given to illustrate the application of the proposed integrated model. The results show that the integrated model outperforms the two stand-alone models with regard to the expected cost per unit time. Finally, a sensitivity analysis is conducted to develop insights into time parameters and cost parameters that influence the integration efforts.     

Confession Governance and Verdicts in China’s Criminal Justice

The objective phenomenon of reliance on confessions is the source of many problems in the practice of criminal justice in China. Although successive amendments to the law have endeavored to improve the handling of confessions, they have been unable to resolve the dilemma of a rigid management mode, incompatible management strategies, and conflicting attitudes towards management. This paper has made a multi-dimensional reinterpretation of it from the judges’ perspective, including functional analysis, cultural interpretation and immediate considerations. It reveals that behind its functions of offering moral comfort or allowing evasion of responsibility, confession, as a medium of physical and intellectual management, has the more important function of supplementing and improving the legitimacy of criminal verdicts. The remolding of the legitimacy of criminal verdicts with the coordinated improvement of the management of confessions will change the mode of association between those who currently govern confession and those they govern and will clash with the current closed mode of criminal justice management, thereby promoting a transformation in the governance of criminal justice.