Probability density estimation with data missing at random when covariables are present

This paper addresses the problem of the probability density estimation in the presence of covariates when data are missing at random (MAR). The inverse probability weighted method is used to define a nonparametric and a semiparametric weighted probability density estimators. A regression calibration technique is also used to define an imputed estimator. It is shown that all the estimators are asymptotically normal with the same asymptotic variance as that of the inverse probability weighted estimator with known selection probability function and weights. Also, we establish the mean squared error (MSE) bounds and obtain the MSE convergence rates. A simulation is carried out to assess the proposed estimators in terms of the bias and standard error.     

Lifetesting and estimation with arbitrary distribution function

We consider the problem of estimating the life–distribution F from censored lifetimes. The observation scheme is renewal testing over a long time horizon although the results can apply to survival testing with repetitions. We exhibit a product–limit estimator of F which is shown to be consistent and to converge weakly to a GAUSsian process. To do this we first extend these properties of the NELSON-AALEN martingale estimator to the family of PoissoN–type counting processes. Our proof of weak convergence is based on the general functional central limit theorems for semimartingales as developed by .JACOB, SHIRYAYEV and others     

On the detection of changes in autoregressive time series,II. Resampling procedures

We study an autoregressive time series model with a possible change in the regression parameters. Approximations to the critical values for change-point tests are obtained through various bootstrapping methods. Theoretical results show that the bootstrapping procedures have the same limiting behavior as their asymptotic counterparts discussed in Hušková et al. [2007. On the detection of changes in autoregressive time series, I. Asymptotics. J. Statist. Plann. Inference 137, 1243–1259]. In fact, a small simulation study illustrates that the bootstrap tests behave better than the original asymptotic tests if performance is measured by the α$\alpha$- and β$\beta$-errors, respectively.     

Confidence bands for ROC curves

We develop two methods to construct confidence bands for the receiver operating characteristic (ROC) curve without estimating the densities of the underlying distributions. The first method is based on the smoothed bootstrap while the second method uses the Bonferroni inequality. As an illustration, we provide confidence bands for the ROC curve using data on Duchanne Muscular Dystrophy.     

Nonparametric estimation of conditional expectation

Denote the integer lattice points in the N  -dimensional Euclidean space by Z^N${\mathbb{Z}}^N$ and assume that (_Xi,_Yi)$(X_{\mathbf{i}},Y_{\mathbf{i}})$, i∈Z^N$\mathbf{i}\in {\mathbb{Z}}^N$ is a mixing random field. Estimators of the conditional expectation r(x)=E[_Yi|_Xi=x]$r(\mathbf{x})=\mathrm{E}[Y_{\mathbf{i}}|X_{\mathbf{i}}=\mathbf{x}]$ by nearest neighbor methods are established and investigated. The main analytical result of this study is that, under general mixing assumptions, the estimators considered are asymptotically normal. Many difficulties arise since points in higher dimensional space N?2$N?2$ cannot be linearly ordered. Our result applies to many situations where parametric methods cannot be adopted with confidence.     

Asymptotic properties of nonparametric M-estimation for mixing functional data

We investigate the asymptotic behavior of a nonparametric M-estimator of a regression function for stationary dependent processes, where the explanatory variables take values in some abstract functional space. Under some regularity conditions, we give the weak and strong consistency of the estimator as well as its asymptotic normality. We also give two examples of functional processes that satisfy the mixing conditions assumed in this paper. Furthermore, a simulated example is presented to examine the finite sample performance of the proposed estimator.     

Improved reduced-bias tail index and quantile estimators

In this paper, we deal with bias reduction techniques for heavy tails, trying to improve mainly upon the performance of classical high quantile estimators. High quantiles depend strongly on the tail index γ$\gamma$, for which new classes of reduced-bias estimators have recently been introduced, where the second-order parameters in the bias are estimated at a level _k1$k_1$ of a larger order than the level k at which the tail index is estimated. Doing this, it was seen that the asymptotic variance of the new estimators could be kept equal to the one of the popular Hill estimators. In a similar way, we now introduce new classes of tail index and associated high quantile estimators, with an asymptotic mean squared error smaller than that of the classical ones for all k in a large class of heavy-tailed models. We derive their asymptotic distributional properties and compare them with those of alternative estimators. Next to that, an illustration of the finite sample behavior of the estimators is also provided through a Monte Carlo simulation study and the application to a set of real data in the field of insurance.     

Semiparametric estimation for count data through weighted distributions

This paper is concerned with semiparametric discrete kernel estimators when the unknown count distribution can be considered to have a general weighted Poisson form. The estimator is constructed by multiplying the Poisson estimate with a nonparametric discrete kernel-type estimate of the Poisson weight function. Comparisons are then carried out with the ordinary discrete kernel probability mass function estimators. The Poisson weight function is thus a local multiplicative correction factor, and is considered as the uniform measure to detect departures from the equidispersed Poisson distribution. In this way, the effects of dispersion and zero-proportion with respect to the standard Poisson distribution are also minimized. This method of estimation is also applied to the weighted binomial form for the count distribution having a finite support. The proposed estimators, in addition to being simple, easy-to-implement and effective, also outperform the competing nonparametric and parametric estimators in finite-sample situations. Two examples illustrate this new semiparametric estimation.     

On statistical models for multimodal optimization

A new approach for the construction of statistical models for multimodal optimization is proposed; the examples of such models are given. The results of the psychological experiment show that the proposed approach is intuitively acceptable.     

Wilcoxon–Mann–Whitney test for stratified samples and Efron's paradox dice

Two-treatment multi-center clinical trials are the most common type of clinical trials in practice. The aim of this paper is to discuss a curious property of certain standard nonparametric procedures used in the analysis of such clinical trials. Different analyses of a simulated data example are presented, which lead to contrasting and surprising results. The source of the potentially misleading outcome is then explored while relating the simulated data with the concept of Efron's paradox dice and the notion of nontransitivity. With the root of the problem established, an alternate nonparametric method from the literature is shown to address the problem. Finally, pointing out an interpretational concern of using the alternate procedure, a modification to this procedure is also suggested and corresponding theoretical results are presented.     

Multivariate inverse Gaussian distribution as a limit of multivariate waiting time distributions

Multivariate inverse Gaussian distribution proposed by Minami [2003. A multivariate extension of inverse Gaussian distribution derived from inverse relationship. Commun. Statist. Theory Methods 32(12), 2285–2304] was derived through multivariate inverse relationship with multivariate Gaussian distributions and characterized as the distribution of the location at a certain stopping time of a multivariate Brownian motion. In this paper, we show that the multivariate inverse Gaussian distribution is also a limiting distribution of multivariate Lagrange distributions, which is a family of waiting time distributions, under certain conditions.     

Adaptation over parametric families of symmetric linear estimators

This paper treats an abstract parametric family of symmetric linear estimators for the mean vector of a standard linear model. The estimator in this family that has smallest estimated quadratic risk is shown to attain, asymptotically, the smallest risk achievable over all candidate estimators in the family. The asymptotic analysis is carried out under a strong Gauss–Markov form of the linear model in which the dimension of the regression space tends to infinity. Leading examples to which the results apply include: (a) penalized least squares fits constrained by multiple, weighted, quadratic penalties; and (b) running, symmetrically weighted, means. In both instances, the weights define a parameter vector whose natural domain is a continuum.     

Simultaneous control and estimation of linear stochastic systems with unknown parameters of disturbances

In the present paper methods of the decision theory are applied to determine minimax policies of simultaneous control and estimation for stochastic system (1). There are solved the following cases:

a)when disturbances of the system have the binomial distribution with unknown parameter

b)when disturbances of the system have distribution belonging to an exponential family dependent on natural unknown parameter and the class of prior distributions of parameter is restricted by fixing the second moment

In both cases open analytical forms of minimax policies are given     

Projection estimation capacity of Hadamard designs

In a screening design, often only a few factors among a large number of potential factors are significantly important. Usually, it is not known which factors will be important ones. Thus, it is of practical interest to know if each projection of a design onto a small subset of factors is able to entertain and estimate all two-factor-interactions along with its main effects, assuming higher order interactions are negligible. In this paper, we investigate the estimation capacity of projections of Hadamard designs with run size up to 60. Possible applications of our results to robust parameter designs are also discussed.     

Design selection criteria for discrimination/estimation for nested models and a binomial response

The aim of an experiment is often to enable discrimination between competing forms for a response model. We investigate the selection of a continuous design for a non-sequential strategy when there are two competing generalized linear models for a binomial response, with a common link function and the linear predictor of one model nested within that of the other. A new criterion, _TE$T_E$-optimality, is defined, based on the difference in the deviances from the two models, and comparisons are made with T$T$-, _Ds$D_s$- and D$D$-optimality. Issues are raised through the study of two examples in which designs are assessed using simulation studies of the power to reject the null hypothesis of the smaller model being correct, when the data are generated from the larger model. Parameter estimation for discrimination designs is also discussed and a simple method is investigated of combining designs to form a hybrid design in order to achieve both model discrimination and estimation. This method has a computational advantage over the use of a compound criterion and the similar performance of the designs obtained from the two approaches is illustrated in an example.     

A bias correction and acceleration approach for the problem of regions

For testing the problem of regions in the space of distribution functions, this paper considers approaches to modify the bootstrap probability to be a second-order accurate p$p$-value based on the familiar bias correction and acceleration method. It is shown that Shimodaira's [2004a. Approximately unbiased tests of regions using multistep-multiscale bootstrap resampling. Ann. Statist. 32, 2616–2641] twostep-multiscale bootstrap method works even in the problem of regions in functional space. In this paper the bias correction quantity is estimated by his onestep-multiscale bootstrap method. Instead of using the twostep-multiscale bootstrap method, the acceleration constant is estimated by a newly proposed jackknife method which requires first-level bootstrap resamplings only. Some numerical examples are illustrated, in which an application to testing significance in model selection is included.     

Die waldsche fundamentalidentität und ein sequentieller quotiententest für eine zufäallige irrfahrtüber einer homogenen irreduziblen markovschen kette mit endlichem zustandsraum

For a random walk {R _n ≥0} given on a homogeneous irreducible finite MARKOV chain {X _n ≥0} the identity (8) is obtained. Generalizations (14)-(16) of WALD's Fundamental Identity and WALD's first and second equations for the two-dimensional process {(R _n ,X _n ), n≥0} are proved. The Average Sample Number (21)-(22) and the Operating Characteristic Function (24)-(25) of a Sequential Probability Ratio Test follow. With this test a decision about two simple hypotheses on the unknown transition probability matrix of {X _n , n≥0} and the unknown parameters of the probability distributions for the increments of {X _n , n≥0} can be made. For a special case these results were proved by PHATARFOD [6] and KüCHLER [3] with other methods.     

Eine bemerkung zur stetigkeit zufälliger felder

In the paper a sequence of bounded regions containing n independent identically and uniformly on D_n distributed points is considered. It is assumed that the d–dimensional volume v(D_n) is asymptotically proportional to n. Under these conditions it is shown that the number of pairs of points within a distance r>0 of each other is asymptotically normally distributed. For proving this among other things a lemma of BOLTHAUSEN is used, whereas even strong estimates for U–statistics are insufficient. The obtained result is applied for testing the hypothesis of randomness