Accent on Teaching Materials

Estimation of covariance components in the multivariate random-effect model with nested covariance structure is discussed. There are two covariance matrices to be estimated, namely, the between-group and the within-group covariance matrices. These two covariance matrices are most often estimated by forming a multivariate analysis of variance and equating mean square matrices to their expectations. Such a procedure involves taking the difference between the between-group mean square and the within-group mean square matrices, and often produces an estimated between-group covariance matrix that is not nonnegative definite. We present estimators of the two covariance matrices that are always proper covariance matrices. The estimators are the restricted maximum likelihood estimators if the random effects are normally distributed. The estimation procedure is extended to more complicated models, including the twofold nested and the mixed-effect models. A numerical example is presented to illustrate the use of the estimation procedure.     

C337. Misspecified location and scale. 11. independence of and s2

A hierarchical Bayesian approach to the problem of estimating the largest normal mean is considered. Calculation of the posterior mean and the posterior variance involves, at worst, 3-dimensional numerical integration, for which an efficient Monte Carlo method of evaluation is given. An example is presented to illustrate the methodology. In the two populations case, computation of the posterior estimates can be substantially simplified and in special cases can actually be performed using closed form solutions. A simulation study has been done to compare mean square errors of some hierarchical Bayesian estimators that are expressed in closed forms and several existing estimators of the larger mean.     

On the Estimation of the Probability Density by Trigonometric Series

We present in this article an estimator based on a new orthogonal trigonometric series. We give its statistical properties (bias, variance, mean square error, and mean integrated square error) and the asymptotic properties (convergence of variance, convergence of the mean square error, convergence of the mean integrated square error, uniform convergence in probability, and the rate of convergence of the mean integrated square error). The comparison by simulation on a test density between the estimator obtained from a new trigonometric series with Fejer estimator also based on orthogonal trigonometric series, shows that our estimator is more performant in the sense of the mean integrated square error.     

The three-fold nested random effects model

Estimation of the variance components and the mean of the balanced and unbalanced threefold nested design is considered. The relative merits of the following procedures are evaluated: Analysis of variance (ANOVA), maximum likelihood (ML), restricted maximum likelihood (REML), and minimum variance quadratic unbiased estimator (MIVQUE). A new procedure called the weighted analysis of means (WAM) estimator which utilizes prior information on the variance components is proposed. It is found to have optimum properties similar to the REML and MIVQUE, and it is also computationally simpler. For the mean, the overall sample average, grand mean, unweighted mean, and generalized least-squares (GLS) estimator with its weights obtained from the above estimators for the variance components are considered. Comparisons of the above procedures for the variance components and the mean are made from exact expressions for the biases and mean square errors (MSEs) of the estimators and from empirical investigations.     

PITMAN NEARNESS COMPARISONS OF ESTIMATES OF TWO ORDERED NORMAL MEANS

Maximum likelihood estimates of ordered means of two normal distributions having common variance have been shown to be better than the usual maximum likelihood estimates (i.e. corresponding sample means) with respect to Pitman Nearness criterion. The maximum likelihood estimate of common variance taking into consideration the order restriction of the means is shown to have smaller mean square error than the unrestricted maximum likelihood estimate of the common variance. These two estimators have also been compared with respect to Pitman Nearness criterion.     

Interpenetrating subsampling regression estimation with or without jackknifing

In this paper four regression estimators are considered for a finite population total based on interpenetrating subsamples, two of which are with jackknifing and the other two are without jackknifing. Both theoretical and empirical comparisons of the four proposed estimators are done with respect to bias, variance and mean square error.     

Selection between models through multi-step-ahead forecasting

We develop and show applications of two new test statistics for deciding if one ARIMA model provides significantly better h-step-ahead forecasts than another, as measured by the difference of approximations to their asymptotic mean square forecast errors. The two statistics differ in the variance estimates used for normalization. Both variance estimates are consistent even when the models considered are incorrect. Our main variance estimate is further distinguished by accounting for parameter estimation, while the simpler variance estimate treats parameters as fixed. Their broad consistency properties offer improvements to what are known as tests of Diebold and Mariano (1995) type, which are tests that treat parameters as fixed and use variance estimates that are generally not consistent in our context. We show how these statistics can be calculated for any pair of ARIMA models with the same differencing operator.     

The total bootstrap median: a robust and efficient estimator of location and scale for small samples

We propose the total bootstrap median (TBM) as a robust and efficient estimator of location and scale for small samples. We demonstrate its performance by estimating the mean and variance of a variety of distributions. We also show that, if the underlying distribution is unknown and there is either no contamination or low to moderate contamination, the TBM provides a better estimate of the mean, in mean square terms, than the sample mean or the sample median. In addition, the TBM is a better estimator of the variance of the underlying distribution than the sample variance or the square of the bias-corrected median absolute deviation from the median estimator. We also show that the TBM is an explicit L-estimator, which allows a direct study of its properties.     

Several nonparametric and semiparametric approaches to linear mixed model regression

Mixed models are powerful tools for the analysis of clustered data and many extensions of the classical linear mixed model with normally distributed response have been established. As with all parametric (P) models, correctness of the assumed model is critical for the validity of the ensuing inference. An incorrectly specified P means model may be improved by using a local, or nonparametric (NP), model. Two local models are proposed by a pointwise weighting of the marginal and conditional variance–covariance matrices. However, NP models tend to fit to irregularities in the data and may provide fits with high variance. Model robust regression techniques estimate mean response as a convex combination of a P and a NP model fit to the data. It is a semiparametric method by which incomplete or incorrectly specified P models can be improved by adding an appropriate amount of the NP fit. We compare the approximate integrated mean square error of the P, NP, and mixed model robust methods via a simulation study and apply these methods to two real data sets: the monthly wind speed data from countries in Ireland and the engine speed data.     

Analysis of experiments using the asymptotic quasi-likelihood approach

Comparison of treatment effects in an experiment is usually done through analysis of variance under the assumption that the errors are normally and independently distributed with zero mean and constant variance. The traditional approach in dealing with non-constant variance is to apply a variance stabilizing transformation and then run the analysis on the transformed data. In this approach, the conclusions of analysis of variance apply only to the transformed population. In this paper, the asymptotic quasi-likelihood method is introduced to the analysis of experimental designs. The weak assumptions of the asymptotic quasi-likelihood method make it possible to draw conclusions on heterogeneous populations without transforming them. This paper demonstrates how to apply the asymptotic quasi-likelihood technique to three commonly used models. This gives a possible way to analyse data given a complex experimental design.     

Estimation of reliability for a multicomponent survival stress-strength model based on exponential distributions

Uniformly minimum variance unbiased estimator (UMVUE) of reliability in stress-strength model (known stress) is obtained for a multicomponent survival model based on exponential distributions for parallel system. The variance of this estimator is compared with Cramer-Rao lower bound (CRB) for the variance of unbiased estimator of reliability, and the mean square error (MSE) of maximum likelihood estimator of reliability in case of two component system.     

The correlation structure of longitudinal measurements of vision in patients with macular degeneration

Background: In age‐related macular degeneration (ARMD) trials, the FDA‐approved endpoint is the loss (or gain) of at least three lines of vision as compared to baseline. The use of such a response endpoint entails a potentially severe loss of information. A more efficient strategy could be obtained by using longitudinal measures of the change in visual acuity. In this paper we investigate, by using data from two randomized clinical trials, the mean and variance–covariance structures of the longitudinal measurements of the change in visual acuity. Methods: Individual patient data were collected in 234 patients in a randomized trial comparing interferon‐ α with placebo and in 1181 patients in a randomized trial comparing three active doses of pegaptanib with sham. A linear model for longitudinal data was used to analyze the repeated measurements of the change in visual acuity. Results: For both trials, the data were adequately summarized by a model that assumed a quadratic trend for the mean change in visual acuity over time, a power variance function, and an antedependence correlation structure. The power variance function was remarkably similar for the two datasets and involved the square root of the measurement time. Conclusions: The similarity of the estimated variance functions and correlation structures for both datasets indicates that these aspects may be a genuine feature of the measurements of changes in visual acuity in patients with ARMD. The feature can be used in the planning and analysis of trials that use visual acuity as the clinical endpoint of interest. Copyright © 2010 John Wiley & Sons, Ltd.     

Estimating the variances of ratio estimates in complex sample surveys with two primary sampling units per stratum—a comparison of balanced replication and jackknife techniques

The balanced half-sample and jackknife variance estimation techniques are used to estimate the variance of the combined ratio estimate. An empirical sampling study is conducted using computer-generated populations to investigate the variance, bias and mean square error of these variance estimators and results are compared to theoretical results derived elsewhere for the linear case. Results indicate that either the balanced half-sample or jackknife method may be used effectively for estimating the variance of the combined ratio estimate.     

Maximum likelihood estimation of variance components in a completely random bib design

ABSTRACT

The purposes of this paper are to abstract from a number of articles variance component estimation procedures which can be used for completely random balanced incomplete block designs, to develop an iterated least squares (ITLS) computing algorithm for calculating maximum likelihood estimates, and to compare these procedures by use of simulated experiments. Based on the simulated experiments, the estimated mean square errors of the ITLS estimates are generally less than*those for previously proposed analysis of variance and symmetric sums estimators.     

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.     

Random effects models for count data

Random effects models are considered for count data obtained in a cross or nested classification. The main feature of the proposed models is the use of the additive effects on the original scale in contrast to the commonly used log scale. The rationale behind this approach is given. The estimation of variance components is based on the usual mean square approach. Directly analogous results to those from the analysis of variance models for continuous data are obtained. The usual Poisson dispersion test procedure can be used not only to test for no overall random effects but also to assess the adequacy of the model. Individual variance component can be tested by using the usual F-test. To get a reliable estimate, a large number of factor levels seem to be required.     

Unbiased Confidence Intervals

A method is given for generating different data sets for homework exercises in simple regression, such that not only are the estimated slopes, intercepts and correlations the same over all data sets, but so also are the analysis of variance tables, with error mean square that is a perfect square. An example is given. A computer program REGDATA is available for generating data sets of this nature.     

A new class of unbiased estimators of the variance of the systematic sample mean

We propose a class of estimators of the variance of the systematic sample mean, which is unbiased under the assumption that the population follows a superpopulation model that satisfies some mild conditions. The approach is based on the separate estimation of the portion of the variance due to the systematic component of the model and that due to the stochastic component. In particular, we deal with two estimators belonging to the proposed class that are based on moving averages and local polynomials to estimate the systematic component of the model. The latter estimators are unbiased under the assumption that the population follows a linear trend and the errors are homoscedastic and uncorrelated. Through a simulation study we show that these estimators generally outperform, in terms of bias and mean square error, the usual estimator based on the first differences also when the superpopulation model departs significantly from linearity and the errors are heteroscedastic.     

Variance components of the linear regression model with a random intercept

Several estimators for the variance components of the above model are derived. Biases and mean square errors of the estimators for small samples are examined. Results on the skewness and kurtosis coefficients and the large sample biases and mean square errors of these estimators are presented in detail.     

A new perspective in functional EIV linear model: Part I

Simple linear regression in the functional errors-in-variables (EIV) model is revisited from a different perspective, where the problem is addressed by using the small-sigma model instead of large sample theory. A general analysis is developed to study the slope’s estimator that minimizes a family of objective functions, of which the least-squares fit and the maximum likelihood estimator are minimizers of such special functions. General formulas for the higher-order terms of the bias, the variance, and the mean square error are derived. Accordingly, two efficient estimators are proposed after implementing the pre- and the post-bias elimination techniques. Numerical tests confirm the superiority of the proposed estimators over others.