共查询到20条相似文献,搜索用时 343 毫秒
1.
In hierarchical mixture models the Dirichlet process is used to specify latent patterns of heterogeneity, particularly when the distribution of latent parameters is thought to be clustered (multimodal). The parameters of a Dirichlet process include a precision parameter α and a base probability measure G0. In problems where α is unknown and must be estimated, inferences about the level of clustering can be sensitive to the choice of prior assumed for α. In this paper an approach is developed for computing a prior for the precision parameter α that can be used in the presence or absence of prior information about the level of clustering. This approach is illustrated in an analysis of counts of stream fishes. The results of this fully Bayesian analysis are compared with an empirical Bayes analysis of the same data and with a Bayesian analysis based on an alternative commonly used prior. 相似文献
2.
We consider the problem of estimating the mean θ of an Np(θ,Ip) distribution with squared error loss ∥δ−θ∥2 and under the constraint ∥θ∥≤m, for some constant m>0. Using Stein's identity to obtain unbiased estimates of risk, Karlin's sign change arguments, and conditional risk analysis, we compare the risk performance of truncated linear estimators with that of the maximum likelihood estimator δmle. We obtain for fixed (m,p) sufficient conditions for dominance. An asymptotic framework is developed, where we demonstrate that the truncated linear minimax estimator dominates δmle, and where we obtain simple and accurate measures of relative improvement in risk. Numerical evaluations illustrate the effectiveness of the asymptotic framework for approximating the risks for moderate or large values of p. 相似文献
3.
In common with other non-linear models, the optimal design for a limiting dilution assay (LDA) depends on the value of the unknown parameter, θ, in the model. Consequently optimal designs cannot be specified unless some assumptions are made about the possible values of θ. If a prior distribution can be specified then a Bayesian approach can be adopted. A proper specification of the Bayesian approach requires the aim of the experiment to be described and quantified through an appropriate utility function. This paper addresses the problem of finding optimal designs for LDAs when the aim is to determine whether θ is above or below a specified threshold, θ0. 相似文献
4.
We determine a credible set A that is the “best” with respect to the variation of the prior distribution in a neighborhood Γ of the starting prior π0(θ). Among the class of sets with credibility γ under π0, the “optimally robust” set will be the one which maximizes the minimum probability of including θ as the prior varies over Γ. This procedure is also Γ-minimax with respect to the risk function, probability of non-inclusion. We find the optimally robust credible set for three neighborhood classes Γ, the ε-contamination class, the density ratio class and the density bounded class. A consequence of this investigation is that the maximum likelihood set is seen to be an optimal credible set from a robustness perspective. 相似文献
5.
The probability function and binomial moments of the number Nn of (upper) records up to time (index) n in a geometrically increasing population are obtained in terms of the signless q-Stirling numbers of the first kind, with q being the inverse of the proportion λ of the geometric progression. Further, a strong law of large numbers and a central limit theorem for the sequence of random variables Nn, n=1,2,…, are deduced. As a corollary the probability function of the time Tk of the kth record is also expressed in terms of the signless q -Stirling numbers of the first kind. The mean of Tk is obtained as a q -series with terms of alternating sign. Finally, the probability function of the inter-record time Wk=Tk-Tk-1 is obtained as a sum of a finite number of terms of q -numbers. The mean of Wk is expressed by a q-series. As k increases to infinity the distribution of Wk converges to a geometric distribution with failure probability q. Additional properties of the q-Stirling numbers of the first kind, which facilitate the present study, are derived. 相似文献
6.
7.
In this paper, we study a random field U?(t,x) governed by some type of stochastic partial differential equations with an unknown parameter θ and a small noise ?. We construct an estimator of θ based on the continuous observation of N Fourier coefficients of U?(t,x), and prove the strong convergence and asymptotic normality of the estimator when the noise ? tends to zero. 相似文献
8.
9.
We consider a linear regression model with regression parameter β=(β1,…,βp) and independent and identically N(0,σ2) distributed errors. Suppose that the parameter of interest is θ=aTβ where a is a specified vector. Define the parameter τ=cTβ-t where the vector c and the number t are specified and a and c are linearly independent. Also suppose that we have uncertain prior information that τ=0. We present a new frequentist 1-α confidence interval for θ that utilizes this prior information. We require this confidence interval to (a) have endpoints that are continuous functions of the data and (b) coincide with the standard 1-α confidence interval when the data strongly contradict this prior information. This interval is optimal in the sense that it has minimum weighted average expected length where the largest weight is given to this expected length when τ=0. This minimization leads to an interval that has the following desirable properties. This interval has expected length that (a) is relatively small when the prior information about τ is correct and (b) has a maximum value that is not too large. The following problem will be used to illustrate the application of this new confidence interval. Consider a 2×2 factorial experiment with 20 replicates. Suppose that the parameter of interest θ is a specified simple effect and that we have uncertain prior information that the two-factor interaction is zero. Our aim is to find a frequentist 0.95 confidence interval for θ that utilizes this prior information. 相似文献
10.
For a loss distribution belonging to a location–scale family, Fμ,σ, the risk measures, Value-at-Risk and Expected Shortfall are linear functions of the parameters: μ+τσ where τ is the corresponding risk measure of the mean-zero and unit-variance member of the family. For each risk measure, we consider a natural estimator by replacing the unknown parameters μ and σ by the sample mean and (bias corrected) sample standard deviation, respectively. The large-sample parametric confidence intervals for the risk measures are derived, relying on the asymptotic joint distribution of the sample mean and sample standard deviation. Simulation studies with the Normal, Laplace and Gumbel families illustrate that the derived asymptotic confidence intervals for Value-at-Risk and Expected Shortfall outperform those of Bahadur (1966) and Brazauskas et al. (2008), respectively. The method can also be effectively applied to Log-location-scale families whose supports are positive reals; an illustrative example is given in the area of financial credit risk. 相似文献
11.
Consider the partially balanced one-way layout for comparing k treatments μi,1?i?k, with a control μ0. We propose a new test which is similar to the test statistics of Marcus [1976. The powers of some tests of the equality of normal means against an ordered alternative. Biometrika 63, 177–183]. By simulation we find that the proposed test has a good power performance when compared with other tests. Moreover, it can produce confidence intervals for μi-μ0,1?i?k. 相似文献
12.
We consider m×m covariance matrices, Σ1 and Σ2, which satisfy Σ2-Σ1=Δ, where Δ has a specified rank. Maximum likelihood estimators of Σ1 and Σ2 are obtained when sample covariance matrices having Wishart distributions are available and rank(Δ) is known. The likelihood ratio statistic for a test about the value of rank(Δ) is also given and some properties of its null distribution are obtained. The methods developed in this paper are illustrated through an example. 相似文献
13.
We consider density estimation for a smooth stationary process Xt, t∈R, based on a discrete sample Yi=XΔi, i=0,…,n=T/Δ. By a suitable interpolation scheme of order p , we augment data to form an approximation Xp,t, t∈[0,T], of the continuous-time process and base our density estimate on the augmented sample path. Our results show that this can improve the rate of convergence (measured in terms of n) of the density estimate. Among other things, this implies that recording n observations using a small Δ can be more efficient than recording n independent observations. 相似文献
14.
15.
For a random sample of size n from an absolutely continuous random vector (X,Y), let Yi:n be ith Y-order statistic and Y[j:n] be the Y-concomitant of Xj:n. We determine the joint pdf of Yi:n and Y[j:n] for all i,j=1 to n, and establish some symmetry properties of the joint distribution for symmetric populations. We discuss the uses of the joint distribution in the computation of moments and probabilities of various ranks for Y[j:n]. We also show how our results can be used to determine the expected cost of mismatch in broken bivariate samples and approximate the first two moments of the ratios of linear functions of Yi:n and Y[j:n]. For the bivariate normal case, we compute the expectations of the product of Yi:n and Y[i:n] for n=2 to 8 for selected values of the correlation coefficient and illustrate their uses. 相似文献
16.
Consider a mixture problem consisting of k classes. Suppose we observe an s-dimensional random vector X whose distribution is specified by the relations P(X∈A|Y=i)=Pi(A), where Y is an unobserved class identifier defined on {1,…,k}, having distribution P(Y=i)=pi. Assuming the distributions Pi having a common covariance matrix, elegant identities are presented that connect the matrix of Fisher information in Y on the parameters p1,…,pk, the matrix of linear information in X, and the Mahalanobis distances between the pairs of P 's. Since the parameters are not free, the information matrices are singular and the technique of generalized inverses is used. A matrix extension of the Mahalanobis distance and its invariant forms are introduced that are of interest in their own right. In terms of parameter estimation, the results provide an independent of the parameter upper bound for the loss of accuracy by esimating p1,…,pk from a sample of X′s, as compared with the ideal estimator based on a random sample of Y′s. 相似文献
17.
Consider the model where there are I independent multivariate normal treatment populations with p×1 mean vectors μi, i=1,…,I, and covariance matrix Σ. Independently the (I+1)st population corresponds to a control and it too is multivariate normal with mean vector μI+1 and covariance matrix Σ. Now consider the following two multiple testing problems. 相似文献
18.
Let X={Xn}n?1 be a nonstationary random field satisfying a long range weak dependence for each coordinate at a time and a local dependence condition that avoids clustering of exceedances of high values. For these random fields, the probability of no exceedances of high values can be approximated by exp(−τ), where τ is the limiting mean number of exceedances. We present a class of nonstationary normal random fields for which this result can be applied. 相似文献
19.
In this paper a new multivariate regression estimate is introduced. It is based on ideas derived in the context of wavelet estimates and is constructed by hard thresholding of estimates of coefficients of a series expansion of the regression function. Multivariate functions constructed analogously to the classical Haar wavelets are used for the series expansion. These functions are orthogonal in L2(μn), where μn denotes the empirical design measure. The construction can be considered as designing adapted Haar wavelets. 相似文献
20.
In this article, we extended the empirical distribution function based test statistic Ik of Skaug and Tjostheim [1993. Nonparametric test of serial independence based on the empirical distribution function. Biometrika 80, 591–602] in the time series setting to Dn for spatial lattice data and derived the asymptotic distribution of the proposed test statistic Dn under the null hypothesis of spatial independence. The size and power of the proposed test statistic under conditional autoregressive model (CAR) were simulated. We applied Dn, Moran's I and Geary's c to the transformed and well-studied sudden infant death syndrome data from North Carolina and found that Dn produced a much smaller p-value in testing spatial independence. 相似文献