A BETA-BINOMIAL MODEL FOR ESTIMATING THE SIZE OF A HETEROGENEOUS POPULATION

This paper compares the properties of various estimators for a beta‐binomial model for estimating the size of a heterogeneous population. It is found that maximum likelihood and conditional maximum likelihood estimators perform well for a large population with a large capture proportion. The jackknife and the sample coverage estimators are biased for low capture probabilities. The performance of the martingale estimator is satisfactory, but it requires full capture histories. The Gibbs sampler and Metropolis‐Hastings algorithm provide reasonable posterior estimates for informative priors.     

MODEL-BASED VARIANCE ESTIMATION IN SURVEYS WITH STRATIFIED CLUSTERED DESIGN

A model-based method for estimating the sampling variances of estimators of (sub-)population means, proportions, quantiles, and regression parameters in surveys with stratified clustered design is described and applied to a survey of US secondary education. The method is compared with the jackknife by a simulation study. The model-based estimators of the sampling variances have much smaller mean squared errors than their jackknife counterparts. In addition, they can be improved by incorporating information about the unknown parameters (variances) from external sources. A regression-based smoothing method for estimating the sampling variances of the estimators for a large number of subpopulation means is proposed. Such smoothing may be invaluable when subpopulations are represented in the sample by only few subjects.     

A comparison of alternative variance estimation strategies for estimating the slope of a linear regression in sample surveys

The balanced half-sample, jackknife and linearization methods are used to estimate the variance of the slope of a linear regression under a variety of computer generated situations. The basic sampling design is one in which two PSU's are selected from each of a number of strata . The variance estimation techniques are compared with a Monte Carlo experiment. Results show that variance estimates may be highly biased and variable unless sizeable numbers of observations are available from each stratum. The jackknife and linearization estimates appear superior to the balanced half sample method - particularly when the number of strata or the number of available observations from each stratum is small.     

Estimation of the parameters of a general student’s t distribution

We have observations for a t distribution with unknown mean, variance, and degrees of freedom, each of which we wish to estimate. The major problem lies in the estimate of the degrees of freedom. We show that a relatively efficient yet very simple estimator is a given function of the ratio of percentile estimates. We derive the appropriate estimator, provide equations for transformation and standard errors, contrast this with other estimators, and give examples.     

SMALL SAMPLE VARIANCE ESTIMATORS FOR U-STATISTICS

New unbiased estimators are presented for the dominant and lower-order terms of the variance expansion for U-statistics. In small samples these provide important corrections to the usual estimate of asymptotic standard error which is based on the leading term in the expansion. The new estimators for the first term cannot be recommended. The ordinary jackknife estimator is found to be more effective than the direct estimates of the separate terms.     

Adjusting the degrees of freedom for the jackknife

The rule of thumb given by Mosteller and Tukey for obtaining the degrees of freedom which are appropriate for the approximate t-statistic associated with the jackknife is discussed.Examples of possible misapplication of this rule are given.     

Robustness of the student t based M-estimator

This paper considers the maximum likelihood type (M) estimator based on Student's t distribution for the location/scale model. The Student t M-estimator is generally thought to be robust to outliers. This paper shows that this is only true if the degrees of freedom parameter is kept fixed. By contrast, if the degrees of freedom parameter is also estimated from the data, the influence functions for the scale and degrees of freedom parameter become unbounded. Moreover, the influence function of the location parameter remains bounded, but its change-of-variance function is unboi~nded. The intuitioil behind these results is explained in the paper. The rates at which both the influence functions and the change-of-variance function diverge to infinity, are very slow. Tliis implies that outliers have to be extremely large in order to become detrimental to the performance of the Student t based M-estimator with estimated degrees of freedom. The theoretical results are illustrated in a a simulation experiment using several related competing estimators and several distributions for the error process.     

Link-tracing sampling with an initial sequential sample of sites: Estimating the size of a hidden human population

We proposed a modification to the variant of link-tracing sampling suggested by Félix-Medina and Thompson [M.H. Félix-Medina, S.K. Thompson, Combining cluster sampling and link-tracing sampling to estimate the size of hidden populations, Journal of Official Statistics 20 (2004) 19–38] that allows the researcher to have certain control of the final sample size, precision of the estimates or other characteristics of the sample that the researcher is interested in controlling. We achieve this goal by selecting an initial sequential sample of sites instead of an initial simple random sample of sites as those authors suggested. We estimate the population size by means of the maximum likelihood estimators suggested by the above-mentioned authors or by the Bayesian estimators proposed by Félix-Medina and Monjardin [M.H. Félix-Medina, P.E. Monjardin, Combining link-tracing sampling and cluster sampling to estimate the size of hidden populations: A Bayesian-assisted approach, Survey Methodology 32 (2006) 187–195]. Variances are estimated by means of jackknife and bootstrap estimators as well as by the delta estimators proposed in the two above-mentioned papers. Interval estimates of the population size are obtained by means of Wald and bootstrap confidence intervals. The results of an exploratory simulation study indicate good performance of the proposed sampling strategy.     

A jackknife variance estimator for unequal probability sampling

Summary.  The jackknife method is often used for variance estimation in sample surveys but has only been developed for a limited class of sampling designs. We propose a jackknife variance estimator which is defined for any without-replacement unequal probability sampling design. We demonstrate design consistency of this estimator for a broad class of point estimators. A Monte Carlo study shows how the proposed estimator may improve on existing estimators.     

A note on biases of jackknife estimators of third central moments

In this paper we study the biases of jackknife estimators of central third moments which play an important role in improving the accuracy of the normal approximation. It has been found in simulation studies that the jackknife estimator of the skewness coefficient, into which the jackknife variance and third moment estimators are substituted, have downward biases. For the jackknife variance estimators, their asymptotic properties are precisely studied and their biases are discussed theoretically, Here we study the biases of the jackknife estimators of the central third moments for U-statistics theoretically, The results show that the biases are not always downward.     

Consistency of jackknife estimators of the variances of simple quantiles

Let σ² be the asymptotic variance of the sample p-quantile (0<p<1). Consistency of the delete-d jackknife estimators of σ² with d being a fraction of n is proved under very weak conditions. Some other results, such as the asymptotic orders of the moments of the jackknife histograms and an analog of the generalized Helly's theorem, are also established.     

Jackknife variances and confidence intervals of the Mantel-Haenszel estimator

In the situation of stratified 2×2 tables, consitency of two different jackknife variances of the Mantel-Haenszel estimator is discussed in the case of increasing sample sizes, but a fixed number of strata. Different principles for constructing confidence limits for the common odds ratio are investigated from a theoretical point of view with regard to the position and the length of the resulting intervals. Monte Carlo experiments compare the finite sample performance of the consistent jackknife variance with that of other noniterative variance estimators. In addition, the properties of these variance estimators are investigated when used for confidence interval estimation.     

Small sample properties of ridge estimators with normal and non-normal disturbances

In this paper we consider five well known and widely used ridge estimators when the convenient assumption of normality of the disturbances is abandoned and report on a Monte Carlo study of their small sample properties. The Monte Carlo experiment is applied to four different data sets with artificially varied degrees of multicollinearity, while the disturbances follow normal, lognormal, uniform and Laplace distributions with small and large variances. The results show that the best estimates are obtained for all ridge estimators when the disturbances follow the lognormal distribution. Also, none of the examined ridge estimators shows a consistent behavior under the different settings considered.     

On the Bayesian Estimation for the Uniform Scale Parameter via a Functional Equation

Bayes uniform model under the squared error loss function is shown to be completely identifiable by the form of the Bayes estimates of the scale parameter. This results in solving a specific functional equation. A complete characterization of differentiable Bayes estimators (BE) and generalized Bayes estimators (GBE) is given as well as relations between degrees of smoothness of the estimators and the priors. Characterizations of strong (generalized Bayes) Bayes sequence (SBS or SGBS) are also investigated. A SBS is a sequence of estimators (one for each sample size) where all its components are BE generated by the same prior measure. A complete solution is given for polynomial Bayesian estimation.     

Jackknife Estimation for Smooth Functions of the Parametric Component in Partially Linear Regression Models

Abstract

It is known that due to the existence of the nonparametric component, the usual estimators for the parametric component or its function in partially linear regression models are biased. Sometimes this bias is severe. To reduce the bias, we propose two jackknife estimators and compare them with the naive estimator. All three estimators are shown to be asymptotically equivalent and asymptotically normally distributed under some regularity conditions. However, through simulation we demonstrate that the jackknife estimators perform better than the naive estimator in terms of bias when the sample size is small to moderate. To make our results more useful, we also construct consistent estimators of the asymptotic variance, which are robust against heterogeneity of the error variances.     

Variance estimation for the finite population distribution function with complete auxiliary information

The authors develop jackknife and analytical variance estimators for the estimator of Chambers & Dunstan (1986) and Rao, Kovar & Mantel (1990) of the finite population distribution function, using complete auxiliary information. They also describe the associated model and show the design consistency of the variance estimators, whose small‐sample performance is examined through a limited simulation study. They highlight the operational advantages of the jackknife in the model‐based setting of Chambers & Dunstan (1986) and its better conditional performance in the design‐based setting of Rao, Kovar & Mantel (1990).     

Estimation of odds ratios under order restrictions

A new method for estimating a set of odds ratios under an order restriction based on estimating equations is proposed. The method is applied to those of the conditional maximum likelihood estimators and the Mantel-Haenszel estimators. The estimators derived from the conditional likelihood estimating equations are shown to maximize the conditional likelihoods. It is also seen that the restricted estimators converge almost surely to the respective odds ratios when the respective sample sizes become large regularly. The restricted estimators are compared with the unrestricted maximum likelihood estimators by a Monte Carlo simulation. The simulation studies show that the restricted estimates improve the mean squared errors remarkably, while the Mantel-Haenszel type estimates are competitive with the conditional maximum likelihood estimates, being slightly worse.     

On the stability of the jackknife variance estimator in ratio estimation

The stability of a slightly modified version of the usual jackknife variance estimator is evaluated exactly in small samples under a suitable linear regression model and compared with that of two different linearization variance estimators. Depending on the degree of heteroscedasticity of the error variance in the model, the stability of the jackknife variance estimator is found to be somewhat comparable to that of one or the other of the linearization variance estimators under conditions especially favorable to ratio estimation (i.e., regression approximately through the origin with a relatively small coefficient of variation in the x population). When these conditions do not hold, however, the jackknife variance estimator is found to be less stable than either of the linearization variance estimators.     

On the Bayesian analysis of the mixture of power function distribution using the complete and the censored sample

The power function distribution is often used to study the electrical component reliability. In this paper, we model a heterogeneous population using the two-component mixture of the power function distribution. A comprehensive simulation scheme including a large number of parameter points is followed to highlight the properties and behavior of the estimates in terms of sample size, censoring rate, parameters size and the proportion of the components of the mixture. The parameters of the power function mixture are estimated and compared using the Bayes estimates. A simulated mixture data with censored observations is generated by probabilistic mixing for the computational purposes. Elegant closed form expressions for the Bayes estimators and their variances are derived for the censored sample as well as for the complete sample. Some interesting comparison and properties of the estimates are observed and presented. The system of three non-linear equations, required to be solved iteratively for the computations of maximum likelihood (ML) estimates, is derived. The complete sample expressions for the ML estimates and for their variances are also given. The components of the information matrix are constructed as well. Uninformative as well as informative priors are assumed for the derivation of the Bayes estimators. A real-life mixture data example has also been discussed. The posterior predictive distribution with the informative Gamma prior is derived, and the equations required to find the lower and upper limits of the predictive intervals are constructed. The Bayes estimates are evaluated under the squared error loss function.     

New modified moment estimators for the two-parameter weighted Lindley distribution

We propose a modification of the moment estimators for the two-parameter weighted Lindley distribution. The modification replaces the second sample moment (or equivalently the sample variance) by a certain sample average which is bounded on the unit interval for all values in the sample space. In this method, the estimates always exist uniquely over the entire parameter space and have consistency and asymptotic normality over the entire parameter space. The bias and mean squared error of the estimators are also examined by means of a Monte Carlo simulation study, and the empirical results show the small-sample superiority in addition to the desirable large sample properties. Monte Carlo simulation study showed that the proposed modified moment estimators have smaller biases and smaller mean-square errors than the existing moment estimators and are compared favourably with the maximum likelihood estimators in terms of bias and mean-square error. Three illustrative examples are finally presented.