Maximum Likelihood Methods for Fitting the Burr Type XII Distribution to Life Test Data

This paper describes mathematical and computational methodology for estimating the parameters of the Burr Type XII distribution by the method of maximum likelihood. Expressions for the asymptotic variances and covariances of the parameter estimates are given, and the modality of the log-likelihood and conditional log-likelihood functions is analyzed. As a result of this analysis for various a priori known and unknown parameter combinations, conditions are given which guarantee that the parameter estimates obtained will, indeed, be maximum likelihood estimates. An efficient numerical method for maximizing the conditional log-likelihood function is described, and mathematical expressions are given for the various numerical approximations needed to evaluate the expressions given for the asymptotic variances and covariances of the parameter estimates. The methodology discussed is applied in a numerical example to life test data arising in a clinical setting.     

Rank Covariance Methods for the Analysis of Survival Data

A procedure for comparing survival times between several groups of patients through rank analysis of covariance was introduced by WOOLSON and LACHENBRUCH (1983). It is a modification of Quade' rank analysis of covariance procedure (1967) and can be used for the analysis of right-censored data. In this paper, two additional modifications of Quade' original test statistic are proposed and compared to the original modification introduced by Woolson and Lachenbruch. These statistics are compared to one another and to the score test from Cox' proportional hazards model by way of a limited Monte Carlo study. One of the statistics, Q_R2, is recommended for general use for the rank analysis of covariance of right-censored survivorship data.     

Bayes Estimator of Survival Probability from Randomly Censored Observations with Weibull Distribution of Time to Death

This paper deals with Bayes estimation of survival probability when the data are randomly censored. Such a situation arises in case of a clinical trial which extends for a limited period T. A fixed number of patients (n) are observed whose times to death have identical Weibull distribution with parameters β and θ. The maximum times of observation for different patients are also independent uniform variables as the patients arrive randomly throughout the trial. For the joint prior distribution of (β, θ) as suggested by Sinha and Kale (1980, page 137) Bayes estimator of survival probability at time t (0<t<T) has been obtained. Considering squared error loss function it is the mean of the survival probability with respect to the posterior distribution of (β, θ). This estimator is then compared with the maximum likelihood estimator, by simulation, for various values of β, θ and censoring percentage. The proposed estimator is found to be better under certain conditions.     

A Model for Censored Survival Data When There is an Intervening Event

A model is discussed for incorporating information from a time-dependent covariable (an intervening event) and covariables independent of time into the analysis of survival data. In the model, it is assumed that individuals are potentially subject to two paths to failure, one including the intervening event and the other not. Additional assumptions are that failure times associated with the two paths are independent and that the time to failure subsequent to the intervening event is dependent on the intervening event time. Allowing the underlying hazard rates for the model to follow a WEIBULL form, use of the model and methods for fitting and hypothesis testing are illustrated by application to a follow-up study involving industrial workers where disability retirement was the intervening event. Extensions of the model to accommodate grouped survival data are presented.     

A Power Law Model for Survival Data Using GLIM

This paper proposes a regression model for the Weibull survival distribution of which the scale parameter is a power function of covariates. The estimation of parameters for partially censored data is pursued by using a statistical package called GLIM. Two sets of carcinogenic data are used to illustrate this procedure.     

On Generalized Log-Logistic Model for Censored Survival Data

In the analysis of survival data with parametric models, it is well known that the Weibull model is not suitable for modeling cases where the hazard rate is non-monotonic. For such cases, log-logistic model is frequently used. However, due to the symmetric property of the log-logistic model, it may be poor for the cases where the hazard rate is skewed or heavily tailed. In this paper, we suggest a generalization of the log-logistic model by introducing a shape parameter. This generalized model is then applied to fit the lung cancer data of Prentice (1973). The results seem to improve over those obtained by using the log-logistic model.     

The Planning of Sample Sizes for Maximum Likelihood Estimation of the Parameters of the Weibull Distribution

Based on maximum likelihood estimators the problem of setting confidence limits for the scale and shape parameters of the Weibull distribution with two parameters is considered. Approximation formulae for the determination of sample sizes for the parameter estimation are given for complete and type II censored samples.     

Modelling Survival Data to Examine Length of Treatment Effectiveness in a Clinical Trial

In some clinical trials where the experimental treatment is found to be effective in increasing survival, an important question is how long should the patient remain under treatment. Although the trial may not be designed to specifically answer this question, comparisons of the hazard curves among the treatment groups can yield some useful information. The survival data may be modelled using a flexible set of hazard functions and specific models are then chosen for further examination. This paper illustrates the approach using data from the Beta-Blocker Heart Attack Trial. Parametric and semi-parametric models are fitted and likelihood methods are used to assess length of treatment effectiveness.     

辐照剂量场分布的两种拟合方法比较研究

为比较不同方法拟合辐照剂量分布的相关度,本文基于空气比释动能计算公式和γ射线指数减弱规律两种方法分别拟合辐照室地面的剂量率分布;计算结果表明,两种方法都能较好反映测量区域的辐射剂量场分布,基于空气比释动能计算公式的拟合结果相关度要稍高于指数减弱规律的拟合结果。     

A Class of Multivariate Survival Functions Characterized by Minimum and Marginal Distribution

This paper characterizes a class of multivariate survival functions in terms of the minimum and marginal distributions.     

Some Methods for Summarizing Survivorship Data in Nonstandard Situations

One difficulty in summarising biological survivorship data is that the hazard rates are often neither constant nor increasing with time or decreasing with time in the entire life span. The promising Weibull model does not work here. The paper demonstrates how bath tub shaped quadratic models may be used in such a case. Further, sometimes due to a paucity of data actual lifetimes are not as certainable. It is shown how a concept from queuing theory namely first in first out (FIFO) can be profitably used here. Another nonstandard situation considered is one in which lifespan of the individual entity is too long compared to duration of the experiment. This situation is dealt with, by using ancilliary information. In each case the methodology is illustrated with numerical examples.     

Estimating Haplotype Effects for Survival Data

Summary Genetic association studies often investigate the effect of haplotypes on an outcome of interest. Haplotypes are not observed directly, and this complicates the inclusion of such effects in survival models. We describe a new estimating equations approach for Cox's regression model to assess haplotype effects for survival data. These estimating equations are simple to implement and avoid the use of the EM algorithm, which may be slow in the context of the semiparametric Cox model with incomplete covariate information. These estimating equations also lead to easily computable, direct estimators of standard errors, and thus overcome some of the difficulty in obtaining variance estimators based on the EM algorithm in this setting. We also develop an easily implemented goodness‐of‐fit procedure for Cox's regression model including haplotype effects. Finally, we apply the procedures presented in this article to investigate possible haplotype effects of the PAF‐receptor on cardiovascular events in patients with coronary artery disease, and compare our results to those based on the EM algorithm.     

Asymptotic Distribution of Score Statistics for Spatial Cluster Detection with Censored Data

Summary Cook, Gold, and Li (2007, Biometrics 63, 540–549) extended the Kulldorff (1997, Communications in Statistics 26, 1481–1496) scan statistic for spatial cluster detection to survival‐type observations. Their approach was based on the score statistic and they proposed a permutation distribution for the maximum of score tests. The score statistic makes it possible to apply the scan statistic idea to models including explanatory variables. However, we show that the permutation distribution requires strong assumptions of independence between potential cluster and both censoring and explanatory variables. In contrast, we present an approach using the asymptotic distribution of the maximum of score statistics in a manner not requiring these assumptions.     

Joint Modelling of Repeated Measures and Survival Time Data

In many clinical trials both repeated measures data and event history data are simultaneously observed from the same subject. These two types of responses are usually correlated, because they are from the same subject. In this article, we propose a joint model for the combined analysis of repeated measures data and event history data in the framework of hierarchical generalized linear models. The correlation between repeated measures and event time is modelled by introducing a shared random effect. The model parameters are estimated using the hierarchical‐likelihood approach. The proposed model is illustrated using a real data set for the renal transplant patients.     

Fitting Regression Models to Ordinal Data

This paper deals with the analysis of ordinal data by means of a threshold model. Maximum likelihood estimation is discussed and two examples are used to illustrate the methods.