On Mixture Failure Rates Ordering

Mixtures of increasing failure rate distributions (IFR) can decrease at least in some intervals of time. Usually, this property can be observed asymptotically as t → ∞. This is due to the fact that the mixture failure rate is “bent down” compared with the corresponding unconditional expectation of the baseline failure rate, which was proved previously for some specific cases. We generalize this result and discuss the “weakest populations are dying first” property, which leads to the change in the failure rate shape. We also consider the problem of mixture failure rate ordering for the ordered mixing distributions. Two types of stochastic ordering are analyzed: ordering in the likelihood ratio sense and ordering in variances when the means are equal.     

On Uniformly Optimal Networks: A Reversal of Fortune?

In this article, the general problem of comparing the performance of two communication networks is examined. The standard approach, using stochastic ordering as a metric, is reviewed, as are the mixed results on the existence of uniformly optimal networks (UONs) which have emerged from this approach. While UONs have been shown to exist for certain classes of networks, it has also been shown that no UON network exists for other classes. Results to date beg the question: Is the problem of identifying a Uniformly Optimal Network (UON) of a given size dead or alive? We reframe the investigation into UONs in terms of network signatures and the alternative metric of stochastic precedence. While the endeavor has been dead, or at least dormant, for some 20 years, the findings in the present article suggest that the question above is by no means settled. Specifically, we examine a class of networks of a particular size for which it was shown that no individual network was uniformly optimal relative to the standard metric (the uniform ordering of reliability polynomials), and we show, using the aforementioned alternative metric, that this class is totally ordered and that a uniformly optimal network exists after all. Optimality with respect to “performance per unit cost” type metrics is also discussed.     

Extensions of empirical likelihood and chi-squared-based tests for ordered alternatives

Several methods for comparing k populations have been proposed in the literature. These methods assess the same null hypothesis of equal distributions but differ in the alternative hypothesis they consider. We focus on two important alternative hypotheses: monotone and umbrella ordering. Two new families of test statistics are proposed, including two known tests, as well as two new powerful tests under monotone ordering. Furthermore, these families are adapted for testing umbrella ordering. We compare some members of the families with respect to power and Type I errors under different simulation scenarios. Finally, the methods are illustrated in several applications to real data.     

An Automatic Test for the Umbrella Alternatives

The paper proposes a new test for detecting the umbrella pattern under a general non‐parametric scheme. The alternative asserts that the umbrella ordering holds while the hypothesis is its complement. The main focus is put on controlling the power function of the test outside the alternative. As a result, the asymptotic error of the first kind of the constructed solution is smaller than or equal to the fixed significance level α on the whole set where the umbrella ordering does not hold. Also, under finite sample sizes, this error is controlled to a satisfactory extent. A simulation study shows, among other things, that the new test improves upon the solution widely recommended in the literature of the subject. A routine, written in R, is attached as the Supporting Information file.     

The Maximum Probability 2 × c Contingency Tables and the Maximum Probability Points of the Multivariate Hypergeometric Distribution

Abstract

The problem of obtaining the maximum probability 2 × c contingency table with fixed marginal sums, R  = (R ₁, R ₂) and C  = (C ₁, … , C _c ), and row and column independence is equivalent to the problem of obtaining the maximum probability points (mode) of the multivariate hypergeometric distribution MH(R ₁; C ₁, … , C _c ). The most simple and general method for these problems is Joe's (Joe, H. (1988 Joe, H. 1988. Extreme probabilities for contingency tables under row and column independence with application to Fisher's exact test. Commun. Statist. Theory Meth., 17(11): 3677–3685. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Extreme probabilities for contingency tables under row and column independence with application to Fisher's exact test. Commun. Statist. Theory Meth. 17(11):3677–3685.) In this article we study a family of MH's in which a connection relationship is defined between its elements. Based on this family and on a characterization of the mode described in Requena and Martín (Requena, F., Martín, N. (2000 Requena, F. and Martín, N. 2000. Characterization of maximum probability points in the multivariate hypergeometric distribution. Statist. Probab. Lett., 50: 39–47.  [Google Scholar]). Characterization of maximum probability points in the multivariate hypergeometric distribution. Statist. Probab. Lett. 50:39–47.), we develop a new method for the above problems, which is completely general, non recursive, very simple in practice and more efficient than the Joe's method. Also, under weak conditions (which almost always hold), the proposed method provides a simple explicit solution to these problems. In addition, the well-known expression for the mode of a hypergeometric distribution is just a particular case of the method in this article.     

A Revised Version of McNemar's Test for Paired Binary Data

In the formula of the McNemar test, a test on 2 × 2 classification tables with pairs of data, only the two categories A and D, which represent changes, are included; the retained parts B and C, which represent concordant responses, are not considered. Generally, it would be more reasonable for the significance of the changes to depend not only on A and D, but also on B and C, or on the sample size, n. To develop the test, two formulae, based on A, D, and n, and on A, D, B, C, and n, respectively, are proposed.     

Modification of Pickands' Dependence Function for Ordered Bivariate Extreme Distribution

We study the characteristics of the Pickands' dependence function for bivariate extreme distribution for minima, BEVM, when considering the stochastics ordering of the two variables, X < Y. The existing Pickand's dependence function terminologies and theories are modified to suit the dependence functions of extreme minimum cases. The main result is the introduction of the restricted logistic dependence function, A _RL , and the restricted exponential function, V _RL (x, y).     

A permutation test for umbrella alternatives

There is a wide variety of stochastic ordering problems where K groups (typically ordered with respect to time) are observed along with a (continuous) response. The interest of the study may be on finding the change-point group, i.e. the group where an inversion of trend of the variable under study is observed. A change point is not merely a maximum (or a minimum) of the time-series function, but a further requirement is that the trend of the time-series is monotonically increasing before that point, and monotonically decreasing afterwards. A suitable solution can be provided within a conditional approach, i.e. by considering some suitable nonparametric combination of dependent tests for simple stochastic ordering problems. The proposed procedure is very flexible and can be extended to trend and/or repeated measure problems. Some comparisons through simulations and examples with the well known Mack & Wolfe test for umbrella alternative and with Page’s test for trend problems with correlated data are investigated.     

Small Sample Confidence Intervals for the Odds Ratio

Abstract

Numerous methods—based on exact and asymptotic distributions—can be used to obtain confidence intervals for the odds ratio in 2 × 2 tables. We examine ten methods for generating these intervals based on coverage probability, closeness of coverage probability to target, and length of confidence intervals. Based on these criteria, Cornfield’s method, without the continuity correction, performed the best of the methods examined here. A drawback to use of this method is the significant possibility that the attained coverage probability will not meet the nominal confidence level. Use of a mid-P value greatly improves methods based on the “exact” distribution. When combined with the Wilson rule for selection of a rejection set, the resulting method is a procedure that performed very well. Crow’s method, with use of a mid-P, performed well, although it was only a slight improvement over the Wilson mid-P method. Its cumbersome calculations preclude its general acceptance. Woolf's (logit) method—with the Haldane–Anscombe correction— performed well, especially with regard to length of confidence intervals, and is recommended based on ease of computation.     

The one-way analysis of variance with ordered alternatives: A modification of bartholomew's Ē2 Test

Tests of homogeneity of normal means with the alternative restricted by an ordering on the means are considered. The simply ordered case, μ₁ ≤ μ₂ ≤ ··· ≤ μ_k, and the simple tree ordering, μ₁ ≤ μ_j, for; j= 2, 3,…, k, are emphasized. A modification of the likelihood-ratio test is proposed which is asymptotically equivalent to it but is more robust to violations of the hypothesized orderings. The new test has power at the points satisfying the hypothesized ordering which is similar to that of the likelihood-ratio test provided the degrees of freedom are not too small. The modified test is shown to be unbiased and consistent.     

Superiority of the r–k Class Estimator Over Some Estimators In A Linear Model

In regression analysis, to overcome the problem of multicollinearity, the r ? k class estimator is proposed as an alternative to the ordinary least squares estimator which is a general estimator including the ordinary ridge regression estimator, the principal components regression estimator and the ordinary least squares estimator. In this article, we derive the necessary and sufficient conditions for the superiority of the r ? k class estimator over each of these estimators under the Mahalanobis loss function by the average loss criterion. Then, we compare these estimators with each other using the same criterion. Also, we suggest to test to verify if these conditions are indeed satisfied. Finally, a numerical example and a Monte Carlo simulation are done to illustrate the theoretical results.     

A STOCHASTIC ORDERING FOR RANDOM VARIABLES WITH APPLICATIONS

A new univariate stochastic ordering is introduced. Some characterization results for such an ordering are stated. It is proved that the ordering is an integral stochastic ordering, obtaining a maximal generator. By means of this generator, the main properties of the ordering are deduced. A method for introducing univariate stochastic orderings, suggested by the new ordering, is analysed. Relationships with other stochastic orderings are also developed. To conclude, an example of an application of the new ordering to the field of medicine is proposed.     

Exact inference for multivariate ordered alternatives

In this paper we discuss three types of ordered alternatives ordered location, stochastic ordering and quadrant dependence. We prove that quadrant dependence is the more general among the three. Then we consider a conditional tests for the equality of c distributions against quadrant dependence in a multivariate setup. An exact simultaneous testing procedure based on dependent conditional tests is presented. Two applications to real data are also given.     

A note on the mean residual life of a coherent system with a cold standby component

In this paper, we investigate the effect of a cold standby component on the mean residual life (MRL) of a system. When the system fails, a cold standby component is immediately put in operation. We particularly focus on the coherent systems in which, after putting the standby component into operation, the failure of the system is due to the next component failure. For these systems, we define MRL functions and obtain their explicit expressions. Also some stochastic ordering results are provided. Such systems include k-out-of-n systems. Hence, our results extend some results in literature.     

Cumulative or adjacent logits: Which choice for an ordinal logistic latent variable model?

With ordinal response items, a graded response model (GRM) is of cumulative logits type, while the polytomous Rasch model (PRM) is based on adjacent logits. In this work, we compare the two approaches. We show that the PRM is superior to the GRM, with interesting properties that we prove. Note S_ν the sum of item responses of individual ν and Θ_ν its latent parameter; we show i) S_ν is a sufficient statistic for θ_ν and ii) a property of “stochastic ordering” of the conditional distributions G_θ/S. The second property, less known, is, to our knowledge, nowhere satisfactorily demonstrated.     

On Upshifted Reversed Mean Residual Life Order

Recently, the concept of reversed mean residual life order based on the mean of the random variable X _t  = (t ? X | X ≤ t), t > 0, called the reversed residual life, defined for the nonnegative random variable X, has been introduced in the literature. In this paper, a stochastic order based on the shifted version of the reversed mean residual life is proposed, based on the reversed mean residual life function for a random variable X with support (l _X , ∞), where l _X may be negative infinity, and its properties are studied. Closure under the Poisson shock model and properties for spare allocation are also discussed.     

DISTRIBUTION-FREE CONFIDENCE PROCEDURE FOR UMBRELLA ORDERINGS

Responses in a one-factor experiment with A; ordered treatments follow an umbrella, ordering if they consist of two piecewise monotone segments, i.e. increasing and then decreasing, or the converse. This paper proposes a non-parametric distribution-free confidence procedure for umbrella orderings, the aim being to identify the treatments that correspond to the optimal effects. It uses a method that joins the seemingly unrelated theories of U-statistics and isotonic regression. A random confidence subset of the ordered treatments is constructed, such that it contains all the unknown peaks (optimal treatments) of an umbrella ordering with any prespecified confidence level. The paper demonstrates that the proposed confidence procedure is nonparametric distribution-free and, further, that the proposed procedure naturally implies a test for umbrella alternatives. Since the proposed confidence procedure is always more informative than tests for umbrella alternatives, it should be used in their place in practice. An example illustrates the proposed procedure.     

An Integral Bhattacharyya Type Bound for the Bayes Risk

In this article, some results are derived on stochastic comparisons of the residual and past lifetimes of an (n ? k + 1)-out-of-n system with dependent components. These findings generalize some recent results obtained on systems with independent components and provide some interesting results for a system with dependent components following an Archimedean copula. An illustrative example is provided.     

Conditional inference under simultaneous stochastic ordering constraints

Testing for stochastic ordering is of considerable importance when increasing does of a treatment are being compared, but in applications involving multivariate responses has received much less attention. We propose a permutation test for testing against multivariate stochastic ordering. This test is distribution-free and no assumption is made about the dependence relations among variables. A comparative simulation study shows that the proposed solution exhibits a good overall performance when compared with existing tests that can be used for the same problem.     

Inference for Adaptive Time Series Models: Stochastic Volatility and Conditionally Gaussian State Space Form

In this paper we model the Gaussian errors in the standard Gaussian linear state space model as stochastic volatility processes. We show that conventional MCMC algorithms for this class of models are ineffective, but that the problem can be alleviated by reparameterizing the model. Instead of sampling the unobserved variance series directly, we sample in the space of the disturbances, which proves to lower correlation in the sampler and thus increases the quality of the Markov chain.

Using our reparameterized MCMC sampler, it is possible to estimate an unobserved factor model for exchange rates between a group of n countries. The underlying n + 1 country-specific currency strength factors and the n + 1 currency volatility factors can be extracted using the new methodology. With the factors, a more detailed image of the events around the 1992 EMS crisis is obtained.

We assess the fit of competitive models on the panels of exchange rates with an effective particle filter and find that indeed the factor model is strongly preferred by the data.