On the Effect of Unbalancedness and Heteroscedasticity on the ANOVA Estimator of Group Variance Component in One-Way Random Model

The expression for rth cumulant of ANOVA estimator of group variance component is derived in the One-way unbalanced random model under heteroscedasticity. The expression is used to study the effect of unbalancedness and heteroscedasticity on the mean and variance of the estimator, numerically. The computed results reveal that the unbalancedness and heteroscedasticity have a combined effect on the mean and variance of the estimator. For certain situations of unequal group sizes and error variances, the mean and variance of the estimator are increased and for certain other situations the values are decreased.     

Estimation of Heritability in Heteroscedastic One‐Way Unbalanced Random Model

A modified estimator of heritability is proposed under heteroscedastic one way unbalanced random model. The distribution, moments and probability of permissible values (PPV) for conventional and modified estimators are derived. The behaviour of two estimators has been investigated, numerically, to devise a suitable estimator of heritability under variance heterogeneity. The numerical results reveal that under balanced case the heteroscedasticity affects the bias, MSE and PPV of conventional estimator, marginally. In case of unbalanced situations, the conventional estimator underestimates the parameter when more variable group has more observations and overestimates when more variable group has less observations, MSE of the conventional estimator decreases when more variable group has more observations and increases when more variable group has less observations and PPV is marginally decreased. The MSE and PPV are comparable for two estimators while the bias of modified estimator is less than the conventional estimator particularly for small and medium values of the parameter. These results suggest the use of modified estimator with equal or more observations for more variable group in presence of variance heterogeneity.     

The dilemma of negative analysis of variance estimators of intraclass correlation

Summary At least two common practices exist when a negative variance component estimate is obtained, either setting it to zero or not reporting the estimate. The consequences of these practices are investigated in the context of the intraclass correlation estimation in terms of bias, variance and mean squared error (MSE). For the one-way analysis of variance random effects model and its extension to the common correlation model, we compare five estimators: analysis of variance (ANOVA), concentrated ANOVA, truncated ANOVA and two maximum likelihood-like (ML) estimators. For the balanced case, the exact bias and MSE are calculated via numerical integration of the exact sample distributions, while a Monte Carlo simulation study is conducted for the unbalanced case. The results indicate that the ANOVA estimator performs well except for designs with family size n = 2. The two ML estimators are generally poor, and the concentrated and truncated ANOVA estimators have some advantages over the ANOVA in terms of MSE. However, the large biases may make the concentrated and truncated ANOVA estimators objectionable when intraclass correlation () is small. Bias should be a concern when a pooled estimate is obtained from the literature since <0.05 in many genetic studies.     

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%.     

Unbalancedness and Efficiency in Estimating Components of Variance: MINQUE and ANOVA Procedure

We consider the one-way classification random model. Using a measure of unbalancedness some MINQ-Estimation procedures and ANOVA-procedure have been investigated with respect to their sensibility to unbalancedness. Further if we define the efficiency of a design N as the ratio of the MSE (.)'s of the balanced and the unbalanced case we obtain for three MINQE-procedures upper and lower bounds of the minimum efficiency.     

Estimating effective population size or mutation rate with microsatellites

Microsatellites are short tandem repeats that are widely dispersed among eukaryotic genomes. Many of them are highly polymorphic; they have been used widely in genetic studies. Statistical properties of all measures of genetic variation at microsatellites critically depend upon the composite parameter theta = 4Nmicro, where N is the effective population size and micro is mutation rate per locus per generation. Since mutation leads to expansion or contraction of a repeat number in a stepwise fashion, the stepwise mutation model has been widely used to study the dynamics of these loci. We developed an estimator of theta, theta; (F), on the basis of sample homozygosity under the single-step stepwise mutation model. The estimator is unbiased and is much more efficient than the variance-based estimator under the single-step stepwise mutation model. It also has smaller bias and mean square error (MSE) than the variance-based estimator when the mutation follows the multistep generalized stepwise mutation model. Compared with the maximum-likelihood estimator theta; (L) by, theta; (F) has less bias and smaller MSE in general. theta; (L) has a slight advantage when theta is small, but in such a situation the bias in theta; (L) may be more of a concern.     

On the small-sample performance of Efron's and of Gill's version of the product limit estimator under nonproportional hazards

Various proposals have been made to extend the product limit estimator to survival times beyond the largest observation in case that observation is censored. Two extreme extensions are examined with respect to bias and mean squared error (MSE). Their quality depends considerably on the "censoring constellation." The MSE of one extension appears to be robust against a wide variety of nonproportionalities of the hazard rates of the distributions of lifelength and time to censoring.     

A comparison of several point estimators of the odds ratio in a single 2 x 2 contingency table.

The relative performance of the unconditioned maximum likelihood estimators (UMLEs), conditional MLEs (CMLEs), and Jewell-type estimators of the odds ratio (OR) and its logarithm were investigated in sets of single 2 x 2 contingency tables. The tables were generated by complete enumeration of all possible cell frequencies consistent with a single fixed margin. The bias, mean squared error (MSE), and average absolute error (AAE) were computed for all estimators using the individual table probabilities as weights. The results showed that, for the OR, Jewell's estimator usually had smaller bias, MSE, and AAE than either of the MLEs. While the differences were often slight for MSE and AAE, for bias it was sometimes substantial. For the log(OR), the UMLE usually had the lowest bias, and its MSE and AAE were only slightly greater than those for the other estimators. Overall, we recommend estimation on the log scale using the UMLE. If OR is to be estimated, Jewell's method had strong merit, although it is nonsymmetric with respect to the table orientation. In view of this, the UMLE may again be favoured in some situations.     

Incorporating Estimates of Group Size in Sightability Models for Wildlife

ABSTRACT Sightability models have been used to estimate population size of many wildlife species; however, a limitation of these models is an assumption that groups of animals observed and counted during aerial surveys are enumerated completely. Replacing these unknown counts with maximum observed counts, as is typically done, produces population size estimates that are negatively biased. This bias can be substantial depending on the degree of undercounting occurring. We first investigated a method-of-moments estimator of group sizes. We then defined a population size estimator using the method-of-moments estimator of group sizes in place of maximum counts in the traditional sightability models, thereby correcting for bias associated with undercounting group size. We also provide associated equations for calculating the variance of our estimator. This estimator is an improvement over existing sightability model techniques because it significantly reduces bias, and variance estimates provide near nominal confidence interval coverage. The data needed for this estimator can be easily collected and implemented by wildlife managers with a field crew of only 3 individuals and little additional flight or personnel time beyond the normal requirements for developing sightability models.     

On Two Measures of Unbalancedness in a One-Way Model and Their Relation to Efficiency

This paper discusses two measures of unbalancedness in a one-way model and shows that for a given statistical procedure they may serve as measures of efficiency of a design. They also allow to compare for example estimation methods for variance components in designs with a fixed level of unbalancedness.     

Estimation of a Secondary Parameter in a Group Sequential Clinical Trial

In this article, we investigate estimation of a secondary parameter in group sequential tests. We study the model in which the secondary parameter is the mean of the normal distribution in a subgroup of the subjects. The bias of the naive secondary parameter estimator is studied. It is shown that the sampling proportions of the subgroup have a crucial effect on the bias: As the sampling proportion of the subgroup at or just before the stopping time increases, the bias of the naive subgroup parameter estimator increases as well. An unbiased estimator for the subgroup parameter and an unbiased estimator for its variance are derived. Using simulations, we compare the mean squared error of the unbiased estimator to that of the naive estimator, and we show that the differences are negligible. As an example, the methods of estimation are applied to an actual group sequential clinical trial, The Beta-Blocker Heart Attack Trial.     

On estimating the proportion of variance in a phenotypic trait attributable to a measured locus

The measured genotype approach can be used to estimate the variance contributions of specific candidate loci to quantitative traits of interest. We show here that both the naive estimate of measured-locus heritability, obtained by invoking infinite-sample theory, and an estimate obtained from a bias-corrected variance estimate based on finite-sample theory, produce biased estimates of heritability. We identify the sources of bias, and quantify their effects. The two sources of bias are: (1) the estimation of heritability from population samples as the ratio of two variances, and (2) the existence of sampling error. We show that neither heritability estimator is less biased (in absolute value) than the other in all situations, and the choice of an ideal estimator is therefore a function of the sample size and magnitude of the locus-specific contribution to the overall phenotypic variance. In most cases the bias is small, so that the practical implications of using either estimator are expected to be minimal.     

Combining multiple imputation and inverse-probability weighting

Two approaches commonly used to deal with missing data are multiple imputation (MI) and inverse-probability weighting (IPW). IPW is also used to adjust for unequal sampling fractions. MI is generally more efficient than IPW but more complex. Whereas IPW requires only a model for the probability that an individual has complete data (a univariate outcome), MI needs a model for the joint distribution of the missing data (a multivariate outcome) given the observed data. Inadequacies in either model may lead to important bias if large amounts of data are missing. A third approach combines MI and IPW to give a doubly robust estimator. A fourth approach (IPW/MI) combines MI and IPW but, unlike doubly robust methods, imputes only isolated missing values and uses weights to account for remaining larger blocks of unimputed missing data, such as would arise, e.g., in a cohort study subject to sample attrition, and/or unequal sampling fractions. In this article, we examine the performance, in terms of bias and efficiency, of IPW/MI relative to MI and IPW alone and investigate whether the Rubin's rules variance estimator is valid for IPW/MI. We prove that the Rubin's rules variance estimator is valid for IPW/MI for linear regression with an imputed outcome, we present simulations supporting the use of this variance estimator in more general settings, and we demonstrate that IPW/MI can have advantages over alternatives. IPW/MI is applied to data from the National Child Development Study.     

Estimating the encounter rate variance in distance sampling

Summary .  The dominant source of variance in line transect sampling is usually the encounter rate variance. Systematic survey designs are often used to reduce the true variability among different realizations of the design, but estimating the variance is difficult and estimators typically approximate the variance by treating the design as a simple random sample of lines. We explore the properties of different encounter rate variance estimators under random and systematic designs. We show that a design-based variance estimator improves upon the model-based estimator of Buckland et al. (2001, Introduction to Distance Sampling. Oxford: Oxford University Press, p. 79) when transects are positioned at random. However, if populations exhibit strong spatial trends, both estimators can have substantial positive bias under systematic designs. We show that poststratification is effective in reducing this bias.     

A random effects model for multivariate failure time data from multicenter clinical trials

A random effects model for analyzing multivariate failure time data is proposed. The work is motivated by the need for assessing the mean treatment effect in a multicenter clinical trial study, assuming that the centers are a random sample from an underlying population. An estimating equation for the mean hazard ratio parameter is proposed. The proposed estimator is shown to be consistent and asymptotically normally distributed. A variance estimator, based on large sample theory, is proposed. Simulation results indicate that the proposed estimator performs well in finite samples. The proposed variance estimator effectively corrects the bias of the naive variance estimator, which assumes independence of individuals within a group. The methodology is illustrated with a clinical trial data set from the Studies of Left Ventricular Dysfunction. This shows that the variability of the treatment effect is higher than found by means of simpler models.     

A Sometimes-Pool Estimator of the Mean Life

The paper deals with the problem of estimating the mean life θ in an exponential model \documentclass{article}\pagestyle{empty}\begin{document}$f(x/\theta ) = \frac{1}{\theta }e - x/\theta$\end{document}. It is assumed that in addition to the current ordered sample, we have a sample collected sometime in the recent past when the mean life might have been β. We have proposed a Sometimes—Pool procedure which is based on the outcome of a preliminary test of H₀: θ=β and obtained the expressions of the bias and MSE. An attempt has been made to locate that region in the parameter space in which the proposed estimator does better (in MSE sense) than the usual estimator based only on the current sample.     

Variance component estimation techniques compared for two mating designs with forest genetic architecture through computer simulation

Computer simulation was used to compare minimum variance quadratic estimation (MIVQUE), minimum norm quadratic unbiased estimation (MINQUE), restricted maximum likelihood (REML), maximum likelihood (ML), and Henderson's Method 3 (HM3) on the basis of variance among estimates, mean square error (MSE), bias and probability of nearness for estimation of both individual variance components and three ratios of variance components. The investigation also compared three procedures for dealing with negative estimates and included the use of both individual observations and plot means as the experimental unit of the analysis. The structure of data simulated (field design, mating designs, genetic architecture and imbalance) represented typical analysis problems in quantitative forest genetics. Results of comparing the estimation techniques demonstrated that: estimates of probability of nearness did not discriminate among techniques; bias was discriminatory among procedures for dealing with negative estimates but not among estimation techniques (except ML); sampling variance among estimates was discriminatory among procedures for dealing with negative estimates, estimation techniques and unit of observation; and MSE provided no additional information to variance of the estimates. HM3 and REML were the closest competitors under these criteria; however, REML demonstrated greater robustness to imbalance. Of the three negative estimate procedures, two are of practical significance and guidelines for their application are presented. Estimates from individual observations were always preferable to those from plot means over the experimental levels of this study.This is Journal Series NO. R-03768 of the Institute of Food and Agricultural Sciences     

A comparison of two methods for making statistical inferences on Nei's measure of genetic distance

The delta and jackknife methods can be used to estimate Nei's measure of genetic distance and calculate confidence intervals for this estimate. Computer stimulations were used to study the bias and variance of each estimator and the accuracy of the corresponding approximate 95% confidence intervals. The simulations were conducted using 3 sets of data and several sample sizes. The results showed: (1) the jackknife reduced bias; (2) in 8 out of 9 cases the variance and mean square error of the jackknife estimator were less; (3) a second order jackknife reduced the bias the most but suffered a corresponding increase in variance; (4) both the first order jackknife and delta methods yielded intervals whose confidence levels were approximately equal but less than 95%.     

Least Squares and Minimum MSE Estimators of Variance Components in Mixed Linear Models

The paper deals with the quadratic invariant estimators of the linear functions of variance components in mixed linear model. The estimator with locally minimal mean square error with respect to a parameter ? is derived. Under the condition of normality of the vector Y the theoretical values of MSE of several types of estimators are compared in two different mixed models; under a different types of distributions a simulation study is carried out for the behaviour of derived estimators.     

Using maximum likelihood to estimate population size from temporal changes in allele frequencies.

We develop a maximum-likelihood framework for using temporal changes in allele frequencies to estimate the number of breeding individuals in a population. We use simulations to compare the performance of this estimator to an F-statistic estimator of variance effective population size. The maximum-likelihood estimator had a lower variance and smaller bias. Taking advantage of the likelihood framework, we extend the model to include exponential growth and show that temporal allele frequency data from three or more sampling events can be used to test for population growth.