Constrained estimation using judgment post-stratification

In ranked-set sampling (RSS) and judgment post-stratification (JPS), more efficient inference is obtained by creating a stratification based on ranking information. Using this stratification exactly as is done in stratified sampling or standard post-stratification leads to the standard nonparametric estimators for RSS and JPS. However, we show that strata obtained from ranking information satisfy additional constraints that need not be met by ordinary strata. Specifically, the in-stratum cumulative distribution functions (CDFs) can be no more extreme, in a certain sense, than the CDFs for order statistics from the overall distribution. The additional constraints can be used to obtain better small-sample estimates of the in-stratum CDFs using either RSS or JPS. In the JPS case, the constraints also lead to better small-sample estimates of the overall CDF and the population mean.     

Semi-parametric estimation of covariate effects using the positive stable frailty model

Many biological and medical studies have as a response of interest the time to occurrence of some event, such as the occurrence of a particular symptom or disease, remission, relapse, death due to some specific disease, or simply death. In this paper we study the problem of assessing the effect of potential risk factors on the outcome event of interest through a parametric or semi-parametric frailty model where the lifetimes have a reason to be considered dependent. This dependence may arise because of multiple endpoints within the same individual or because, when studying a single endpoint, there are natural groupings between study subjects. The objective of this paper is to extend both parametric and semi-parametric approaches to regression analysis in which the lifetimes of individuals in a group are effected by the same random frailty which follows a positive stable distribution. Some comparisons of the properties of this frailty distribution with other frailty distributions are made and an example which assesses the effect of a treatment in a litter-matched tumorigenesis study is presented.     

A predictive estimator of finite population mean using nonparametric regression

The paper considers the problem of estimating the population mean using auxiliary information. We propose a new model-based estimator of the population mean, based on local polynomial regression. This estimator exhibits several attractive properties under the model-based approach. The estimator is compared to a number of methods which have been proposed in the literature via a simulation study based on several populations.     

Multi-stage nonparametric estimation of density function using orthonormal systems

In this note, the authors propose a new nonparametric method of estimation of density using orthonormal systems iteratively. The asymptotic mean integrated square error of the estimate at each stage is less than or equal to that of the preceding stage. The new estimate is better, in some cases, than the traditional estimate based upon orthonormal functions from the point of view of the mean integrated square error in the limit.     

Bivariate nonparametric estimation of the Pickands dependence function using Bernstein copula with kernel regression approach

In this study, a new nonparametric approach using Bernstein copula approximation is proposed to estimate Pickands dependence function. New data points obtained with Bernstein copula approximation serve to estimate the unknown Pickands dependence function. Kernel regression method is then used to derive an intrinsic estimator satisfying the convexity. Some extreme-value copula models are used to measure the performance of the estimator by a comprehensive simulation study. Also, a real-data example is illustrated. The proposed Pickands estimator provides a flexible way to have a better fit and has a better performance than the conventional estimators.     

A note on nonparametric density estimation for dependent variables using a delta sequence

Summary A general method based on “delta sequences” due to Walter and Blum [12] is extended to sequences of strictly stationary mixing random variables having the same marginal distribution admitting a Lebesgue probability density function. It is proved that, under certain conditions, the rate of mean square convergence obtained in the i.i.d. case by Walter and Blum, continues to hold. University of Petroleum and Minerals     

Adaptive nonparametric estimation of a multivariate regression function

We consider the kernel estimation of a multivariate regression function at a point. Theoretical choices of the bandwidth are possible for attaining minimum mean squared error or for local scaling, in the sense of asymptotic distribution. However, these choices are not available in practice. We follow the approach of Krieger and Pickands (Ann. Statist.9 (1981) 1066–1078) and Abramson (J. Multivariate Anal.12 (1982), 562–567) in constructing adaptive estimates after demonstrating the weak convergence of some error process. As consequences, efficient data-driven consistent estimation is feasible, and data-driven local scaling is also feasible. In the latter instance, nearest-neighbor-type estimates and variance-stabilizing estimates are obtained as special cases.     

A nonparametric estimation of the infection curve

Predicting the future course of an epidemic depends on being able to estimate the current numbers of infected individuals.However,while back-projection techniques allow reliable estimation of the numbers of infected individuals in the more distant past,they are less reliable in the recent past.We propose two new nonparametric methods to estimate the unobserved numbers of infected individuals in the recent past in an epidemic.The proposed methods are noniterative,easily computed and asymptotically normal wit...     

On the estimation of a restricted normal mean

Let X ∼ N(θ,1), where θ ϵ [−m, m], for some m > 0, and consider the problem of estimating θ with quadratic loss. We show that the Bayes estimator δ_m, corresponding to the uniform prior on [−m, m], dominates δ⁰ (x) = x on [−m, m] and it also dominates the MLE over a large part of the parameter interval. We further offer numerical evidence to suggest that δ_m has quite satisfactory risk performance when compared with the minimax estimators proposed by Casella and Strawderman (1981) and the estimators proposed by Bickel (1981).     

On the estimation of a monotone conditional variance in nonparametric regression

A monotone estimate of the conditional variance function in a heteroscedastic, nonparametric regression model is proposed. The method is based on the application of a kernel density estimate to an unconstrained estimate of the variance function and yields an estimate of the inverse variance function. The final monotone estimate of the variance function is obtained by an inversion of this function. The method is applicable to a broad class of nonparametric estimates of the conditional variance and particularly attractive to users of conventional kernel methods, because it does not require constrained optimization techniques. The approach is also illustrated by means of a simulation study.     

Nonlinear weighted least squares estimation of a three-parameter Weibull density with a nonparametric start

This paper is concerned with the parameter estimation problem for the three-parameter Weibull density which is widely employed as a model in reliability and lifetime studies. Our approach is a combination of nonparametric and parametric methods. The basic idea is to start with an initial nonparametric density estimate which needs to be as good as possible, and then apply the nonlinear least squares method to estimate the unknown parameters. As a main result, a theorem on the existence of the least squares estimate is obtained. Some simulations are given to show that our approach is satisfactory if the initial density is of good enough quality.     

Optimal nonparametric estimation of the density of regression errors with finite support

Knowledge of the probability distribution of error in a regression problem plays an important role in verification of an assumed regression model, making inference about predictions, finding optimal regression estimates, suggesting confidence bands and goodness of fit tests as well as in many other issues of the regression analysis. This article is devoted to an optimal estimation of the error probability density in a general heteroscedastic regression model with possibly dependent predictors and regression errors. Neither the design density nor regression function nor scale function is assumed to be known, but they are suppose to be differentiable and an estimated error density is suppose to have a finite support and to be at least twice differentiable. Under this assumption the article proves, for the first time in the literature, that it is possible to estimate the regression error density with the accuracy of an oracle that knows “true” underlying regression errors. Real and simulated examples illustrate importance of the error density estimation as well as the suggested oracle methodology and the method of estimation.     

Quantifying the risk using copulae with nonparametric marginals

We show that copulae and kernel estimation can be mixed to estimate the risk of an economic loss. We analyze the properties of the Sarmanov copula. We find that the maximum pseudo-likelihood estimation of the dependence parameter associated with the copula with double transformed kernel estimation to estimate marginal cumulative distribution functions is a useful method for approximating the risk of extreme dependent losses when we have large data sets. We use a bivariate sample of losses from a real database of auto insurance claims.     

Abrupt change in mean using block bootstrap and avoiding variance estimation

We deal with sequences of weakly dependent observations that are naturally ordered in time. Their constant mean is possibly subject to change at most once at some unknown time point. The aim is to test whether such an unknown change has occurred or not. The change point methods presented here rely on ratio type test statistics based on maxima of the cumulative sums. These detection procedures for the abrupt change in mean are also robustified by considering a general score function. The main advantage of the proposed approach is that the variance of the observations neither has to be known nor estimated. The asymptotic distribution of the test statistic under the no change null hypothesis is derived. Moreover, we prove the consistency of the test under the alternatives. A block bootstrap method is developed in order to obtain better approximations for the test’s critical values. The validity of the bootstrap algorithm is shown. The results are illustrated through a simulation study, which demonstrates computational efficiency of the procedures. A practical application to real data is presented as well.     

Ratio methods to the mean estimation with known quantiles

We address the problem of estimating the finite population mean in survey sampling, by exploiting any available auxiliary information in order to increase the precision of classical estimators. The idea is to use any population quantiles of the available auxiliary variables which are known in many real situation from census, administrative files, etc. This is achieved using these known quantities in the construction of the estimators, by modifying the usual ratio estimation methods and afterwards defining a general class of exponentiation ratio estimators. The advantages of the proposed estimators are demonstrated using theoretical asymptotic tools and through a simulation study.     

Bayes empirical Bayes estimation of a poisson mean

In the empirical Bayes (EB) decision problem consisting of squared error estimation of a Poisson mean, a prior distribution λ is placed on the gamma family of prior distributions to produce Bayes EB estimators which are admissible. A subclass of such estimators is shown to be asymptotically optimal (a.o.). The results of a Monte Carlo study are presented to demonstrate the favorable a.o. property of the Bayes EB estimators in comparison with other competitors.     

One-armed bandit process with a covariate

We generalize the bandit process with a covariate introduced by Woodroofe in several significant directions: a linear regression model characterizing the unknown arm, an unknown variance for regression residuals and general discounting sequence for a non-stationary model. With the Bayesian regression approach, we assume a normal-gamma conjugate prior distribution of the unknown parameters. It is shown that the optimal strategy is determined by a sequence of index values which are monotonic and determined by the observed value of the covariate and updated posterior distributions. We further show that the myopic strategy is not optimal in general. Such structural properties help to understand the tradeoff between information gathering and immediate expected payoff and may provide certain insight for covariate adjusted response adaptive design of clinical trials.