A note on the mean past and the mean residual life of a (n − k + 1)-out-of-n system under multi monitoring

In the study of the reliability of technical systems, k-out-of-n systems play an important role. In the present paper, we consider a (n − k + 1)-out-of-n system consisting of n identical components such that the lifetimes of components are independent and have a common distribution function F. It is assumed that the number of monitoring is l and the total number of failures of the components at time t _i is m _i , i = 1, . . . , l − 1. Also at time t _l (t ₁ < . . . < t _l ) the system have failed or the system is still working. Under these conditions, the mean past lifetime, the mean residual lifetime of system and their properties are investigated.     

Number of observations near the maximum in an F α-scheme

Let K _n (a) be the number of observations in the interval (M _n ,?a, M _n ), where M _n is the maximum value in a sequence of size n. We study the asymptotic properties of K _n (a) under the F ^α-scheme and discuss the influence of the associated sequence α _n on the limit behaviour of this random variable.     

Precedence tests and Lehmann alternatives

G = F ^k (k > 1); G = 1 − (1−F) ^k (k < 1); G = F ^k (k < 1); and G = 1 − (1−F) ^k (k > 1), where F and G are two continuous cumulative distribution functions. If an optimal precedence test (one with the maximal power) is determined for one of these four classes, the optimal tests for the other classes of alternatives can be derived. Application of this is given using the results of Lin and Sukhatme (1992) who derived the best precedence test for testing the null hypothesis that the lifetimes of two types of items on test have the same distibution. The test has maximum power for fixed κ in the class of alternatives G = 1 − (1−F) ^k , with k < 1. Best precedence tests for the other three classes of Lehmann-type alternatives are derived using their results. Finally, a comparison of precedence tests with Wilcoxon's two-sample test is presented. Received: February 22, 1999; revised version: June 7, 2000     

Approximations of the Distributions of Test Statistics for Homogeneity of a Product Multinomial Model

Statistics R ^a based on power divergence can be used for testing the homogeneity of a product multinomial model. All R ^a have the same chi-square limiting distribution under the null hypothesis of homogeneity. R ⁰ is the log likelihood ratio statistic and R ¹ is Pearson's X ² statistic. In this article, we consider improvement of approximation of the distribution of R ^a under the homogeneity hypothesis. The expression of the asymptotic expansion of distribution of R ^a under the homogeneity hypothesis is investigated. The expression consists of continuous and discontinuous terms. Using the continuous term of the expression, a new approximation of the distribution of R ^a is proposed. A moment-corrected type of chi-square approximation is also derived. By numerical comparison, we show that both of the approximations perform much better than that of usual chi-square approximation for the statistics R ^a when a ≤ 0, which include the log likelihood ratio statistic.     

Note on the bias in the estimation of the serial correlation coefficient of AR(1) processes

We derive approximating formulas for the mean and the variance of an autocorrelation estimator which are of practical use over the entire range of the autocorrelation coefficient ρ. The least-squares estimator ∑ ^{n −1} _{i =1}ε _i ε _{i +1} / ∑ ^{n −1} _{i =1}ε² _i is studied for a stationary AR(1) process with known mean. We use the second order Taylor expansion of a ratio, and employ the arithmetic-geometric series instead of replacing partial Cesàro sums. In case of the mean we derive Marriott and Pope's (1954) formula, with (n− 1)⁻¹ instead of (n)⁻¹, and an additional term α (n− 1)⁻². This new formula produces the expected decline to zero negative bias as ρ approaches unity. In case of the variance Bartlett's (1946) formula results, with (n− 1)⁻¹ instead of (n)⁻¹. The theoretical expressions are corroborated with a simulation experiment. A comparison shows that our formula for the mean is more accurate than the higher-order approximation of White (1961), for |ρ| > 0.88 and n≥ 20. In principal, the presented method can be used to derive approximating formulas for other estimators and processes. Received: November 30, 1999; revised version: July 3, 2000     

On residual lifetimes in sequential (n ? k?+?1)-out-of-n systems

Sequential order statistics is an extension of ordinary order statistics. They model the successive failure times in sequential k-out-of-n systems, where the failures of components possibly affect the residual lifetimes of the remaining ones. In this paper, we consider the residual lifetime of the components after the kth failure in the sequential (n − k + 1)-out-of-n system. We extend some results on the joint distribution of the residual lifetimes of the remaining components in an ordinary (n − k + 1)-out-of-n system presented in Bairamov and Arnold (Stat Probab Lett 78(8):945–952, 2008) to the case of the sequential (n − k + 1)-out-of-n system.     

Estimating Average Worth of the Selected Subset from Two-Parameter Exponential Populations

ABSTRACT

Suppose independent random samples are available from k(k ≥ 2) exponential populations ∏₁,…,∏ _k with a common location θ and scale parameters σ₁,…,σ _k , respectively. Let X _i and Y _i denote the minimum and the mean, respectively, of the ith sample, and further let X = min{X ₁,…, X _k } and T _i  = Y _i  ? X; i = 1,…, k. For selecting a nonempty subset of {∏₁,…,∏ _k } containing the best population (the one associated with max{σ₁,…,σ _k }), we use the decision rule which selects ∏ _i if T _i  ≥ c max{T ₁,…,T _k }, i = 1,…, k. Here 0 < c ≤ 1 is chosen so that the probability of including the best population in the selected subset is at least P* (1/k ≤ P* < 1), a pre-assigned level. The problem is to estimate the average worth W of the selected subset, the arithmetic average of means of selected populations. In this article, we derive the uniformly minimum variance unbiased estimator (UMVUE) of W. The bias and risk function of the UMVUE are compared numerically with those of analogs of the best affine equivariant estimator (BAEE) and the maximum likelihood estimator (MLE).     

Improved approaches to interval estimation of effective doses for a logistic dose-response curve

Asymptotic confidence (delta) intervals and intervals based upon the use of Fieller's theorem are alternative methods for constructing intervals for the <$>\gamma<$>% effective doses (ED<$>_\gamma<$>). Sitter and Wu (1993) provided a comparison of the two approaches for the ED<$>_{50}<$>, for the case in which a logistic dose response curve is assumed. They showed that the Fieller intervals are generally superior. In this paper, we introduce two new families of intervals, both of which include the delta and Fieller intervals as special cases. In addition we consider interval estimation of the ED<$>_{90}<$> as well as the ED<$>_{50}<$>. We provide a comparison of the various methods for the problem of constructing a confidence interval for the ED<$>_\gamma<$>.     

Distribution of record statistics in a geometrically increasing population

The probability function and binomial moments of the number _Nn$N_n$ of (upper) records up to time (index) n in a geometrically increasing population are obtained in terms of the signless q-Stirling numbers of the first kind, with q   being the inverse of the proportion λ$\lambda$ of the geometric progression. Further, a strong law of large numbers and a central limit theorem for the sequence of random variables _Nn$N_n$, n=1,2,…,$n=1,2,\dots ,$ are deduced. As a corollary the probability function of the time _Tk$T_k$ of the kth record is also expressed in terms of the signless q  -Stirling numbers of the first kind. The mean of _Tk$T_k$ is obtained as a q  -series with terms of alternating sign. Finally, the probability function of the inter-record time _Wk=_Tk-_Tk-1$W_k=T_k-T_{k-1}$ is obtained as a sum of a finite number of terms of q  -numbers. The mean of _Wk$W_k$ is expressed by a q-series. As k   increases to infinity the distribution of _Wk$W_k$ converges to a geometric distribution with failure probability q. Additional properties of the q-Stirling numbers of the first kind, which facilitate the present study, are derived.     

Improved estimators for the selected location parameters

_i , i = 1, 2, ..., k be k independent exponential populations with different unknown location parameters θ _i , i = 1, 2, ..., k and common known scale parameter σ. Let Y _i denote the smallest observation based on a random sample of size n from the i-th population. Suppose a subset of the given k population is selected using the subset selection procedure according to which the population π _i is selected iff Y _i ≥Y ₍₁₎−d, where Y ₍₁₎ is the largest of the Y _i 's and d is some suitable constant. The estimation of the location parameters associated with the selected populations is considered for the squared error loss. It is observed that the natural estimator dominates the unbiased estimator. It is also shown that the natural estimator itself is inadmissible and a class of improved estimators that dominate the natural estimator is obtained. The improved estimators are consistent and their risks are shown to be O(kn ⁻²). As a special case, we obtain the coresponding results for the estimation of θ₍₁₎, the parameter associated with Y ₍₁₎. Received: January 6, 1998; revised version: July 11, 2000     

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.     

Generalized polya eggenberger family of distributions and its relation to lagrangian katz family

Janardan (1973) introduced the generalized Polya Eggenberger family of distributions (GPED) as a limiting distribution of the generalized Markov-Polya distribution (GMPD). Janardan and Rao (1982) gave a number of characterizing properties of the generalized Markov-Polya and generalized Polya Eggenberger distributions. Here, the GPED family characterized by four parameters, is formally defined and studied. The probability generating function, its moments, and certain recurrence relations with the moments are provided. The Lagrangian Katz family of distributions (Consul and Famoye (1996)) is shown to be a sub-class of the family of GPED (or GPED ₁ ) as it is called in this paper). A generalized Polya Eggenberger distribution of the second kind (GPED ₂ ) is also introduced and some of it's properties are given. Recurrence relations for the probabilities of GPED ₁ and GPED ₂ are given. A number of other structural and characteristic properties of the GPED ₁ are provided, from which the properties of Lagrangian Katz family follow. The parameters of GMPD ₁ are estimated by the method of moments and the maximum likelihood method. An application is provided.     

Non-parametric recursive estimates of a probability density function and its derivatives

Recursive estimates of a probability density function (pdf) are known. This paper presents recursive estimates of a derivative of any desired order of a pdf. Let f be a pdf on the real line and p?0 be any desired integer. Based on a random sample of size n from f, estimators f^(p)_n of f^(p), the pth order derivatives of f, are exhibited. These estimators are of the form $\text{n}{}^{\mathrm{?1}}\text{∑n}\text{j=1}\text{δ}{}_{\mathrm{jp}}$, where δ_jp depends only on p and the jth observation in the sample, and hence can be computed recursively as the sample size increases. These estimators are shown to be asymptotically unbiased, mean square consistent and strongly consistent, both at a point and uniformly on the real line. For pointwise properties, the conditions on f^(p) have been weakened with a little stronger assumption on the kernel function.     

On an application of concomitants of order statistics

ABSTRACT

Let (X_i, Y_i), i = 1, …, n be a pair where the first coordinate X_i represents the lifetime of a component, and the second coordinate Y_i denotes the utility of the component during its lifetime. Then the random variable Y_{[r: n]} which is known to be the concomitant of the rth order statistic defines the utility of the component which has the rth smallest lifetime. In this paper, we present a dynamic analysis for an n component system under the above-mentioned concomitant setup.     

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     

Bayesian decision theoretic approach to hypothesis problems with skewed alternatives

Many hypothesis problems in practice require the selection of the left side or the right side alternative when the null is rejected. For parametric models, this problem can be stated as _H0:θ=_θ0$H_0:\theta ={\theta }_0$vs.  _H−:θ<_θ0$H_{-}:\theta <{\theta }_0$ or _H+:θ>_θ0$H_{+}:\theta >{\theta }_0$. Frequentists use Type-III error (directional error) to develop statistical methodologies. This approach and other approaches considered in the literature do not take into account the situations where the selection of one side may be more important or when one side may be more probable than the other. This problem can be tackled by specifying a loss function and/or by specifying a hierarchical prior structure with allowing the skewness in the alternatives. Based on this, we develop a Bayesian decision theoretic methodology and show that the resulted Bayes rule perform better in the side of the alternatives which is more probable. The methodology can be also used in a frequentist's framework when it is desired to discover an alternative that is more important. We also consider the multiple hypotheses problem and develop new false discovery rates for the selection of the left and the right sides of alternatives. These discovery rates would be useful in the situations when one side of the alternatives are more important or more probable than the other.