Practical considerations when analyzing discrete survival times using the grouped relative risk model

The grouped relative risk model (GRRM) is a popular semi-parametric model for analyzing discrete survival time data. The maximum likelihood estimators (MLEs) of the regression coefficients in this model are often asymptotically efficient relative to those based on a more restrictive, parametric model. However, in settings with a small number of sampling units, the usual properties of the MLEs are not assured. In this paper, we discuss computational issues that can arise when fitting a GRRM to small samples, and describe conditions under which the MLEs can be ill-behaved. We find that, overall, estimators based on a penalized score function behave substantially better than the MLEs in this setting and, in particular, can be far more efficient. We also provide methods of assessing the fit of a GRRM to small samples.     

Exact Inference for a Simple Step-stress Model with Generalized Type-I Hybrid Censored Data from the Exponential Distribution

In this article, the simple step-stress model is considered based on generalized Type-I hybrid censored data from the exponential distribution. The maximum likelihood estimators (MLEs) of the unknown parameters are derived assuming a cumulative exposure model. We then derive the exact distributions of the MLEs of the parameters using conditional moment generating functions. The Bayesian estimators of the parameters are derived and then compared with the MLEs. We also derive confidence intervals for the parameters using these exact distributions, asymptotic distributions of the MLEs, Bayesian, and the parametric bootstrap methods. The problem of determining the optimal stress-changing point is discussed and the MLEs of the pth quantile and reliability functions at the use condition are obtained. Finally, Monte Carlo simulation and some numerical results are presented for illustrating all the inferential methods developed here.     

Bivariate log Birnbaum–Saunders distribution

Univariate Birnbaum–Saunders distribution has received a considerable amount of attention in recent years. Rieck and Nedelman (A log-linear model for the Birnbaum–Saunders distribution. Technometrics, 1991;33:51–60) introduced a log Birnbaum–Saunders distribution. The main aim of this paper is to introduce bivariate log Birnbaum–Saunders distribution. The proposed model is symmetric and it has five parameters. It can be obtained using Gaussian copula. Different properties can be obtained using copula structure. It is observed that the maximum likelihood estimators (MLEs) cannot be obtained explicitly. Two-dimensional profile likelihood approach may be adopted to compute the MLEs. We propose some alternative estimators also, which can be obtained quite conveniently. The analysis of one data set is performed for illustrative purposes. Finally, it is observed that this model can be used as a bivariate log-linear model, and its multivariate generalization is also quite straight forward.     

On MLEs of the parameters of a modified Weibull distribution based on record values

The modified Weibull distribution can be used quite effectively to model complex data from mechanical engineering or survival analysis studies that posses a monotonic or a bathtub-shape hazard rate. In this paper, we study the MLEs of the parameters of a modified Weibull distribution model in the presence of upper kth record values. The existence and uniqueness of the MLEs are proven in this case. Real data analysis is performed for illustrative purposes.     

Analysis of hybrid censored competing risks data

In this paper, we consider the analysis of hybrid censored competing risks data, based on Cox's latent failure time model assumptions. It is assumed that lifetime distributions of latent causes of failure follow Weibull distribution with the same shape parameter, but different scale parameters. Maximum likelihood estimators (MLEs) of the unknown parameters can be obtained by solving a one-dimensional optimization problem, and we propose a fixed-point type algorithm to solve this optimization problem. Approximate MLEs have been proposed based on Taylor series expansion, and they have explicit expressions. Bayesian inference of the unknown parameters are obtained based on the assumption that the shape parameter has a log-concave prior density function, and for the given shape parameter, the scale parameters have Beta–Gamma priors. We propose to use Markov Chain Monte Carlo samples to compute Bayes estimates and also to construct highest posterior density credible intervals. Monte Carlo simulations are performed to investigate the performances of the different estimators, and two data sets have been analysed for illustrative purposes.     

Weighted Marshall–Olkin bivariate exponential distribution

Recently, Gupta and Kundu [R.D. Gupta and D. Kundu, A new class of weighted exponential distributions, Statistics 43 (2009), pp. 621–634] have introduced a new class of weighted exponential (WE) distributions, and this can be used quite effectively to model lifetime data. In this paper, we introduce a new class of weighted Marshall–Olkin bivariate exponential distributions. This new singular distribution has univariate WE marginals. We study different properties of the proposed model. There are four parameters in this model and the maximum-likelihood estimators (MLEs) of the unknown parameters cannot be obtained in explicit forms. We need to solve a four-dimensional optimization problem to compute the MLEs. One data set has been analysed for illustrative purposes and finally we propose some generalization of the proposed model.     

Missing covariates in generalized linear models when the missing data mechanism is non-ignorable

We propose a method for estimating parameters in generalized linear models with missing covariates and a non-ignorable missing data mechanism. We use a multinomial model for the missing data indicators and propose a joint distribution for them which can be written as a sequence of one-dimensional conditional distributions, with each one-dimensional conditional distribution consisting of a logistic regression. We allow the covariates to be either categorical or continuous. The joint covariate distribution is also modelled via a sequence of one-dimensional conditional distributions, and the response variable is assumed to be completely observed. We derive the E- and M-steps of the EM algorithm with non-ignorable missing covariate data. For categorical covariates, we derive a closed form expression for the E- and M-steps of the EM algorithm for obtaining the maximum likelihood estimates (MLEs). For continuous covariates, we use a Monte Carlo version of the EM algorithm to obtain the MLEs via the Gibbs sampler. Computational techniques for Gibbs sampling are proposed and implemented. The parametric form of the assumed missing data mechanism itself is not `testable' from the data, and thus the non-ignorable modelling considered here can be viewed as a sensitivity analysis concerning a more complicated model. Therefore, although a model may have `passed' the tests for a certain missing data mechanism, this does not mean that we have captured, even approximately, the correct missing data mechanism. Hence, model checking for the missing data mechanism and sensitivity analyses play an important role in this problem and are discussed in detail. Several simulations are given to demonstrate the methodology. In addition, a real data set from a melanoma cancer clinical trial is presented to illustrate the methods proposed.     

Estimating the parameters of a doubly truncated normal distribution

This paper deals with the maximum likelihood estimation of parameters for a doubly truncated normal distribution when the truncation points are known. We prove, in this case, that the MLEs are nonexistent (become infinite) with positive probability. For estimators that exist with probability one, the class of Bayes modal estimators or modified maximum likelihood estimators is explored. Another useful estimating procedure, called mixed estimation, is proposed. Simulations compare the behavior of the MLEs, the modified MLEs, and the mixed estimators which reveal that the MLE, in addition to being nonexistent with positive probability, behaves poorly near the upper boundary of the interval of its existence. The modified MLEs and the mixed estimators are seen to be remarkably better than the MLE     

Estimating the Pareto parameters under progressive censoring data for constant-partially accelerated life tests

Accelerated life testing is widely used in product life testing experiments since it provides significant reduction in time and cost of testing. In this paper, assuming that the lifetime of items under use condition follow the two-parameter Pareto distribution of the second kind, partially accelerated life tests based on progressively Type-II censored samples are considered. The likelihood equations of the model parameters and the acceleration factor are reduced to a single nonlinear equation to be solved numerically to obtain the maximum-likelihood estimates (MLEs). Based on normal approximation to the asymptotic distribution of MLEs, the approximate confidence intervals (ACIs) for the parameters are derived. Two bootstrap CIs are also proposed. The classical Bayes estimates cannot be obtained in explicit form, so we propose to apply Markov chain Monte Carlo method to tackle this problem, which allows us to construct the credible interval of the involved parameters. Analysis of a simulated data set has also been presented for illustrative purposes. Finally, a Monte Carlo simulation study is carried out to investigate the precision of the Bayes estimates with MLEs and to compare the performance of different corresponding CIs considered.     

On classical and Bayesian order restricted inference for multiple exponential step stress model

In this article, we consider the multiple step stress model based on the cumulative exposure model assumption. Here, it is assumed that for a given stress level, the lifetime of the experimental units follows exponential distribution and the expected lifetime decreases as the stress level increases. We mainly focus on the order restricted inference of the unknown parameters of the lifetime distributions. First we consider the order restricted maximum likelihood estimators (MLEs) of the model parameters. It is well known that the order restricted MLEs cannot be obtained in explicit forms. We propose an algorithm that stops in finite number of steps and it provides the MLEs. We further consider the Bayes estimates and the associated credible intervals under the squared error loss function. Due to the absence of explicit form of the Bayes estimates, we propose to use the importance sampling technique to compute Bayes estimates. We provide an extensive simulation study in case of three stress levels mainly to see the performance of the proposed methods. Finally the analysis of one real data set has been provided for illustrative purposes.     

Regular mles for nonregular distributions

Distributions whose extremity values of the support depend on unknown pa¬rameters are usually known as nonregular distributions. In most cases, the MLEs for these parameters cannot be obtained by differentiation. Familiar examples are the uniform distribution on the interval (0,0) and the truncated exponential distribution with truncation parameter 0. However, there exist distributions whose extremity points of the support depend on unknown pa¬rameters, which nevertheless are regular in the sense that the MLEs can be obtained by differentiation. This note provides a method of constructing such nonregular distributions with regular MLEs.     

Exact inference on multiple exponential populations under a joint type-II progressive censoring scheme

ABSTRACT

Recently Mondal and Kundu [Mondal S, Kundu D. A new two sample type-II progressive censoring scheme. Commun Stat Theory Methods. 2018. doi:10.1080/03610926.2018.1472781] introduced a Type-II progressive censoring scheme for two populations. In this article, we extend the above scheme for more than two populations. The aim of this paper is to study the statistical inference under the multi-sample Type-II progressive censoring scheme, when the underlying distributions are exponential. We derive the maximum likelihood estimators (MLEs) of the unknown parameters when they exist and find out their exact distributions. The stochastic monotonicity of the MLEs has been established and this property can be used to construct exact confidence intervals of the parameters via pivoting the cumulative distribution functions of the MLEs. The distributional properties of the ordered failure times are also obtained. The Bayesian analysis of the unknown model parameters has been provided. The performances of the different methods have been examined by extensive Monte Carlo simulations. We analyse two data sets for illustrative purposes.     

Consistency of Normal Distribution Based Pseudo Maximum Likelihood Estimates When Data Are Missing at Random

This paper shows that, when variables with missing values are linearly related to observed variables, the normal-distribution-based pseudo MLEs are still consistent. The population distribution may be unknown while the missing data process can follow an arbitrary missing at random mechanism. Enough details are provided for the bivariate case so that readers having taken a course in statistics/probability can fully understand the development. Sufficient conditions for the consistency of the MLEs in higher dimensions are also stated, while the details are omitted.     

Censored Exponential Data: Large Deviations for MLEs and Posterior Distributions

In this article, we consider several statistical models for censored exponential data. We prove a large deviation result for the maximum likelihood estimators (MLEs) of each model, and a unique result for the posterior distributions which works well for all the cases. Finally, comparing the large deviation rate functions for MLEs and posterior distributions, we show that a typical feature fails for one model; moreover, we illustrate the relation between this fact and a well-known result for curved exponential models.     

Statistical inference for the extreme value distribution under adaptive Type-II progressive censoring schemes

Adaptive Type-II progressive censoring schemes have been shown to be useful in striking a balance between statistical estimation efficiency and the time spent on a life-testing experiment. In this article, some general statistical properties of an adaptive Type-II progressive censoring scheme are first investigated. A bias correction procedure is proposed to reduce the bias of the maximum likelihood estimators (MLEs). We then focus on the extreme value distributed lifetimes and derive the Fisher information matrix for the MLEs based on these properties. Four different approaches are proposed to construct confidence intervals for the parameters of the extreme value distribution. Performance of these methods is compared through an extensive Monte Carlo simulation.     

Exact inference for a simple step-stress model with Type-I hybrid censored data from the exponential distribution

In reliability and life-testing experiments, the researcher is often interested in the effects of extreme or varying stress factors such as temperature, voltage and load on the lifetimes of experimental units. Step-stress test, which is a special class of accelerated life-tests, allows the experimenter to increase the stress levels at fixed times during the experiment in order to obtain information on the parameters of the life distributions more quickly than under normal operating conditions. In this paper, we consider the simple step-stress model under the exponential distribution when the available data are Type-I hybrid censored. We derive the maximum likelihood estimators (MLEs) of the parameters assuming a cumulative exposure model with lifetimes being exponentially distributed. The exact distributions of the MLEs of parameters are obtained through the use of conditional moment generating functions. We also derive confidence intervals for the parameters using these exact distributions, asymptotic distributions of the MLEs and the parametric bootstrap methods, and assess their performance through a Monte Carlo simulation study. Finally, we present two examples to illustrate all the methods of inference discussed here.     

Bayesian Estimation of the Parameters of Bivariate Exponential Distributions

In this article, we develop an empirical Bayesian approach for the Bayesian estimation of parameters in four bivariate exponential (BVE) distributions. We have opted for gamma distribution as a prior for the parameters of the model in which the hyper parameters have been estimated based on the method of moments and maximum likelihood estimates (MLEs). A simulation study was conducted to compute empirical Bayesian estimates of the parameters and their standard errors. We use moment estimators or MLEs to estimate the hyper parameters of the prior distributions. Furthermore, we compare the posterior mode of parameters obtained by different prior distributions and the Bayesian estimates based on gamma priors are very close to the true values as compared to improper priors. We use MCMC method to obtain the posterior mean and compared the same using the improper priors and the classical estimates, MLEs.     

Inference for the extreme value distribution under progressive Type-II censoring

The extreme value distribution has been extensively used to model natural phenomena such as rainfall and floods, and also in modeling lifetimes and material strengths. Maximum likelihood estimation (MLE) for the parameters of the extreme value distribution leads to likelihood equations that have to be solved numerically, even when the complete sample is available. In this paper, we discuss point and interval estimation based on progressively Type-II censored samples. Through an approximation in the likelihood equations, we obtain explicit estimators which are approximations to the MLEs. Using these approximate estimators as starting values, we obtain the MLEs using an iterative method and examine numerically their bias and mean squared error. The approximate estimators compare quite favorably to the MLEs in terms of both bias and efficiency. Results of the simulation study, however, show that the probability coverages of the pivotal quantities (for location and scale parameters) based on asymptotic normality are unsatisfactory for both these estimators and particularly so when the effective sample size is small. We, therefore, suggest the use of unconditional simulated percentage points of these pivotal quantities for the construction of confidence intervals. The results are presented for a wide range of sample sizes and different progressive censoring schemes. We conclude with an illustrative example.     

On the existence and uniqueness of the maximum likelihood estimates of the parameters of Birnbaum–Saunders distribution based on Type-I,Type-II and hybrid censored samples

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution and is a common model for analysing lifetime data. In this paper, we discuss the existence and uniqueness of the maximum likelihood estimates (MLEs) of the parameters of BS distribution based on Type-I, Type-II and hybrid censored samples. The line of proof is based on the monotonicity property of the likelihood function. We then describe the numerical iterative procedure for determining the MLEs of the parameters, and point out briefly some recently developed simple methods of estimation in the case of Type-II censoring. Some graphical illustrations of the approach are given for three real data from the reliability literature. Finally, for illustrative purpose, we also present an example in which the MLEs do not exist.     

Lattice conditional independence models for seemingly unrelated regressions

Seemingly unrelated regressions (SUR) models appear frequently in econometrics and in the analyses of repeated measures designs and longitudinal data. It is known that iterative algorithms are generally required to obtain the MLEs of the regression parameters. Under a minimal set of lattice conditional independence (LCI) restrictions imposed on the covariance structure, however, closed-form MLEs can be obtained by standard linear regression techniques (Andersson and Perlman, 1993, 1994, 1998). In this paper, simulation is used to study the efficiency of these LCI model-based estimators. We also propose two possible improvements of the usual two-stage estimators for the regression parameters.