Two-Stage Least Squares Analysis of Hierarchical Linear Models

This paper has two major objectives. The first is to present a two-stage least squares procedure for estimation of the parameters in a linear model whose parameters are in themselves linear functions of some hyperparameters. The second, and perhaps more important point, is that the new estimator can be shown to be generally more precise than either the Bayesian or the generalized single-stage least squares estimator reported by LINDLEY and SMITH (1072).     

On the Optimality of F-test on Variance Components in Mixed Linear Models

The analysis-of-variance tests for hypotheses on random effects in regular linear models are considered. Conditions are given for these tests to be uniformly most powerful unbiased or uniformly most powerful invariant unbiased. An example shows that the difference between these conditions can be serious.     

Unbiasedness of the Estimator of the Function of Expected Value in the Mixed Linear Model

The traditional method for estimating the linear function of fixed parameters in mixed linear model is a two-stage procedure. In the first stage of this procedure the variance components estimators are calculated and next in the second stage these estimators are taken as true values of variance components to estimating the linear function of fixed parameters according to generalized least squares method. In this paper the general mixed linear model is considered in which a matrix related to fixed parameters and or/a dispersion matrix of observation vector may be deficient in rank. It is shown that the estimators of a set of functions of fixed parameters obtained in second stage are unbiased if only the observation vector is symmetrically distributed about its expected value and the estimators of variance components from first stage are translation-invariant and are even functions of the observation vector.     

On the mean square error associated with adaptive generalized ridge regression

Some numerical results are presented for generalized ridge regression where the additive constants are based on the data. The adaptive estimator so obtained is compared with the least-squares estimator on the basis of mean square error (MSE). It is shown that the MSE of each component of the vector of ridge estimators may be as low as 47.1% of the variance of the corresponding component of the least squares vector or as high as 125.2%.     

具有亲缘关系矩阵的混合模型的方差组分迭代估计——应用于计算BLUP的方差比值和遗传参数

论述的是来自非均街资料的混合模型中具有亲缘关系矩阵时利用迭代法估计方差组分问题。这篇文章表明计算程序是可行的,只要能够按照混合模型中固定效应的结构矩阵和Henderson方法3的固定效应的假设条件正确地计算二次型约化平方和,就可获得较为精确的方差组分估计值;而且表明方差初始比值k偏高或偏低,不影响迭代求解的最后结果,这是因为在迭代过程中可以通过结构矩阵x'x和x'x的控制而自行调整。这些方差组分不仅可应用于选种种畜用的BLUP计算,还可用来估计遗传参数。     

Non-Iterative Least Squares Estimation of Missing Values in Hyper-Graeco-Latin Square Designs

In this paper an attempt has been made to obtain explicit expressions for the estimators of the several missing values in hyper-graeco-latin square designs. Further it has been shown that the estimates of the missing values in latin square designs and graeco-latin square designs are obtained as a particular case of the estimates of the missing values in hyper-graeco-latin square designs.     

A New Algorithm for Explicit Determination of Variance Component Estimations in Mixed Models and Mixed Classifications of Balanced ANOVA

This paper deals with the balanced case of the analysis of variance. The use of a classification function leads to an easy determination of all possible sources of variation of any mixed classification. For mixed models a new method is derived, which allows to represent explicit the ANOVA-estimations of the variance components respectively the estimation of the mean sum of squares of the fixed effects for all sources of variation. Thereby the corresponding F-quotients and the approximate confidence intervals of variance components are received in a simple way.     

Estimation of a Conditional Mean in a Linear Regression Model after a Preliminary Test on Regression Coefficient

Consider the two linear regression models of Y_ij on X_ij, namely Y_ij = β_io + β_ij, X_ij + E_ij = 1, 2,…, n_i, i = 1, 2, where E_ij are assumed to be normally distributed with zero mean and common unknown variance σ². The problem of estimating the conditional mean of Y₁ for a given value of X₁ is considered when it is a priori suspected that β₁₀ = β₂₀ and β₁₁ = β₂₁. The preliminary test estimator is proposed. The exact expressions for the bias and the mean square error of the estimator are derived. The relative efficiency of the new estimator to the usual least square estimator based on the first regression alone is computed and is used to determine the appropriate value of the significance level of the preliminary test β₁₀ = β₂₀ and β₁₁ = β₂₁.     

Class of Estimators of Mean and Variance Using Auxiliary Information when Correlation Coefficient is also Known

For estimating the mean of a finite population using information on an auxiliary variable, a class of estimators which also uses the value of the correlation coefficient between the two variables which is assumed known, is defined. Expression for its asymptotic mean squared error and its minimum value is obtained. An expression by which the minimum mean squared error of this class is smaller than those which use only the sample mean and the sample variance of the auxiliary variable is obtained. A similar class of estimators is considered for the estimation of the population variance. The gain in efficiency is illustrated for six populations considered in literature.     

A Note on Balanced Sampling and Robust Genetic Predictors

The concept of balanced sampling is applied to prediction in finite samples using model based inference procedures. Necessary and sufficient conditions are derived for a general linear model with arbitrary covariance structure to yield the expansion estimator as the best linear unbiased predictor for the mean. The analysis is extended to produce a robust estimator for the mean squared error under balanced sampling and the results are discussed in the context of statistical genetics where appropriate sampling produces simple efficient and robust genetic predictors free from unnecessary genetic assumptions.     

The Applicability of Ordinary Least Squares to Consistently Short Distances between Taxa in Phylogenetic Tree Construction and the Normal Distribution Test Consequences

Short phylogenetic distances between taxa occur, for example, in studies on ribosomal RNA-genes with slow substitution rates. For consistently short distances, it is proved that in the completely singular limit of the covariance matrix ordinary least squares (OLS) estimates are minimum variance or best linear unbiased (BLU) estimates of phylogenetic tree branch lengths. Although OLS estimates are in this situation equal to generalized least squares (GLS) estimates, the GLS chi-square likelihood ratio test will be inapplicable as it is associated with zero degrees of freedom. Consequently, an OLS normal distribution test or an analogous bootstrap approach will provide optimal branch length tests of significance for consistently short phylogenetic distances. As the asymptotic covariances between branch lengths will be equal to zero, it follows that the product rule can be used in tree evaluation to calculate an approximate simultaneous confidence probability that all interior branches are positive.     

选择回归方程自变量的条件数法及其在RK手术中的应用

选择合适的自变量是确定线性回归模型的首要问题,本文以消除自变量之间的复共线性为目标,介绍了一种选择回归方程自变量的条件数法,并在ＲＫ手术的结果预测问题中采用了这一方法。     

The Estimation of Skewness and Kurtosis of Random Effects in the Linear Model

Generalising the ANOVA method of estimating variance components in mixed linear models a simple procedure is presented to estimate skewness and kurtosis of the distributions of the random effects of the model. For the model II of a one-way classification this procedure is demonstrated explicitly.     

S形增长模型之比较、组合预测及应用

对５种Ｓ形增长模型进行比较,以组合模型的误差平方和最小为目标建立组合预测模型,并以实例说明其应用。     

Use of replicates in statistical analyses in papers submitted for publication in Animal Feed Science and Technology

This editorial is an annex, essentially an amplification of earlier Animal Feed Science and Technology (AFST) editorials which focused on analytical and statistical issues of papers submitted for publication in AFST. This amplification is needed because, since publication of those editorials, there has been a sharp increase in papers submitted to AFST (particularly in vitro gas, ruminal in sacco, continuous culture fermenter and mini-silo studies), in which there have been substantive disagreements among authors, reviewers and editors as to what, exactly, constitutes acceptable statistical replicates as opposed to pseudo replicates. In this editorial, the Co-Editors in Chief of AFST provide clarification of their view relative to use of analytical observations in statistical analyses in papers submitted for consideration for publication in AFST. If the objective is to compare feeds and from the results, make inferences to populations, only multiple samples of each feed represent an acceptable feed base. This suggests that means comparisons based on repeated assays of the same sample, whether by chemical, physical or microbiological methods, will almost certainly be rejected by the CEIC.     

Estimation and Verification of Hypotheses in Some Zyskind-Martin Models with Missing Values

Several theorems on estimation and verification of linear hypotheses in some Zyskind-Martin (ZM) models are given. The assumptions are as follows. Let y = Xβ + e or (y, Xβ, σ²V) be a fixed model where y is a vector of n observations, X is a known matrix nXp with rank r(X) = r ≦ p < n, where p is a number of coordinates of the unknown parameter vector β, e is a random vector of errors with covariance matrix σ²V, where σ² is unknown scalar parameter, V is a known non-negative definite matrix such that R(X) ? R(V). Symbol R(A) denotes a vector space generated by columns of matrix A. The expected value of y is Xβ. In this paper four following Zyskind-Martin (ZM) models are considered: ZMd, ZMa, ZMc and ZMqd (definitions in sec. 1) when vector y y₁ y₂ involves a vector y₁ of m missing values and a vector y₂ with (n — m) observed values. A special transformation of ZM model gives again ZM model (cf. theorem 2.1). Ten properties of actual (ZMa) and complete (ZMc) Zyskind-Martin models with missing values (cf. theorem 2.2) test functions F are given in (2.11)) are presented. The third propriety constitutes a generalization of R. A. Fisher's rule from standard model (y, Xβ, σ²I) to ZM model. Estimation of vector y₁ (cf. 3.3) of vector β (cf. th. 3.2) and of scalar σ² (cf. th. 3.4) in actual ZMa model and in diagonal quasi-ZM model (ZMqd) are presented. Relation between y? ₁ and β is given in theorem 3.1. The results of section 2 are illustrated by numerical example in section 4.