On the LBI criterion for the multivariate one-way random effects model under non-normality

The asymptotic null distribution of the locally best invariant (LBI) test criterion for testing the random effect in the one-way multivariable analysis of variance model is derived under normality and non-normality. The error of the approximation is characterized as O(1/n). The non-null asymptotic distribution is also discussed. In addition to providing a way of obtaining percentage points and p-values, the results of this paper are useful in assessing the robustness of the LBI criterion. Numerical results are presented to illustrate the accuracy of the approximation.     

Effects of non-normality on the power function in a one-way random model

The effects of non-normality on type-I and type-II errors in a one-way random model are investigated for moderate departures from normality. It is found that the probabilities of both errors are more sensitive to the kurtosis of between group effects than that of within group effects.     

Estimation of variance components in an unbalanced one-way classification

In the unbalanced one-way random effects model the weighted least squares approach with estimated weights is used to develop a relatively simple estimator of variance components. As the number of classes increases, the proposed estimator is seen not only to be best asymptotically normal but also to be asymptotically equivalent to the maximum likelihood estimator.     

Confidence intervals for the between group variance in the unbalanced one-way random effects model of anaylsis of variance

A confidence interval for the between group variance is proposed which is deduced from Wald'sexact confidence interval for the rtio of the two variance components in the one-way random effects model and the exact confidence interval for the error variance resp.an unbiased estimator of the error variance. In a simulation study the confidence coeffecients for these two intervals are compared with the confidence coefficients of two other commonly used confidence intervals. There the confidence interval derived here yields confidence coefficiends which are always greater than the prescriped level.     

Identification of the variance components in the general two-variance linear model

Bayesian analyses frequently employ two-stage hierarchical models involving two-variance parameters: one controlling measurement error and the other controlling the degree of smoothing implied by the model's higher level. These analyses can be hampered by poorly identified variances which may lead to difficulty in computing and in choosing reference priors for these parameters. In this paper, we introduce the class of two-variance hierarchical linear models and characterize the aspects of these models that lead to well-identified or poorly identified variances. These ideas are illustrated with a spatial analysis of a periodontal data set and examined in some generality for specific two-variance models including the conditionally autoregressive (CAR) and one-way random effect models. We also connect this theory with other constrained regression methods and suggest a diagnostic that can be used to search for missing spatially varying fixed effects in the CAR model.     

A likelihood ratio test for bimodality in two-component mixtures with application to regional income distribution in the EU

We propose a parametric test for bimodality based on the likelihood principle by using two-component mixtures. The test uses explicit characterizations of the modal structure of such mixtures in terms of their parameters. Examples include the univariate and multivariate normal distributions and the von Mises distribution. We present the asymptotic distribution of the proposed test and analyze its finite sample performance in a simulation study. To illustrate our method, we use mixtures to investigate the modal structure of the cross-sectional distribution of per capita log GDP across EU regions from 1977 to 1993. Although these mixtures clearly have two components over the whole time period, the resulting distributions evolve from bimodality toward unimodality at the end of the 1970s.     

Identification of outliers in a one-way random effects model

We distinguish between three types of outliers in a one-way random effects model. These are formally described in terms of their position relative to the main part of the observations. We propose simple rules for identifying such outliers and give an example which involves median-based statistics.     

A Bayesian compartmental model for the evaluation of 1,3-butadiene metabolism

Summary. We propose a Bayesian model for physiologically based pharmacokinetics of 1,3-butadiene (BD). BD is classified as a suspected human carcinogen and exposure to it is common, especially through cigarette smoke as well as in urban settings. The main aim of the methodology and analysis that are presented here is to quantify variability in the rates of BD metabolism by human subjects. A three-compartmental model is described, together with informative prior distributions for the population parameters, all of which represent real physiological variables. The model is described in detail along with the meanings and interpretations of the associated parameters. A four-compartment model is also given for comparison. Markov chain Monte Carlo methods are described for fitting the model proposed. The model is fitted to toxicokinetic data obtained from 133 healthy subjects (males and females) from the four major racial groups in the USA, with ages ranging from 19 to 62 years. Subjects were exposed to 2 parts per million of BD for 20 min through a face mask by using a computer-controlled exposure and respiratory monitoring system. Stratification by ethnic group results in major changes in the physiological parameters. Sex and age were also tested but not found to have a significant effect.     

Improved estimation of variance components in balanced hierarchical mixed models

The problem of simultaneous estimation of variance components is considered for a balanced hierarchical mixed model under a sum of squared error loss. A new class of estimators is suggested which dominate the usual sensible estimators. These estimators shrink towards the geometric mean of the component mean squares that appear in the ANOVA table. Numerical results are tabled to exhibit the improvement in risk under a simple model.     

The Superiorities of Empirical Bayes Estimation of Variance Components in Random Effects Model

In this article, the Bayes estimators of variance components are derived and the parametric empirical Bayes estimators (PEBE) for the balanced one-way classification random effects model are constructed. The superiorities of the PEBE over the analysis of variance (ANOVA) estimators are investigated under the mean square error (MSE) criterion, some simulation results for the PEBE are obtained. Finally, a remark for the main results is given.     

Between group variance component interval estimation for the unbalanced heteroscedastic one-way random effects model

The studied topic is motivated by the problem of interlaboratory comparisons. This paper focuses on the confidence interval estimation of the between group variance in the unbalanced heteroscedastic one-way random effects model. Several interval estimators are proposed and compared by means of the simulation study. The most recommended (safest) is the confidence interval based on Bonferroni's inequality.     

Independent Quadratic Forms in 3-Variance-Component Models

This article addresses derivation and existence of quadratic forms that were suggested by Burch (2007 Burch , B. D. ( 2007 ). Generalized confidence intervals for proportions of total variance in mixed linear models . J. Statist. Plann. Infer. 137 : 2394 – 2404 .[Crossref], [Web of Science ®] , [Google Scholar]) for procedures for inference on variance components in mixed linear models in combination with generalized fiducial inference. A relatively simple algorithm leading to the required quadratic forms in a general 3-variance-component model is stated and designs for two-way ANOVA models without interactions that permit Burch's procedure are characterized. This complements developments in the original article by Burch.     

The equivalence between likelihood ratio test and F-test for testing variance component in a balanced one-way random effects model

In mixed linear models, it is frequently of interest to test hypotheses on the variance components. F-test and likelihood ratio test (LRT) are commonly used for such purposes. Current LRTs available in literature are based on limiting distribution theory. With the development of finite sample distribution theory, it becomes possible to derive the exact test for likelihood ratio statistic. In this paper, we consider the problem of testing null hypotheses on the variance component in a one-way balanced random effects model. We use the exact test for the likelihood ratio statistic and compare the performance of F-test and LRT. Simulations provide strong support of the equivalence between these two tests. Furthermore, we prove the equivalence between these two tests mathematically.     

A multilevel Bayesian model for contextual effect of material deprivation

The relationship between socioeconomic factors and health has been studied in many circumstances. Whether the association takes place at individual level only, or also at population level (contextual effect) is still unclear. We present a multilevel hierarchical Bayesian model to investigate the joint contribution of individual and population-based socioeconomic factors to mortality, using data from the census cohort of the general population of the city of Florence, Italy (Tuscany Longitudinal Study, 1991-1995). Evidence supporting a contextual effect of deprivation on mortality at the very fine level of aggregation is found. Inappropriate modelling of individual and aggregate variables could strongly bias effect estimates.Received: 10 January 2002, Revised: 23 June 2003, The research on Tuscany Longitudinal Study (Studio Longitudinale Toscano, SLTo) was supported by the Regione Toscana Servizio Statistica.     

Variance components estimation for the balanced random effects model under mixed prior distributions

For the balanced random effects models, when the variance components are correlated either naturally or through common prior structures, by assuming a mixed prior distribution for the variance components, we propose some new Bayesian estimators. To contrast and compare the new estimators with the minimum variance unbiased (MVUE) and restricted maximum likelihood estimators (RMLE), some simulation studies are also carried out. It turns out that the proposed estimators have smaller mean squared errors than the MVUE and RMLE.     

Improved inference for a linear mixed-effects model when the subpopulation effects are clustered

We provide Bayesian methodology to relax the assumption that all subpopulation effects in a linear mixed-effects model have, after adjustment for covariates, a common mean. We expand the model specification by assuming that the m subpopulation effects are allowed to cluster into d groups where the value of d, 1?d?m, and the composition of the d groups are unknown, a priori. Specifically, for each partition of the m effects into d groups we only assume that the subpopulation effects in each group are exchangeable and are independent across the groups. We show that failure to take account of this clustering, as with the customary method, will lead to serious errors in inference about the variances and subpopulation effects, but the proposed, expanded, model leads to appropriate inferences. The efficacy of the proposed method is evaluated by contrasting it with both the customary method and use of a Dirichlet process prior. We use data from small area estimation to illustrate our method.     

Posterior Predictive Comparisons for the Two-sample Problem

The two-sample problem of inferring whether two random samples have equal underlying distributions is formulated within the Bayesian framework as a comparison of two posterior predictive inferences rather than as a problem of model selection. The suggested approach is argued to be particularly advantageous in problems where the objective is to evaluate evidence in support of equality, along with being robust to the priors used and being capable of handling improper priors. Our approach is contrasted with the Bayes factor in a normal setting and finally, an additional example is considered where the observed samples are realizations of Markov chains.     

Information analysis of local suppression scheme based on a spatial-temporal model

In a wireless sensor network, data collection is relatively cheap whereas data transmission is relatively expensive. Thus, preserving battery life is critical. If the process of interest is sufficiently predictable, the suppression in transmission can be adopted to improve efficiency of sensor networks because the loss of information is not great. The prime interest lies in finding an inference-efficient way to support suppressed data collection application. In this paper, we present a suppression scheme for a multiple nodes setting with spatio-temporal processes, especially when process knowledge is insufficient. We also explore the impact of suppression schemes on the inference of the regional processes under various suppression levels. Finally, we formalize the hierarchical Bayesian model for these schemes.     

A bayesian analysis of the mixed linear model

The mixed model is defined. The exact posterior distribution for the fixed effect vector is obtained. The exact posterior distribution for the error variance is obtained. The exact posterior mean and variance of a Bayesian estimator for the variances of random effects is also derived. All computations are non-iterative and avoid numerical integrations.     

Nonparametric Bayesian analysis for multi-site hidden Markov model

The hidden Markov model (HMM) provides an attractive framework for modeling long-term persistence in a variety of applications including pattern recognition. Unlike typical mixture models, hidden Markov states can represent the heterogeneity in data and it can be extended to a multivariate case using a hierarchical Bayesian approach. This article provides a nonparametric Bayesian modeling approach to the multi-site HMM by considering stick-breaking priors for each row of an infinite state transition matrix. This extension has many advantages over a parametric HMM. For example, it can provide more flexible information for identifying the structure of the HMM than parametric HMM analysis, such as the number of states in HMM. We exploit a simulation example and a real dataset to evaluate the proposed approach.