Estimation of the Density of G Given Observations of Y=G+E

We study the problem of estimating the density of a random variable G, given observations of a random variable Y = G + E. The random variable E is independent of G and its probability distribution function is considered as known. We build a family of estimators of the density of G using characteristic functions. We then derive a family of estimators of the density of Y based on the model for Y. The estimators are shown to be asymptotically unbiased and consistent. Simulations show that these estimators are better, as measured by integrated squared error, than the standard kernel estimators. Finally, we give an example of the use of this method for the detection of major genes in animal populations.     

Outlier robust model‐assisted small area estimation

Small area estimation with M‐quantile models was proposed by Chambers and Tzavidis ( 2006 ). The key target of this approach to small area estimation is to obtain reliable and outlier robust estimates avoiding at the same time the need for strong parametric assumptions. This approach, however, does not allow for the use of unit level survey weights, making questionable the design consistency of the estimators unless the sampling design is self‐weighting within small areas. In this paper, we adopt a model‐assisted approach and construct design consistent small area estimators that are based on the M‐quantile small area model. Analytic and bootstrap estimators of the design‐based variance are discussed. The proposed estimators are empirically evaluated in the presence of complex sampling designs.     

Heteroscedastic Nonlinear Regression Models with Random Effects and Their Application to Enzyme Kinetic Data

We present a new modification of nonlinear regression models for repeated measures data with heteroscedastic error structures by combining the transform-both-sides and weighting model from Caroll and Ruppert (1988) with the nonlinear random effects model from Lindstrom and Bates (1990). The proposed parameter estimators are a combination of pseudo maximum likelihood estimators for the transform-both-sides and weighting model and maximum likelihood (ML) or restricted maximum likelihood (REML) estimators for linear mixed effects models. The new method is investigated by analyzing simulated enzyme kinetic data published by Jones (1993).     

Estimating the kernel parameters of premises-based stochastic models of farmed animal infectious disease epidemics using limited, incomplete, or ongoing data

Three different estimators are presented for the types of parameters present in mathematical models of animal epidemics. The estimators make use of the data collected during an epidemic, which may be limited, incomplete, or under collection on an ongoing basis. When data are being collected on an ongoing basis, the estimated parameters can be used to evaluate putative control strategies. These estimators were tested using simulated epidemics based on a spatial, discrete-time, gravity-type, stochastic mathematical model containing two parameters. Target epidemics were simulated with the model and the three estimators were implemented using various combinations of collected data to independently determine the two parameters.     

Weighted rank regression for clustered data analysis

Summary .   We consider ranked-based regression models for clustered data analysis. A weighted Wilcoxon rank method is proposed to take account of within-cluster correlations and varying cluster sizes. The asymptotic normality of the resulting estimators is established. A method to estimate covariance of the estimators is also given, which can bypass estimation of the density function. Simulation studies are carried out to compare different estimators for a number of scenarios on the correlation structure, presence/absence of outliers and different correlation values. The proposed methods appear to perform well, in particular, the one incorporating the correlation in the weighting achieves the highest efficiency and robustness against misspecification of correlation structure and outliers. A real example is provided for illustration.     

Nonparametric inverse-probability-weighted estimators based on the highly adaptive lasso

Inverse-probability-weighted estimators are the oldest and potentially most commonly used class of procedures for the estimation of causal effects. By adjusting for selection biases via a weighting mechanism, these procedures estimate an effect of interest by constructing a pseudopopulation in which selection biases are eliminated. Despite their ease of use, these estimators require the correct specification of a model for the weighting mechanism, are known to be inefficient, and suffer from the curse of dimensionality. We propose a class of nonparametric inverse-probability-weighted estimators in which the weighting mechanism is estimated via undersmoothing of the highly adaptive lasso, a nonparametric regression function proven to converge at nearly $n^{-1/3}$-rate to the true weighting mechanism. We demonstrate that our estimators are asymptotically linear with variance converging to the nonparametric efficiency bound. Unlike doubly robust estimators, our procedures require neither derivation of the efficient influence function nor specification of the conditional outcome model. Our theoretical developments have broad implications for the construction of efficient inverse-probability-weighted estimators in large statistical models and a variety of problem settings. We assess the practical performance of our estimators in simulation studies and demonstrate use of our proposed methodology with data from a large-scale epidemiologic study.     

Statistical Analysis of Nonrectangular Family Data

MINQUE (Minimum Norm Quadratic Unbiased Estimators) theory is applied to the problem of estimation of variance components in family data (siblings) with variable family size. Using this approach, the traditional iterative maximum likelihood estimators are shown to be asymptotically normal, even though the data come from non-identical parent distributions. Asymptotic expressions are also obtained for the variance of the MINQUE estimators which hold even if the data are decidedly non-normal (e.g. a mixture of normals). In the case of normal data, exact small-sample variance estimates are derived. Simulations demonstrate the fast rate of convergence to asymptotic properties as the number of families increases. These desirable qualities suggest that the easy to compute MINQUE class of estimators may provide a useful alternative method for modelling familial aggregation.     

Statistical Properties and Control Algorithms of Recursive Quantile Estimators

The present study is an extension of the investigations made by Grieszbach and Schack (1993) where the recursive estimators of the quantile were introduced. Attention is focused on statistical properties and on the controlling of these estimators in order to reduce their variance and to improve their capability of adaptation. Using methods of stochastic approximation, several control algorithms have been developed, where both the consistent and the adaptive estimation are considered. Due to the recursive computation formula the estimators are suitable for the analysis of large data sets and for sets whose elements are obtained sequentially. In this study, application examples from the analysis of EEG‐records are presented, where quantiles are used as threshold values.     

Mixed Moment Estimation in Hermite Distribution and its Comparison with Other Estimation Procedures from Sample Size Considerations

In this paper (1) expressions (correct to n^?2 terms) for biases, variances, and covariances of the estimators a and b of Hermite distribution with probability generating function Exp[a(t–1) + b(t–1)] are obtained for two mixed moment estimates; (2) for the biases and variance-covariances, approximate regions of the parameter space (a>0, b>0) have been outlined where a sample of size 100 can be considered as “safe” in the sense that contribution of second order terms in them is 5% of that from the first order term; (3) comparison of the biases and variance-covariances of these two sets of estimators are made with those for the moment estimators, maximum likelihood estimates and the even point estimators for a sample of size 100 using the terms up to order n^?2; (4) the comparisons based on n^?2 terms in (3) have not only provided information on the estimation procedures included in the Hermite distribution, but also demonstrated the importance of higher order terms in the sampling properties of the various alternative techniques for the Hermite distribution.     

Influence functions and robust Bayes and empirical Bayes small area estimation

We introduce new robust small area estimation procedures basedon area-level models. We first find influence functions correspondingto each individual area-level observation by measuring the divergencebetween the posterior density functions of regression coefficientswith and without that observation. Next, based on these influencefunctions, properly standardized, we propose some new robustBayes and empirical Bayes small area estimators. The mean squarederrors and estimated mean squared errors of these estimatorsare also found. A small simulation study compares the performanceof the robust and the regular empirical Bayes estimators. Whenthe model variance is larger than the sample variance, the proposedrobust empirical Bayes estimators are superior.     

Nonparametric estimation for the three-stage irreversible illness-death model

In this paper, we present new nonparametric estimators of the stage-occupation probabilities in the three-stage irreversible illness-death model. These estimators use a fractional risk set and a reweighting approach and are valid under stage-dependent censoring. Using a simulated data set, we compare the behavior of our estimators with previously proposed estimators. We also apply our estimators to data on time to Pneumocystis pneumonia and death obtained from an AIDS cohort study.     

Nonparametric estimation of a Markov 'illness-death' process from interval-censored observations, with application to diabetes survival data

The nonparametric estimation of the cumulative transition intensityfunctions in a threestate time-nonhomogeneous Markov processwith irreversible transitions, an ‘illness-death’model, is considered when times of the intermediate transition,e.g. onset of a disease, are interval-censored. The times of‘death’ are assumed to be known exactly or to beright-censored. In addition the observed process may be left-truncated.Data of this type arise when the process is sampled periodically.For example, when the patients are monitored through periodicexaminations the observations on times of change in their diseasestatus will be interval-censored. Under the sampling schemeconsidered here the Nelson–Aalen estimator (Aalen, 1978)for a cumulative transition intensity is not applicable. Inthe proposed method the maximum likelihood estimators of someof the transition intensities are derived from the estimatorsof the corresponding subdistribution functions. The maximumlikelihood estimators are shown to have a self-consistency property.The self-consistency algorithm is developed for the computationof the estimators. This approach generalises the results fromTurnbull (1976) and Frydman (1992). The methods are illustratedwith diabetes survival data.     

Evaluating the performance of species richness estimators: sensitivity to sample grain size

1. Fifteen species richness estimators (three asymptotic based on species accumulation curves, 11 nonparametric, and one based in the species-area relationship) were compared by examining their performance in estimating the total species richness of epigean arthropods in the Azorean Laurisilva forests. Data obtained with standardized sampling of 78 transects in natural forest remnants of five islands were aggregated in seven different grains (i.e. ways of defining a single sample): islands, natural areas, transects, pairs of traps, traps, database records and individuals to assess the effect of using different sampling units on species richness estimations. 2. Estimated species richness scores depended both on the estimator considered and on the grain size used to aggregate data. However, several estimators (ACE, Chao 1, Jackknifel and 2 and Bootstrap) were precise in spite of grain variations. Weibull and several recent estimators [proposed by Rosenzweig et al. (Conservation Biology, 2003, 17, 864-874), and Ugland et al. (Journal of Animal Ecology, 2003, 72, 888-897)] performed poorly. 3. Estimations developed using the smaller grain sizes (pair of traps, traps, records and individuals) presented similar scores in a number of estimators (the above-mentioned plus ICE, Chao2, Michaelis-Menten, Negative Exponential and Clench). The estimations from those four sample sizes were also highly correlated. 4. Contrary to other studies, we conclude that most species richness estimators may be useful in biodiversity studies. Owing to their inherent formulas, several nonparametric and asymptotic estimators present insensitivity to differences in the way the samples are aggregated. Thus, they could be used to compare species richness scores obtained from different sampling strategies. Our results also point out that species richness estimations coming from small grain sizes can be directly compared and other estimators could give more precise results in those cases. We propose a decision framework based on our results and on the literature to assess which estimator should be used to compare species richness scores of different sites, depending on the grain size of the original data, and of the kind of data available (species occurrence or abundance data).     

Upper bound estimators of the population size based on ordinal models for capture‐recapture experiments

Capture‐recapture studies have attracted a lot of attention over the past few decades, especially in applied disciplines where a direct estimate for the size of a population of interest is not available. Epidemiology, ecology, public health, and biodiversity are just a few examples. The estimation of the number of unseen units has been a challenge for theoretical statisticians, and considerable progress has been made in providing lower bound estimators for the population size. In fact, it is well known that consistent estimators for this cannot be provided in the very general case. Considering a case where capture‐recapture studies are summarized by a frequency of frequencies distribution, we derive a simple upper bound of the population size based on the cumulative distribution function. We introduce two estimators of this bound, without any specific parametric assumption on the distribution of the observed frequency counts. The behavior of the proposed estimators is investigated using several benchmark datasets and a large‐scale simulation experiment based on the scheme discussed by Pledger.     

Robust inference on the average treatment effect using the outcome highly adaptive lasso

Many estimators of the average effect of a treatment on an outcome require estimation of the propensity score, the outcome regression, or both. It is often beneficial to utilize flexible techniques, such as semiparametric regression or machine learning, to estimate these quantities. However, optimal estimation of these regressions does not necessarily lead to optimal estimation of the average treatment effect, particularly in settings with strong instrumental variables. A recent proposal addressed these issues via the outcome-adaptive lasso, a penalized regression technique for estimating the propensity score that seeks to minimize the impact of instrumental variables on treatment effect estimators. However, a notable limitation of this approach is that its application is restricted to parametric models. We propose a more flexible alternative that we call the outcome highly adaptive lasso. We discuss the large sample theory for this estimator and propose closed-form confidence intervals based on the proposed estimator. We show via simulation that our method offers benefits over several popular approaches.     

Estimation of multiple biomass growth rates and biomass concentration in a class of bioprocesses

In this paper, an approach to the estimation of multiple biomass growth rates and biomass concentration is proposed for a class of aerobic bioprocesses characterized by on-line measurements of dissolved oxygen and carbon dioxide concentrations, as well as off-line measurements of biomass concentration. The approach is based on adaptive observer theory and includes two steps. In the first step, an adaptive estimator of two out of three biomass growth rates is designed. In the second step, the third biomass growth rate and the biomass concentration are estimated, using two different adaptive estimators. One of them is based on on-line measurements of dissolved oxygen concentration and off-line measurement of biomass concentrations, while the other needs only on-line measurements of the carbon dioxide concentration. Simulations demonstrated good performance of the proposed estimators under continuous and batch-fed conditions.     

Dose finding using the biased coin up-and-down design and isotonic regression

We are interested in finding a dose that has a prespecified toxicity rate in the target population. In this article, we investigate five estimators of the target dose to be used with the up-and-down biased coin design (BCD) introduced by Durham and Flournoy (1994, Statistical Decision Theory and Related Topics). These estimators are derived using maximum likelihood, weighted least squares, sample averages, and isotonic regression. A linearly interpolated isotonic regression estimate is shown to be simple to derive and to perform as well as or better than the other target dose estimators in terms of mean square error and average number of subjects needed for convergence in most scenarios studied.     

Adjustment for Missingness Using Auxiliary Information in Semiparametric Regression

Summary .  In this article, we study the estimation of mean response and regression coefficient in semiparametric regression problems when response variable is subject to nonrandom missingness. When the missingness is independent of the response conditional on high-dimensional auxiliary information, the parametric approach may misspecify the relationship between covariates and response while the nonparametric approach is infeasible because of the curse of dimensionality. To overcome this, we study a model-based approach to condense the auxiliary information and estimate the parameters of interest nonparametrically on the condensed covariate space. Our estimators possess the double robustness property, i.e., they are consistent whenever the model for the response given auxiliary covariates or the model for the missingness given auxiliary covariate is correct. We conduct a number of simulations to compare the numerical performance between our estimators and other existing estimators in the current missing data literature, including the propensity score approach and the inverse probability weighted estimating equation. A set of real data is used to illustrate our approach.     

A general approach for two-stage analysis of multilevel clustered non-Gaussian data

In this article, we propose a two-stage approach to modeling multilevel clustered non-Gaussian data with sufficiently large numbers of continuous measures per cluster. Such data are common in biological and medical studies utilizing monitoring or image-processing equipment. We consider a general class of hierarchical models that generalizes the model in the global two-stage (GTS) method for nonlinear mixed effects models by using any square-root-n-consistent and asymptotically normal estimators from stage 1 as pseudodata in the stage 2 model, and by extending the stage 2 model to accommodate random effects from multiple levels of clustering. The second-stage model is a standard linear mixed effects model with normal random effects, but the cluster-specific distributions, conditional on random effects, can be non-Gaussian. This methodology provides a flexible framework for modeling not only a location parameter but also other characteristics of conditional distributions that may be of specific interest. For estimation of the population parameters, we propose a conditional restricted maximum likelihood (CREML) approach and establish the asymptotic properties of the CREML estimators. The proposed general approach is illustrated using quartiles as cluster-specific parameters estimated in the first stage, and applied to the data example from a collagen fibril development study. We demonstrate using simulations that in samples with small numbers of independent clusters, the CREML estimators may perform better than conditional maximum likelihood estimators, which are a direct extension of the estimators from the GTS method.     

Confidence Intervals of the Attributable Risk under Cross‐Sectional Sampling with Confounders

Since it can account for both the strength of the association between exposure to a risk factor and the underlying disease of interest and the prevalence of the risk factor, the attributable risk (AR) is probably the most commonly used epidemiologic measure for public health administrators to locate important risk factors. This paper discusses interval estimation of the AR in the presence of confounders under cross‐sectional sampling. This paper considers four asymptotic interval estimators which are direct generalizations of those originally proposed for the case of no confounders, and employs Monte Carlo simulation to evaluate the finite‐sample performance of these estimators in a variety of situations. This paper finds that interval estimators using Wald's test statistic and a quadratic equation suggested here can consistently perform reasonably well with respect to the coverage probability in all the situations considered here. This paper notes that the interval estimator using the logarithmic transformation, that is previously found to consistently perform well for the case of no confounders, may have the coverage probability less than the desired confidence level when the underlying common prevalence rate ratio (RR) across strata between the exposure and the non‐exposure is large (≥4). This paper further notes that the interval estimator using the logit transformation is inappropriate for use when the underlying common RR ≐ 1. On the other hand, when the underlying common RR is large (≥4), this interval estimator is probably preferable to all the other three estimators. When the sample size is large (≥400) and the RR ≥ 2 in the situations considered here, this paper finds that all the four interval estimators developed here are essentially equivalent with respect to both the coverage probability and the average length.