Several Major Distributions of Survival Function

September 3, 2017 Medical Statistics, Oncology, Research No comments , , , , , , , , , , , , , ,

Exponential Distribution

The simplest and most important distribution in survival studies is the exponential distribution. In the late 1940s, researchers began to choose the exponential distribution to describe the life pattern of electronic systems. The exponential distribution has since continued to play a role in lifetime studies analogous to that of the normal distribution in other areas of statistics. The exponential distribution is often referred to as a purely random failure pattern. It is famous for its unique "lack of memory," which requires that the age of the animal or person does not affect future survival. Although many survival data cannot be described adequately by the exponential distribution, an understanding of it facilitates the treatment of more general situations.

The exponential distribution is characterized by a constant hazard rate 𝜆, its only parameter. A high ðœ† value indicates high risk and short survival; a low ðœ† value indicates low risk and long survival. When the surival time T follows the exponential distribution with a parameter ðœ†, the probability density function is defined as

The cumulative distribution function is

and the survivorship function is then

and the hazard function is

Note that the hazard function is a constant, ðœ†, independent of t. Because the exponential distribution is characterized by a constant hazard rate, independent of the age of the person, there is no aging or wearing out, and failure or death is a random event indepdendent of time. When natural logarithms of the survivorship function are taken, log S(t) = -𝜆t, which is a linear function of t.

Weibull Distribution

The Weibull distribution is a generalization of the exponential distribution. However, unlike the exponential distribution, it does not assume a constant hazard rate and therefore has broader application. The Weibull distribution is characterized by two parameters, 𝛾 and 𝜆. The value of ð›¾ determines the shape of the distribution curve and the value of ðœ† determines its scaling. Consequently, ð›¾ and 𝜆 are called the shape and scale parameters, respectively. When ð›¾ = 1, the hazard rate remains constant as time increases; this is the exponential case. The hazard rate increases when ð›¾ >1 and decrease when ð›¾ <1 as t increases. Thus, the Weibull distribution may be used to model the survival distribution of a population with increasing, decreasing, or constant risk.

The probability density function, cumulative distribution functions are, survivorship function, and hazard function are:

Weibull distribution is named after Swedish mathematician Waloddi Weibull, who described it in detail in 1951, although it was first identified by Frechet and first applied by Rosin & Rammler to describe a particle size distribution.

Lognormal Distribution

In its simplest form the lognormal distribution can be defined as the distribution of a variable whose logarithm follows the normal distribution. Its origin may be traced as far back as 1879, when McAlister described explicitly a theory of the distribution. Most of its aspects have since been under study. Gaddum gave a review of its application in biology, followed by Boag's applications in cancer research. Its history, properties, estimation problem, and uses in economics have been discussed in detail by AItchison and Brown. Later, other investigators also observed that the age at onset of Alzheimer's disease and the distribution of survival time of several diseases such as Hodgkin's disease and chronic leukemia could be rather closely approximated by a lognormal distribution since they are markedly skewed to the right and the logarithms of survival times are approxiamtely normally distributed.

Consider the survival time T such that log T is normally distributed with mean 𝜇 and variance 𝜎2. We then say that T is lognormally distributed and write T as 𝛬(𝜇, ðœŽ2). It should be noted that ðœ‡ and 𝜎2 are not the mean and variance of the lognormal distribution. The hazard function of the lognormal distrition increases initially to a maximum and then decreases (almost as soon as the median is passed) to zero as time approaches infinity. Therefore, the lognormal distribution is suitable for survival patterns with an initially increasing and then decreasing hazard rate. By a central limit theorem, it can be shown that the distribution of the product of n independent positive variates approaches a lognormal distribution under very general conditions: for example, the distribution of the size of an organism whose growth is subject to many small impulses, the effect of each of which is proportional to the momentary size of the organism.

Gamma Distributions

The gamma distribution, which includes the exponential and chi-square distribution, was used a long time ago by Brown and Flood to describe the life of glass tumblers circulating in a cafeteria and by Birnbaum and Saunders as a statistical model for life length of materials. Since then, this distribution has been used frequently as a model for industrial reliability problems and human survival.

Suppose that failure or death takes place in n stages or as soon as n subfailures have happened. At the end of the first stage, after time T1, the first subfailure occurs; after that the second stage begins and the second subfailure occurs after time T2; and so on. Total failure or death occurs at the end of the nth stage, when the nth subfailure happens. The survival time, T, is then T1 + T2 + … + Tn. The times T1, T2, …, Tn spent in each stage are assumed to be independently exponentially distributed with probability density function 𝜆exp(-𝜆ti), i = 1, …, n. That is, the subfailures occur independently at a constant rate ðœ†. The distribution of T is then called the Erlangian distribution. There is no need for the stages to have physical significance since we can always assume that death occurs in the n-stage process just described. This idea, introduced by A. K. Erlang in his study of congestion in telephone systems, has been used widely in queuing theory and life processes.

The gamma distribution is characterized by two parameters, 𝛾 and ðœ†. When 0 < ð›¾ < 1, there is negative aging and the hazard rate decreases monotonically from infinity to ðœ† as time increases from 0 to infinity. When ð›¾ > 1, there is positive aging and the hazard rate increases monotonically from 0 to ðœ† as time increases from 0 to infinity. When ð›¾ = 1, the hazard rate equals ðœ†, a constant, as in the exponential case.

Log-logistic Distribution

The survival time T has a log-logistic distribution if log(T) has a logistic distribution. The density, survivorship, hazard, and cumulative hazard functions of the log-logistic distribution are, respectively,

The log-logistic distribution is characterized by two parameters, 𝛼, and 𝛾. The median of the log-logistic distribution is ð›¼-1/𝛾. When ð›¾ > 1, the log-logistic hazard has the value 0 at time 0, increases to a peak at a specific t, and then declines, which is similar to the lognormal hazard. When ð›¾ = 1, the hazard starts at ð›¼1/𝛾 and then declines monotonically. When ð›¾ < 1, the hazard starts at infinity and then declines, which is similar to the Weibull distribution. The hazard function declines toward 0 as t approaches infinity. Thus, the log-logistic distribution may be used to describe a first increasing and then decreasing hazard or a monotonically decreasing hazard.

Other Survival Distributions

Many other distributions can be used as models of survival time, three of which we discuss briefly in this section: the linear exponential, the Gompertz, and a distribution whose hazard rate is a step function. The linear-exponential model and the Gompertz distribution are extensions of the exponential distribution. Both describe survival patterns that have a constant initial hazard rate. The hazard rate varies as a linear function of time or age in the linear-exponential model and as an exponential function of time or age in the Gompertz distribution.

In demonstrating the use of the linear-exponential model, Broadbent, uses as an example the serivce of milk bottles that are filled in a dairy, circulated to customers, and returned empty to the dairy. The model was also used by Carbone et al. to describe the survival pattern of patients with plasmacytic myeloma. The hazard function of the linear-exponential distribution is

where ðœ† and ð›¾ can be values such that h(t) is nonnegative. The hazard rate increases from 𝜆 with time if ð›¾ > 0, decrease if ð›¾ < 0, and remains constant (an exponential case) if ð›¾ = 0. The probability density function and the survivorship function are, respectively,

The Gompertz distribution is also characterized by two parameters, ðœ† and 𝛾. The hazard function, survival function, and the probability density function are below, respectively,

Finally, we consider a distribution where the hazard rate is a step function. The hazard rate, survival function, and probability density function are below, respective,


One application of this distribution is the life-table analysis. In a life-table analysis, time is divided into intervals and the harzard rate is assumed to be constant in each interval. However, the overall hazard rate is not necessarily consrtant.


The nine distributions described above are, among others, reasonable model for survival time distribution. All have been designed by considering a biological failure, a death process, or an aging property. They may or may not be appropriate for many practical situations, but the objective here is to illustrate the various possible techniques, assumptions, and arguments that can be used to choose the most appropriate model. If none of these distributions fits the data, investigators might have to derive an original model to suit the particular data, perhaps by using some of the ideas presented here.

Specific Immunosuppressive Therapy

July 20, 2016 Hematology, Immunology, Infectious Diseases, Oncology, Pharmacology, Transplantation No comments , , , , , , , , , , , , , , , ,

The ideal immunosuppressant would be antigen-specific, inhibiting the immune response to the alloantigens present in the graft (or vice versa alloantigens present in recipient in GVHD) while preserving the recipient's ability to respond to other foreign antigens. Although this goal has not yet been achieved, several more targeted immunosuppressive agents have been developed. Most involve the use of monoclonal antibodies (mAbs) or soluble ligands that bind specific cell-surface molecules. On limitation of most first-generation of mAbs came from their origin in animals. Recipients of these frequently developed an immune response to the nonhuman epitopes, rapidly clearing the mAbs from the body. This limitation has been overcome by the construction of humanized mAbs and mouse-human chimeric antibodies.

Many different mAbs have been tested in transplantation settings, and the majority work by either depleting the recipient of a particular cell population or by blocking a key step in immune signaling. Antithymocyte globulin (ATG), prepared from animals exposed to human lymphocytes, can be used to deplete lymphocytes in patients prior to transplantation, but has significant side effects. A more subset-specific strategy uses a mAb to the CD3 molecule of the TCR, called OKT3, and rapidly depletes mature T cells from the circulation. This depletion appears to be caused by binding of antibody-coated T cells to Fc receptors on phagocytic cells, which then phagocytose and clear the T cells from the circulation. In a further refinement of this strategy, a cytotoxic agent such as diphtheria toxin is coupled with the mAb. Antibody-bound cells then internalize the toxin and die. Another technique uses mAbs specific for the high-affinity IL-2 receptor CD25. Since this receptor is expressed only on activated T cells, this treatment specifically blocks proliferation of T cells activated in response to the alloantigens of the graft. However, since TREG cells also express CD25 and may aid in alloantigen tolerance, this strategy may have drawbacks. More recently, a mAb against CD20 has been used to deplete mature B cells and is aimed at suppressing AMR (antibody-mediated rejection) responses. Finally, in cases of bone marrow transplantation, mAbs against T-cell-specific markers have been used to pretreat the donor's bone marrow to destory immunocompetent T cells that may react with the recipient tissues, causing GVHD.

Because cytokines appear to play an important role in allograft rejection, these compounds can also be specifically targeted. Animal studies have explored the use of mAbs specific for the cytokines implicated in transplant rejection, particularly TNF-alpha, IFN-gamma, and IL-2. In mice, anti-TNF-alpha mAbs prolong bone marrow transplants and reduce the incidence of GVHD. Antibodies to IFN-gamma and to IL-2 have each been reported in some cases to prolong cardiac transplants in rats.

TH-cell activation requires a costimulatory signal in addition to the signal mediated by the TCR. The interaction between CD80/86 on the membrane of APCs and the CD28 or CTLA-4 molecule on T cells provides one such signal. Without this costimulatory signal, antigen-activated T cells become anergic. CD28 is expressed on both resting and activated T cells, while CTLA-4 is expressed only on activated T cells and binds CD80/86 with a 20-fold-higher affinity. In mice, D. J. Lenschow, J. A. Bluestone, and colleagues demonstrated prolonged graft survival by blocking CD80/86 signaling with a soluble fusion protein consisting of the extracellular domain of CTLA-4 fused to human IgG1 heavy chain. This new drug, belatacept, was shown to induce anergy in T cells directed against the graft tissue and has been approved by the FDA for prevention of organ rejection in adult kidney transplant pateints.