Let $X_1, \cdots, X_n$ be $n \geqslant 2$ positive random variables and $G(\mathbf{X})$ a positive variable satisfying $G(a\mathbf{X}) = aG(\mathbf{X})$ for all $a > 0$. Then $G$ is a size variable, and $\mathbf{X}/G$ is a shape vector. If $X_1, \cdots, X_n$ are independent, then the independence of shape and the size variable $G(\mathbf{X})$ characterizes (i) the lognormal distribution if $G(\mathbf{X}) = \Pi X^{1/n}_i$, (ii) the generalized gamma distribution if $G(\mathbf{X}) = (\sum X^b_i)^{1/b}$…
Ordinarily a Bayesian estimation procedure uses one prior distribution to obtain a unique estimation rule (its Bayes rule). From the decision theoretical point of view, this procedure can be regarded as a convenient way to obtain admissible decision rules. However, many intuitively appealing, admissible estimation rules cannot be obtained directly in this way. We propose a new mechanism, called the Stepwise Bayesian Procedure (SBP). When the parameter space contains only…
Let $X$ be an observation from a $p$-variate normal distribution $(p \geqslant 3)$ with mean vector $\theta$ and unknown positive definite covariance matrix $\sum$. We wish to estimate $\theta$ under the quadratic loss $L(\delta; \theta, \sum) = \lbrack\mathrm{tr}(Q\sum)\rbrack^{-1}(\delta - \theta)'Q(\delta - \theta)$, where $Q$ is a known positive definite matrix. Estimators of the following form are considered: $\delta_{k, h}(X, W) = \lbrack I - kh(X'W^{-1}X)\lambda_1(QW/n^\ast)Q^{-1}W^{-1}\rbrack X,$…
The multivariate normal patterned mean and covariance matrix testing problem is studied for general one and $k$-population hypotheses. T. W. Anderson's iterative algorithm for finding the maximum likelihood estimates, the forms of the likelihood ratio tests, and asymptotic chi-square distributions of these tests under the null hypothesis are given. The nonnull asymptotic normal distribution is derived using the standard delta method. This derivation involves using several…
Let $T \subseteq I$ be sets of real numbers. Let $\{Y(t): t\in I\}$ be a real time series whose covariance kernel is assumed known and positive definite. The mean is assumed either to be known or to be an unknown member of a known class of functions on $I$. For each fixed $s\in I, Y(s)$ is predicted by a minimum mean square error unbiased linear predictor $\hat{Y}(s)$ based on $\{Y(t): t\in T\}$. If $\hat{y}(s)$ is the evaluation of $\hat{Y}(s)$ given that the sample path…
It is shown that under fairly mild conditions linear combinations of independent normally distributed random variables with random coefficients tend to zero almost everywhere. The result is applied to parameter estimation in linear regression models.
Statistics of the form $H_n = \sum d_{ijn}h_n(X_i, X_j)$ are considered, where $X_1, X_2, \cdots$, are independent and identically distributed random variables, the diagonal terms, $d_{iin}$, are equal to zero, and $h_n(x, y)$ is a symmetric real valued function. The asymptotic normality of such statistics is proven and the result then combined with work of Jogdeo on statistics that are weighted sums of bivariate functions of ranks to find sufficient conditions for…
A vector $\mu = (\mu_1, \cdots, \mu_n)$ is said to be upper [lower] starshaped if $\mu_{m + 1} \geqslant 0 \lbrack \leqslant \mu_{m + 1} \leqslant \bar{\mu}_m\rbrack m = 1, \cdots, n - 1$, where $\bar{\mu}_m$ is a weighted average of $\mu_1, \cdots, \mu_m$. Obtained is the maximum likelihood estimate of $\mu$ when the $\mu_i$'s are the means of $n$ Poisson or normal populations and $\mu$ is known to be starshaped. The method is applied to obtain estimators of IHRA…
The likelihood function is the common basis of all parametric inference. However, with the exception of an ad hoc definition by Fisher, there has been no such unifying basis for prediction of future events, given past observations. This article proposes a definition of predictive likelihood which can help to remove some nonuniqueness problems in sampling-theory predictive inference, and which can produce a simple prediction analog of the Bayesian parametric result, posterior…
Multivariate generalizations of the Wald-Wolfowitz runs statistic and the Smirnov maximum deviation statistic for the two-sample problem are presented. They are based on the minimal spanning tree of the pooled sample points. Some null distribution results are derived and a simulation study of power is reported.
Estimators are given for the two main parameters of the branching process with random environments, namely, the mean of the environmental mean (the average number of offspring per individual) and the mean of the logarithm of the environmental mean. The latter parameter indicates the certain extinction or possible explosion of the process. The consistency and asymptotic normality of the estimators are shown in the context of Smith and Wilkinson. Consistent estimators are also…
Let $\{X_n\}$ be a Markov chain having a stationary transition function and assume that the state set is an arbitrary set in a Euclidean space. The state transition law of the chain is given by a function $F(y|x) = P\lbrack X_{n+1} \leqslant y|X_n = x\rbrack$, which is assumed defined and continuous for all $x$. In this paper we give a statistical procedure for determining a function $F_n(y\mid x)$ on the basis of the sample $\{X_j\}^n_{j=1}, n = 1, 2,\cdots,$ and prove that…
We consider the class of regression functions $\mathscr{M}(F, G) = \{m(x) = E\lbrack Y\mid X = x\rbrack, (X, Y) \in \Pi(F, G)\}$ where $\Pi(F, G)$ denotes the set of random vectors with marginal distributions $F$ and $G$. A characterization of $\mathscr{M}(F, G)$ is given together with a representation for the projection operator it induces in an appropriate Hilbert space.
A population consisting of an unknown number of distinct species is searched by selecting one member at a time. No a priori information is available concerning the probability that an object selected from this population will represent a particular species. Based on the information available after an $n$-stage search it is desired to predict the conditional probability that the next selection will represent a species not represented in the $n$-stage sample. Properties of a…
T. Kowalczyk · A. Kowalski · A. Matuszewski · E. Pleszczynska·刊期 1979-05-01 · 卷 7(3)·The Annals of Statistics已编入正刊
多元情形下的筛查函数与单调依赖函数AI 解读
An approach to simultaneous treatment of dependence and screening problems is presented. New characterizations of dependence of a random variable $X$ on a random vector $Y$ are obtained by functions $\nu_{X, Y}: (0, 1)\rightarrow \lbrack 0, 1\rbrack$ and $\mu_{X, Y} : (0, 1) \rightarrow \lbrack -1, 1\rbrack$ called respectively screening and monotonic dependence functions. These functions are shown to be appropriate measures of the intensity of connection and concordance of $X$…
It has been conjectured by Khatri and Mardia that with probability one MLEs for the parameters of the von Mises-Fisher matrix distribution exist and are unique. We prove that, except for small sample sizes, this conjecture is true, both in the case where the parameter matrix has known rank and in the unrestricted case. The corresponding result for the matrix Bingham distribution is proven also.
In multivariate analysis under normality assumptions, many likelihood ratio criteria $(\lambda^{(n)})$ are distributed as $k\prod^m_{i=1} Z^a_{li}(1 - Z_{li})^{b_i}\prod^{m'}_{j=1} Z^{c_j}_{2j}$ for some constants, $k, m, m', a_i, b_i,$ and $c_j$ when their associated null hypotheses are true, where $Z_{ij}$ are independently distributed beta variates. Let $T^{(n)} = -n^{-1} \ln \lambda^{(n)}$. This paper shows that a sequence $\{T^{(n)}\}$ of this kind is asymptotically…
This paper deals with approximations to Bayesian sequential tests of composite hypotheses. If the distributions of the data form an exponential or truncation family, then such tests may be described by a continuation region in the space of $n$, the sample size, and $M_n$, the sufficient statistics, which are of fixed dimension. In this case Schwarz has been able to describe the asymptotic shape of the continuation region as the sampling cost $c$ approaches zero. We have…
The theory of Dirichlet processes is applied to the empirical Bayes estimation problem in the binomial case. The approach is Bayesian rather than being empirical Bayesian. When the prior is a Dirichlet process the posterior is a mixture of Dirichlet processes. Explicit estimators are given for the case of 2 and 3 parameters and compared with other empirical Bayes estimators by way of examples. Since the number of calculations become enormous when the number of parameters…
Let $(X_1, Y_1),\cdots, (X_N, Y_N)$ be i.i.d. rv's where the $X$'s are nonnegative integer-valued. Conditional on $\sigma X_k$ the asymptotic distribution of $\sigma Y_k$ and $\sigma a_k X_k$ are derived by general methods. Some applications are briefly discussed: sampling without replacement, the classical occupancy problem, the Wilcoxon statistic, the Poisson index of dispersion, testing geometric versus Poisson distribution.
R. L. Anderson proposed a $\chi^2$-type rank statistic for the nonparametric analysis of a randomized blocks design. In this paper the asymptotic distribution of the test statistic is derived under a sequence of alternatives contiguous to the null hypothesis. Using Bahadur's concept of local approximate slope, it is shown that the test is optimal within a large class of rank tests including the tests proposed by Friedman and by Page. The results are extended to BIB designs.…
Estimation as a reduction of data is usually accompanied by some loss of information. This paper theoretically compares asymptotically efficient estimation methods for parameters in Gaussian linear processes. By means of the concept of "asymptotic information loss" suitably defined, estimates equivalent to the order of $N^{-\frac{1}{2}}$ are differentiated. This problem was studied by C. R. Rao for multinomial distributions and by K. Takeuchi for the exponential family of…
An infinite distributed lag system is approximated by a truncated one, which can be consistently estimated. We investigate the rate at which the parameter space in the approximate model can be increased as sample size increases, in such a way that these estimates will provide consistent estimates of the underlying system.
A general finite parameter model for stationary ergodic nondeterministic vector time series is considered. A central limit theorem for parameter estimates, obtained by maximising frequency domain approximations to the Gaussian likelihood, is established. The treatment given extends the central limit theorem of Dunsmuir and Hannan in that the innovations covariance matrix and the linear transfer function need not be separately parameterised. Models for a stationary vector…
Let X1,⋯,Xn be n⩾2 positive random variables and G(X) a positive variable satisfying G(aX)=aG(X) for all a>0. Then G is a size variable, and X/G is a shape vector. If X1,⋯,Xn are independent, then the independence of shape and the size variable G(X) characterizes (i) the lognormal distribution if G(X)=ΠXi1/n, (ii) the generalized gamma distribution if G(X)=(∑Xib)1/b…
Let X be an observation from a p-variate normal distribution (p⩾3) with mean vector θ and unknown positive definite covariance matrix ∑. We wish to estimate θ under the quadratic loss L(δ;θ,∑)=[tr(Q∑)]−1(δ−θ)′Q(δ−θ), where Q is a known positive definite matrix. Estimators of the following form are considered: δk,h(X,W)=[I−kh(X′W−1X)λ1(QW/n∗)Q−1W−1]X,…
The multivariate normal patterned mean and covariance matrix testing problem is studied for general one and k-population hypotheses. T. W. Anderson's iterative algorithm for finding the maximum likelihood estimates, the forms of the likelihood ratio tests, and asymptotic chi-square distributions of these tests under the null hypothesis are given. The nonnull asymptotic normal distribution is derived using the standard delta method. This derivation involves using several…
Let T⊆I be sets of real numbers. Let {Y(t):t∈I} be a real time series whose covariance kernel is assumed known and positive definite. The mean is assumed either to be known or to be an unknown member of a known class of functions on I. For each fixed s∈I,Y(s) is predicted by a minimum mean square error unbiased linear predictor Y^(s) based on {Y(t):t∈T}. If y^(s) is the evaluation of Y^(s) given that the sample path…
Statistics of the form Hn=∑dijnhn(Xi,Xj) are considered, where X1,X2,⋯, are independent and identically distributed random variables, the diagonal terms, diin, are equal to zero, and hn(x,y) is a symmetric real valued function. The asymptotic normality of such statistics is proven and the result then combined with work of Jogdeo on statistics that are weighted sums of bivariate functions of ranks to find sufficient conditions for…
A vector μ=(μ1,⋯,μn) is said to be upper [lower] starshaped if μm+1⩾0[⩽μm+1⩽μˉm]m=1,⋯,n−1, where μˉm is a weighted average of μ1,⋯,μm. Obtained is the maximum likelihood estimate of μ when the μi's are the means of n Poisson or normal populations and μ is known to be starshaped. The method is applied to obtain estimators of IHRA…
Let {Xn} be a Markov chain having a stationary transition function and assume that the state set is an arbitrary set in a Euclidean space. The state transition law of the chain is given by a function F(y∣x)=P[Xn+1⩽y∣Xn=x], which is assumed defined and continuous for all x. In this paper we give a statistical procedure for determining a function Fn(y∣x) on the basis of the sample {Xj}j=1n,n=1,2,⋯, and prove that…
We consider the class of regression functions M(F,G)={m(x)=E[Y∣X=x],(X,Y)∈Π(F,G)} where Π(F,G) denotes the set of random vectors with marginal distributions F and G. A characterization of M(F,G) is given together with a representation for the projection operator it induces in an appropriate Hilbert space.
A population consisting of an unknown number of distinct species is searched by selecting one member at a time. No a priori information is available concerning the probability that an object selected from this population will represent a particular species. Based on the information available after an n-stage search it is desired to predict the conditional probability that the next selection will represent a species not represented in the n-stage sample. Properties of a…
An approach to simultaneous treatment of dependence and screening problems is presented. New characterizations of dependence of a random variable X on a random vector Y are obtained by functions νX,Y:(0,1)→[0,1] and μX,Y:(0,1)→[−1,1] called respectively screening and monotonic dependence functions. These functions are shown to be appropriate measures of the intensity of connection and concordance of X…
In multivariate analysis under normality assumptions, many likelihood ratio criteria (λ(n)) are distributed as k∏i=1mZlia(1−Zli)bi∏j=1m′Z2jcj for some constants, k,m,m′,ai,bi, and cj when their associated null hypotheses are true, where Zij are independently distributed beta variates. Let T(n)=−n−1lnλ(n). This paper shows that a sequence {T(n)} of this kind is asymptotically…
This paper deals with approximations to Bayesian sequential tests of composite hypotheses. If the distributions of the data form an exponential or truncation family, then such tests may be described by a continuation region in the space of n, the sample size, and Mn, the sufficient statistics, which are of fixed dimension. In this case Schwarz has been able to describe the asymptotic shape of the continuation region as the sampling cost c approaches zero. We have…
Let (X1,Y1),⋯,(XN,YN) be i.i.d. rv's where the X's are nonnegative integer-valued. Conditional on σXk the asymptotic distribution of σYk and σakXk are derived by general methods. Some applications are briefly discussed: sampling without replacement, the classical occupancy problem, the Wilcoxon statistic, the Poisson index of dispersion, testing geometric versus Poisson distribution.
R. L. Anderson proposed a χ2-type rank statistic for the nonparametric analysis of a randomized blocks design. In this paper the asymptotic distribution of the test statistic is derived under a sequence of alternatives contiguous to the null hypothesis. Using Bahadur's concept of local approximate slope, it is shown that the test is optimal within a large class of rank tests including the tests proposed by Friedman and by Page. The results are extended to BIB designs.…
Estimation as a reduction of data is usually accompanied by some loss of information. This paper theoretically compares asymptotically efficient estimation methods for parameters in Gaussian linear processes. By means of the concept of "asymptotic information loss" suitably defined, estimates equivalent to the order of N−21 are differentiated. This problem was studied by C. R. Rao for multinomial distributions and by K. Takeuchi for the exponential family of…