Rate of convergence for asymptotic variance of the Horvitz–Thompson estimator

Drawing distinct units without replacement and with unequal probabilities from a population is a problem often considered in the literature (e.g. Hanif and Brewer, 1980, Int. Statist. Rev. 48, 317–355). In such a case, the sample mean is a biased estimator of the population mean. For this reason, we use the unbiased Horvitz–Thompson estimator (1951). In this work, we focus our interest on the variance of this estimator. The variance is cumbersome to compute because it requires the calculation of a large number of second-order inclusion probabilities. It would be helpful to use an approximation that does not need heavy calculations. The Hájek (1964) variance approximation provides this advantage as it is free of second-order inclusion probabilities. Hájek (1964) proved that this approximation is valid under restrictive conditions that are usually not fulfilled in practice. In this paper, we give more general conditions and we show that this approximation remains acceptable for most practical problems.     

A Monte Carlo comparison of some variance estimators of the Horvitz–Thompson estimator

At the design and estimation stage of a survey, large survey organization often uses auxiliary information. This article discusses various procedures for improving variance estimation of the Horvitz–Thompson estimator of a finite population total with the aid of auxiliary information. To study the design-based properties of the proposed variance estimators relative to the standard one, a small scale Monte Carlo study is performed.     

On the comparison of Horvitz–Thompson and Murthy's sampling strategies for estimating finite population proportions in direct and randomized response surveys

We consider the problem of unbiased estimation of a finite population proportion and compare the relative efficiency of the unequal probability sampling strategies due to Horvitz and Thompson (1952 Horvitz, D.G., Thompson, D.J. (1952). A generalization of sampling without replacement. J Am Stat Assoc. 47:663–685.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Murthy (1957 Murthy, M.N. (1957). Ordered and unordered estimators in sampling without replacement. Sankhya 18:379–390. [Google Scholar]) under a super-population model. It is shown that the model expected variance is smaller for the Murthy's (1957 Murthy, M.N. (1957). Ordered and unordered estimators in sampling without replacement. Sankhya 18:379–390. [Google Scholar]) strategy both when these two sampling strategies are based on data obtained from (i) a direct survey, and (ii) a randomized response (RR) survey employing some RR technique following a general RR model.     

An extension of the Blinder–Oaxaca decomposition to nonlinear models

This paper develops a general Blinder–Oaxaca decomposition that allows the differences in an outcome variable between two groups to be decomposed into (i) a part that is explained by differences in observed characteristics and (ii) a part attributable to differences in the estimated coefficients also for nonlinear regression models. Based on this general model, we show how it can be applied to different models with discrete and limited dependent variables.     

The uniformly minimum variance unbiased estimator of odds ratio in case–control studies under inverse sampling

The stated goal of this paper is to propose the uniformly minimum variance unbiased estimator of odds ratio in case–control studies under inverse sampling design. The problem of estimating odds ratio plays a central role in case–control studies. However, the traditional sampling schemes appear inadequate when the expected frequencies of not exposed cases and exposed controls can be very low. In such a case, it is convenient to use the inverse sampling design, which requires that random drawings shall be continued until a given number of relevant events has emerged. In this paper we prove that a uniformly minimum variance unbiased estimator of odds ratio does not exist under usual binomial sampling, while the standard odds ratio estimator is uniformly minimum variance unbiased under inverse sampling. In addition, we compare these two sampling schemes by means of large-sample theory and small-sample simulation.     

An extension of the Anderson–Darling k-sample test to arbitrary sample space partition sizes

In this paper we first show that the k-sample Anderson–Darling test is basically an average of Pearson statistics in 2?×?k contingency tables that are induced by observation-based partitions of the sample space. As an extension, we construct a family of rank test statistics, indexed by c?∈??, which is based on similarly constructed c?×?k partitions. An extensive simulation study, in which we compare the new test with others, suggests that generally very high powers are obtained with the new tests. Finally we propose a decomposition of the test statistic in interpretable components.     

Analytic approach to coset–leader decoding leads to extension of the Welch–Berlekamp theorem

This paper provides a complete proof of the Welch–Berlekamp theorem on which the Welch–Berlekamp algorithm was founded. By introducing an analytic approach to coset–leader decoders for Reed–Solomon codes, the Welch–Berlekamp key-equation of error corrections is enlarged and a complete proof of the Welch–Berlekamp theorem is derived in a natural way, and the theorem is extended such that the BCH-bound constraint is moved.     

Optimality of the quasi-score estimator in a mean–variance model with applications to measurement error models

We consider a regression of y$y$ on x$x$ given by a pair of mean and variance functions with a parameter vector θ$\theta$ to be estimated that also appears in the distribution of the regressor variable x$x$. The estimation of θ$\theta$ is based on an extended quasi-score (QS) function. We show that the QS estimator is optimal within a wide class of estimators based on linear-in-y$y$ unbiased estimating functions. Of special interest is the case where the distribution of x$x$ depends only on a subvector α$\alpha$ of θ$\theta$, which may be considered a nuisance parameter. In general, α$\alpha$ must be estimated simultaneously together with the rest of θ$\theta$, but there are cases where α$\alpha$ can be pre-estimated. A major application of this model is the classical measurement error model, where the corrected score (CS) estimator is an alternative to the QS estimator. We derive conditions under which the QS estimator is strictly more efficient than the CS estimator.     

Distribution of the mean reversion estimator in the Ornstein–Uhlenbeck process

ABSTRACT

We derive the exact distribution of the maximum likelihood estimator of the mean reversion parameter (κ) in the Ornstein–Uhlenbeck process using numerical integration through analytical evaluation of a joint characteristic function. Different scenarios are considered: known or unknown drift term, fixed or random start-up value, and zero or positive κ. Monte Carlo results demonstrate the remarkably reliable performance of our exact approach across all the scenarios. In comparison, misleading results may arise under the asymptotic distributions, including the advocated infill asymptotic distribution, which performs poorly in the tails when there is no intercept in the regression and the starting value of the process is nonzero.     

A consistent estimator of the smoothing operator in the functional Hodrick–Prescott filter

In this article, we consider a version of the functional Hodrick–Prescott filter for functional time series. We show that the associated optimal smoothing operator preserves the “noise-to-signal ratio” structure. Moreover, as the main result, we propose a consistent estimator of this optimal smoothing operator.     

Ridge Regression and Generalized Maximum Entropy: An improved version of the Ridge–GME parameter estimator

In this article, the Ridge–GME parameter estimator, which combines Ridge Regression and Generalized Maximum Entropy, is improved in order to eliminate the subjectivity in the analysis of the ridge trace. A serious concern with the visual inspection of the ridge trace to define the supports for the parameters in the Ridge–GME parameter estimator is the misinterpretation of some ridge traces, in particular where some of them are very close to the axes. A simulation study and two empirical applications are used to illustrate the performance of the improved estimator. A MATLAB code is provided as supplementary material.     

A new estimator of Kullback–Leibler information and its application in goodness of fit tests*

We propose here a general statistic for the goodness of fit test of statistical distributions. The proposed statistic is constructed based on an estimate of Kullback–Leibler information. The proposed test is consistent and the limiting distribution of the test statistic is derived. Then, the established results are used to introduce goodness of fit tests for the normal, exponential, Laplace and Weibull distributions. A simulation study is carried out for examining the power of the proposed test and to compare it with those of some existing procedures. Finally, some illustrative examples are presented and analysed, and concluding comments are made.     

Nonparametric estimator for the survival function of quality adjusted lifetime (QAL) in a three-state illness–death model

In this work, we consider nonparametric estimation of QAL distribution in a three-state illness–death model. In our approach, we first write down the expression for the distribution of QAL in terms of the joint distribution of the sojourn times in the three states. The estimate of the QAL distribution is obtained by substituting the estimates of sojourn time distributions in the expression of the QAL distribution. The proposed nonparametric estimate, assuming independence between time to illness and sojourn time in the state of illness, is uniformly consistent. Asymptotic normality has also been established. An estimate of asymptotic variance has been obtained. The performance of the proposed estimator is investigated by simulation. A data set of the Stanford Heart Transplant program has been analyzed for illustration.     

Robustness of the EWMA and the combined X¯–EWMA control charts with variable sampling intervals to non-normality

The exponentially weighted moving average (EWMA) control charts with variable sampling intervals (VSIs) have been shown to be substantially quicker than the fixed sampling intervals (FSI) EWMA control charts in detecting process mean shifts. The usual assumption for designing a control chart is that the data or measurements are normally distributed. However, this assumption may not be true for some processes. In the present paper, the performances of the EWMA and combined X¯–EWMA control charts with VSIs are evaluated under non-normality. It is shown that adding the VSI feature to the EWMA control charts results in very substantial decreases in the expected time to detect shifts in process mean under both normality and non-normality. However, the combined X¯–EWMA chart has its false alarm rate and its detection ability is affected if the process data are not normally distributed.     

Statistical considerations in bioequivalence of two area under the concentration–time curves obtained from serial sampling data

In this paper, we study the bioequivalence (BE) inference problem motivated by pharmacokinetic data that were collected using the serial sampling technique. In serial sampling designs, subjects are independently assigned to one of the two drugs; each subject can be sampled only once, and data are collected at K distinct timepoints from multiple subjects. We consider design and hypothesis testing for the parameter of interest: the area under the concentration–time curve (AUC). Decision rules in demonstrating BE were established using an equivalence test for either the ratio or logarithmic difference of two AUCs. The proposed t-test can deal with cases where two AUCs have unequal variances. To control for the type I error rate, the involved degrees-of-freedom were adjusted using Satterthwaite's approximation. A power formula was derived to allow the determination of necessary sample sizes. Simulation results show that, when the two AUCs have unequal variances, the type I error rate is better controlled by the proposed method compared with a method that only handles equal variances. We also propose an unequal subject allocation method that improves the power relative to that of the equal and symmetric allocation. The methods are illustrated using practical examples.     

Optimal design of repetitive group sampling plans for Weibull and gamma distributions with applications and comparison to the Birnbaum–Saunders distribution

In this paper, we propose a multiple deferred state repetitive group sampling plan which is a new sampling plan developed by incorporating the features of both multiple deferred state sampling plan and repetitive group sampling plan, for assuring Weibull or gamma distributed mean life of the products. The quality of the product is represented by the ratio of true mean life and specified mean life of the products. Two points on the operating characteristic curve approach is used to determine the optimal parameters of the proposed plan. The plan parameters are determined by formulating an optimization problem for various combinations of producer's risk and consumer's risk for both distributions. The sensitivity analysis of the proposed plan is discussed. The implementation of the proposed plan is explained using real-life data and simulated data. The proposed plan under Weibull distribution is compared with the existing sampling plans. The average sample number (ASN) of the proposed plan and failure probability of the product are obtained under Weibull, gamma and Birnbaum–Saunders distributions for a specified value of shape parameter and compared with each other. In addition, a comparative study is made between the ASN of the proposed plan under Weibull and gamma distributions.     

An investigation of the spatial correlations for relative purchasing power in Baden–Württemberg

The relative purchasing power—i.e., the purchasing power per inhabitant—is one of the key characteristics for businesses deciding on site selection. Apart from that it also plays a major role in regional planning, pricing policy and market research. In this study we investigate the spatial correlations for relative purchasing power of the townships in Baden–Württemberg. In particular, changes in relative purchasing power are analysed for three different time intervals, 1987–1993, 1993–1998 and 1998–2004, by means of distance-dependent characteristics like the mark–correlation function, the Simpson indices α(r) and β(r) and by tests on random labelling. It is shown that there are positive correlations for small distances between different townships but that these positive correlations become weaker over the years until they are almost nonexistent (in the sense that hypotheses of random labelling are no longer rejected). A conclusion from this loss of spatial correlations with time is that the relative purchasing power might become more and more purely random. This means that the relative purchasing power in a township is less and less influenced by the relative purchasing power of townships nearby. We further analysed these changes in the Bodensee–Oberschwaben and Stuttgart regions to compare the development of the relative purchasing power in both urban and rural environments. This analysis was done in close cooperation with W. Brachat–Schwarz and W. Walla of the Statistical Office of Baden–Württemberg. S.E. is supported by a grant of the graduate college 1100 at Ulm University.