Nonparametric model checking for k-out-of-n systems

It is an important problem in reliability analysis to decide whether for a given k-out-of-n system the static or the sequential k-out-of-n model is appropriate. Often components are redundantly added to a system to protect against failure of the system. If the failure of any component of the system induces a higher rate of failure of the remaining components due to increased load, the sequential k-out-of-n model is appropriate. The increase of the failure rate of the remaining components after a failure of some component implies that the effects of the component redundancy are diminished. On the other hand, if all the components have the same failure distribution and whenever a failure occurs, the remaining components are not affected, the static k-out-of-n model is adequate. In this paper, we consider nonparametric hypothesis tests to make a decision between these two models. We analyze test statistics based on the profile score process as well as test statistics based on a multivariate intensity ratio and derive their asymptotic distribution. Finally, we compare the different test statistics.     

Discrete q-distributions on Bernoulli trials with a geometrically varying success probability

Consider a sequence of independent Bernoulli trials and assume that the odds of success (or failure) or the probability of success (or failure) at the ith trial varies (increases or decreases) geometrically with rate (proportion) q, for increasing i=1,2,…. Introducing the notion of a geometric sequence of trials as a sequence of Bernoulli trials, with constant probability, that is terminated with the occurrence of the first success, a useful stochastic model is constructed. Specifically, consider a sequence of independent geometric sequences of trials and assume that the probability of success at the jth geometric sequence varies (increases or decreases) geometrically with rate (proportion) q, for increasing j=1,2,…. On both models, let X_n be the number of successes up the nth trial and T_k (or W_k) be the number of trials (or failures) until the occurrence of the kth success. The distributions of these random variables turned out to be q-analogues of the binomial and Pascal (or negative binomial) distributions. The distributions of X_n, for n→∞$n\to \infty$, and the distributions of W_k, for k→∞$k\to \infty$, can be approximated by a q  -Poisson distribution. Also, as k→0$k\to 0$, a zero truncated negative q  -binomial distribution _Uk=_Wk|_Wk>0$U_k=W_k|W_k>0$ can be approximated by a q-logarithmic distribution. These discrete q-distributions and their applications are reviewed, with critical comments and additions. Finally, consider a sequence of independent Bernoulli trials and assume that the probability of success (or failure) is a product of two sequences of probabilities with one of these sequences depending only the number of trials and the other depending only on the number of successes (or failures). The q-distributions of the number X_n of successes up to the nth trial and the number T_k of trials until the occurrence of the kth success are similarly reviewed.     

Concomitants of order statistics from multivariate elliptical distributions

In this paper, by considering a 2n-dimensional elliptically contoured random vector (X^T,Y^T)^T=(X₁,…,X_n,Y₁,…,Y_n)^T, we derive the exact joint distribution of linear combinations of concomitants of order statistics arising from X. Specifically, we establish a mixture representation for the distribution of the rth concomitant order statistic, and also for the joint distribution of the rth order statistic and its concomitant. We show that these distributions are indeed mixtures of multivariate unified skew-elliptical distributions. The two most important special cases of multivariate normal and multivariate t distributions are then discussed in detail. Finally, an application of the established results in an inferential problem is outlined.     

Concomitants of multivariate order statistics from multivariate elliptical distributions

ABSTRACT

In this article, we consider a (k + 1)n-dimensional elliptically contoured random vector (X^T₁, X₂^T, …, X^T_k, Z^T)^T = (X₁₁, …, X_1n, …, X_k1, …, X_kn, Z₁, …, Z_n)^T and derive the distribution of concomitant of multivariate order statistics arising from X₁, X₂, …, X_k. Specially, we derive a mixture representation for concomitant of bivariate order statistics. The joint distribution of the concomitant of bivariate order statistics is also obtained. Finally, the usefulness of our result is illustrated by a real-life data.     

On a system of nonhomogeneous components sharing a common frailty

The components of a reliability system subjected to a common random environment usually have dependent lifetimes. This paper studies the stochastic properties of such a system with lifetimes of the components following multivariate frailty models and multivariate mixed proportional reversed hazard rate (PRHR) models, respectively. Through doing stochastic comparison, we devote to throwing a new light on how the random environment affects the number of working components of a reliability system and on assessing the performance of a k-out-of-n system.     

Circular success and failure runs in a sequence of exchangeable binary trials

Let X₁,…,X_n be an exchangeable sequence of binary trials arranged on a circle with possible values “1” (success) or “0” (failure). In an exchangeable sequence, the joint distribution of X₁,X₂,…,X_n is invariant under the permutation of its arguments. For the circular sequence, general expressions for the joint distributions of run statistics based on the joint distribution of success and failure run lengths are obtained. As a special case, we present our results for Bernoulli trials. The results presented consist of combinatorial terms and therefore provide easier calculations. For illustration purposes, some numerical examples are given and the reliability of the circular combined k-out-of-n:G and consecutive k_c-out-of-n:G system under stress–strength setup is evaluated.     

The q-Bernstein basis as a q-binomial distribution

The q-Bernstein basis, used in the definition of the q-Bernstein polynomials, is shown to be the probability mass function of a q-binomial distribution. This distribution is defined on a sequence of zero–one Bernoulli trials with probability of failure at any trial increasing geometrically with the number of previous failures. A modification of this model, with the probability of failure at any trial decreasing geometrically with the number of previous failures, leads to a second q-binomial distribution that is also connected to the q-Bernstein polynomials. The q-factorial moments as well as the usual factorial moments of these distributions are derived. Further, the q-Bernstein polynomial B_n(f(t),q;x) is expressed as the expected value of the function f([X_n]_q/[n]_q) of the random variable X_n obeying the q-binomial distribution. Also, using the expression of the q-moments of X_n, an explicit expression of the q-Bernstein polynomial B_n(f_r(t),q;x), for f_r(t) a polynomial, is obtained.     

A note on the mean residual life function of a coherent system with exchangeable or nonidentical components

This article investigates some properties of the mean residual life function of (n−k+1)-out-of-n systems, when the lifetimes of the system components are independent random variables but not necessarily identically distributed and when the joint distribution of the component lifetimes is exchangeable, extending the results of Asadi and Goliforushani (2008) [On the mean residual life function of coherent systems. IEEE Transactions on Reliability 57 (4) 574-580] for the case of independent and identically distributed components. The extension to a coherent system with exchangeable components is also given.     

A q-Pólya urn model and the q-Pólya and inverse q-Pólya distributions

A q-Pólya urn model is introduced by assuming that the probability of drawing a white ball at a drawing varies geometrically, with rate q, both with the number of drawings and the number of white balls drawn in the previous drawings. Then, the probability mass functions and moments of (a) the number of white balls drawn in a specific number of drawings and (b) the number of black balls drawn until a specific number of white balls are drawn are derived. These two distributions turned out to be q-analogs of the Pólya and the inverse Pólya distributions, respectively. Also, the limiting distributions of the q-Pólya and the inverse q-Pólya distributions, as the number of balls in the urn tends to infinity, are shown to be a q-binomial and a negative q-binomial distribution, respectively. In addition, the positive or negative q-hypergeometric distribution is obtained as conditional distribution of a positive or negative q-binomial distribution, given its sum with another positive or negative q-binomial distribution, independent of it.     

Optimal k-circulant supersaturated designs

Optimal k-circulant supersaturated designs have been constructed in literature using computer intensive methods. A systematic method of construction for multi-level experiments based on balanced incomplete block designs is presented in this paper. The method is also applicable to two-level experiments. Illustrative examples are also given.     

An optimal response adaptive biased coin design with k heteroscedastic treatments

For comparing treatments in clinical trials, Atkinson (1982) introduced optimal biased coins for balancing patients across treatment assignments by using D-optimality under the assumption of homoscedastic responses of different treatments. However, this assumption can be violated in many real applications. In this paper, we relax the homoscedasticity assumption in the k treatments setting with k>2. A general family of optimal response adaptive biased coin designs are proposed following Atkinson's procedure. Asymptotic properties of the proposed designs are obtained. Some advantages of the proposed design are discussed.     

Dependence Structure of Spacings of Order Statistics

Let X_1:n ≤ X_2:n ≤···≤ X_n:n denote the order statistics of a sample of n independent random variables X₁, X₂,…, X_n, all identically distributed as some X. It is shown that if X has a log-convex [log-concave] density function, then the general spacing vector (X_k1:n, X_k2:n ? X_k1:n,…, X_kr:n ? X_kr?1:n) is MTP₂ [S-MRR₂] whenever 1 ≤ k₁ < k₂ <···< k_r ≤ n and 1 ≤ r ≤ n. Multivariate likelihood ratio ordering of such general spacing vectors corresponding to two random samples is also considered. These extend some of the results in the literature for usual spacing vectors.     

Regularized parameter estimation of high dimensional t distribution

We propose penalized-likelihood methods for parameter estimation of high dimensional t distribution. First, we show that a general class of commonly used shrinkage covariance matrix estimators for multivariate normal can be obtained as penalized-likelihood estimator with a penalty that is proportional to the entropy loss between the estimate and an appropriately chosen shrinkage target. Motivated by this fact, we then consider applying this penalty to multivariate t distribution. The penalized estimate can be computed efficiently using EM algorithm for given tuning parameters. It can also be viewed as an empirical Bayes estimator. Taking advantage of its Bayesian interpretation, we propose a variant of the method of moments to effectively elicit the tuning parameters. Simulations and real data analysis demonstrate the competitive performance of the new methods.     

Rank-based empirical likelihood inference on medians of k populations

We propose a nonparametric method, called rank-based empirical likelihood (REL), for making inferences on medians and cumulative distribution functions (CDFs) of k populations. The standard distribution-free approach to testing the equality of k medians requires that the k population distributions have the same shape. Our REL-ratio (RELR) test for this problem requires fewer assumptions and can effectively use the symmetry information when the distributions are symmetric. Furthermore, our RELR statistic does not require estimation of variance, and achieves asymptotic pivotalness implicitly. When the k populations have equal medians we show that the REL method produces valid inferences for the common median and CDFs of k populations. Simulation results show that the REL approach works remarkably well in finite samples. A real data example is used to illustrate the proposed REL method.     

Reliability analysis of consecutive k-out-of-n systems with non-identical components lifetimes

In recent years, the study of reliability properties of consecutive k-out-of-n systems has attracted a great deal of attention from both theoretical and practical perspectives. In this paper we consider linear and circular consecutive k-out-of-n systems. It is assumed that lifetimes of components of the systems are independent but their probability distributions are non-identical. We study the reliability properties of the residual lifetimes of such systems under the condition that at least (n−r+1), r≤n, components of the system are operating. We also investigate the probability that a specific number of components of the above-mentioned system operate at time t, t>0, under the condition that the system is alive at time t.     

Order restricted statistical inference for scale parameters based on sequential order statistics

In the situation of a multi-sample experiment consisting of differently equipped sequential k-out-of-n systems, scale parameters of underlying distributions from a general location-scale family of distributions are estimated under an order restriction. In each sample, the case of missing the smallest observations is included. Moreover, based on a profile likelihood a homogeneity test against an ordered alternative is proposed and analyzed. This work extends an approach of Bhattacharya [2007. Testing for ordered failure rates under general progressive censoring. J. Statist. Plann. Inference 137, 1775–1786] in the progressive Type-II censoring framework.     

Lower bounds on the expectations of upper record values

For infinite sequences of independent random variables with identical continuous distributions, we establish optimal lower bounds on the deviations of the expectations of record values from population means in units generated by the central absolute moments of various orders. The bounds are non-negative for the classic record values, and non-positive for the other kth records with k?2. We also provide analogous bounds for the record increments.     

On k-th record times,record values and their moments

We give new formula for moments of k-th record values in terms of Stirling numbers of the first kind. In particular, the formulae allow to derive the explicit formulae for moments of k-th lower record values from exponential distribution which have not been known yet. Moreover, some interesting identities involving harmonic numbers are also obtained as corollaries to presented results.     

Reliability Evaluation for a Multi-State System Under Stress-Strength Setup

The two most commonly used reliability models in engineering applications are binary k-out-of-n:G and consecutive k-out-of-n:G systems. Multi-state k-out-of-n:G and multi-state consecutive k-out-of-n:G systems have been proposed as an extension of these systems and they have been found to be more flexible tool for modeling engineering systems. In this article, multi-state systems, in particular, multi-state k-out-of-n:G and multi-state consecutive k-out-of-n:G, are considered in a stress-strength setup. The states of the system are classified considering the number of components whose strengths above (below) the multiple stresses available in an environment. The exact state probabilities are provided and the results are illustrated for various stress-strength distributions. Maximum likelihood estimators of state probabilities are also presented.