A one-sided test for testino homogeneity of scale parameters against simple ordered alternative

In an earlier paper the authors (1997) extended the results of Hayter (1990) to the two parameter exponential probability model. This paper addressee the extention to the scale parameter case under location-scale probability model. Consider k (k≧3) treatments or competing firms such that an observation from with treatment or firm follows a distribution with cumulative distribution function (cdf) F_i(x)=F[(x-μ_i)/Q_i], where F(·) is any absolutely continuous cdf, i=1,…,k. We propose a test to test the null hypothesis H₀:θ₁=…=θ_k against the simple ordered alternative H₁:θ₁≦…≦θ_k, with at least one strict inequality, using the data X_i,j, i=1,…k; j=1,…,n₁. Two methods to compute the critical points of the proposed test have been demonstrated by talking k two parameter exponential distributions. The test procedure also allows us to construct simultaneous one sided confidence intervals (SOCIs) for the ordered pairwise ratios θ_j/θ_i, 1≦i<j≦k. Statistical simulation revealed that: 9i) actual sizes of the critical points are almost conservative and (ii) power of the proposed test relative to some existing tests is higher.     

Confidence intervals for the largest mean from k correlated normal populations

Let X= (X₁,…, X_k)’ be a k-variate (k ≥ 2) normal random vector with unknown population mean vector μ = (μ₁ ,…, μ_k)’ and covariance matrix Σ of order k and let μ_[1] ≤ … ≤ μ_[k] be the ordered values of the μ ’ s. No prior knowledge of the pairing of the μ_[i] with the X_j. (or μ_[i] with the σ_j ²) is assumed for any i and j (1 ≤ i, j ≤ k). Based on a random sample of N independent vector observations on X, this paper considers both upper and lower (one-sided) and two-sided 100γ% (0 < γ < 1) confidence intervals for μ_[k] and μ_[1], the largest and the smallest mean, respectively, when Σ is known and when Σ is equal to σ²R with common unknown variance σ² > 0 and correlation matrix R known, respectively. An optimum two-sided confidence interval via finding the shortest length from this class is also considered. Necessary tables and computer program to actually apply these procedures are provided.     

A class of fixed-width confidence intervals for a ranked normal mean

Suppose that we are given k(≥ 2) independent and normally distributed populations π₁, …, π_k where π_i has unknown mean μ_i and unknown variance σ² _i (i = 1, …, k). Let μ_[i] (i = 1, …, k) denote the i^th smallest one of μ₁, …, μ_k. A two-stage procedure is used to construct lower and upper confidence intervals for μ_[i] and then use these to obtain a class of two-sided confidence intervals on μ_[i] with fixed width. For i = k, the interval given by Chen and Dudewicz (1976) is a special case. Comparison is made between the class of two-sided intervals and a symmetric interval proposed by Chen and Dudewicz (1976) for the largest mean, and it is found that for large values of k at least one of the former intervals requires a smaller total sample size. The tables needed to actually apply the procedure are provided.     

Estimation of the location parameter and the average worth of the selected subset of two parameter exponential populations under LINEX loss function

Suppose a subset of populations is selected from k exponential populations with unknown location parameters θ₁, θ₂, …, θ_k and common known scale parameter σ. We consider the estimation of the location parameter of the selected population and the average worth of the selected subset under an asymmetric LINEX loss function. We show that the natural estimator of these parameters is biased and find the uniformly minimum risk-unbiased (UMRU) estimator of these parameters. In the case of k = 2, we find the minimax estimator of the location parameter of the smallest selected population. Furthermore, we compare numerically the risk of UMRU, minimax, and the natural estimators.     

Subset selection based on likelihood from uniform and related popupulations

Let π_i (i=1,2,…, k) be charceterized by the uniform distribution on (a_i;b_i), where exactly one of a_i and b_i is unknown. With unequal sample sizes, suppose that from the k (>=2) given populations, we wish to select a random-size subset containing the one with the smllest value of θ_i= b_i - a_i. RuleR_i selects π if a likelihood-based k-dimensional confidence region for the unknown (θ₁,… θ_k) contains at least one point having θ_i as its smallest component. A second rule, R , is derived through a likelihood ratio and turns out to be that of Barr and prabhu whenthe sample sizes are equal. Numerical comparisons are made. The results apply to the larger class of densities g ( z ; θ_i) =M(z)Q(θ_i) if a(θ_i) < z <b(θ_i). Extensions to the cases when both a_i and b_i are unknown and when θ_j isof interest are indicated. 1<=j<=k     

Estimation of the Scale Parameter of the Selected Gamma Population Under the Entropy Loss Function

Let X ₁, X ₂,…, X _k be k (≥2) independent random variables from gamma populations Π₁, Π₂,…, Π _k with common known shape parameter α and unknown scale parameter θ _i , i = 1,2,…,k, respectively. Let X _(i) denotes the ith order statistics of X ₁,X ₂,…,X _k . Suppose the population corresponding to largest X _(k) (or the smallest X ₍₁₎) observation is selected. We consider the problem of estimating the scale parameter θ _M (or θ _J ) of the selected population under the entropy loss function. For k ≥ 2, we obtain the Unique Minimum Risk Unbiased (UMRU) estimator of θ _M (and θ _J ). For k = 2, we derive the class of all linear admissible estimators of the form cX ₍₂₎ (and cX ₍₁₎) and show that the UMRU estimator of θ _M is inadmissible. The results are extended to some subclass of exponential family.     

Some subset selection problems

Consider that we have a collection of k populations π₁, π₂…,π_k. The quality of the ith population is characterized by a real parameter θ_i and the population is to be designated as superior or inferior depending on how much the θ_i differs from θ_max = max{θ₁, θ₂,…,θ_k}. From the set {π₁, π₂,…,π_k}, we wish to select the subset of superior populations. In this paper we devise rules of selection which have the property that their selected set excludes all the inferior populations with probability at least 1?α, where a is a specified number.     

Estimating the parameter of selected uniform population under the squared log error loss function

Let π₁, …, π_k be k (? 2) independent populations, where π_i denotes the uniform distribution over the interval (0, θ_i) and θ_i > 0 (i = 1, …, k) is an unknown scale parameter. The population associated with the largest scale parameter is called the best population. For selecting the best population, We use a selection rule based on the natural estimators of θ_i, i = 1, …, k, for the case of unequal sample sizes. Consider the problem of estimating the scale parameter θ_L of the selected uniform population when sample sizes are unequal and the loss is measured by the squared log error (SLE) loss function. We derive the uniformly minimum risk unbiased (UMRU) estimator of θ_L under the SLE loss function and two natural estimators of θ_L are also studied. For k = 2, we derive a sufficient condition for inadmissibility of an estimator of θ_L. Using these condition, we conclude that the UMRU estimator and natural estimator are inadmissible. Finally, the risk functions of various competing estimators of θ_L are compared through simulation.     

ON ASYMPTOTIC OPTIMALITY OF SOME TESTS FOR EXPONENTIAL DISTRIBUTIONS

The probability density function (pdf) of a two parameter exponential distribution is given by f(x; p, s?) =s?^-1 exp {-(x - ρ)/s?} for x≥ρ and 0 elsewhere, where 0 < ρ < ∞ and 0 < s?∞. Suppose we have k independent random samples where the ith sample is drawn from the ith population having the pdf f(x; ρ_i, s?_i), 0 < ρ_i < ∞, 0 < s?_i < s?_i < and f(x; ρ, s?) is as given above. Let X_i1 < X_i2 <… < X_iri denote the first r_i order statistics in a random sample of size n_i, drawn from the ith population with pdf f(x; ρ_i, s?_i), i = 1, 2,…, k. In this paper we show that the well known tests of hypotheses about the parameters ρ_i, s?_i, i = 1, 2,…, k based on the above observations are asymptotically optimal in the sense of Bahadur efficiency. Our results are similar to those for normal distributions.     

SELECTION PROCEDURES FOR SCALE PARAMETERS USING TWO-SAMPLE U-STATISTICS

Let be k independent populations having the same known quantile of order p (0 p 1) and let F(x)=F(x/_i) be the absolutely continuous cumulative distribution function of the ith population indexed by the scale parameter ₁, i = 1,…, k. We propose subset selection procedures based on two-sample U-statistics for selecting a subset of k populations containing the one associated with the smallest scale parameter. These procedures are compared with the subset selection procedures based on two-sample linear rank statistics given by Gill & Mehta (1989) in the sense of Pitman asymptotic relative efficiency, with interesting results.     

A Primer on Visualizations for Comparing Populations,Including the Issue of Overlapping Confidence Intervals

In comparing a collection of K populations, it is common practice to display in one visualization confidence intervals for the corresponding population parameters θ₁, θ₂, …, θ_K. For a pair of confidence intervals that do (or do not) overlap, viewers of the visualization are cognitively compelled to declare that there is not (or there is) a statistically significant difference between the two corresponding population parameters. It is generally well known that the method of examining overlap of pairs of confidence intervals should not be used for formal hypothesis testing. However, use of a single visualization with overlapping and nonoverlapping confidence intervals leads many to draw such conclusions, despite the best efforts of statisticians toward preventing users from reaching such conclusions. In this article, we summarize some alternative visualizations from the literature that can be used to properly test equality between a pair of population parameters. We recommend that these visualizations be used with caution to avoid incorrect statistical inference. The methods presented require only that we have K sample estimates and their associated standard errors. We also assume that the sample estimators are independent, unbiased, and normally distributed.     

Optimal transitive procedures for comparing large numbers of parameters

Independent random samples are selected from each of a set of N independent populations, P₁,…,P_n. Interest centers around comparing N (unknown) scalar parameters θ₁,…,θ_N associated respectively with the N populations P₁,…,P_n. Procedures are constructed for estimating the magnitude of each of the differences δ_t,j = θ_i ? θ_j (1 ≤ i,j ≤ N) between pairs of populations. A loss function which adopts appropriate penalties for magnitude errors in estimation of differences is constructed. Magnitude estimators of differences are called transitive if they give rise to a transitive (i.e., consistent) relationship between pairwise differences of parameters. We show how to construct optimal effcient transitive magnitude–estimation procedures and demonstrate their usefulness through an example involving estimating the magnitude of the differences between disease incidence in paired towns for different pairs. Optimal transitive pairwise–comparison procedures are optimum (i.e., have the smallest posterior Bayes risks) in the class of all transitive pairwise–comparison procedures; as such they replace classical Bayes procedures which are usually not transitive when the number N of parameters compared is large. The posterior Bayes risk of optimal transitive pairwise comparison procedures are compared with that for alternative ‘adapted’ procedures, constructed from optimal simultaneous estimators and adapted for the purpose of pairwise comparisons. It is shown that the optimal transitive pairwise comparison procedures dominated the adapted procedures (in posterior Bayes risk) and typically represent only a small increase in posterior risk over the classical Bayes procedures which generally fail to be consistent. Optimal Bayes procedures are shown, for large numbers of parameters to be reasonably easy to construct using the algorithms outlined in this paper     

A RENEWAL THEOREM IN MULTIDIMENSIONAL TIME

Let Y_l, Y₂,… be i.i.d., positive, integer-valued random variables with means, μ. Let the sequences {Y_ij, j= 1,2,…}, i= 1,…, r be independent copies of {Y₁, Y₂,…}. For n={n₁,…, n_r.}, n₁≥1, let S_n=S?ⁿ¹_k1=1= 1 …S?^nr_kr=1 Y_ik1… Y_rkr. We show that S?^N_k=1S?^∞_k1=1…S?^∞_nr=1 P[[S_n= k] ? [μ^-r N log^r-1 (N)/(r-1)!] as N →∞.     

Simultaneous selection of extreme populations:a bayesian approach

The problem of simultaneously selecting two non-empty subsets, S_Land S_U, of k populations which contain the lower extreme population (LEP) and the upper extreme population (UEP), respectively, is considered. Unknown parameters θ₁,…,θ_kcharacterize the populations π₁,…,π_kand the populations associated with θ_[1]=min θ_i. and θ_[k]= max θ_i. are called the LEP and the UEP, respectively. It is assumed that the underlying distributions possess the monotone likelihood ratio property and that the prior distribution of θ= (θ₁,…,θ_k) is exchangeable. The Bayes rule with respect to a general loss function is obtained. Bayes rule with respect to a semi-additive and non-negative loss function is also determined and it is shown that it is minimax and admissible. When the selected subsets are required to be disjoint, it shown that the Bayes rule with respect to a specific loss function can be obtained by comparing certain computable integrals, Application to normal distributions with unknown means θ₁,…,θ_kand a common known variance is also considered.     

DISTRIBUTION FREE MULTIPLE COMPARISONS FOR MONOTONICALLY ORDERED TREATMENT EFFECTS

A method of calculating simultaneous one-sided confidence intervals for all ordered pairwise differences of the treatment effects_j–_i, 1 i < j k, in a one-way model without any distributional assumptions is discussed. When it is known a priori that the treatment effects satisfy the simple ordering_1k, these simultaneous confidence intervals offer the experimenter a simple way of determining which treatment effects may be declared to be unequal, and is more powerful than the usual two-sided Steel-Dwass procedure. Some exact critical points required by the confidence intervals are presented for k= 3 and small sample sizes, and other methods of critical point determination such as asymptotic approximation and simulation are discussed.     

Non-parametric estimation of an acceleration parameter

Let X₁,… X_m be a random sample of m failure times under normal conditions with the underlying distribution F(x) and Y₁,…,Y_n a random sample of n failure times under accelerated condititons with underlying distribution G(x);G(x)=1?[1?F(x)]^θ with θ being the unknown parameter under study.Define:U_ij=1 otherwise.The joint distribution of _ijdoes not involve the distribution F and thus can be used to estimate the acceleration parameter θ.The second approach for estimating θ is to use the ranks of the Y-observations in the combined X- and Y-samples.In this paper we establish that the rank of the Y-observations in the pooled sample form a sufficient statistic for the information contained in the U_ii 's about the parameter θ and that there does not exist an unbiassed estimator for the parameter θ.We also construct several estimators and confidence interavals for the parameter θ.     

Estimation After Selection Under Reflected Normal Loss Function

Let X ₁ and X ₂ be two independent random variables from normal populations Π₁, Π₂ with different unknown location parameters θ₁ and θ₂, respectively and common known scale parameter σ. Let X ₍₂₎ = max (X ₁, X ₂) and X ₍₁₎ = min (X ₁, X ₂). We consider the problem of estimating the location parameter θ _M (or θ _J ) of the selected population under the reflected normal loss function. We obtain minimax estimators of θ _M and θ _J . Also, we provide sufficient conditions for the inadmissibility of invariant estimators of θ _M and θ _J .     

Testing retrospective hypotheses

This paper develops a conditional approach to testing hypotheses set up after viewing the data. For example, suppose X_i are estimates of location parameters θ_i, i = 1,…n. We show how to compute p-values for testing whether θ₁ is one of the three largest θ_i after observing that X₁ is one of the three largest X_i, or for testing whether θ₁ > θ₂ > … > θ_n after observing X₁ >X₂> … >X_n.