Flowgraph Models for Generalized Phase Type Distributions Having Non-Exponential Waiting Times

ABSTRACT. Aalen (1995) introduced phase type distributions based on Markov processes for modelling disease progression in survival analysis. For tractability and to maintain the Markov property, these use exponential waiting times for transitions between states. This article extends the work of Aalen (1995) by generalizing these models to semi-Markov processes with non-exponential waiting times. The generalization allows more realistic modelling of the stages of a disease where the Markov property and exponential waiting times may not hold. Flowgraph models are introduced to provide a closed form for the distributions in situations involving non-exponential waiting times. Flowgraph models work where traditional methods of stochastic processes are intractable. Saddlepoint approximations are used in the analysis. Together, generalized phase type distributions, flowgraphs, and saddlepoint approximations create exciting and innovative prospects for the analysis of survival data.     

An alternative competing risk model to the Weibull distribution for modelling aging in lifetime data analysis

A simple competing risk distribution as a possible alternative to the Weibull distribution in lifetime analysis is proposed. This distribution corresponds to the minimum between exponential and Weibull distributions. Our motivation is to take account of both accidental and aging failures in lifetime data analysis. First, the main characteristics of this distribution are presented. Then, the estimation of its parameters are considered through maximum likelihood and Bayesian inference. In particular, the existence of a unique consistent root of the likelihood equations is proved. Decision tests to choose between an exponential, Weibull and this competing risk distribution are presented. And this alternative model is compared to the Weibull model from numerical experiments on both real and simulated data sets, especially in an industrial context.     

A Bayesian piecewise exponential phase II design for monitoring a time-to-event endpoint

A robust Bayesian design is presented for a single-arm phase II trial with an early stopping rule to monitor a time to event endpoint. The assumed model is a piecewise exponential distribution with non-informative gamma priors on the hazard parameters in subintervals of a fixed follow up interval. As an additional comparator, we also define and evaluate a version of the design based on an assumed Weibull distribution. Except for the assumed models, the piecewise exponential and Weibull model based designs are identical to an established design that assumes an exponential event time distribution with an inverse gamma prior on the mean event time. The three designs are compared by simulation under several log-logistic and Weibull distributions having different shape parameters, and for different monitoring schedules. The simulations show that, compared to the exponential inverse gamma model based design, the piecewise exponential design has substantially better performance, with much higher probabilities of correctly stopping the trial early, and shorter and less variable trial duration, when the assumed median event time is unacceptably low. Compared to the Weibull model based design, the piecewise exponential design does a much better job of maintaining small incorrect stopping probabilities in cases where the true median survival time is desirably large.     

A semi-Markov multistate model for estimation of the mean quality-adjusted survival for non-progressive processes

We discuss the estimation of the expected value of the quality-adjusted survival, based on multistate models. We generalize an earlier work, considering the sojourn times in health states are not identically distributed, for a given vector of covariates. Approaches based on semiparametric and parametric (exponential and Weibull distributions) methodologies are considered. A simulation study is conducted to evaluate the performance of the proposed estimator and the jackknife resampling method is used to estimate the variance of such estimator. An application to a real data set is also included.     

A class of absolutely continuous bivariate distributions

Block and Basu bivariate exponential distribution is one of the most popular absolutely continuous bivariate distributions. Extensive work has been done on the Block and Basu bivariate exponential model over the past several decades. Interestingly it is observed that the Block and Basu bivariate exponential model can be extended to the Weibull model also. We call this new model as the Block and Basu bivariate Weibull model. We consider different properties of the Block and Basu bivariate Weibull model. The Block and Basu bivariate Weibull model has four unknown parameters and the maximum likelihood estimators cannot be obtained in closed form. To compute the maximum likelihood estimators directly, one needs to solve a four dimensional optimization problem. We propose to use the EM algorithm for computing the maximum likelihood estimators of the unknown parameters. The proposed EM algorithm can be carried out by solving one non-linear equation at each EM step. Our method can be also used to compute the maximum likelihood estimators for the Block and Basu bivariate exponential model. One data analysis has been preformed for illustrative purpose.     

A review of model selection procedures relevant to the weibull distribution

A number of goodness-of-fit and model selection procedures related to the Weibull distribution are reviewed. These procedures include probability plotting, correlation type goodness-of-fit tests, and chi-square goodness-of-fit tests. Also the Kolmogorow-Smirniv, Kuiper, and Cramer-Von Mises test statistics for completely specified hypothesis based on censored data are reviewed, and these test statistics based on complete samples for the unspecified parameters case are considered. Goodness-of-fit tests based on sample spacings, and a goodness-of-fit test for the Weibull process, is also discussed.

Model selection procedures for selecting between a Weibull and gamma model, a Weibull and lognormal model, and for selecting from among all three models are considered. Also tests of exponential versus Weibull and Weibull versus generalized gamma are mentioned.     

Dynamic latent trait models with mixed hidden Markov structure for mixed longitudinal outcomes

We propose a general Bayesian joint modeling approach to model mixed longitudinal outcomes from the exponential family for taking into account any differential misclassification that may exist among categorical outcomes. Under this framework, outcomes observed without measurement error are related to latent trait variables through generalized linear mixed effect models. The misclassified outcomes are related to the latent class variables, which represent unobserved real states, using mixed hidden Markov models (MHMMs). In addition to enabling the estimation of parameters in prevalence, transition and misclassification probabilities, MHMMs capture cluster level heterogeneity. A transition modeling structure allows the latent trait and latent class variables to depend on observed predictors at the same time period and also on latent trait and latent class variables at previous time periods for each individual. Simulation studies are conducted to make comparisons with traditional models in order to illustrate the gains from the proposed approach. The new approach is applied to data from the Southern California Children Health Study to jointly model questionnaire-based asthma state and multiple lung function measurements in order to gain better insight about the underlying biological mechanism that governs the inter-relationship between asthma state and lung function development.     

An improper form of Weibull distribution for competing risks analysis with Bayesian approach

ABSTRACT

In survival analysis, individuals may fail due to multiple causes of failure called competing risks setting. Parametric models such as Weibull model are not improper that ignore the assumption of multiple failure times. In this study, a novel extension of Weibull distribution is proposed which is improper and then can incorporate to the competing risks framework. This model includes the original Weibull model before a pre-specified time point and an exponential form for the tail of the time axis. A Bayesian approach is used for parameter estimation. A simulation study is performed to evaluate the proposed model. The conducted simulation study showed identifiability and appropriate convergence of the proposed model. The proposed model and the 3-parameter Gompertz model, another improper parametric distribution, are fitted to the acute lymphoblastic leukemia dataset.     

Report of the ASA Technical Panel on the Census Undercount

A new model is proposed for the joint distribution of paired survival times generated from clinical trials and certain reliability settings. The new model can be considered an extension to the bivariate exponential models studied in the literature. Here, a more flexible bivariate Weibull model will be derived, and two exact parametric tests for testing the equality of marginal survival distributions are developed.     

Expanded renal transplantation: a competing risk model approach

Multi-state models (MSMs) are useful to analyze survival data when, besides the event of main interest, one or more intermediate states of the individual are identified. These models take the several existing states and the possible transitions among them into account. At the same time, covariate effects on each transition intensity may be investigated separately and, therefore, MSMs are more flexible than the standard Cox proportional hazards model. In this work, we use MSMs to investigate the impact of the quality of a transplanted kidney for a group of patients at the Hospital Universitario Central de Asturias. Specifically, we use an illness-death model to study the evolution of patients with kidney disease who received a renal transplant after a dialysis period. The intermediate state is defined as the failure of the received organ, while the terminating state is the death of the patient. In order to increase the potential number of organs available for transplant, the standards of quality for the transplanted kidneys were relaxed (the new criteria are labeled expanded criteria), and these ‘expanded kidneys’ were transplanted in appropriate candidates (older patients, with higher prevalence of diabetes mellitus). Results suggest that the expanded kidneys have a minor effect on survival, while both the kidney mortality and the risk of death increase with the patient's age and the serum creatinine and serum hemoglobin levels.     

Weibull prediction of event times in clinical trials

In clinical trials with interim analyses planned at pre-specified event counts, one may wish to predict the times of these landmark events as a tool for logistical planning. Currently available methods use either a parametric approach based on an exponential model for survival (Bagiella and Heitjan, Statistics in Medicine 2001; 20:2055) or a non-parametric approach based on the Kaplan-Meier estimate (Ying et al., Clinical Trials 2004; 1:352). Ying et al. (2004) demonstrated the trade-off between bias and variance in these models; the exponential method is highly efficient when its assumptions hold but potentially biased when they do not, whereas the non-parametric method has minimal bias and is well calibrated under a range of survival models but typically gives wider prediction intervals and may fail to produce useful predictions early in the trial. As a potential compromise, we propose here to make predictions under a Weibull survival model. Computations are somewhat more difficult than with the simpler exponential model, but Monte Carlo studies show that predictions are robust under a broader range of assumptions. We demonstrate the method using data from a trial of immunotherapy for chronic granulomatous disease.     

Testing Symmetry of a NIG Distribution

Exponential and Weibull models are commonly used models with former being the special case of the latter. In their most general forms, the exponential model involves both threshold and scale parameters whereas the Weibull model involves threshold, scale and shape parameters. The article analyzes the two models in a Bayesian framework and examines the feasibility of generality versus particularity in the sense that it tests for the possibility of (not) having a threshold and/or a shape parameter in the data arising from exponential (Weibull) model. The results are illustrated based on both complete and censored datasets from the models.     

A Parametric Model for the Interval Censored Survival Times of Acacia Mangium Plantation in a Spacing Trial

Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.     

Generalized exponential distribution: Existing results and some recent developments

Mudholkar and Srivastava [1993. Exponentiated Weibull family for analyzing bathtub failure data. IEEE Trans. Reliability 42, 299–302] introduced three-parameter exponentiated Weibull distribution. Two-parameter exponentiated exponential or generalized exponential distribution is a particular member of the exponentiated Weibull distribution. Generalized exponential distribution has a right skewed unimodal density function and monotone hazard function similar to the density functions and hazard functions of the gamma and Weibull distributions. It is observed that it can be used quite effectively to analyze lifetime data in place of gamma, Weibull and log-normal distributions. The genesis of this model, several properties, different estimation procedures and their properties, estimation of the stress-strength parameter, closeness of this distribution to some of the well-known distribution functions are discussed in this article.     

A new class of weighted exponential distributions

Introducing a shape parameter to an exponential model is nothing new. There are many ways to introduce a shape parameter to an exponential distribution. The different methods may result in variety of weighted exponential (WE) distributions. In this article, we have introduced a shape parameter to an exponential model using the idea of Azzalini, which results in a new class of WE distributions. This new WE model has the probability density function (PDF) whose shape is very close to the shape of the PDFS of Weibull, gamma or generalized exponential distributions. Therefore, this model can be used as an alternative to any of these distributions. It is observed that this model can also be obtained as a hidden truncation model. Different properties of this new model have been discussed and compared with the corresponding properties of well-known distributions. Two data sets have been analysed for illustrative purposes and it is observed that in both the cases it fits better than Weibull, gamma or generalized exponential distributions.     

Effects of censoring on the robustness of exponential-based confidence intervals for median lifetime

Statistical procedures for constructing confidence intervals for median lifetime often rest on a distributional assumption for failure times.This paper explores the interplay between censoring levels and robustness for two construction procedures based on exponential lifetime, subject to general right-censoring. Data are simulated from nearby Weibull distributions. As expected, the simulations indicate that when the exponential assumption is not satisfied, observed coverage by the confidence intervals may differ substantially from the specified coverage level. The marked improvement in the robustness properties of the intervals as the level of censoring increases suggests questions for future research.     

Analysis of survival data by a Weibull–Bessel distribution

In survival analysis and reliability studies, problems with random sample size arise quite frequently. More specifically, in cancer studies, the number of clonogens is unknown and the time to relapse of the cancer is defined by the minimum of the incubation times of the various clonogenic cells. In this article, we have proposed a new model where the distribution of the incubation time is taken as Weibull and the distribution of the random sample size as Bessel, giving rise to a Weibull–Bessel distribution. The maximum likelihood estimation of the model parameters is studied and a score test is developed to compare it with its special submodel, namely, exponential–Bessel distribution. To illustrate the model, two real datasets are examined, and it is shown that the proposed model, presented here, fits better than several other existing models in the literature. Extensive simulation studies are also carried out to examine the performance of the estimates.     

A generalization of discrete Weibull distribution

In this paper, we introduce a new family of discrete distributions and study its properties. It is shown that the new family is a generalization of discrete Marshall-Olkin family of distributions. In particular, we study generalized discrete Weibull distribution in detail. Discrete Marshall-Olkin Weibull distribution, exponentiated discrete Weibull distribution, discrete Weibull distribution, discrete Marshall-Olkin generalized exponential distribution, exponentiated geometric distribution, generalized discrete exponential distribution, discrete Marshall-Olkin Rayleigh distribution and exponentiated discrete Rayleigh distribution are sub-models of generalized discrete Weibull distribution. We derive some basic distributional properties such as probability generating function, moments, hazard rate and quantiles of the generalized discrete Weibull distribution. We can see that the hazard rate function can be decreasing, increasing, bathtub and upside-down bathtub shape. Estimation of the parameters are done using maximum likelihood method. A real data set is analyzed to illustrate the suitability of the proposed model.     

The beta exponentiated Weibull distribution

The Weibull distribution is one of the most important distributions in reliability. For the first time, we introduce the beta exponentiated Weibull distribution which extends recent models by Lee et al. [Beta-Weibull distribution: some properties and applications to censored data, J. Mod. Appl. Statist. Meth. 6 (2007), pp. 173–186] and Barreto-Souza et al. [The beta generalized exponential distribution, J. Statist. Comput. Simul. 80 (2010), pp. 159–172]. The new distribution is an important competitive model to the Weibull, exponentiated exponential, exponentiated Weibull, beta exponential and beta Weibull distributions since it contains all these models as special cases. We demonstrate that the density of the new distribution can be expressed as a linear combination of Weibull densities. We provide the moments and two closed-form expressions for the moment-generating function. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The density of the order statistics can also be expressed as a linear combination of Weibull densities. We obtain the moments of the order statistics. The expected information matrix is derived. We define a log-beta exponentiated Weibull regression model to analyse censored data. The estimation of the parameters is approached by the method of maximum likelihood. The usefulness of the new distribution to analyse positive data is illustrated in two real data sets.     

General results for the beta Weibull distribution

In this paper, we study some mathematical properties of the beta Weibull (BW) distribution, which is a quite flexible model in analysing positive data. It contains the Weibull, exponentiated exponential, exponentiated Weibull and beta exponential distributions as special sub-models. We demonstrate that the BW density can be expressed as a mixture of Weibull densities. We provide their moments and two closed-form expressions for their moment-generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and two entropies. The density of the BW-order statistics is a mixture of Weibull densities and two closed-form expressions are derived for their moments. The estimation of the parameters is approached by two methods: moments and maximum likelihood. We compare the performances of the estimates obtained from both the methods by simulation. The expected information matrix is derived. For the first time, we introduce a log-BW regression model to analyse censored data. The usefulness of the BW distribution is illustrated in the analysis of three real data sets.