The Residual Likelihood Ratio Test for the Variance Ratio in a Linear Model with Two Variance Components

Residual maximum likelihood has proved to be a successful approach to the estimation of variance components. In this paper, its counterpart in testing, the residual likelihood ratio test, is applied to testing the ratio of two variance components. The test is compared with the Wald test and the locally most powerful invariant test.     

A Note on Maximum Likelihood Ratio Test of No Outliers in Regression Models

A derivation of the maximum likelihood ratio test for testing no outliers in regression models is given using the method of WETHERILL (1981, pp. 106–107) for estimating the regression parameters. This method is essentially similar to the one outlined in BARNETT and LEWIS (1978, p. 263), although by our detailed derivation it is easier to see that the maximum likelihood estimate of θ of model (3) under the hypothesis that the i^th observation in an outlier is the same as that obtained from model (1) when the i^th observation is removed.     

Finite Sample Properties of the Maximum Likelihood Estimator and of Likelihood Ratio Tests in Hidden Markov Models

Hidden Markov models were successfully applied in various fields of time series analysis, especially for analyzing ion channel recordings. The maximum likelihood estimator (MLE) has recently been proven to be asymptotically normally distributed. Here, we investigate finite sample properties of the MLE and of different types of likelihood ratio tests (LRTs) by means of simulation studies. The MLE is shown to reach the asymptotic behavior within sample sizes that are common for various applications. Thus, reliable estimates and confidence intervals can be obtained. We give an approximative scaling function for the estimation error for finite samples, and investigate the power of different LRTs suitable for applications to ion channels, including tests for superimposed hidden Markov processes. Our results are applied to physiological sodium channel data.     

Accuracy and Power of the Likelihood Ratio Test for Comparing Evolutionary Rates Among Genes

Sequences for multiple protein-coding genes are now commonly available from several, often closely related species. These data sets offer intriguing opportunities to test hypotheses regarding whether different types of genes evolve under different selective pressures. Although maximum likelihood (ML) models of codon substitution that are suitable for such analyses have been developed, little is known about the statistical properties of these tests. We use a previously developed fixed-sites model and computer simulations to examine the accuracy and power of the likelihood ratio test (LRT) in comparing the nonsynonymous-to-synonymous substitution rate ratio (=dN/dS) between two genes. Our results show that the LRT applied to fixed-sites models may be inaccurate in some cases when setting significance thresholds using a ² approximation. Instead, we use a parametric bootstrap to describe the distribution of the LRT statistic for fixed-sites models and examine the power of the test as a function of sampling variables and properties of the genes under study. We find that the power of the test is high (>80%) even when sampling few taxa (e.g., six species) if sequences are sufficiently diverged and the test is largely unaffected by the tree topology used to simulate data. Our simulations show fixed-sites models are suitable for comparing substitution parameters among genes evolving under even strong evolutionary constraint ( 0.05), although relative rate differences of 25% or less may be difficult to detect.Reviewing Editor: Dr. Rosmus Nielsen     

Use of the Likelihood Ratio Test on the Uniform Distribution

In this article we give a simple procedure to determine the exact distribution of the likelihood ratio test of a statistical hypothesis regarding the parameter of the uniform distribution. The resulting distribution will be shown to serve as an approximation to the distribution of the likelihood ratio statistic for testing the equality of scale parameters of k independent Exponential populations.     

Bootstrap Restricted Likelihood Ratio Test for the Detection of Rare Variants

In this paper the detection of rare variants association with continuous phenotypes of interest is investigated via the likelihood-ratio based variance component test under the framework of linear mixed models. The hypothesis testing is challenging and nonstandard, since under the null the variance component is located on the boundary of its parameter space. In this situation the usual asymptotic chisquare distribution of the likelihood ratio statistic does not necessarily hold. To circumvent the derivation of the null distribution we resort to the bootstrap method due to its generic applicability and being easy to implement. Both parametric and nonparametric bootstrap likelihood ratio tests are studied. Numerical studies are implemented to evaluate the performance of the proposed bootstrap likelihood ratio test and compare to some existing methods for the identification of rare variants. To reduce the computational time of the bootstrap likelihood ratio test we propose an effective approximation mixture for the bootstrap null distribution. The GAW17 data is used to illustrate the proposed test.     

The Wald Statistic in Proportional Hazards Hypothesis Testing

In survivorship modelling using the proportional hazards model of Cox (1972, Journal of the Royal Statistical Society, Series B, 34, 187–220), it is often desired to test a subset of the vector of unknown regression parameters β in the expression for the hazard rate at time t. The likelihood ratio test statistic is well behaved in most situations but may be expensive to calculate. The Wald (1943, Transactions of the American Mathematical Society 54, 426–482) test statistic is easier to calculate, but has some drawbacks. In testing a single parameter in a binomial logit model, Hauck and Donner (1977, Journal of the American Statistical Association 72, 851–853) show that the Wald statistic decreases to zero the further the parameter estimate is from the null and that the asymptotic power of the test decreases to the significance level. The Wald statistic is extensively used in statistical software packages for survivorship modelling and it is therefore important to understand its behavior. The present work examines empirically the behavior of the Wald statistic under various departures from the null hypothesis and under the presence of Type I censoring and covariates in the model. It is shown via examples that the Wald statistic's behavior is not as aberrant as found for the logistic model. For the single parameter case, the asymptotic non-null distribution of the Wald statistic is examined.     

Testing the Equality of the Intra‐Subject Treatment Covariance Matrices in a Crossover Study

In a 2 × 2 crossover bioavailability study, the sets of estimates of the pharmacokinetic parameters quite often have a symmetric covariance structure between the two treatments. For testing the equality of the intra‐subject covariance matrices for the two treatments in such studies, we suggest in this paper some statistical tests. When the response vectors are bivariate, we propose an exact test. Since the statistical procedures depend on the assumption of a symmetric covariance structure between the two treatments, we put forth some statistical tests for this assumption. We then apply the discussed tests to real data from a crossover bioavailability trial.     

Testing for group effect in a 2 x k heteroscedastic ANOVA model

It is natural to want to relax the assumption of homoscedasticity and Gaussian error in ANOVA models. For a two-way ANOVA model with 2 x k cells, one can derive tests of main effect for the factor with two levels (referred to as group) without assuming homoscedasticity or Gaussian error. Empirical likelihood can be used to derive testing procedures. An approximate empirical likelihood ratio test (AELRT) is derived for the test of group main effect. To approximate the distributions of the test statistics under the null hypothesis, simulation from the approximate empirical maximum likelihood estimate (AEMLE) restricted by the null hypothesis is used. The homoscedastic ANOVA F -test and a Box-type approximation to the distribution of the heteroscedastic ANOVA F -test are compared to the AELRT in level and power. The AELRT procedure is shown by simulation to have appropriate type I error control (although possibly conservative) when the distribution of the test statistics are approximated by simulation from the constrained AEMLE. The methodology is motivated and illustrated by an analysis of folate levels in the blood among two alcohol intake groups while accounting for gender.     

On the use of historical control information for trend test in carcinogenesis

In this article a new non-model-based significance test for detecting dose-response relationship with the incorporation of historical control data is proposed. This non-model-based test is considered simpler from a regulatory perspective because it does not require validating any modeling assumptions. Moreover, our test is especially appropriate to those studies in which the intravenous doses for the investigational chemical are labeled as, e.g., low, medium and high or the dose labels do not suggest any obvious choices of dose scores. This test can be easily adopted for detecting general dose-response shape, such as an umbrella pattern. Simple adjustments will be proposed for better control of the actual Type I error. Data sets from two carcinogenesis studies will be used to illustrate our method. We also evaluate the performance of the proposed test and the famous model-based Tarone's trend test with respect to size and power.