A Frailty Model in Crossover Studies with Time-to-Event Response Variable

In this article, we develop a model to study treatment, period, carryover, and other applicable effects in a crossover design with a time-to-event response variable. Because time-to-event outcomes on different treatment regimens within the crossover design are correlated for an individual, we adopt a proportional hazards frailty model. If the frailty is assumed to have a gamma distribution, and the hazard rates are piecewise constant, then the likelihood function can be determined via closed-form expressions. We illustrate the methodology via an application to a data set from an asthma clinical trial and run simulations that investigate sensitivity of the model to data generated from different distributions.     

Bayes Estimation Based on Joint Progressive Type II Censored Data Under LINEX Loss Function

In the lifetime experiments, the joint censoring scheme is useful for planning comparative purposes of two identical products manufactured coming from different lines. In this article, we will confine ourselves to the data obtained by conducting a joint progressive Type II censoring scheme on the basis of the two combined samples selected from the two lines. Moreover, it is supposed that the distributions of lifetimes of the two products satisfy in a proportional hazard model. A general form for the distributions is considered, and we tackle the problem of obtaining Bayes estimates under the squared error and linear-exponential (LINEX) loss functions. As a special case, the Weibull distribution is discussed in more detail. Finally, the estimated risks of the various estimators obtained are compared using the Monte Carlo method.     

Survival models with two phases and length biased sampling

A survival model is presented in which all patients go through a first phase of disease; some then die and the remainder progress to a second phase of disease. The data observed are the path taken and the total sojourn time in the system but not the time, if ever, at which the second phase is entered. The sojourn time in each phase is assumed to be exponentially distributed with possibly different rates for the two phases. Themodel describes serious diseases that progress through one or two phases, and can be extended to multiple phases. The model is extended to account for several length-biased sampling situations. Censoring is considered in all models. Maximum like lihood estimates for the parameters involved exist, are consistent and are a symptotically normal. One of the proposed models is applied to data from the Veterans Administration involving a study of coronary arterial occlusive disease.     

Lower Confidence Limits on the Generalized Exponential Distribution Percentiles Under Progressive Type-I Interval Censoring

In industrial life test and survival analysis, the percentile estimation is always a practical issue with lower confidence bound required for maintenance purpose. Sampling distributions for the maximum likelihood estimators of percentiles are usually unknown. Bootstrap procedures are common ways to estimate the unknown sampling distributions. Five parametric bootstrap procedures are proposed to estimate the confidence lower bounds on maximum likelihood estimators for the generalized exponential (GE) distribution percentiles under progressive type-I interval censoring. An intensive simulation is conducted to evaluate the performances of proposed procedures. Finally, an example of 112 patients with plasma cell myeloma is given for illustration.     

Maximum likelihood estimation in the proportional hazards model of random censorship

The maximum likelihood estimator (MLE) for the survival function S_Tunder the proportional hazards model of censorship is derived and shown to differ from the Abdushukurov-Cheng-Lin estimator when the class of allowable distributions includes all continuous and discrete distributions. The estimators are compared via an example. The MLE is calculated using a Newton-Raphson iterative procedure and implemented via a FORTRAN algorithm.     

Computing maximum likelihood estimates from type II doubly censored exponential data

It is well-known that, under Type II double censoring, the maximum likelihood (ML) estimators of the location and scale parameters, θ and δ, of a twoparameter exponential distribution are linear functions of the order statistics. In contrast, when θ is known, theML estimator of δ does not admit a closed form expression. It is shown, however, that theML estimator of the scale parameter exists and is unique. Moreover, it has good large-sample properties. In addition, sharp lower and upper bounds for this estimator are provided, which can serve as starting points for iterative interpolation methods such as regula falsi. Explicit expressions for the expected Fisher information and Cramér-Rao lower bound are also derived. In the Bayesian context, assuming an inverted gamma prior on δ, the uniqueness, boundedness and asymptotics of the highest posterior density estimator of δ can be deduced in a similar way. Finally, an illustrative example is included.     

Inference of stress-strength for the Type-II generalized logistic distribution under progressively Type-II censored samples

This article studies the estimation of the reliability R = P[Y < X] when X and Y come from two independent generalized logistic distributions of Type-II with different parameters, based on progressively Type-II censored samples. When the common scale parameter is unknown, the maximum likelihood estimator and its asymptotic distribution are proposed. The asymptotic distribution is used to construct an asymptotic confidence interval of R. Bayes estimator of R and the corresponding credible interval using the Gibbs sampling technique have been proposed too. Assuming that the common scale parameter is known, the maximum likelihood estimator, uniformly minimum variance unbiased estimator, Bayes estimation, and confidence interval of R are extracted. Monte Carlo simulations are performed to compare the different proposed methods. Analysis of a real dataset is given for illustrative purposes. Finally, methods are extended for proportional hazard rate models.     

A new flexible hazard rate distribution: Application and estimation of its parameters

We introduce a new class of flexible hazard rate distributions which have constant, increasing, decreasing, and bathtub-shaped hazard function. This class of distributions obtained by compounding the power and exponential hazard rate functions, which is called the power-exponential hazard rate distribution and contains several important lifetime distributions. We obtain some distributional properties of the new family of distributions. The estimation of parameters is obtained by using the maximum likelihood and the Bayesian methods under squared error, linear-exponential, and Stein’s loss functions. Also, approximate confidence intervals and HPD credible intervals of parameters are presented. An application to real dataset is provided to show that the new hazard rate distribution has a better fit than the other existing hazard rate distributions and some four-parameter distributions. Finally , to compare the performance of proposed estimators and confidence intervals, an extensive Monte Carlo simulation study is conducted.     

Revisiting Beal's Confidence Intervals for the Difference of Two Binomial Proportions

Confidence interval construction for the difference of two independent binomial proportions is a well-known problem with a full panoply of proposed solutions. In this paper, we focus largely on the family of intervals proposed by Beal (1987 Beal , S. ( 1987 ). Asymptotic confidence intervals for the difference between two binomial parameters for use with small samples . Biometrics 43 : 941 – 950 . [CSA] [CROSSREF] [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). This family, which includes the Haldane and Jeffreys–Perks intervals as special cases, assumes a symmetric prior distribution for the population proportions p ₁ and p ₂. We propose new methods that allow the currently observed data to set the prior distribution by taking a parametric empirical-Bayes approach; in addition, we also provide an investigation of the new interval' behaviors in small-sample situations. Unlike other solutions, our intervals can be used adaptively for experiments conducted in multiple stages over time. We illustrate this notion using data from an Argentinean study involving the Mal Rio Cuarto virus and its transmission to susceptible maize crops.     

Estimation of Survival Function in Cox Regression Model Under Random Censoring From Both Sides

In this article we present three types of parametric–non parametric estimators for conditional survival function in Cox proportional hazards regression model when the lifetime of interest is subjected to random censorship from both sides. We prove consistency and asymptotic normality of estimators.     

A Note on Nonunique MLE's and Sufficient Statistics

A simple classroom example is presented that illustrates nonunique maximum likelihood estimates and the fact that such estimates need not be functions of sufficient statistics.     

Non parametric hazard estimation with dependent censoring using penalized likelihood and an assumed copula

This article introduces a novel non parametric penalized likelihood hazard estimation when the censoring time is dependent on the failure time for each subject under observation. More specifically, we model this dependence using a copula, and the method of maximum penalized likelihood (MPL) is adopted to estimate the hazard function. We do not consider covariates in this article. The non negatively constrained MPL hazard estimation is obtained using a multiplicative iterative algorithm. The consistency results and the asymptotic properties of the proposed hazard estimator are derived. The simulation studies show that our MPL estimator under dependent censoring with an assumed copula model provides a better accuracy than the MPL estimator under independent censoring if the sign of dependence is correctly specified in the copula function. The proposed method is applied to a real dataset, with a sensitivity analysis performed over various values of correlation between failure and censoring times.     

Maximum Likelihood Estimation in Compound Geometric Failure Model with Changing Parameters From Type-I Two-Stage Progressively Censored and Group Censored Samples

Boardman and Kendell (1970 Boardman , T. J. , Kendell , P. J. ( 1970 ). Estimation in compound failure models . Technometrics 12 : 891 – 908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimation with respect to Type-I censoring when an item is subjected to only one of the two causes of failure assuming exponential model. Patel and Gajjar (1992 Patel , M. N. , Gajjar , A. V. ( 1992 ). Maximum likelihood estimation in compound exponential failure model with changing failure rates from Type-I progressively censored and group censored samples . Commun. Statist. Theor. Meth. 21 ( 10 ): 2899 – 2908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered extension of the Boardman and Kendell's results in case of two-stage progressive censoring. Here we have considered geometric competing risk failure model with two independent causes of failures. Maximum likelihood estimation of the parameters is carried out using Type-I two-stage progressively censored and group censored samples. Asymptotic standard errors of the estimators are obtained for both the cases. Two illustrative examples are cited for ungroup and group competing risk models.     

Two-stage informative cluster sampling—estimation and prediction with applications for small-area models

This paper considers the effects of informative two-stage cluster sampling on estimation and prediction. The aims of this article are twofold: first to estimate the parameters of the superpopulation model for two-stage cluster sampling from a finite population, when the sampling design for both stages is informative, using maximum likelihood estimation methods based on the sample-likelihood function; secondly to predict the finite population total and to predict the cluster-specific effects and the cluster totals for clusters in the sample and for clusters not in the sample. To achieve this we derive the sample and sample-complement distributions and the moments of the first and second stage measurements. Also we derive the conditional sample and conditional sample-complement distributions and the moments of the cluster-specific effects given the cluster measurements. It should be noted that classical design-based inference that consists of weighting the sample observations by the inverse of sample selection probabilities cannot be applied for the prediction of the cluster-specific effects for clusters not in the sample. Also we give an alternative justification of the Royall [1976. The linear least squares prediction approach to two-stage sampling. Journal of the American Statistical Association 71, 657–664] predictor of the finite population total under two-stage cluster population. Furthermore, small-area models are studied under informative sampling.     

Extended Hazard Regression Model for Reliability and Survival Analysis

We propose an extended hazard regression model which allows the spread parameter to be dependent on covariates. This allows a broad class of models which includes the most common hazard models, such as the proportional hazards model, the accelerated failure time model and a proportional hazards/accelerated failure time hybrid model with constant spread parameter. Simulations based on sub-classes of this model suggest that maximum likelihood performs well even when only small or moderate-size data sets are available and the censoring pattern is heavy. The methodology provides a broad framework for analysis of reliability and survival data. Two numerical examples illustrate the results.     

Estimation of a Change Point in a Hazard Function Based on Censored Data

The hazard function plays an important role in reliability or survival studies since it describes the instantaneous risk of failure of items at a time point, given that they have not failed before. In some real life applications, abrupt changes in the hazard function are observed due to overhauls, major operations or specific maintenance activities. In such situations it is of interest to detect the location where such a change occurs and estimate the size of the change. In this paper we consider the problem of estimating a single change point in a piecewise constant hazard function when the observed variables are subject to random censoring. We suggest an estimation procedure that is based on certain structural properties and on least squares ideas. A simulation study is carried out to compare the performance of this estimator with two estimators available in the literature: an estimator based on a functional of the Nelson-Aalen estimator and a maximum likelihood estimator. The proposed least squares estimator tums out to be less biased than the other two estimators, but has a larger variance. We illustrate the estimation method on some real data sets.     

Stress-Strength Reliability of a Two-Parameter Bathtub-shaped Lifetime Distribution Based on Progressively Censored Samples

Based on progressively Type-II censored samples, this article deals with inference for the stress-strength reliability R = P(Y < X) when X and Y are two independent two-parameter bathtub-shape lifetime distributions with different scale parameters, but having the same shape parameter. Different methods for estimating the reliability are applied. The maximum likelihood estimate of R is derived. Also, its asymptotic distribution is used to construct an asymptotic confidence interval for R. Assuming that the shape parameter is known, the maximum likelihood estimator of R is obtained. Based on the exact distribution of the maximum likelihood estimator of R an exact confidence interval of that has been obtained. The uniformly minimum variance unbiased estimator are calculated for R. Bayes estimate of R and the associated credible interval are also got under the assumption of independent gamma priors. Monte Carlo simulations are performed to compare the performances of the proposed estimators. One data analysis has been performed for illustrative purpose. Finally, we will generalize this distribution to the proportional hazard family with two parameters and derive various estimators in this family.     

A complementary generalized linear failure rate-geometric distribution

In this article, we shall attempt to introduce a new class of lifetime distributions, which enfolds several known distributions such as the generalized linear failure rate distribution and covers both positive as well as negative skewed data. This new four-parameter distribution allows for flexible hazard rate behavior. Indeed, the hazard rate function here can be increasing, decreasing, bathtub-shaped, or upside-down bathtub-shaped. We shall first study some basic distributional properties of the new model such as the cumulative distribution function, the density of the order statistics, their moments, and Rényi entropy. Estimation of the stress-strength parameter as an important reliability property is also studied. The maximum likelihood estimation procedure for complete and censored data and Bayesian method are used for estimating the parameters involved. Finally, application of the new model to three real datasets is illustrated to show the flexibility and potential of the new model compared to rival models.     

Exact Bayesian Variable Sampling Plans for the Exponential Distribution Based on Type-I and Type-II Hybrid Censored Samples

A sampling plan with a polynomial loss function for the exponential distribution is considered. From the distribution of the maximum likelihood estimator of the mean of an exponential distribution based on Type-I and Type-II hybrid censored samples, we obtain an explicit expression for the Bayes risk of a sampling plan with a quadratic loss function. Some numerical examples and comparisons are given to illustrate the effectiveness of the proposed method, and a robustness study reveals that the proposed optimal sampling plans are quite robust.     

Estimation using hybrid censored data from a generalized inverted exponential distribution

ABSTRACT

We consider point and interval estimation of the unknown parameters of a generalized inverted exponential distribution in the presence of hybrid censoring. The maximum likelihood estimates are obtained using EM algorithm. We then compute Fisher information matrix using the missing value principle. Bayes estimates are derived under squared error and general entropy loss functions. Furthermore, approximate Bayes estimates are obtained using Tierney and Kadane method as well as using importance sampling approach. Asymptotic and highest posterior density intervals are also constructed. Proposed estimates are compared numerically using Monte Carlo simulations and a real data set is analyzed for illustrative purposes.