Inference for domains under imputation for missing survey data

The authors study the estimation of domain totals and means under survey‐weighted regression imputation for missing items. They use two different approaches to inference: (i) design‐based with uniform response within classes; (ii) model‐assisted with ignorable response and an imputation model. They show that the imputed domain estimators are biased under (i) but approximately unbiased under (ii). They obtain a bias‐adjusted estimator that is approximately unbiased under (i) or (ii). They also derive linearization variance estimators. They report the results of a simulation study on the bias ratio and efficiency of alternative estimators, including a complete case estimator that requires the knowledge of response indicators.     

Point estimation of the stress–strength reliability parameter for parallel system with independent and non-identical components

In this article, the estimation problem of the multicomponent stress–strength reliability parameter is considered where the stress and the strength systems have arbitrary fixed numbers of independent and non-identical parallel components. It is assumed that the distribution functions of the stress and the strength components satisfy the proportional reversed hazard rate model. The study is done in more details when the baseline distributions are exponential. Maximum likelihood and uniformly minimum variance unbiased estimators are obtained and compared. Also, Bayes and empirical Bayes estimators are discussed and Monte Carlo simulations are carried out to compare their performances.     

The Posterior Distribution of the Parameters of Component Lifetimes Based on Autopsy Data in a Shock Model

In this paper we consider a binary, monotone system whose component states are dependent through the possible occurrence of independent common shocks, i.e. shocks that destroy several components at once. The individual failure of a component is also thought of as a shock. Such systems can be used to model common cause failures in reliability analysis. The system may be a technological one, or a human being. It is observed until it fails or dies. At this instant, the set of failed components and the failure time of the system are noted. The failure times of the components are not known. These are the so-called autopsy data of the system. For the case of independent components, i.e. no common shocks, Meilijson (1981), Nowik (1990), Antoine et al . (1993) and GTsemyr (1998) discuss the corresponding identifiability problem, i.e. whether the component life distributions can be determined from the distribution of the observed data. Assuming a model where autopsy data is known to be enough for identifia bility, Meilijson (1994) goes beyond the identifiability question and into maximum likelihood estimation of the parameters of the component lifetime distributions based on empirical autopsy data from a sample of several systems. He also considers life-monitoring of some components and conditional life-monitoring of some other. Here a corresponding Bayesian approach is presented for the shock model. Due to prior information one advantage of this approach is that the identifiability problem represents no obstacle. The motivation for introducing the shock model is that the autopsy model is of special importance when components can not be tested separately because it is difficult to reproduce the conditions prevailing in the functioning system. In Gåsemyr & Natvig (1997) we treat the Bayesian approach to life-monitoring and conditional life- monitoring of components     

Topp–Leone normal distribution with application to increasing failure rate data

In this article, we proposed a new three-parameter probability distribution, called Topp–Leone normal, for modelling increasing failure rate data. The distribution is obtained by using Topp–Leone-X family of distributions with normal as a baseline model. The basic properties including moments, quantile function, stochastic ordering and order statistics are derived here. The estimation of unknown parameters is approached by the method of maximum likelihood, least squares, weighted least squares and maximum product spacings. An extensive simulation study is carried out to compare the long-run performance of the estimators. Applicability of the distribution is illustrated by means of three real data analyses over existing distributions.     

Simultaneous equivariant estimation of the parameters of matrix scale and matrix location-scale models

Eaton and Olkin (1987) discussed the problem of best equivariant estimator of the matrix scale parameter with respect to different scalar loss functions. Edwin Prabakaran and Chandrasekar (1994) developed simultaneous equivariant estimation approach and illustrated the method with examples. The problems considered in this paper are simultaneous equivariant estimation of the parameters of (i) a matrix scale model and (ii) a multivariate location-scale model. By considering matrix loss function (Klebanov, Linnik and Ruhin, 1971) a characterization of matrix minimum risk equivariant (MMRE) estimator of the matrix parameter is obtained in each case. Illustrative examples are provided in which MMRE estimators are obtained with respect to two matrix loss functions.     

Constant Stress Accelerated Life Test on a Multiple-Component Series System under Weibull Lifetime Distributions

We will discuss the reliability analysis of the constant stress accelerated life test on a series system connected with multiple components under independent Weibull lifetime distributions whose scale parameters are log-linear in the level of the stress variable. The system lifetimes are collected under Type I censoring but the components that cause the systems to fail may or may not be observed. The data are so called masked for the latter case. Maximum likelihood approach and the Bayesian method are considered when the data are masked. Statistical inference on the estimation of the underlying model parameters as well as the mean time to failure and the reliability function will be addressed. Simulation study for a three-component case shows that Bayesian analysis outperforms the maximum likelihood approach especially when the data are highly masked.     

MSE performance of the weighted average estimators consisting of shrinkage estimators

In this paper, we consider a regression model and propose estimators which are the weighted averages of two estimators among three estimators; the Stein-rule (SR), the minimum mean squared error (MMSE), and the adjusted minimum mean-squared error (AMMSE) estimators. It is shown that one of the proposed estimators has smaller mean-squared error (MSE) than the positive-part Stein-rule (PSR) estimator over a moderate region of parameter space when the number of the regression coefficients is small (i.e., 3), and its MSE performance is comparable to the PSR estimator even when the number of the regression coefficients is not so small.     

Two-stage hierarchical modeling for analysis of subpopulations in conditional distributions

In this work, we develop modeling and estimation approach for the analysis of cross-sectional clustered data with multimodal conditional distributions where the main interest is in analysis of subpopulations. It is proposed to model such data in a hierarchical model with conditional distributions viewed as finite mixtures of normal components. With a large number of observations in the lowest level clusters, a two-stage estimation approach is used. In the first stage, the normal mixture parameters in each lowest level cluster are estimated using robust methods. Robust alternatives to the maximum likelihood estimation are used to provide stable results even for data with conditional distributions such that their components may not quite meet normality assumptions. Then the lowest level cluster-specific means and standard deviations are modeled in a mixed effects model in the second stage. A small simulation study was conducted to compare performance of finite normal mixture population parameter estimates based on robust and maximum likelihood estimation in stage 1. The proposed modeling approach is illustrated through the analysis of mice tendon fibril diameters data. Analyses results address genotype differences between corresponding components in the mixtures and demonstrate advantages of robust estimation in stage 1.     

Estimation of cumulative incidence functions in competing risks studies under an order restriction

In the competing risks problem an important role is played by the cumulative incidence function (CIF), whose value at time t is the probability of failure by time t for a particular type of failure in the presence of other risks. Its estimation and asymptotic distribution theory have been studied by many. In some cases there are reasons to believe that the CIFs due to two types of failure are order restricted. Several procedures have appeared in the literature for testing for such orders. In this paper we initiate the study of estimation of two CIFs subject to a type of stochastic ordering, both when there are just two causes of failure and when there are more than two causes of failure, treating those other than the two of interest as a censoring mechanism. We do not assume independence of the two types of failure of interest; however, these are assumed to be independent of the other causes in the censored case. Weak convergence results for the estimators have been derived. It is shown that when the order restriction is strict, the asymptotic distributions are the same as those for the empirical estimators without the order restriction. Thus we get the restricted estimators “free of charge”, at least in the asymptotic sense. When the two CIFs are equal, the asymptotic MSE is reduced by using the order restriction. For finite sample sizes simulations seem to indicate that the restricted estimators have uniformly smaller MSEs than the unrestricted ones in all cases.     

Weighted exponential distribution: properties and different methods of estimation

In this article, we investigate various properties and methods of estimation of the Weighted Exponential distribution. Although, our main focus is on estimation (from both frequentist and Bayesian point of view) yet, the stochastic ordering, the Bonferroni and the Lorenz curves, various entropies and order statistics are derived first time for the said distribution. Different types of loss functions are considered for Bayesian estimation. Furthermore, the Bayes estimators and their respective posterior risks are computed and compared using Gibbs sampling. The different reliability characteristics including hazard function, stress and strength analysis, and mean residual life function are also derived. Monte Carlo simulations are performed to compare the performances of the proposed methods of estimation and two real data sets have been analysed for illustrative purposes.     

Estimation of R=P[Y<X] for three-parameter generalized Rayleigh distribution

Surles and Padgett [Inference for reliability and stress–strength for a scaled Burr type X distribution. Lifetime Data Anal. 2001;7:187–200] introduced a two-parameter Burr-type X distribution, which can be described as a generalized Rayleigh distribution. In this paper, we consider the estimation of the stress–strength parameter R=P[Y<X], when X and Y are both three-parameter generalized Rayleigh distributions with the same scale and locations parameters but different shape parameters. It is assumed that they are independently distributed. It is observed that the maximum-likelihood estimators (MLEs) do not exist, and we propose a modified MLE of R. We obtain the asymptotic distribution of the modified MLE of R, and it can be used to construct the asymptotic confidence interval of R. We also propose the Bayes estimate of R and the construction of the associated credible interval based on importance sampling technique. Analysis of two real data sets, (i) simulated and (ii) real, have been performed for illustrative purposes.     

Nonparametric Estimation of the Average Availability

The average availability of a repairable system is the expected proportion of time that the system is operating in the interval [0, t]. The present article discusses the nonparametric estimation of the average availability when (i) the data on ‘n’ complete cycles of system operation are available, (ii) the data are subject to right censorship, and (iii) the process is observed upto a specified time ‘T’. In each case, a nonparametric confidence interval for the average availability is also constructed. Simulations are conducted to assess the performance of the estimators.     

Estimation of p{y p> max(y1, y2, …, yp-1)}in theexponential case

This paper deals with the derivation of (i) the MLE (ii) the MVUE (iii) a Bayes estimator of the probability in the title, for the case p = 2. Simulation studies are carried out to compare these estimators. The results suggest that the MLE and the Bayes estimator are biased and the Bayes estimator have the smallest MSE. In the general case, explicit expression for the probability in the title is derived and the MLE and Bayes estimator are obtained. A general method of deriving the MVUE is pointed out. Because of the simulation studies for p = 2 it is recommended that the Bayes or predictive estimator should be used.     

A sufficient condition for the MSE dominance of the positive-part shrinkage estimator when each individual regression coefficient is estimated in a misspecified linear regression model

In this paper, assuming that there exist omitted explanatory variables in the specified model, we derive the exact formula for the mean squared error (MSE) of a general family of shrinkage estimators for each individual regression coefficient. It is shown analytically that when our concern is to estimate each individual regression coefficient, the positive-part shrinkage estimators have smaller MSE than the original shrinkage estimators under some conditions even when the relevant regressors are omitted. Also, by numerical evaluations, we showed the effects of our theorem for several specific cases. It is shown that the positive-part shrinkage estimators have smaller MSE than the original shrinkage estimators for wide region of parameter space even when there exist omitted variables in the specified model.     

Simultaneous estimation of the multivariate normal mean under balanced loss function

This paper considers simultaneous estimation of multivariate normal mean vector using Zellner's(1994) balanced loss function which is defined as follows:

where 0 < w < 1 and for i = 1,…,p and j = 1,…,n, X_ij is distributed as normal with mean θi and variance 1. It is shown that the sample mean, X, is admissible when p <3. For p ≥3, we obtain that James-Stein type estimator which has uniformly smaller risk than that of sample mean X.     

On improved estimation of a gamma shape parameter

This paper deals with improved estimation of a gamma shape parameter from a decision-theoretic point of view. First we study the second-order properties of three estimators – (i) the maximum-likelihood estimator (MLE), (ii) a bias corrected version of the MLE, and (iii) an improved version (in terms of mean squared error) of the MLE. It is shown that all the three estimators mentioned above are second-order inadmissible. Next, we obtain superior estimators which are second order better than the above three estimators. Simulation results are provided to study the relative risk improvement of each improved estimator over the MLE.     

Estimation of poisson means under weighted squared error loss

In estimating the means of several independent Poisson distributions, we show that the maximum likelihood estimator is inadmissible when general weighted squared error loss is the criterion. Using this result, we extend the known results on estimation of several Poisson means (Peng 1975, Hudson 1978) to the case where possibly more than one observation is taken from each Poisson distribution and the samples are not necessarily of the same size.     

Ordering Properties of Systems with Two Dependent Components

The study of systems with dependent components from a reliability point of view is a very important topic. However, the majority of the articles study the case of independent components. In this article, we study how the dependency influences the performance of the system. We extend some comparison results obtained in the case of independent components to the case of two dependent components. We show that the more diverse the exponential parameters of the two components, the stronger (weaker) the parallel (series) system in the stochastic ordering. We apply our general results to some common bivariate models in the reliability theory.     

A new estimator combining the ridge regression and the restricted least squares methods of estimation

It is well-known in the literature on multicollinearity that one of the major consequences of multicollinearity on the ordinary least squares estimator is that the estimator produces large sampling variances, which in turn might inappropriately lead to exclusion of otherwise significant coefficients from the model. To circumvent this problem, two accepted estimation procedures which are often suggested are the restricted least squares method and the ridge regression method. While the former leads to a reduction in the sampling variance of the estimator, the later ensures a smaller mean square error value for the estimator. In this paper we have proposed a new estimator which is based on a criterion that combines the ideas underlying these two estimators. The standard properties of this new estimator have been studied in the paper. It has also been shown that this estimator is superior to both the restricted least squares as well as the ordinary ridge regression estimators by the criterion of mean sauare error of the estimator of the regression coefficients when the restrictions are indeed correct. The conditions for superiority of this estimator over the other two have also been derived for the situation when the restrictions are not correct.     

Additional comparisons of designs to estimate variance components in a two-way classification model

In this paper, the research of Muse and Anderson is extended to include additional comparisons of designs, featuring planned unbalance, for the estimation of variance components in a two-way cross classification model. Their results are extended to Include the following: (i) a small sample study of the original off-diagonal (OD) design and (ii) an asymptotic maximum likelihood investigation of three modifica-tions of the balanced diagonal rectangles (BD) design and one modification of the 01) design to permit the estimation of row, column, interaction and error variance components. Also a general iterative least.