A general algorithm for computing simultaneous prediction intervals for the (log)-location-scale family of distributions

Making predictions of future realized values of random variables based on currently available data is a frequent task in statistical applications. In some applications, the interest is to obtain a two-sided simultaneous prediction interval (SPI) to contain at least k out of m future observations with a certain confidence level based on n previous observations from the same distribution. A closely related problem is to obtain a one-sided upper (or lower) simultaneous prediction bound (SPB) to exceed (or be exceeded) by at least k out of m future observations. In this paper, we provide a general approach for computing SPIs and SPBs based on data from a particular member of the (log)-location-scale family of distributions with complete or right censored data. The proposed simulation-based procedure can provide exact coverage probability for complete and Type II censored data. For Type I censored data, our simulation results show that our procedure provides satisfactory results in small samples. We use three applications to illustrate the proposed simultaneous prediction intervals and bounds.     

NONPARAMETRIC RELIABILITY ESTIMATION BASED ON A FEW ORDERED OBSERVATIONS

The usual one-sided Kolmogorov-Smirnov distance is generalized to obtain an improved lower confidence region for the extreme left tail of the reliability function based on k observations in a “k out of n censored” plan. Finite sample and asymptotic critical values necessary for implementation are given. The two numerical comparisons with existing parametric procedures for the case of complete or censored samples demonstrate the applicability of the proposed nonparametric procedure.     

Parameter estimation for ihe birnbaum-saunders distribution based on symmetrically censored samples

Ihe Bimbaum-Saunders distribution was derived to model fatigue life. Frequently, it becomes necessary to stop a life testing process before all the test items have failed. Therefore, estimation procedures need to be developed for use when censoring occurs. In this article, we have derived estimators for the parameters of this distribution which may be used for complete samples or Type II symmetrically censored samples A simulation study was also conducted to examine the performance of these estimators.     

A product-limit estimator for use with length-biased data

The following life-testing situation is considered. At some time in the distant past, n objects, from a population with life distribution F, were put in use; whenever an object failed, it was promptly replaced. At some time τ, long after the start of the process, a statistician starts observing the n objects in use at that time; he knows the age of each of those n objects, and observes each of them for a fixed length of time? ∞, or until failure, whichever occurs first. In the case where T is finite, some of the observations may be censored; in the case where T =∞, there is no censoring. The total life of an object in use at time ∞ is a length-biased observation from F. A nonparametric estimator of the (cumulative) hazard function is proposed, and is used to construct an estimator of F which is of the product-limit type. Strong uniform consistency results (for n → ∞) are obtained. An “Aalen-Johansen” identity, satisfied by any pair of life distributions and their (cumulative) hazard functions, is used in obtaining rate-of-convergence results.     

Influence diagnostics for censored regression models with autoregressive errors

Observations collected over time are often autocorrelated rather than independent, and sometimes include observations below or above detection limits (i.e. censored values reported as less or more than a level of detection) and/or missing data. Practitioners commonly disregard censored data cases or replace these observations with some function of the limit of detection, which often results in biased estimates. Moreover, parameter estimation can be greatly affected by the presence of influential observations in the data. In this paper we derive local influence diagnostic measures for censored regression models with autoregressive errors of order p (hereafter, AR(p)‐CR models) on the basis of the Q‐function under three useful perturbation schemes. In order to account for censoring in a likelihood‐based estimation procedure for AR(p)‐CR models, we used a stochastic approximation version of the expectation‐maximisation algorithm. The accuracy of the local influence diagnostic measure in detecting influential observations is explored through the analysis of empirical studies. The proposed methods are illustrated using data, from a study of total phosphorus concentration, that contain left‐censored observations. These methods are implemented in the R package ARCensReg.     

ON TESTING THE EQUALITY OF LOCATION PARAMETERS IN k CENSORED EXPONENTIAL DISTRIBUTIONS

Kumar and Patel (1971) have considered the problem of testing the equality of location parameters of two exponential distributions on the basis of samples censored from above, when the scale parameters are the same and unknown. The test proposed by them is shown to be biased for n₁≠n₂, while for n₁=n₂ the test possesses the property of monotonicity and is equivalent to the likelihood ratio test, which is considered by Epstein and Tsao (1953) and Dubey (1963a, 1963b). Epstein and Tsao state that the test is unbiased. We may note that when the scale parameters of k exponential distributions are unknown the problem of testing the equality of location parameters is reducible to that of testing the equality of parameters in k rectangular populations for which a test and its power function were given by Khatri (1960, 1965); Jaiswal (1969) considered similar problems in his thesis. Here we extend the problem of testing the equality of k exponential distributions on the basis of samples censored from above when the scale parameters are equal and unknown, and we establish the likelihood ratio test (LET) and the union-intersection test (UIT) procedures. Using the results previously derived by Jaiswal (1969), we obtain the power function for the LET and for k= 2 show that the test possesses the property of monotonicity. The power function of the UIT is also given.     

Control chart for monitoring the Weibull shape parameter under two competing risks

ABSTRACT

In this paper, we propose a control chart to monitor the Weibull shape parameter where the observations are censored due to competing risks. We assume that the failure occurs due to two competing risks that are independent and follow Weibull distribution with different shape and scale parameters. The control charts are proposed to monitor one or both of the shape parameters of competing risk distributions and established based on the conditional expected values. The proposed control chart for both shape parameters is used in certain situations and allows to monitor both shape parameters in only one chart. The control limits depend on the sample size, number of failures due to each risk and the desired stable average run length (ARL). We also consider the estimation problem of the target parameters when the Phase I sample is incomplete. We assumed that some of the products that fail during the life testing have a cause of failure that is only known to belong to a certain subset of all possible failures. This case is known as masking. In the presence of masking, the expectation-maximization (EM) algorithm is proposed to estimate the parameters. For both cases, with and without masking, the behaviour of ARLs of charts is studied through the numerical methods. The influence of masking on the performance of proposed charts is also studied through a simulation study. An example illustrates the applicability of the proposed charts.     

Control percentile test procedures for censored data

In life testing and survival analyses which involve the use of expensive equipment the cost of continuing an experiment until all the items on test have failed can be quite high. In these situations it is reasonable to make a statistical test when a pre-specified percentile, e.g. median of the control group has been observed. This article adapts some existing procedures for complete samples to randomly censored data. The results of Lo and Singh (1985) who extended the Bahadur representation of quantiles to the censored case enable us to use the methods of Gastwirth (1968) and Hettmansperger (1973) which were based on Bahadur's result to extend the procedures of Mathisen (1943), Gart (1963) and Slivka (1970).The large sample efficiency of the control median test is the same as that of Brookmeyer and Crowley's (1982) extension of the usual median test. For the two-sample shift problem with observations following the double-exponential law, the median remains the optimum percentile to use until the censoring becomes quite heavy. On the other hand, in the two-sample scale parameter problem for data from an exponential distribution the percentile (80th in the uncensored case) yielding the asymptotically most powerful test in the family of control percentile tests no longer is optimum. The effect becomes noticeable when 25% or more of the data is censored.     

Nonparametric predictive comparison of lifetime data under progressive censoring

In reliability and lifetime testing, comparison of two groups of data is a common problem. It is often attractive, or even necessary, to make a quick and efficient decision in order to save time and costs. This paper presents a nonparametric predictive inference (NPI) approach to compare two groups, say X and Y, when one (or both) is (are) progressively censored. NPI can easily be applied to different types of progressive censoring schemes. NPI is a statistical approach based on few assumptions, with inferences strongly based on data and with uncertainty quantified via lower and upper probabilities. These inferences consider the event that the lifetime of a future unit from Y is greater than the lifetime of a future unit from X.     

Recurrence relations for moments of progressively censored order statistics from logistic distribution with applications to inference

In this paper, we establish several recurrence relations for the single and product moments of progressively Type-II right censored order statistics from a logistic distribution. The use of these relations in a systematic manner allows us to compute all the means, variances and covariances of progressively Type-II right censored order statistics from the logistic distribution for all sample sizes n, effective sample sizes m, and all progressive censoring schemes (R₁,…,R_m). The results established here generalize the corresponding results for the usual order statistics due to [Shah, 1966] and [Shah, 1970]. These moments are then utilized to derive best linear unbiased estimators of the location and scale parameters of the logistic distribution. A comparison of these estimators with the maximum likelihood estimations is then made. The best linear unbiased predictors of censored failure times are briefly discussed. Finally, an illustrative example is presented.     

A pseudo-complete sample technique for estimation from censored samples

An iterative procedure is presented whereby singly right censored samples are transformed into pseudo-complete samples. This involves the use of order statistics to calculate expected “complete” values of the censored observations. Estimates of the distribution parameters are then calculated by employing standard complete sample estimators. This procedure is applicable to various types of distributions.     

Multiple imputation of censored survival data in the presence of missing covariates using restricted mean survival time

Missing covariates data with censored outcomes put a challenge in the analysis of clinical data especially in small sample settings. Multiple imputation (MI) techniques are popularly used to impute missing covariates and the data are then analyzed through methods that can handle censoring. However, techniques based on MI are available to impute censored data also but they are not much in practice. In the present study, we applied a method based on multiple imputation by chained equations to impute missing values of covariates and also to impute censored outcomes using restricted survival time in small sample settings. The complete data were then analyzed using linear regression models. Simulation studies and a real example of CHD data show that the present method produced better estimates and lower standard errors when applied on the data having missing covariate values and censored outcomes than the analysis of the data having censored outcome but excluding cases with missing covariates or the analysis when cases with missing covariate values and censored outcomes were excluded from the data (complete case analysis).     

Hypothesis Testing in Functional Comparative Calibration Models

The purpose of this article is to investigate hypothesis testing in functional comparative calibration models. Wald type statistics are considered which are asymptotically distributed according to the chi-square distribution. The statistics are based on maximum likelihood, corrected score approach, and method of moment estimators of the model parameters, which are shown to be consistent and asymptotically normally distributed. Results of analytical and simulation studies seem to indicate that the Wald statistics based on the method of moment estimators and the corrected score estimators are, as expected, less efficient than the Wald type statistic based on the maximum likelihood estimators for small n. Wald statistic based on moment estimators are simpler to compute than the other Wald statistics tests and their performance improves significantly as n increases. Comparisons with an alternative F statistics proposed in the literature are also reported.     

Une application des méthodes séquentielles aux tests de permutations

To carry out a permutation test we have to examine the n! permutations of the observations. In order to make the permutation test feasible, Dwass (1957) proposed to examine only a sample of these permutations. With the help of sequential methods, we obtain a test which is never less efficient than that proposed by Dwass or the permutation test itself, in the sense that it is as powerful and never requires more permutations to make a decision. In practice, we can expect to gain much efficiency.     

Note on the Cramér-Rao inequality in the nonregular case: the family of uniform distributions

We consider the family of uniform distributions with range of unit length. The main result of this note asserts that the average variance of any unbiased estimator of the midpoint of the range is not less than (2(n+1))(n+2))^-1 and this lower bound is sharp. The proof is based upon a nonregular version of the Cramér-Rao inequality.     

Monitoring processes with data censored owing to competing risks by using exponentially weighted moving average control charts

In industry, process monitoring is widely employed to detect process changes rapidly. However, in some industrial applications observations are censored. For example, when testing breaking strengths and failure times often a limited stress test is performed. With censored observations, a direct application of traditional monitoring procedures is not appropriate. When the censoring occurs due to competing risks, we propose a control chart based on conditional expected values to detect changes in the mean strength. To protect against possible confounding caused by changes in the mean of the censoring mechanism we also suggest a similar chart to detect changes in the mean censoring level. We provide an example of monitoring bond strength to illustrate the application of this methodology.     

Comparing the Fisher information in record data and random observations

We compare the Fisher information (FI) contained in the firstn record values and record times with the FI inn i. i. d. observations. General results are established for exponential family and Weibull type setups, and a summary table is provided listing several common distributions. We show that the FI in record data improves notably once the record times are included, often changing from being less to being equal or greater than the FI in a random sample of the same size. The behavior in the Weibull case is surprising. There it depends onn, whether the record or the i.i. d. observations have more FI. We propose new estimators based on record data. The results may be of interest in some life testing situations. Supported in part by Fondo Nacional de Desarrollo Cientifico y Tecnologico (FONDECYT) grant # 1010222 of Chile.     

Inference about parameters in Binomial-Poisson distribution with additional incomplete data

ABSTRACT

This article considers the distribution of Binomial-Poisson random vector which has two components and includes two parameters: one is the rate of a Poisson distribution, the other is the proportion in a Binomial distribution. The inference about the two parameters is usually made based on only paired observations. However, the number of paired observations is, in general, not large enough because of either technical difficulty or budget limitation, and so one can not make efficient inferences with only paired data. Instead, it is often much easier and not too costly to have incomplete observation on only one component independently. In this article we will combine both the paired complete data and unpaired incomplete data for estimating the two parameters. The performances of various estimators are compared both analytically and numerically. It is observed that fully using the unpaired incomplete data can always improve the inference, and the improvement is very significant in the case when there are only a few paired complete observations.     

Bayes and Empirical Bayes Two-Stage Procedures for Sampling Inspection

A batch of M items is inspected for defectives. Suppose there are d defective items in the batch. Let d ₀ be a given standard used to evaluate the quality of the population where 0 < d ₀ < M. The problem of testing H ₀: d < d ₀ versus H ₁: d ≥ d ₀ is considered. It is assumed that past observations are available when the current testing problem is considered. Accordingly, the empirical Bayes approach is employed. By using information obtained from the past data, an empirical Bayes two-stage testing procedure is developed. The associated asymptotic optimality is investigated. It is proved that the rate of convergence of the empirical Bayes two-stage testing procedure is of order O (exp(? c? n)), for some constant c? > 0, where n is the number of past observations at hand.