Evaluation of P(Xt>Yt ) when both Xt and Yt are residual lifetimes of two systems

Let the random variables X and Y denote the lifetimes of two systems. In reliability theory to compare between the lifetimes of X and Y there are several approaches. Among the most popular methods of comparing the lifetimes are to compare the survival functions, the failure rates and the mean residual lifetime functions of X and Y. Assume that both systems are operating at time t > 0. Then the residual lifetimes of them are X_t=X?t | X>t and Y_t=Y?t | Y>t, respectively. In this paper, we introduce, by taking into account the age of systems, a time‐dependent criterion to compare the residual lifetimes of them. In other words, we concentrate on function R(t ):=P(X_t>Y_t) which enables one to obtain, at time t, the probability that the residual lifetime X_t is greater than the residual lifetime Y_t. It is mentioned, in Brown and Rutemiller (IEEE Transactions on Reliability, 22 , 1973) that the probability of type R(t) is important for designing as long‐lived a product as possible. Several properties of R(t) and its connection with well‐known reliability measures are investigated. The estimation of R(t) based on samples from X and Y is also discussed.     

Dispersive ordering of fail-safe systems with heterogeneous exponential components

Let X ₁, . . . , X _n be independent exponential random variables with respective hazard rates λ₁, . . . , λ _n , and Y ₁, . . . , Y _n be independent and identically distributed random variables from an exponential distribution with hazard rate λ. Then, we prove that X _2:n , the second order statistic from X ₁, . . . , X _n , is larger than Y _2:n , the second order statistic from Y ₁, . . . , Y _n , in terms of the dispersive order if and only if

$\lambda\geq \sqrt{\frac{1}{{n\choose 2}}\sum_{1\leq i < j\leq n}\lambda_i\lambda_j}.$

We also show that X _2:n is smaller than Y _2:n in terms of the dispersive order if and only if

$ \lambda\le\frac{\sum^{n}_{i=1} \lambda_i-{\rm max}_{1\leq i\leq n} \lambda_i}{n-1}. $

Moreover, we extend the above two results to the proportional hazard rates model. These two results established here form nice extensions of the corresponding results on hazard rate, likelihood ratio, and MRL orderings established recently by Pǎltǎnea (J Stat Plan Inference 138:1993–1997, 2008), Zhao et al. (J Multivar Anal 100:952–962, 2009), and Zhao and Balakrishnan (J Stat Plan Inference 139:3027–3037, 2009), respectively.

     

Sulle operazioni associative tra variabili casuali

Si mostra che, sotto condizioni di regolarità, seo è un’operazione associativa tra variabili casuali reali e indipendenti, è definibile una trasformata integrale ξ delle loro funzioni di ripartizione con la proprietà: ξx ₀ Y (t)=ξx(t)·ξ _y (t). Si indicano alcune proprietà di tale trasformata e si tratta della possibilità di estendere a un’operazione associativa risultati noti per l’addizione tra variabili casuali. In particolare ci si occupa dell’« infinita divisibilità » fornendo condizioni perché una variabile casualeX ammetta la rappresentazioneX=X _1O X _2O … _o X _n per ognin naturale con leX _i indipendenti e identicamente ditribuite.     

Characterization of the multivariate normal distribution by conditional normal distributions

Summary LetX andY be two random vectors with values in ℝ ^k and ℝ∝, respectively. IfZ=(X ^T,Y ^T) ^T is multivariate normal thenX givenY=y andY givenX=x are (multivariate) normal; the converse is wrong. In this paper simple additional conditions are stated such that the converse is true, too. Furthermore, the case is treated that the random vectorZ=(X ₁ ^T , …,X _t ^T ) ^T is splitted intot≥3 partsX ₁, …,X _t.     

Characterization of the exponential distribution function by Properties of the differenceX k+s∶n−X k∶n of order statistics

We considern independent and identically distributed random variables with common continuous distribution functionF concentrated on (0, ∞). LetX _1∶n≤X_2∶n...≤X_n∶n be the corresponding order statistics. Put $$d_s \left( x \right) = P\left( {X_{k + s:n} - X_{k:n} \geqslant x} \right) - P\left( {X_{s:n - k} \geqslant x} \right), x \geqslant 0,$$ and $$\delta _s \left( {x, \rho } \right) = P\left( {X_{k + s:n} - X_{k:n} \geqslant x} \right) - e^{ - \rho \left( {n - k} \right)x} ,\rho > 0,x \geqslant 0.$$ Fors=1 it is well known that each of the conditions d₁(x)=O ?x≥0 and δ₁ (x, p) = O ?x≥0 implies thatF is exponential; but the analytic tools in the proofs of these two statements are radically different. In contrast to this in the present paper we present a rather elementary method which permits us to derive the above conclusions for somes, 1≤n —k, using only asymptotic assumptions (either forx→0 orx→∞) ond _s(x) and δ₁ (x, p), respectively.     

Correcting for truncation bias caused by a latent truncation variable

We discuss estimation of the model Y_i=X_ib_Y+e_Yi, T_i=X_ib_T+ e_Ti, when data on the continuous dependent variable Y and on the independent variables X are observed iff the ‘truncation variable’ T>0 and when T is latent. This case is distinct from both (i) the‘censored sample’ case, in which Y data are available iff T>0, T is latent and X data are available for all observations, and (ii) the ‘observed truncation variable’ case, in which both Y and X are observed iff T>0 and in which the actual value of T is observed whenever T>0. We derive a maximum-likelihood procedure for estimating this model and discuss identification and estimation.     

On generalized transforms of distribution functions

In this note we discuss the following problem. LetX andY to be two real valued independent r.v.'s with d.f.'sF and ?. Consider the d.f.F*? of the r.v.X oY, being o a binary operation among real numbers. We deal with the following equation: $$\mathcal{G}^1 (F * \phi ,s) = \mathcal{G}^2 (F,s)\square \mathcal{G}^3 (\phi ,s)\forall s \in S$$ where \(\mathcal{G}^1 ,\mathcal{G}^2 ,\mathcal{G}^3 \) are real or complex functionals, т another binary operation ands a parameter. We give a solution, that under stronger assumptions (Aczél 1966), is the only one, of the problem. Such a solution is obtained in two steps. First of all we give a solution in the very special case in whichX andY are degenerate r.v.'s. Secondly we extend the result to the general case under the following additional assumption: $$\begin{gathered} \mathcal{G}^1 (\alpha F + (1 - \alpha )\phi ,s) = H[\mathcal{G}^i (F,s),\mathcal{G}^i (\phi ,s);\alpha ] \hfill \\ \forall \alpha \in [0,1]i = 1,2,3 \hfill \\ \end{gathered} $$ .     

Identification and estimation of nonlinear dynamic panel data models with unobserved covariates

This paper considers nonparametric identification of nonlinear dynamic models for panel data with unobserved covariates. Including such unobserved covariates may control for both the individual-specific unobserved heterogeneity and the endogeneity of the explanatory variables. Without specifying the distribution of the initial condition with the unobserved variables, we show that the models are nonparametrically identified from two periods of the dependent variable _Yit$Y_{it}$ and three periods of the covariate _Xit$X_{it}$. The main identifying assumptions include high-level injectivity restrictions and require that the evolution of the observed covariates depends on the unobserved covariates but not on the lagged dependent variable. We also propose a sieve maximum likelihood estimator (MLE) and focus on two classes of nonlinear dynamic panel data models, i.e., dynamic discrete choice models and dynamic censored models. We present the asymptotic properties of the sieve MLE and investigate the finite sample properties of these sieve-based estimators through a Monte Carlo study. An intertemporal female labor force participation model is estimated as an empirical illustration using a sample from the Panel Study of Income Dynamics (PSID).     

The average number of events per time unit and its estimation

Let ζ _t be the number of events which will be observed in the time interval [0;t] and define as the average number of events per time unit if this limit exists. In the case of i.i.d. waiting-times between the events,E[ζ _t ] is the renewal function and it follows from well-known results of renewal theory thatA exists and is equal to 1/τ, if τ>0 is the expectation of the waiting-times. This holds true also when τ = ∞.A may be estimate by ζ _t /t or where is the mean of the firstn waiting-timesX ₁,X ₂, ...,X _n . Both estimators converage with probability 1 to 1/τ if theX _i are i.i.d.; but the expectation of may be infinite for alln and also if it is finite, is in general a positively biased estimator ofA. For a stationary renewal process, ζ _t /t is unbiased for eacht; if theX _i are i.i.d. with densityf(x), then ζ _t /t has this property only iff(x) is of the exponential type and only for this type the numbers of events in consecutive time intervals [0,t], [t, 2t], ... are i.i.d. random variables for arbitraryt > 0.     

Distributions determined by conditioning on a single order statistic

LetX ₁,X ₂, …,X _n(n ? 2) be a random sample on a random variablex with a continuous distribution functionF which is strictly increasing over (a, b), ?∞ ?a <b ? ∞, the support ofF andX _1:n ?X _2:n ? … ?X _n:n the corresponding order statistics. Letg be a nonconstant continuous function over (a, b) with finiteg(a ⁺) andE {g(X)}. Then for some positive integers, 1 <s ?n $$E\left\{ {\frac{1}{{s - 1}}\sum\limits_{i - 1}^{s - 1} {g(X_{i:n} )|X_{s:n} } = x} \right\} = 1/2(g(x) + g(a^ + )), \forall x \in (a,b)$$ iffg is bounded, monotonic and \(F(x) = \frac{{g(x) - g(a^ + )}}{{g(b^ - ) - g(a^ + )}},\forall x \in (a,b)\) . This leads to characterization of several distribution functions. A general form of this result is also stated.     

Almost sure central limit theorem for the products of <Emphasis Type="Italic">U</Emphasis>-statistics

Let (X _n ) be a sequence of i.i.d random variables and U _n a U-statistic corresponding to a symmetric kernel function h, where h ₁(x ₁) = Eh(x ₁, X ₂, X ₃, . . . , X _m ), μ = E(h(X ₁, X ₂, . . . , X _m )) and ? ₁ = Var(h ₁(X ₁)). Denote \({\gamma=\sqrt{\varsigma_{1}}/\mu}\), the coefficient of variation. Assume that P(h(X ₁, X ₂, . . . , X _m ) > 0) = 1, ? ₁ > 0 and E|h(X ₁, X ₂, . . . , X _m )|³ < ∞. We give herein the conditions under which

$\lim_{N\rightarrow\infty}\frac{1}{\log N}\sum_{n=1}^{N}\frac{1}{n}g\left(\left(\prod_{k=m}^{n}\frac{U_{k}}{\mu}\right)^{\frac{1}{m\gamma\sqrt{n}}}\right) =\int\limits_{-\infty}^{\infty}g(x)dF(x)\quad {\rm a.s.}$

for a certain family of unbounded measurable functions g, where F(·) is the distribution function of the random variable \({\exp(\sqrt{2} \xi)}\) and ξ is a standard normal random variable.

     

Optimal process parameters under LINEX loss function with general input quality characteristic

In this paper we generalize the quality and cost trade-off problem of Chang and Hung (Qual Quant 41: 291–301, 2007) under the LINEX loss function. We consider the general input characteristic given by the random variable X with moment generating function m _X (t) and output characteristic given by the deterministic transformation Y  =  g(X). The two cases we consider are when g(X) is an affine function of X and X follows (1) the gamma distribution, and (2) the double exponential distribution.     

Alarm the quickly changing of performance evaluation decisive indexes—qualitative to quantitative?

This paper research about a special phenomenon and what is the main fators which is decided the evaluation total scores of local government in diferent years of China. The contents of the paper included introduction, hypothesis and research methodology, data and quantitative model, editing data, analyzing data, data showed and conclusion. The paper tests the new hypothesis of the quantitative evaluation, and it is right that new public policy and its tools of public management are effectively in China, so the situation is changed that the quantitative evaluation becomes the decisive factor in the performance total score by the model of original evaluation data analysis. The following public policy is effective in the year of China local government: the new public policy increases strength of evaluation; the new public policy deal with the weakness of evaluation before 2009; make new public policy to fight threats of evaluation; and make new policy for opportunities of evaluation. The paper also tests the hypothesis of the quantitative evaluations and which fator is positive correlated with results of performance, the quantitative evaluation is the decisive factors of evaluation scores in local governments of China by the years of changing. This paper research about the ratio of quantitative scores and the relationship between different fators of Chinese-style performance evaluation, and the data showing that the caves of related fators are almost the same level, it is proved that the reason of the qualitative scores become the negative correlation factors of performance score after 2009. At last, the paper builds the evaluation model by econometrics and thinking about performance evaluation practice of Chinese local governments as following: $${\begin{array}{ll}Y=0.23+0.81_{1} Y_{1} -0.18Y_{2} +0.44Y_{3} +0.01Y_{4} +e_{t2}\\ \hspace*{2.7pc}(0.35) \hspace*{1.2pc}(0.07) \hspace*{1.2pc} (0.05)\hspace*{1.2pc} (0.06) \hspace*{1.2pc} (0.0003)\\ \hspace*{6pc}{N}={12\quad R^2}=0.99\end{array}}$$ The value of constant 0.23 means Y = exp 0.23 × 10⁹ while every variable goes back to zero. The value of Y will be increase 0.81%, ?0.18%, 0.44%, 0.01% while Y ₁, Y ₂, Y ₃, Y ₄ increase 1%, and the prefer condition is that the other variables are not change. The article shows a changing map of performance evaluation in China.     

A new non-linear AR(1) time series model having approximate beta marginals

We consider the mixed AR(1) time series model $$X_t=\left\{\begin{array}{ll}\alpha X_{t-1}+ \xi_t \quad {\rm w.p.} \qquad \frac{\alpha^p}{\alpha^p-\beta ^p},\\ \beta X_{t-1} + \xi_{t} \quad {\rm w.p.} \quad -\frac{\beta^p}{\alpha^p-\beta ^p} \end{array}\right.$$ for ?1 < β ^p ≤ 0 ≤ α ^p  < 1 and α ^p ? β ^p  > 0 when X _t has the two-parameter beta distribution B₂(p, q) with parameters q > 1 and ${p \in \mathcal P(u,v)}$ , where $$\mathcal P(u,v) = \left\{u/v : u < v,\,u,v\,{\rm odd\,positive\,integers} \right\}.$$ Special attention is given to the case p = 1. Using Laplace transform and suitable approximation procedures, we prove that the distribution of innovation sequence for p = 1 can be approximated by the uniform discrete distribution and that for ${p \in \mathcal P(u,v)}$ can be approximated by a continuous distribution. We also consider estimation issues of the model.     

Beware of averages!

Summary Let X₁,.,., X_m, and Y₁, Y_n, be two independent samples from the same distribution and let X and Y be the means of these samples. What is the maximal value of P(X < Y)?     

Consistency of a nonparametric estimation of a density functional

LetX be a random variable with distribution functionF and density functionf. Let ? and ψ be known measurable functions defined on the real lineR and the closed interval [0, 1], respectively. This paper proposes a smooth nonparametric estimate of the density functional \(\theta = \int\limits_R \phi (x) \psi \left[ {F (x)} \right]f^2 (x) dx\) based on a random sampleX ₁, ...,X _n fromF using a kernel functionk. The proposed estimate is given by \(\hat \theta = (n^2 a_n )^{ - 1} \mathop \sum \limits_{i = 1}^n \mathop \sum \limits_{j = 1}^n \phi (X_i ) \psi \left[ {\hat F (X_i )} \right]k\left[ {(X_i - X_j )/a_n } \right]\) , where \(\hat F(x) = n^{ - 1} \mathop \sum \limits_{i = 1}^n K\left[ {(x - X_i )/a_n } \right]\) with \(K (w) = \int\limits_{ - \infty }^w {k (u) } du\) . The estimate \(\hat \theta \) is shown to be consistent both in the weak and strong sense and is used to estimate the asymptotic relative efficiency of various nonparametric tests, with particular reference to those using the Chernoff-Savage statistic.     

Estimating features of a distribution from binomial data

We propose estimators of features of the distribution of an unobserved random variable W. What is observed is a sample of Y,V,X where a binary Y equals one when W exceeds a threshold V determined by experimental design, and X are covariates. Potential applications include bioassay and destructive duration analysis. Our empirical application is referendum contingent valuation in resource economics, where one is interested in features of the distribution of values W (willingness to pay) placed by consumers on a public good such as endangered species. Sample consumers with characteristics X are asked whether they favor (with Y=1 if yes and zero otherwise) a referendum that would provide the good at a cost V specified by experimental design. This paper provides estimators for quantiles and conditional on X moments of W under both nonparametric and semiparametric specifications.     

A new measure of association between random variables

We propose a new measure of association between two continuous random variables X and Y based on the covariance between X and the log-odds rate associated to Y. The proposed index of correlation lies in the range [\(-1\), 1]. We show that the extremes of the range, i.e., \(-1\) and 1, are attainable by the Fr\(\acute{\mathrm{e}}\)chet bivariate minimal and maximal distributions, respectively. It is also shown that if X and Y have bivariate normal distribution, the resulting measure of correlation equals the Pearson correlation coefficient \(\rho \). Some interpretations and relationships to other variability measures are presented. Among others, it is shown that for non-negative random variables the proposed association measure can be represented in terms of the mean residual and mean inactivity functions. Some illustrative examples are also provided.     

Asymptotic normality of M-estimators in a semiparametric model with longitudinal data

Suppose that the longitudinal observations (Y _ij , X _ij , t _ij ) for i = 1, . . . ,n; j = 1, . . . ,m _i are modeled by the semiparamtric model where β ₀ is a k × 1 vector of unknown parameters, g(·) is an unknown estimated function and e _ij are unobserved disturbances. This article consider M-type regressions which include mean, median and quantile regressions. The M-estimator of the slope parameter β ₀ is obtained through piecewise local polynomial approximation of the nonparametric component. The local M-estimator of g(·) is also obtained by replacing β ₀ in model with its M-estimator and using local linear approximation. The asymptotic distribution of the estimator of β ₀ is derived. The asymptotic distributions of the local M-estimators of g(·) at both interior and boundary points are also established. Various applications of our main results are given. The research is supported in part by National Natural Science Foundation of China (Grant No. 10671089).