Note on estimation of reliability under bivariate pareto stress-strength model

In this paper, we discuss the problem of estimating reliability (R) of a component based on maximum likelihood estimators (MLEs). The reliability of a component is given byR=P[Y<X]. Here X is a random strength of a component subjected to a random stress(Y) and (X,Y) follow a bivariate pareto(BVP) distribution. We obtain an asymptotic normal(AN) distribution of MLE of the reliability(R).     

Markov Chain Markov Field dynamics: Models and statistics

This study deals with time dynamics of Markov fields defined on a finite set of sites with state space <$>E<$>, focussing on Markov Chain Markov Field (MCMF) evolution. Such a model is characterized by two families of potentials: the instantaneous interaction potentials, and the time delay potentials. Four models are specified: auto-exponential dynamics (<$>E = {\of R}^+<$>), auto-normal dynamics (<$>E = {\of R}<$>), auto-Poissonian dynamics (<$>E = {\of N}<$>) and auto-logistic dynamics ( E qualitative and finite). Sufficient conditions ensuring ergodicity and strong law of large numbers are given by using a Lyapunov criterion of stability, and the conditional pseudo-likelihood statistics are summarized. We discuss the identification procedure of the two Markovian graphs and look for validation tests using martingale central limit theorems. An application to meteorological data illustrates such a modelling.     

Unbiastld estimation of parametric functions in sampling from two two-truncation parameter families

We consider the problem of UMVU estimation of a U-estimable function of four unknown truncation parameters based on two independent random samples from two two-truncation parameter families. In particular, we obtain the UMVU estimator of functional, P (Y > X). Also the confidence intervals for some parametric functions are obtained.     

Estimation of the Reliability of a Stress-Strength System from Power Lindley Distributions

In this paper, we are interested in the estimation of the reliability parameter R = P(X > Y) where X, a component strength, and Y, a component stress, are independent power Lindley random variables. The point and interval estimation of R, based on maximum likelihood, nonparametric and parametric bootstrap methods, are developed. The performance of the point estimate and confidence interval of R under the considered estimation methods is studied through extensive simulation. A numerical example, based on a real data, is presented to illustrate the proposed procedure.     

The Robustness of Several Estimators of the Survivorship Function with Randomly Censored Data

The problem of estimating the survivorship function, R(t) = P(T > t), arises frequently in both engineering and biomedical sciences. In many applications the data one sees are censored due to the occurrence of some competing cause of failure such as withdrawal from the study, failure from some cause not under study, etc. In the biomedical sciences the distribution free estimator suggested by Kaplan and Meier (JASA 1958) is routinely used, while in the engineering sciences a parametric approach is more commonly used. In this report we study the efficiency of these two techniques when a particular parametric model such as the exponential, Weibull, normal, log normal, exponential power, Pareto, Gompertz, gamma, or bathtub shaped hazard distribution is assumed under a variety of censoring schemes and underlying failure models. We conclude that in most cases the parametric estimators outperform the distribution free estimator. The results are particularly striking if the Weibull forms of these estimators are used routinely.     

Estimating population size with truncated sampling

Let X₁, X₂,…,X_n be independent, indentically distributed random variables with density f(x,θ) with respect to a σ-finite measure μ. Let R be a measurable set in the sample space X. The value of X is observable if X ? (X?R) and not otherwise. The number J of observable X’s is binomial, N, Q, Q = 1?P(X ? R). On the basis of J observations, it is desired to estimate N and θ. Estimators considered are conditional and unconditional maximum likelihood and modified maximum likelihood using a prior weight function to modify the likelihood before maximizing. Asymptotic expansions are developed for the [Ncirc]’s of the form [Ncirc] = N + α√N + β + o_p(1), where α and β are random variables. All estimators have the same α, which has mean 0, variance σ² (a function of θ) and is asymptotically normal. Hence all are asymptotically equivalent by the usual limit distributional theory. The β’s differ and Eβ can be considered an “asymptotic bias”. Formulas are developed to compare the asymptotic biases of the various estimators. For a scale parameter family of absolutely continuous distributions with X = (0,∞) and R = (T,∞), special formuli are developed and a best estimator is found.     

UMVU estimation of the reliability function of the generalized life distributions

The problems of estimating the reliability function and P=P_rX > Y are considered for the generalized life distributions. Uniformly minimum variance unbiased estimators (UMVUES) of the powers of the parameter involved in the probabilistic model and the probability density function (pdf) at a specified point are derived. The UMVUE of the pdf is utilized to obtain the UMVUE of the reliability function and ‘P’. Our method of obtaining these estimators is quite simple than the traditional approaches. A theoretical method of studying the behaviour of the hazard-rate is provided.     

Tests for the mean of a multivariate normal population

We study the problem of testing: H₀ : μ ∈ P against H₁ : μ ? P, based on a random sample of N observations from a p-dimensional normal distribution N_p(μ, Σ) with Σ > 0 and P a closed convex positively homogeneous set. We develop the likelihood-ratio test (LRT) for this problem. We show that the union-intersection principle leads to a test equivalent to the LRT. It also gives a large class of tests which are shown to be admissible by Stein's theorem (1956). Finally, we give the α-level cutoff points for the LRT.     

Robust Liu estimator for regression based on an M-estimator

Consider the regression model y = beta 0 1 + Xbeta + epsilon. Recently, the Liu estimator, which is an alternative biased estimator beta L (d) = (X'X + I) -1 (X'X + dI)beta OLS , where 0<d<1 is a parameter, has been proposed to overcome multicollinearity . The advantage of beta L (d) over the ridge estimator beta R (k) is that beta L (d) is a linear function of d. Therefore, it is easier to choose d than to choose k in the ridge estimator. However, beta L (d) is obtained by shrinking the ordinary least squares (OLS) estimator using the matrix (X'X + I) -1 (X'X + dI) so that the presence of outliers in the y direction may affect the beta L (d) estimator. To cope with this combined problem of multicollinearity and outliers, we propose an alternative class of Liu-type M-estimators (LM-estimators) obtained by shrinking an M-estimator beta M , instead of the OLS estimator using the matrix (X'X + I) -1 (X'X + dI).     

Estimation of reliability of a component subjected to bivariate exponential stress

In this paper, we estimate the reliability of a component subjected to two different stresses which are independent of the strength of a component. We assume that the distribution of stresses follow a bivariate exponential (BVE) distribution. If X is the strength of a component subjected to two stresses (Y ₁,Y ₂), then the reliability of a component is given by R=P[Y ₁+Y ₂<X]. We estimate R when (Y ₁,Y ₂) follow different BVE models proposed by Marshall-Olkin (1967), Block-Basu-(1974), Freund (1961) and Proschan-Sullo (1974). The distribution of X is assumed to be exponential. The asymptotic normal (AN) distributions of these estimates of R are obtained.     

Exponential models,brownian motion,and independence

Some examples of steep, reproductive exponential models are considered. These models are shown to possess a τ-parallel foliation in the terminology of Barndorff-Nielsen and Blaesild. The independence of certain functions follows directly from the foliation. Suppose X(t) is a Wiener process with drift where X(t) = W(t) + ct, 0 < t < T. Furthermore let Y = max [X(s), 0 < s < T]. The joint density of Y and X = X(T), the end value, is studied within the framework of an exponential model, and it is shown that Y(Y – X) is independent of X. It is further shown that Y(Y – X) suitably scaled has an exponential distribution. Further examples are considered by randomizing on T.     

Nonparametric estimator for mean residual life and vitality function

In this paper, we study the estimation of the vitality function(v(x)=E(X|X>x) and mean residual life function(e(x)=E(X-x|X>x) from a sample ofX using the empirical estimator and kernel estimator. Under suitable conditions of regularity, the asymptotic normality of the kernel estimator is obtained. Partially supported by Consejeria de Cultura y Ed. (C.A.R.M.), under Grant PIB 95/90.     

Estimation of p(y<x) in the burr case: a comparative study

This paper provides a simulation study which compares three estimators for R = P(Y<X) when Y and X are two independent but not identically distributed Burr random variables. These estimators are the minimum variance unbiased, the maximum likelihood and Bayes estimators. Moreover, the sensitivity of Bayes estimator to the prior parameters is considered.     

Tables of two-sided tolerance factors for a normal distribution

Given a random sample of size N from a normal distribution, we consider tolerance intervals of the form X ? ks to X + ks, where X is the sample mean and s is the sample standard deviation. The value of k is chosen so that the interval covers a given proportion P of the population with confidence γ. Exact values of k, computed from numerical integration, are given for N = 2(1)100; P = 0.75, 0.90, 0.95, 0.975, 0.99, 0.995, 0.999; and γ = 0.5, 0.75, 0.90, 0.95, 0.975, 0.99, 0.995. The exact values are compared with the values obtained from an approximation developed by Wald and Wolfowitz (1946).     

Empirical bayes estimation of functionals of unknown probability measures

The empirical Dayes approach to one and two sal-npie problcrns has beeir considered by Korwar and Hollander (1976), Holiander and Korwar (1976) and Phadia and Susarla (1979). In this article we essen- tially generalize their empirical Bayes results by replacing the inlicaro-functions of. the sets (?∞,x) and {X≦Y} by arbitrary mea5, irable functions h(x) and h(x,y). More speclfically, the ernpiricaion yes estimation of esrimabie paramerers of degree one ani KG,I;ti kliown probability measure Pon (R,R) is considered. The asymptotic optimality of the these estimators, obtaining the exact risk expressions, is established. Also the results of Dalal and Phad (1983) we extended to the estimation of an estimable parametric function of an unknow probability measure P on (R² , B²)     

Selecting the t-best cells in a multinomial distribution

The problem of selecting the t-best cells in a multinomial distribution with t + k cells, k > 1, 2 <= t is considered under the fixed sample-size indifference zone approach. The least favourable configuration is derived for the usual procedure of selection, for large values of N (the sample size). The result settles Conjecture I (for large N) and Conjecture IV of Chen and Hwang (Commun. Statist. - Theory Meth. 13 (10), 1289-1298, 1984) in the affirmative.     

Detecting changes in the mean of functional observations

Summary.  Principal component analysis has become a fundamental tool of functional data analysis. It represents the functional data as X _i ( t )= μ ( t )+Σ_{1≤ l <∞} η _{i ,  l} +  v _l ( t ), where μ is the common mean, v _l are the eigenfunctions of the covariance operator and the η _{i ,  l} are the scores. Inferential procedures assume that the mean function μ ( t ) is the same for all values of i . If, in fact, the observations do not come from one population, but rather their mean changes at some point(s), the results of principal component analysis are confounded by the change(s). It is therefore important to develop a methodology to test the assumption of a common functional mean. We develop such a test using quantities which can be readily computed in the R package fda. The null distribution of the test statistic is asymptotically pivotal with a well-known asymptotic distribution. The asymptotic test has excellent finite sample performance. Its application is illustrated on temperature data from England.     

A Comparison of the Runs Test for Volatility Forecastability and the LM Test for GARCH Using Aggregated Returns

Christoffersen and Diebold (2000 Christoffersen , P. F. , Diebold , F. X. ( 2000 ). How relevant is volatility forecasting for financial risk management? Review of Economics and Statistics 82 : 12 – 22 .[Crossref] , [Google Scholar]) have introduced a runs test for forecastable volatility in aggregated returns. In this note, we compare the size and power of their runs test and the more conventional LM test for GARCH by Monte Carlo simulation. When the true daily process is GARCH, EGARCH, or stochastic volatility, the LM test has better power than the runs test for the moderate-horizon returns considered by Christoffersen and Diebold. For long-horizon returns, however, the tests have very similar power. We also consider a qualitative threshold GARCH model. For this process, we find that the runs test has greater power than the LM test. Theresults support the use of the runs test with aggregated returns.     

Confidence intervals on p(y < x) for small sample sizes

When F = G^a, confidence intervals are derived and presented in graphs for p = P(Y < X), when X and Y are independent and the sample sizes are at most 25. Also, it is demonstrated via a monte carlo simulation that this is a robust procedure, when the distributions of X and Y differ by a location parameter.