Generalized p-values and generalized confidence regions for the multivariate Behrens-Fisher problem and MANOVA

For two multivariate normal populations with unequal covariance matrices, a procedure is developed for testing the equality of the mean vectors based on the concept of generalized p-values. The generalized p-values we have developed are functions of the sufficient statistics. The computation of the generalized p-values is discussed and illustrated with an example. Numerical results show that one of our generalized p-value test has a type I error probability not exceeding the nominal level. A formula involving only a finite number of chi-square random variables is provided for computing this generalized p-value. The formula is useful in a Bayesian solution as well. The problem of constructing a confidence region for the difference between the mean vectors is also addressed using the concept of generalized confidence regions. Finally, using the generalized p-value approach, a solution is developed for the heteroscedastic MANOVA problem.     

Reliability and expectation bounds for coherent systems with exchangeable components

Sharp upper and lower bounds are obtained for the reliability functions and the expectations of lifetimes of coherent systems based on dependent exchangeable absolutely continuous components with a given marginal distribution function, by use of the concept of Samaniego's signature. We first show that the distribution of any coherent system based on exchangeable components with absolutely continuous joint distribution is a convex combination of distributions of order statistics (equivalent to the k-out-of-n systems) with the weights identical with the values of the Samaniego signature of the system. This extends the Samaniego representation valid for the case of independent and identically distributed components. Combining the representation with optimal bounds on linear combinations of distribution functions of order statistics from dependent identically distributed samples, we derive the corresponding reliability and expectation bounds, dependent on the signature of the system and marginal distribution of dependent components. We also present the sequences of exchangeable absolutely continuous joint distributions of components which attain the bounds in limit. As an application, we obtain the reliability bounds for all the coherent systems with three and four exchangeable components, expressed in terms of the parent marginal reliability function and specify the respective expectation bounds for exchangeable exponential components, comparing them with the lifetime expectations of systems with independent and identically distributed exponential components.     

Multivariate dependence of spacings of generalized order statistics

Multivariate dependence of spacings of generalized order statistics is studied. It is shown that spacings of generalized order statistics from DFR (IFR) distributions have the CIS (CDS) property. By restricting the choice of the model parameters and strengthening the assumptions on the underlying distribution, stronger dependence relations are established. For instance, if the model parameters are decreasingly ordered and the underlying distribution has a log-convex decreasing (log-concave) hazard rate, then the spacings satisfy the MTP₂ (S- MRR₂) property. Some consequences of the results are given. In particular, conditions for non-negativity of the best linear unbiased estimator of the scale parameter in a location-scale family are obtained. By applying a result for dual generalized order statistics, we show that in the particular situation of usual order statistics the assumptions can be weakened.     

On weighted interval entropy

Measures of uncertainty in past and residual lifetime distributions have been proposed in the information-theoretic literature. Recently, Di Crescenzo and Longobardi (2006) introduced weighted differential entropy and its dynamic versions. These information-theoretic uncertainty measures are shift-dependent. In this paper, we study the weighted differential information measure for two-sided truncated random variables. This new measure is a generalization of recent dynamic weighted entropy measures. We study various properties of this measure, including its connection with weighted residual and past entropies, and we obtain its upper and lower bounds.     

A unified approach to testing for and against a set of linear inequality constraints in the product multinomial setting

A problem that is frequently encountered in statistics concerns testing for equality of multiple probability vectors corresponding to independent multinomials against an alternative they are not equal. In applications where an assumption of some type of stochastic ordering is reasonable, it is desirable to test for equality against this more restrictive alternative. Similar problems have been considered heretofore using the likelihood ratio approach. This paper aims to generalize the existing results and provide a unified technique for testing for and against a set of linear inequality constraints placed upon on any probability vectors corresponding to r independent multinomials. The paper shows how to compute the maximum likelihood estimates under all hypotheses of interest and obtains the limiting distributions of the likelihood ratio test statistics. These limiting distributions are of chi bar square type and the expression of the weighting values is given. To illustrate our theoretical results, we use a real life data set to test against second-order stochastic ordering.     

On the dependence structure of order statistics

Given a random sample from a continuous variable, it is observed that the copula linking any pair of order statistics is independent of the parent distribution. To compare the degree of association between two such pairs of ordered random variables, a notion of relative monotone regression dependence (or stochastic increasingness) is considered. Using this concept, it is proved that for i<j, the dependence of the jth order statistic on the ith order statistic decreases as i and j draw apart. This extends earlier results of Tukey (Ann. Math. Statist. 29 (1958) 588) and Kim and David (J. Statist. Plann. Inference 24 (1990) 363). The effect of the sample size on this type of dependence is also investigated, and an explicit expression is given for the population value of Kendall's coefficient of concordance between two arbitrary order statistics of a random sample.     

A new dependence ordering with applications

In this paper, we introduce a new copula-based dependence order to compare the relative degree of dependence between two pairs of random variables. Relationship of the new order to the existing dependence orders is investigated. In particular, the new ordering is stronger than the partial ordering, more monotone regression dependence as developed by Avérous et al. [J. Avérous, C. Genest, S.C. Kochar, On dependence structure of order statistics, Journal of Multivariate Analysis 94 (2005) 159-171]. Applications of this partial order to order statistics, k-record values and frailty models are given.     

Multivariate order statistics via multivariate concomitants

Let denote a set of n independent identically distributed k-dimensional absolutely continuous random variables. A general class of complete orderings of such random vectors is supplied by viewing them as concomitants of an auxiliary random variable. The resulting definitions of multivariate order statistics subsume and extend orderings that have been previously proposed such as norm ordering and N-conditional ordering. Analogous concepts of multivariate record values and multivariate generalized order statistics are also described.     

Multivariate analysis of variance with fewer observations than the dimension

In this article, we consider the problem of testing a linear hypothesis in a multivariate linear regression model which includes the case of testing the equality of mean vectors of several multivariate normal populations with common covariance matrix Σ, the so-called multivariate analysis of variance or MANOVA problem. However, we have fewer observations than the dimension of the random vectors. Two tests are proposed and their asymptotic distributions under the hypothesis as well as under the alternatives are given under some mild conditions. A theoretical comparison of these powers is made.     

Likelihood ratio tests of correlated multivariate samples

We develop methods to compare multiple multivariate normally distributed samples which may be correlated. The methods are new in the context that no assumption is made about the correlations among the samples. Three types of null hypotheses are considered: equality of mean vectors, homogeneity of covariance matrices, and equality of both mean vectors and covariance matrices. We demonstrate that the likelihood ratio test statistics have finite-sample distributions that are functions of two independent Wishart variables and dependent on the covariance matrix of the combined multiple populations. Asymptotic calculations show that the likelihood ratio test statistics converge in distribution to central Chi-squared distributions under the null hypotheses regardless of how the populations are correlated. Following these theoretical findings, we propose a resampling procedure for the implementation of the likelihood ratio tests in which no restrictive assumption is imposed on the structures of the covariance matrices. The empirical size and power of the test procedure are investigated for various sample sizes via simulations. Two examples are provided for illustration. The results show good performance of the methods in terms of test validity and power.     

Sequential order statistics with an order statistics prior

In the model of sequential order statistics, prior distributions are considered for the model parameters, which, for example, describe increasing load put on remaining components. Gamma priors are examined as well as priors out of a class of extended truncated Erlang distributions (ETED), which is introduced along with some properties. The choice of independent priors in both set-ups leads to respective independent, conjugate posterior distributions for the model parameters of sequential order statistics. Since, in practical applications, the model parameters will often be increasingly ordered, a multivariate prior is applied being the joint distribution of common ETED-order statistics. Whatever baseline distribution of the sequential order statistics is chosen, the joint posterior distribution turns out to be a Weinman multivariate exponential distribution. Posterior moments are given explicitly, and HPD credible sets for the model parameters are stated.     

Arithmetic equivalence for function fields, the Goss zeta function and a generalisation

A theorem of Tate and Turner says that global function fields have the same zeta function if and only if the Jacobians of the corresponding curves are isogenous. In this note, we investigate what happens if we replace the usual (characteristic zero) zeta function by the positive characteristic zeta function introduced by Goss. We prove that for function fields whose characteristic exceeds their degree, equality of the Goss zeta function is the same as Gaßmann equivalence (a purely group theoretical property), but this statement can fail if the degree exceeds the characteristic. We introduce a ‘Teichmüller lift’ of the Goss zeta function and show that equality of such is always the same as Gaßmann equivalence.     

Total positivity order and the normal distribution

Unlike the usual stochastic order, total positivity order is closed under conditioning. Here we provide a general formulation of the preservation properties of the order under conditioning; we study certain properties of the order including translation properties and the implications of having equality in the inequality defining the order. Specializing to the multivariate normal distribution, the study of total positivity order leads to new cones defined in terms of covariance M-matrices related to positive dependence, whose properties we study.     

On qualitative robustness of support vector machines

Support vector machines (SVMs) have attracted much attention in theoretical and in applied statistics. The main topics of recent interest are consistency, learning rates and robustness. We address the open problem whether SVMs are qualitatively robust. Our results show that SVMs are qualitatively robust for any fixed regularization parameter λ. However, under extremely mild conditions on the SVM, it turns out that SVMs are not qualitatively robust any more for any null sequence λ_n, which are the classical sequences needed to obtain universal consistency. This lack of qualitative robustness is of a rather theoretical nature because we show that, in any case, SVMs fulfill a finite sample qualitative robustness property.For a fixed regularization parameter, SVMs can be represented by a functional on the set of all probability measures. Qualitative robustness is proven by showing that this functional is continuous with respect to the topology generated by weak convergence of probability measures. Combined with the existence and uniqueness of SVMs, our results show that SVMs are the solutions of a well-posed mathematical problem in Hadamard’s sense.     

TESTS OF COVARIANCE MATRIX BY USING PROJECTION PURSUIT AND BOOTSTRAP METHOD

Testing equality of covariance matrix has long been an interesting issue in statistics inference, To overcome the sparseness of data points in high-dimensional space and deal with the general cases, the author suggests several projection pursuit type statistics. Some results on the limiting distidbutions of the statistics are obtained. Some properties of bootstrap approximation are investigated. Furthermore, for computational reasons an approximation for the statistics based on number-theoretic roethod is applied. Several simulation experiments are performed.     

A note on the maximum correlation for Baker’s bivariate distributions with fixed marginals

We investigate Baker’s bivariate distributions with fixed marginals which are based on order statistics, and find conditions under which the correlation converges to the maximum for Fréchet-Hoeffding upper bound as the sample size tends to infinity. The convergence rate of the correlation is also investigated for some specific cases.     

The fluted fragment with transitive relations

The fluted fragment is a fragment of first-order logic (without equality) in which, roughly speaking, the order of quantification of variables coincides with the order in which those variables appear as arguments of predicates. It is known that this fragment has the finite model property. We consider extensions of the fluted fragment with various numbers of transitive relations, as well as the equality predicate. In the presence of one transitive relation (together with equality), the finite model property is lost; nevertheless, we show that the satisfiability and finite satisfiability problems for this extension remain decidable. We also show that the corresponding problems in the presence of two transitive relations (with equality) or three transitive relations (without equality) are undecidable, even for the two-variable sub-fragment.     

Estimation of location and scale parameters for the Burr XII distribution using generalized order statistics

The minimum variance linear unbiased estimators (MVLUE), the best linear invariant estimators (BLIE) and the maximum likelihood estimators (MLE) based on n-selected generalized order statistics are presented for the parameters of the Burr XII distribution.     

A graph theoretic upper bound on the permanent of a nonnegative integer matrix. II. The extremal case

Matrices A for which the upper bound per(A)?1+min{Π(c_i?1), Π(r_i?1)} holds with equality are characterized. Cases where the bound is achieved correspond to multigraphs with the property that there exists a unique path from any vertex to any disjoint cycle union. This occurs precisely when some multigraph associated with A has the mastercycle property: all cycles thread all branchpoints in the same circular order. Such multigraphs may also be characterized as a circular concatenation of certain acyclic multigraphs, each having a unique source and sink. This analysis yields two normal forms for the extremal matrices, one based on a nested block decomposition, and another based on an overlapping block decomposition. The extremal cases are invariant under contractions, yielding another characterization.     

Restricted one way analysis of variance using the empirical likelihood ratio test

Empirical likelihood (EL) ratio tests are developed for testing for or against the hypothesis that k-population means μ₁,μ₂,…,μ_k are isotonic with respect to some quasi-order ? on {1,2,…,k}. The null asymptotic distributions are derived and are shown to be of chi-bar squared type. The asymptotic power of the proposed test for testing for equality of these means against the order restriction is derived under contiguous alternatives and a simulation study is carried out to investigate the finite sample behaviors of this test. In addition, an adjusted EL test is used to improve the small size performance of our test and an example is also discussed to illustrate the theoretical results.