Necessary and sufficient conditions for the weak consistency of the sample median of independent, but not identically distributed random variables are given and discussed.
Dominique Fourdrinier · William E. Strawderman · Martin T. Wells·刊期 1998-04-01 · 卷 26(2)·The Annals of Statistics已编入正刊
贝叶斯极小极大估计量的构造AI 解读
Bayes estimation of the mean of a multivariate normal distribution is considered under quadratic loss. We show that, when particular spherical priors are used, the superharmonicity of the square root of the marginal density provides a viable method for constructing (possibly proper) Bayes (and admissible) minimax estimators. Examples illustrate the theory; most notably it is shown that a multivariate Student-$t$ prior yields a proper Bayes minimax estimate.
Selection models arise when the data are selected to enter the sample only if they occur in a certain region of the sample space. When this selection occurs according to some probability distribution, the resulting model is often instead called a weighted distribution model. In either case the "original" density becomes multiplied by a "weight function" $w(x)$. Often there is considerable uncertainty concerning this weight function; for instance, it may be known only that $w$…
We prove that the posterior distribution in a possibly incorrect parametric model a.s. concentrates in a strong sense on the set of pseudotrue parameters determined by the true distribution. As a consequence, we obtain in the case of a unique pseudotrue parameter the strong consistency of pseudo-Bayes estimators w.r.t. general loss functions. Further, we present a simple example based on normal distributions and having two different pseudotrue parameters, where pseudo-Bayes…
Nonlinear principal components for an absolutely continuous random vector X with positive bounded density are defined as the solution of a variational problem in a suitable function space. In this way transformations depending on all the components of X are obtained. Some properties of nonlinear principal components are proved: in particular, it is shown that the set of nonlinear principal transformations of X is an orthonormal basis for the function space associated with…
If $U$ and $V$ are independent random variables which are gamma distributed with the same scale parameter, then there exist $a$ and $b$ in $\\mathbb{R}$ such that $$\\mathbb{E}(U|U + V) = a(U + V)$$ and $$\\mathbb{E}(U^2|U + V) = b(U + V)^2$$. This, in fact, is characteristic of gamma distributions. Our paper extends this property to the Wishart distributions in a suitable way, by replacing the real number $U^2$ by a pair of quadratic functions of the symmetric matrix $U$.…
A class of multivariate normal models with symmetry restrictions given by a finite group and conditional independence restrictions given by a finite distributive lattice is defined and studied. The statistical properties of these models including maximum likelihood inference, invariance and hypothesis testing are discussed.
In this paper interest is focused on some theoretical investigations concerning the comparison of two popular multiple test procedures, so-called step-down and step-up procedures, in terms of their defining critical values. Such procedures can be applied, for example, to multiple comparisons with a control. In the definition of the critical values for these procedures order statistics play a central role. For $k \epsilon N_0$ fixed we consider the joint cumulative…
Some new probability inequalities involving the ordered components of an $MTP_2$ random vector are derived, which provide an analytical proof of an important conjecture in the field of multiple hypothesis testing. This conjecture has been mostly validated so far using simulation.
Samples of directional or axial measurements arise in geophysical, biological and econometric contexts. We represent the rotational difference between two mean directions (or two mean axes) as a direction (or axis). We then construct nonparametric simultaneous confidence sets for all pair-wise rotational differences among the mean directions or mean axes of $s$ samples. By specialization, this methodology yields nonparametric simultaneous tests for pairwise equality of…
Hall and Welsh established the best attainable rate of convergence for estimates of a positive extreme value index $\gamma$ under a certain second order condition implying that the distribution function of the maximum of n random variables converges at an algebraic rate to the pertaining extreme value distribution. As a first generalization, we obtain a surprisingly sharp bound on the estimation error if $\gamma$ is still assumed to be positive, but the rate of convergence…
We begin by showing how piecewise linear bounds may be devised, which bound both above and below any concave log-density in general dimensions. We then show how these bounds may be used to gain an upper bound to the volume in the tails outside the convex hull of the sample path in order to assess how well the sampler has explored the target distribution. This method can be used as a stand-alone diagnostic to determine when the sampler output provides a reliable basis for…
We construct Markov chain algorithms for sampling from discrete exponential families conditional on a sufficient statistic. Examples include contingency tables, logistic regression, and spectral analysis of permutation data. The algorithms involve computations in polynomial rings using Gröbner bases.
We estimate a probability density function p which is related by a linear operator K to a density function q in sequences of regular exponential families based on a random sample from q. In this paper deconvolution and positron emission tomography are considered. The logarithm of the density function is approximated by basis functions consisting of singular functions of K. While direct maximum likelihood (or minimum Kullback-Leibler) density estimation in exponential…
When using nonparametric estimates of the mean curve, surface or image underlying noisy observations, the selection of "smoothing parameters" is generally crucial. This paper gives a theoretical comparison of the performances of generalized cross-validation (GCV) and of its fast randomized version (RGCV), as selection criteria. This is mainly done by studying the asymptotic distribution of the excess error for each selector, that is, the difference between the (data-driven)…
Suppose one observes a process V on the unit interval, where $dV = f_o + dW$ with an unknown parameter $f_o \epsilon L_1[0, 1]$ and standard Brownian motion W. We propose a particular test of one-point hypotheses about $f_o$ which is based on suitably standardized increments of V. This test is shown to have desirable consistency properties if, for instance, $f_o$ is restricted to various Hölder classes of functions. The test is mimicked in the context of nonparametric…
An example is provided to show that the natural asymptotic equivalence does not hold between any pairs of three nonparametric experiments: density problem, white noise with drift and nonparametric regression, when the smoothness index of the unknown nonparametric function class is ½
A data-driven estimate is given that, over a Sobolev space, is simultaneously asymptotically sharp minimax for estimating both the function and its derivatives under integrated squared error loss. It is also shown that linear estimates cannot be simultaneously asymptotically sharp minimax over a given Sobolev space.
Jens P. Nielsen · Oliver Linton · Peter J. Bickel·刊期 1998-02-01 · 卷 26(1)·The Annals of Statistics已编入正刊
具有灵活协变量效应的半参数生存模型AI 解读
A semiparametric hazard model with parametrized time but general covariate dependency is formulated and analyzed inside the framework of counting process theory. A profile likelihood principle is introduced for estimation of the parameters: the resulting estimator is $n^{1/2}$-consistent, asymptotically normal and achieves the semiparametric efficiency bound. An estimation procedure for the nonparametric part is also given and its asymptotic properties are derived. We…
The frailty model is a generalization of Cox's proportional hazard model, where a shared unobserved quantity in the intensity induces a positive correlation among the survival times. Murphy showed consistency and asymptotic normality of the nonparametric maximum likelihood estimator (NPMLE) for the shared gamma-frailty model without covariates. In this paper we extend this result to the correlated gamma-frailty model, and we allow for covariates. We discuss the definition of…
The Kaplan-Meier estimator of a survival function is well known to be asymptotically efficient when cause of failure is always observed. It has been an open problem, however, to find an efficient estimator when failure indicators are missing at random. Lo showed that nonparametric maximum likelihood estimators are inconsistent, and this has led to several proposals of ad hoc estimators, none of which are efficient. We now introduce a sieved nonparametric maximum likelihood…
Jane-Ling Wang · Hans-Georg Müller · William B. Capra·刊期 1998-02-01 · 卷 26(1)·The Annals of Statistics已编入正刊
最高龄死亡率分析:重新审视生命表AI 解读
This paper provides a data analysis and some methodological advances which contribute to an ongoing scientific debate about the patterns of aging. One of the problems we address is how to estimate a hazard function when only aggregated information on the lifetimes in the form of a lifetable is available. This problem affects the the estimation of oldest-old mortality which in turn plays an important role in the quantification of biological lifespan and longevity. We…
This paper investigates the maximum likelihood estimator (MLE) for unstable autoregressive moving-average (ARMA) time series with the noise sequence satisfying a general autoregressive heteroscedastic (GARCH) process. Under some mild conditions, it is shown that the MLE satisfying the likelihood equation exists and is consistent. The limiting distribution of the MLE is derived in a unified manner for all types of characteristic roots on or outside the unit circle and is…
We propose a sieve bootstrap procedure for time series with a deterministic trend. The sieve for constructing the bootstrap is based on nonparametric trend estimation and autoregressive approximation for some noise process. The bootstrap scheme itself does i.i.d. resampling of estimated innovations from fitted autoregressive models. We show the validity and indicate second-order correctness of such sieve bootstrap approximations for the limiting distribution of nonparametric…
Bayes estimation of the mean of a multivariate normal distribution is considered under quadratic loss. We show that, when particular spherical priors are used, the superharmonicity of the square root of the marginal density provides a viable method for constructing (possibly proper) Bayes (and admissible) minimax estimators. Examples illustrate the theory; most notably it is shown that a multivariate Student-t prior yields a proper Bayes minimax estimate.
Selection models arise when the data are selected to enter the sample only if they occur in a certain region of the sample space. When this selection occurs according to some probability distribution, the resulting model is often instead called a weighted distribution model. In either case the "original" density becomes multiplied by a "weight function" w(x). Often there is considerable uncertainty concerning this weight function; for instance, it may be known only that w…
If U and V are independent random variables which are gamma distributed with the same scale parameter, then there exist a and b in mathbbR such that mathbbE(U∣U+V)=a(U+V) and mathbbE(U2∣U+V)=b(U+V)2. This, in fact, is characteristic of gamma distributions. Our paper extends this property to the Wishart distributions in a suitable way, by replacing the real number U2 by a pair of quadratic functions of the symmetric matrix U.…
In this paper interest is focused on some theoretical investigations concerning the comparison of two popular multiple test procedures, so-called step-down and step-up procedures, in terms of their defining critical values. Such procedures can be applied, for example, to multiple comparisons with a control. In the definition of the critical values for these procedures order statistics play a central role. For kϵN0 fixed we consider the joint cumulative…
Some new probability inequalities involving the ordered components of an MTP2 random vector are derived, which provide an analytical proof of an important conjecture in the field of multiple hypothesis testing. This conjecture has been mostly validated so far using simulation.
Samples of directional or axial measurements arise in geophysical, biological and econometric contexts. We represent the rotational difference between two mean directions (or two mean axes) as a direction (or axis). We then construct nonparametric simultaneous confidence sets for all pair-wise rotational differences among the mean directions or mean axes of s samples. By specialization, this methodology yields nonparametric simultaneous tests for pairwise equality of…
Hall and Welsh established the best attainable rate of convergence for estimates of a positive extreme value index γ under a certain second order condition implying that the distribution function of the maximum of n random variables converges at an algebraic rate to the pertaining extreme value distribution. As a first generalization, we obtain a surprisingly sharp bound on the estimation error if γ is still assumed to be positive, but the rate of convergence…
Suppose one observes a process V on the unit interval, where dV=fo+dW with an unknown parameter foϵL1[0,1] and standard Brownian motion W. We propose a particular test of one-point hypotheses about fo which is based on suitably standardized increments of V. This test is shown to have desirable consistency properties if, for instance, fo is restricted to various Hölder classes of functions. The test is mimicked in the context of nonparametric…
A semiparametric hazard model with parametrized time but general covariate dependency is formulated and analyzed inside the framework of counting process theory. A profile likelihood principle is introduced for estimation of the parameters: the resulting estimator is n1/2-consistent, asymptotically normal and achieves the semiparametric efficiency bound. An estimation procedure for the nonparametric part is also given and its asymptotic properties are derived. We…