A class of tests for testing ‘new is better than used’

A class of tests is proposed for testing H₀ F?(x) = e^?λx, λ > 0, x≥0 vs. H₁ F?(x + y) ≤ F?(x)F?(y), x, y≥0, with strict inequality for some x, y ≥ 0 (F = new is better than used). Efficiency comparisons of some tests within the class are made and a new test is proposed on the basis of these comparisons. Consistency and the asymptotic normality of the class of tests is proved under fairly broad conditions on the underlying entities.     

One-sided range test for testing against an ordered alternative under heteroscedasticity

In a one-way fixed effects analysis of variance model, when normal variances are unknown and possibly unequal, a one-sided range test for testing the null hypothesis H ₀ : μ ₁ = … = μ_k against an ordered alternative H_a : μ ₁ ≤ … ≤ μ_k by a single-stage and a two-stage procedure, respectively, is proposed. The critical values under H ₀ and the power under a specific alternative are calculated. Relation between the one-stage and the two-stage test procedures is discussed. A numerical example to illustrate these procedures is given.     

Tests based on equivariant estimators

Let ?⁽¹⁾ and ?⁽²⁾ be location-equivariant estimators of an unknown location parameter μ. It is shown that the test for H₀: μ ≤ μ₀ versus H_A : μ > μ₀ that rejects H₀ if ?⁽¹⁾ is large is uniformly more powerful than the one that rejects H₀ if ?⁽²⁾ is large if and only if ?⁽²⁾ is “more dispersed” than ?⁽¹⁾. A similar result is obtained for tests on scale using the star-shaped ordering. Examples are given.     

Nonparametric predictive inference for diagnostic test thresholds

Abstract

Measuring the accuracy of diagnostic tests is crucial in many application areas including medicine, machine learning and credit scoring. The receiver operating characteristic (ROC) curve and surface are useful tools to assess the ability of diagnostic tests to discriminate between ordered classes or groups. To define these diagnostic tests, selecting the optimal thresholds that maximize the accuracy of these tests is required. One procedure that is commonly used to find the optimal thresholds is by maximizing what is known as Youden’s index. This article presents nonparametric predictive inference (NPI) for selecting the optimal thresholds of a diagnostic test. NPI is a frequentist statistical method that is explicitly aimed at using few modeling assumptions, enabled through the use of lower and upper probabilities to quantify uncertainty. Based on multiple future observations, the NPI approach is presented for selecting the optimal thresholds for two-group and three-group scenarios. In addition, a pairwise approach has also been presented for the three-group scenario. The article ends with an example to illustrate the proposed methods and a simulation study of the predictive performance of the proposed methods along with some classical methods such as Youden index. The NPI-based methods show some interesting results that overcome some of the issues concerning the predictive performance of Youden’s index.     

Asymptotic properties of Bayes risk for one-sided tests

Consider a given sequence {T_n} of estimators for a real-valued parameter θ. This paper studies asymptotic properties of restricted Bayes tests of the following form: reject H₀:θ ≤ θ₀ in favour of the alternative θ > θ₀ if T_n ≤ C_n, where the critical point C_n is determined to minimize among all tests of this form the expected probability of error with respect to the prior distribution. Such tests may or may not be fully Bayes tests, and so are called T_n-Bayes. Under fairly broad conditions it is shown that and the T_n-Bayes risk where a_n is the order of the standard error of T_n, - is the prior density, and μ is the median of F, the limit distribution of (T_n – θ)/a_nb(θ). Several examples are given.     

Lois limites pour les statistiques d'ordre dans le cas non identiquement distribue

Let X₁, X₂,… be a sequence of independent random variables with distribution functions F₁, where 1 ≤ i ≤ n, and for each n ≥ 1 let X_1,n ≤… ≤ X_n,n denote the order statistics of the first n random variables. Under suitable hypotheses about the F₁, we characterize the limit distribution functions H(x) for which P(X_k,n ? a_nx + b_n) → H(x), where a_n > 0 and b_n are real constants. We consider the cases where κ = κ(n) satisfies √n {κ(n)/n — λ} → 0 and √n {κ(n)/n — λ} → ∞ separately.     

Testing the Equality and Exchangeability of Bivariate Distributions

We consider three methods (oments, cut-points, and ranks) for testing the hypotheses of equality of two bivariate distribution functions (H _0a ) and exchangeability (H _0b ). To test H _0a , the asymptotic normality of the vector of mixed moments provides a statistic with an asymptotic chi-square distribution. With every observation, method of cut-points associates three 2 × 2 tables to record the proportions of the X, Y, and the combined samples that fall in the four regions around the observation. We measure the total squared deviations of the proportions in the combined sample from X and Y samples. The two methods are compared with the method of ranks based on the Puri and Sen (1971 Puri , M. L. , Sen , P. K. ( 1971 ). Nonparametric Methods in Multivariate Analysis . New York : John Wiley and Sons . [Google Scholar]) multivariate two-sample rank test for location.

To test H _0b we identify two bivariate distributions, one above and the other below the line of symmetry X = Y, to which a test of H _0a is applied. Under H _0b , matrix of mixed moments is symmetric and a quadratic form in differences of (r,s)-th and (s, r)-th mixed moments provides an asymptotic chi-square distribution. A permutation test is devised to apply the method of cut-points to the observations above and below the line of symmetry after they are folded. We also describe an adaption of the Puri-Sen rank test to assess H _0b . To estimate the power of the above methods under different types of alternatives and compare them to existing tests, we report on a Monte Carlo experiment that evaluates the finite-sample performance of these methods under the Plackett's family of bivariate distributions.     

Tests for the mean of a multivariate normal population

We study the problem of testing: H₀ : μ ∈ P against H₁ : μ ? P, based on a random sample of N observations from a p-dimensional normal distribution N_p(μ, Σ) with Σ > 0 and P a closed convex positively homogeneous set. We develop the likelihood-ratio test (LRT) for this problem. We show that the union-intersection principle leads to a test equivalent to the LRT. It also gives a large class of tests which are shown to be admissible by Stein's theorem (1956). Finally, we give the α-level cutoff points for the LRT.     

Lower percentage points of hartley's extremal quotient statistic and their applications

Consider K(>2) independent populations π₁,..,π _k such that observations obtained from π _k are independent and normally distributed with unknown mean µ _i and unknown variance θ _i i = 1,…,k. In this paper, we provide lower percentage points of Hartley's extremal quotient statistic for testing an interval hypothesisH ₀ θ _[k] θ _[k] > δ vs. H _a : θ _[k] θ _[1] ≤ δ , where δ ≥ 1 is a predetermined constant and θ _[k](θ _[1]) is the max (min) of the θ_i,…,θ _k . The least favorable configuration (LFC) for the test under H ₀ is determined in order to obtain the lower percentage points. These percentage points can also be used to construct an upper confidence bound for θ_[k]/θ_[1].     

A Note on Dirichlet One-Armed Bandits

One of the two independent stochastic processes (or ‘arms’) is selected and observed sequentially at each of n(≤ ∝) stages. Arm 1 yields observations identically distributed with unknown probability measure P with a Dirichlet process prior whereas observations from arm 2 have known probability measure Q. Future observations are discounted and at stage m, the payoff is a _m(≥0) times the observation Z _m at that stage. The objective is to maximize the total expected payoff. Clayton and Berry (1985) consider this problem when a _m equals 1 for m ≤ n and 0 for m > n(< ∝) In this paper, the Clayton and Berry (1985) results are extended to the case of regular discount sequences of horizon n, which may also be infinite. The results are illustrated with numerical examples. In case of geometric discounting, the results apply to a bandit with many independent unknown Dirichlet arms.     

FIRST-EXIT TIMES FOR POISSON SHOT NOISE

For a shot-noise process X(t) with Poisson arrival times and exponentially diminishing shocks of i.i.d. sizes, we consider the first time T _b at which a given level b > 0 is exceeded. An integral equation for the joint density of T _b and X(T _b ) is derived and, for the case of exponential jumps, solved explicitly in terms of Laplace transforms (LTs). In the general case we determine the ordinary LT of the function ? P(T _b > t) in terms of certain LTs derived from the distribution function H(x; t) = P(X(t) ≤ x), considered as a function of both variables x and t. Moreover, for G(t, u) = P(T _b > t, X(t) < u), that is the joint distribution function of sup_{0 ≤ s ≤ t} X(s) and X(t), an integro-differential equation is presented, whose unique solution is G(t, u).     

Three optimal cut-point selection criteria based on sensitivity and specificity with user-defined weights

Methods: Based on the index S (S = SENSITIVITY (SEN) × SPECIFICITY (SPE)), the new weighted product index S_w is defined as S_w = (SEN)^2w × (SPE)^2(1-w), where (0≤w≤1). The S_w is developed to be a new tool to select the optimal cut point in ROC analysis and be compared with the other two commonly used criteria.

Results: Comparing the optimal cut point for the three criteria, the wave range of the optimal cut point for the maximized weighted Youden index criterion is the widest, the weighted closest-to-(0,1) criterion is the narrowest and the weighted product index S_w criterion lays between the ranges of the two criteria.     

Tests for the umbrella alternative under normality

Consider a one-way layout of the analysis of variance assuming independence, normality, and homogeneity of variance. Test the null hypothesis H_o that the means, j., of each of Amp; columns, i = 1,…, k are equal versus the alternative that they follow an umbrella pattern. That is, the alternative is H₁-H₀ where H₁: μ₁> μ₂>… > μk, and m is known. We derive a class of tests which are unbiased and lie in a nontrivial complete class. We recommend specific tests within the class. A simulation of the power functions of some tests is contrasted with the simulated power function of the likelihood ratio test.     

Csiszar’s ϕ-divergences for testing the order in a Markov chain

x ₁, ..., x _n+r can be treated as the sample values of a Markov chain of order r or less (chain in which the dependence extends over r+1 consecutive variables only), and consider the problem of testing the hypothesis H ₀ that a chain of order r− 1 will be sufficient on the basis of the tools given by the Statistical Information Theory: _ϕ-Divergences. More precisely, if p _{a 1 ....., a r}: _{a r +1} denotes the transition probability for a r ^th order Markov chain, the hypothesis to be tested is H ₀:p _{a 1 ....., a r}: _{a r +1} = p _{a 2 ....., a r}: _{a r +1}, a _i ∈{1, ..., s}, i = 1, ..., r + 1 The tests given in this paper, for the first time, will have as a particular case the likelihood ratio test and the test based on the chi-squared statistic. Received: August 3, 1998; revised version: November 25, 1999     

BAHADUR SLOPES FOR COMPARISONS OF TESTS BASED ON RESAMPLING

Let X₁,…, X_n be random variables symmetric about θ from a common unknown distribution F_θ(x) =F(x–θ). To test the null hypothesis H₀:θ= 0 against the alternative H₁:θ > 0, permutation tests can be used at the cost of computational difficulties. This paper investigates alternative tests that are computationally simpler, notably some bootstrap tests which are compared with permutation tests. Of these the symmetrical bootstrap-f test competes very favourably with the permutation test in terms of Bahadur asymptotic efficiency, so it is a very attractive alternative.     

Robustness of the Multiple Correlation Coefficient When Sampling From a Mixture of Two Multivariate Normal Populations

The density of the multiple correlation coefficient is derived by direct integration when the sample covariance matrix has a linear non-central distribution. Using the density, we deduce the null and non-null distribution of the multiple correlation coefficient when sampling from a mixture of two multivariate normal populations with the same covariance matrix. We also compute actual significance levels of the test of the hypothesis H_o : ρ_1·2…p = 0 versus H_a:ρ_1·2…p > 0, given the mixture model.     

Weighted product index and its two-independent-sample comparison based on weighted sensitivity and specificity

Background: On the basis of statistical methods about index S (S = SEN × SPE), we develop a new weighted ways (weighted product index S_w) of combining sensitivity and specificity with user-defined weights. Methods: The new weighted product index S_w is defined as S_w = (SEN) (Youden 1950)^2w × (SPE) (Youden 1950) ^2(1?w) Results: For the large sample, the test statistics Z of two-independent-sample weighted product indices can either be a monotonous increasing/decreasing function or a no-monotonous function of weight w. Type I error of this statistics can be guaranteed close to the nominal level of 5%, which is more conservative than the weighted Youden index from simulation.