Estimation and comparison of changes in the presence of informative right censoring: conditional linear model

A general linear regression model for the usual least squares estimated rate of change (slope) on censoring time is described as an approximation to account for informative right censoring in estimating and comparing changes of a continuous variable in two groups. Two noniterative estimators for the group slope means, the linear minimum variance unbiased (LMVUB) estimator and the linear minimum mean squared error (LMMSE) estimator, are proposed under this conditional model. In realistic situations, we illustrate that the LMVUB and LMMSE estimators, derived under a simple linear regression model, are quite competitive compared to the pseudo maximum likelihood estimator (PMLE) derived by modeling the censoring probabilities. Generalizations to polynomial response curves and general linear models are also described.     

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.     

A new approach to weighting and inference in sample surveys

The validity of design-based inference is not dependent on anymodel assumption. However, it is well known that estimatorsderived through design-based theory may be inefficient for theestimation of population totals when the design weights areweakly related to the variables of interest and have widelydispersed values. We propose estimators that have the potentialto improve the efficiency of any estimator derived under thedesign-based theory. Our main focus is limited to the improvementof the Horvitz–Thompson estimator, but we also discussthe extension to calibration estimators. The new estimatorsare obtained by smoothing design or calibration weights usingan appropriate model. Our approach to inference requires themodelling of only one variable, the weight, and it leads toa single set of smoothed weights in multipurpose surveys. Thisis to be contrasted with other model-based approaches, suchas the prediction approach, in which it is necessary to postulateand validate a model for each variable of interest leading potentiallyto variable-specific sets of weights. Our proposed approachis first justified theoretically and then evaluated througha simulation study.     

On the potential for illogic with logically defined outcomes

Logically defined outcomes are commonly used in medical diagnoses and epidemiological research. When missing values in the original outcomes exist, the method of handling the missingness can have unintended consequences, even if the original outcomes are missing completely at random. In this note, we consider 2 binary original outcomes, which are missing completely at random. For estimating the prevalence of a logically defined "or" outcome, we discuss the properties of 4 estimators: the complete-case estimator, the available-case estimator, the maximum likelihood estimator (MLE), and a moment-based estimator. With the exception of the available-case case estimator, all the estimators are consistent. The MLE exhibits superior performance and should be generally adopted.     

Some Improved Regression Estimators of Heritability

Three new improved regression estimators of heritability viz. modified range restricted estimator, minimum quadratic loss estimator and minimax linear restricted estimator are proposed. In addition, these estimators are illustrated and compared numerically with the existing restricted estimator based on linear stochastic constraint.     

On Combination of Ratio and PPS Estimators

AGARWAL and KUMAR (1980) proposed an estimator, combining ratio and pps estimators of population mean and proved that the proposed estimator would always be better (in minimum mean square error sense) than the pps estimator or the ratio estimator under pps sampling scheme for optimum value of constant k (parameter). The optimum value of k is rarely known in practice, hence the alternative is to replace k from the sample-values. In this paper, an estimator depending on estimated optimum value of k based on sample-values, under pps sampling scheme is proposed and studied.     

Design‐based treatment of unit nonresponse in environmental surveys using calibration weighting

Unit nonresponse is often a problem in sample surveys. It arises when the values of the survey variable cannot be recorded for some sampled units. In this paper, the use of nonresponse calibration weighting to treat nonresponse is considered in a complete design‐based framework. Nonresponse is viewed as a fixed characteristic of the units. The approach is suitable in environmental and forest surveys when sampled sites cannot be reached by field crews. Approximate expressions of design‐based bias and variance of the calibration estimator are derived and design‐based consistency is investigated. Choice of auxiliary variables to perform calibration is discussed. Sen–Yates–Grundy, Horvitz–Thompson, and jackknife estimators of the sampling variance are proposed. Analytical and Monte Carlo results demonstrate the validity of the procedure when the relationship between survey and auxiliary variables is similar in respondent and nonrespondent strata. An application to a forest survey performed in Northeastern Italy is considered.     

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.     

An Efficient Regression Type Estimator in Survey Sampling

For the estimation of population mean in simple random sampling, an efficient regression-type estimator is proposed which is more efficient than the conventional regression estimator and hence than mean per unit estimator, ratio and product estimators and many other estimators proposed by various authors. Some numerical examples are included for illustration.     

Analysis of longitudinal marginal structural models

In this article we construct and study estimators of the causal effect of a time-dependent treatment on survival in longitudinal studies. We employ a particular marginal structural model (MSM), proposed by Robins (2000), and follow a general methodology for constructing estimating functions in censored data models. The inverse probability of treatment weighted (IPTW) estimator of Robins et al. (2000) is used as an initial estimator and forms the basis for an improved, one-step estimator that is consistent and asymptotically linear when the treatment mechanism is consistently estimated. We extend these methods to handle informative censoring. The proposed methodology is employed to estimate the causal effect of exercise on mortality in a longitudinal study of seniors in Sonoma County. A simulation study demonstrates the bias of naive estimators in the presence of time-dependent confounders and also shows the efficiency gain of the IPTW estimator, even in the absence such confounding. The efficiency gain of the improved, one-step estimator is demonstrated through simulation.     

Design‐based inference on Bernstein type estimators for continuous populations

The estimation of the values of a variable at any point of a study area is performed using Bernstein polynomials when the sampling scheme is implemented by selecting a point in each polygon of a regular grid overimposed onto the area. The evaluation of the precision of the resulting estimates is investigated under a completely design‐based framework. Moreover, as the main contribution to the mean squared error of the Bernstein‐type estimator is due to the bias, also a pseudo‐jackknife estimator is proposed. The performance of both estimators is investigated theoretically and by means of a simulation study. An application to a soil survey performed in Berkshire Downs in Oxfordshire (UK) is considered.     

Improved Shrinkage Estimation of Relative Potency

This article considers the asymptotic estimation theory for the log relative potency in a symmetric parallel bioassay when uncertain prior information about the true log relative potency is assumed to be a known quantity. Three classes of point estimation, namely, the unrestricted estimator, the shrinkage restricted estimator and shrinkage preliminary test estimator are proposed. Their asymptotic mean squared errors are derived and compared. The relative dominance picture of the estimators is presented. Interestingly, proposed shrinkage preliminary test estimator dominates the unrestricted estimator in a range that is wider than that of the usual preliminary test estimator. Most importantly, the size of the preliminary test is much appropriate than the usual preliminary test estimator.     

On an Efficient Regression Type Estimator

The regression type estimator proposed by KAUR (1985) is considered. Another expression for the approximated mean square error (AMSE), to a first degree of approximation, is obtained. This AMSE is also minimized with respect to a parameter α. Three numerical examples are included. These numerical examples show that this estimator is not significantly more efficient than regression estimator and with respect to ratio and sample mean estimators, it does not always exhibit a high efficiency, as was contended by KAUR (1985). Moreover, an upper bound for the relative precision of the proposed estimator with respect to linear regression estimator is derived.     

An estimator of number of species from quadrat sampling

We consider the problem of estimating the number of distinct species S in a study area from the recorded presence or absence of species in each of a sample of quadrats. A generalized jackknife estimator of S is derived, along with an estimate of its variance. It is compared with the jackknife estimator for S proposed by Heltshe and Forrester and the empirical Bayes estimator of Mingoti and Meeden. We show that the empirical Bayes estimator has the form of a generalized jackknife estimator under a specific model for species distribution. We compare the new estimators of S to the empirical Bayes estimator via simulation. We characterize circumstances under which each is superior.     

Time-dependent ROC curves for censored survival data and a diagnostic marker

ROC curves are a popular method for displaying sensitivity and specificity of a continuous diagnostic marker, X, for a binary disease variable, D. However, many disease outcomes are time dependent, D(t), and ROC curves that vary as a function of time may be more appropriate. A common example of a time-dependent variable is vital status, where D(t) = 1 if a patient has died prior to time t and zero otherwise. We propose summarizing the discrimination potential of a marker X, measured at baseline (t = 0), by calculating ROC curves for cumulative disease or death incidence by time t, which we denote as ROC(t). A typical complexity with survival data is that observations may be censored. Two ROC curve estimators are proposed that can accommodate censored data. A simple estimator is based on using the Kaplan-Meier estimator for each possible subset X > c. However, this estimator does not guarantee the necessary condition that sensitivity and specificity are monotone in X. An alternative estimator that does guarantee monotonicity is based on a nearest neighbor estimator for the bivariate distribution function of (X, T), where T represents survival time (Akritas, M. J., 1994, Annals of Statistics 22, 1299-1327). We present an example where ROC(t) is used to compare a standard and a modified flow cytometry measurement for predicting survival after detection of breast cancer and an example where the ROC(t) curve displays the impact of modifying eligibility criteria for sample size and power in HIV prevention trials.     

On a Class of Almost Unbiased Ratio Estimators

A class of almost unbiased ratio estimators for population mean σ is derived by weighting sample σ = (1/n) σ y_i, ratio estimators σ and an estimator, σ (y_i/x_i). It is shown that NIETO DE PASCUAL (1961) estimator is a particular member of the class and an optimum estimator in the class (in the minimum variance sense) is identified. The results are illustrated through two numerical examples.     

利用多重抽样估计稀有物种的种类数

利用矩估计和二个稳健估计方法（jackknife估计,bootstrap估计）来处理野外生态学工作者的调查数据,在假定已经发现一些稀有物种的情形下,通过统计推断得到那些未被发现的物种的种类数。利用本文所提出的方法调查水稻水稻田的昆虫群落和林地的在面植被群落的稀有种是十分有效的。     

Robust best linear estimation for regression analysis using surrogate and instrumental variables

We investigate methods for regression analysis when covariates are measured with errors. In a subset of the whole cohort, a surrogate variable is available for the true unobserved exposure variable. The surrogate variable satisfies the classical measurement error model, but it may not have repeated measurements. In addition to the surrogate variables that are available among the subjects in the calibration sample, we assume that there is an instrumental variable (IV) that is available for all study subjects. An IV is correlated with the unobserved true exposure variable and hence can be useful in the estimation of the regression coefficients. We propose a robust best linear estimator that uses all the available data, which is the most efficient among a class of consistent estimators. The proposed estimator is shown to be consistent and asymptotically normal under very weak distributional assumptions. For Poisson or linear regression, the proposed estimator is consistent even if the measurement error from the surrogate or IV is heteroscedastic. Finite-sample performance of the proposed estimator is examined and compared with other estimators via intensive simulation studies. The proposed method and other methods are applied to a bladder cancer case-control study.     

On the inadmissibility of Watterson's estimator

We consider the estimation of the scaled mutation parameter θ, which is one of the parameters of key interest in population genetics. We provide a general result showing when estimators of θ can be improved using shrinkage when taking the mean squared error as the measure of performance. As a consequence, we show that Watterson’s estimator is inadmissible, and propose an alternative shrinkage-based estimator that is easy to calculate and has a smaller mean squared error than Watterson’s estimator for all possible parameter values 0<θ<∞. This estimator is admissible in the class of all linear estimators. We then derive improved versions for other estimators of θ, including the MLE. We also investigate how an improvement can be obtained both when combining information from several independent loci and when explicitly taking into account recombination. A simulation study provides information about the amount of improvement achieved by our alternative estimators.     

Optimizing performance of nonparametric species richness estimators under constrained sampling

Understanding the functional relationship between the sample size and the performance of species richness estimators is necessary to optimize limited sampling resources against estimation error. Nonparametric estimators such as Chao and Jackknife demonstrate strong performances, but consensus is lacking as to which estimator performs better under constrained sampling. We explore a method to improve the estimators under such scenario. The method we propose involves randomly splitting species‐abundance data from a single sample into two equally sized samples, and using an appropriate incidence‐based estimator to estimate richness. To test this method, we assume a lognormal species‐abundance distribution (SAD) with varying coefficients of variation (CV), generate samples using MCMC simulations, and use the expected mean‐squared error as the performance criterion of the estimators. We test this method for Chao, Jackknife, ICE, and ACE estimators. Between abundance‐based estimators with the single sample, and incidence‐based estimators with the split‐in‐two samples, Chao2 performed the best when CV < 0.65, and incidence‐based Jackknife performed the best when CV > 0.65, given that the ratio of sample size to observed species richness is greater than a critical value given by a power function of CV with respect to abundance of the sampled population. The proposed method increases the performance of the estimators substantially and is more effective when more rare species are in an assemblage. We also show that the splitting method works qualitatively similarly well when the SADs are log series, geometric series, and negative binomial. We demonstrate an application of the proposed method by estimating richness of zooplankton communities in samples of ballast water. The proposed splitting method is an alternative to sampling a large number of individuals to increase the accuracy of richness estimations; therefore, it is appropriate for a wide range of resource‐limited sampling scenarios in ecology.