Tutorials › VEE Mathematical Statistics › Deriving the Exponential CDF and Survival Function

VEE Mathematical Statistics · Uniform and Exponential Models · Tutorial 181 of 880

Deriving the Exponential CDF and Survival Function

Derive both exponential functions from the density, then apply them carefully at positive times and at the boundary of the model’s support.

Beginner 10 min read

What You'll Learn

  • Integrate the exponential density to obtain its cumulative distribution function.
  • Derive the survival function both as a complement and as an integral.
  • State the cumulative distribution and survival functions for all real inputs, including boundary cases.
  • Distinguish a density value from a probability and from the derivative of a CDF.
  • Check exponential calculations for signs, support, units, and probability totals.

Why Derive Both Functions?

For a waiting time or lifetime, the cumulative distribution function answers how likely it is that the event has occurred by a specified time. The survival function answers how likely it is that the waiting time exceeds that time. Both functions come from the same exponential density, but the integral limits and the meaning of each result differ. Deriving them makes those distinctions explicit and helps prevent sign and support errors.

Assume throughout that \(T\) is a nonnegative waiting time with an exponential distribution of rate \(\lambda>0\). The rate is measured in inverse time units: if time is measured in months, for example, \(\lambda\) is per month. The density is zero for negative times and is \(\lambda e^{-\lambda t}\) for nonnegative times. This is a continuous model, so it assigns probability zero to any single time, including \(T=0\).

Definition: For a random variable \(T\), its cumulative distribution function is \(F_T(t)=P(T\le t)\). Its survival function is \(S_T(t)=P(T>t)\). Under the exponential model with rate \(\lambda>0\), the density is \(f_T(t)=0\) for \(t<0\), and \(f_T(t)=\lambda e^{-\lambda t}\) for \(t\ge 0\).

Deriving the Cumulative Distribution Function

By definition, \(F_T(t)\) is the probability that \(T\) lies at or below \(t\). Since \(T\) cannot be negative, \(F_T(t)=0\) whenever \(t<0\). For \(t\ge 0\), integrate the density over the possible values from 0 through \(t\):

$$ \begin{aligned} F_T(t) &=P(T\le t)\\ &=\int_0^t \lambda e^{-\lambda u}\,du\\ &=\left[-e^{-\lambda u}\right]_0^t\\ &=-e^{-\lambda t}-(-e^0)\\ &=1-e^{-\lambda t}. \end{aligned} $$

The variable \(u\) is an integration variable; \(t\) is the time at which the CDF is evaluated. The antiderivative has a minus sign because the derivative of \(e^{-\lambda u}\) is \(-\lambda e^{-\lambda u}\). Consequently, the CDF for every real \(t\) is

$$ F_T(t)= \begin{cases} 0, & t<0,\\ 1-e^{-\lambda t}, & t\ge 0. \end{cases} $$

At \(t=0\), the second expression gives \(1-e^0=0\), agreeing with the first branch’s limiting value. Thus the CDF is continuous at zero, as expected for a continuous distribution with no probability mass at that point. For \(t\ge0\), it starts at zero and increases toward one as \(t\) increases.

Derivation: The exponential CDF results from integrating the density over the interval of possible values below \(t\). The negative sign in the antiderivative is essential: it produces \(1-e^{-\lambda t}\), not \(e^{-\lambda t}-1\).

Deriving the Survival Function

For any real-valued random variable, the events \(T\le t\) and \(T>t\) are complementary. Therefore, \(S_T(t)=1-F_T(t)\). For \(t\ge0\), substituting the CDF gives

$$ S_T(t)=1-\left(1-e^{-\lambda t}\right)=e^{-\lambda t}. $$

The same result follows directly by integrating the density above \(t\). For \(t\ge0\), the upper limit is infinity:

$$ \begin{aligned} S_T(t) &=P(T>t)\\ &=\int_t^\infty \lambda e^{-\lambda u}\,du\\ &=\left[-e^{-\lambda u}\right]_t^\infty\\ &=0-(-e^{-\lambda t})\\ &=e^{-\lambda t}. \end{aligned} $$

The term at infinity is zero because \(\lambda>0\), so \(e^{-\lambda u}\) tends to zero as \(u\) tends to infinity. For \(t<0\), all possible values of \(T\) exceed \(t\), so \(S_T(t)=1\). At \(t=0\), \(S_T(0)=P(T>0)=1\), because the model puts no probability mass at zero. Thus

$$ S_T(t)= \begin{cases} 1, & t<0,\\ e^{-\lambda t}, & t\ge 0. \end{cases} $$

For every \(t\), the two functions satisfy \(F_T(t)+S_T(t)=1\). The strict inequality in the survival definition causes no discrepancy here: \(P(T=t)=0\) for every fixed \(t\), so \(P(T>t)=P(T\ge t)\) under this continuous model.

Worked Example: Cumulative Probability and Survival at a Fixed Time

A hypothetical equipment component has operating time \(T\) that follows an exponential distribution with rate \(\lambda=0.4\) per operating year. Find the probability that it operates for at most 2 years and the probability that it operates for more than 2 years. Give probabilities to four decimal places.

Method. Since 2 years is nonnegative, use \(F_T(t)=1-e^{-\lambda t}\) for the first probability and \(S_T(t)=e^{-\lambda t}\) for the second. The exponent is dimensionless: \((0.4\text{ per year})(2\text{ years})=0.8\).

The probability of operating for at most 2 years is

$$ \begin{aligned} P(T\le2) &=F_T(2)\\ &=1-e^{-(0.4)(2)}\\ &=1-e^{-0.8}\\ &=1-0.449328964\\ &=0.550671036\\ &\approx 0.5507. \end{aligned} $$

The probability of operating for more than 2 years is

$$ \begin{aligned} P(T>2) &=S_T(2)\\ &=e^{-(0.4)(2)}\\ &=e^{-0.8}\\ &=0.449328964\\ &\approx 0.4493. \end{aligned} $$

Check. The rounded probabilities sum to \(0.5507+0.4493=1.0000\), consistent with the complementary events. The common mistake to avoid is using \(e^{-0.8}\) for the probability of operating for at most 2 years; that is the survival probability, not the CDF.

Density, CDF, and Boundary Conditions

The CDF is a probability accumulated up to a threshold; the density is not itself a probability. For an absolutely continuous distribution, a density equals the derivative of the CDF almost everywhere. A density may be changed at a single point without changing any probability, so the general statement is not that a density must equal the derivative at every point.

For this exponential model, the displayed density \(\lambda e^{-\lambda t}\) equals the derivative of the displayed CDF at every \(t>0\). At the endpoint \(t=0\), treat the boundary separately: the CDF’s left and right derivatives do not agree there under the piecewise definition, and the value assigned to the density at that single point does not change the distribution. In particular, the density value \(f_T(t)\) has units of inverse time and should not be interpreted as \(P(T=t)\).

Worked Example: Boundary Values and a Density Value

A hypothetical sensor’s time \(T\) until a specified event follows an exponential distribution with rate \(\lambda=0.25\) per month. Find \(F_T(t)\) at \(t=-1\), \(t=0\), and \(t=2\) months. Also find the density at 2 months. Give probabilities to four decimal places and the density to four decimal places per month.

Method. Use the negative-time branch for \(t=-1\), the endpoint value for \(t=0\), and the positive-time formulas for \(t=2\).

Because \(-1<0\), the model assigns no probability to \(T\le-1\):

$$ F_T(-1)=0. $$

At zero, the positive-time expression gives

$$ F_T(0)=1-e^{-(0.25)(0)}=1-1=0. $$

At 2 months, \((0.25\text{ per month})(2\text{ months})=0.5\), so

$$ \begin{aligned} F_T(2) &=1-e^{-0.5}\\ &=1-0.606530660\\ &=0.393469340\\ &\approx 0.3935. \end{aligned} $$

The density at 2 months is

$$ \begin{aligned} f_T(2) &=(0.25)e^{-(0.25)(2)}\\ &=(0.25)(0.606530660)\\ &=0.151632665\\ &\approx 0.1516\text{ per month}. \end{aligned} $$

Check. The CDF values are between zero and one and do not decrease as time increases: \(0=F_T(-1)=F_T(0)<F_T(2)=0.3935\). The density value is not a probability and is not reported as one; its units are per month. The common mistake is to substitute a negative time into \(1-e^{-\lambda t}\), which would give an invalid negative CDF instead of using the zero branch.

Using the Formulas Reliably

Before evaluating either function, identify the threshold’s sign and confirm that the time units match the rate. The exponent \(-\lambda t\) must be dimensionless. If time is in days, for instance, the rate used in the product must be per day. Unit consistency is part of the model calculation, not merely a reporting detail.

A useful check is to compare the result with the expected direction. For \(t\ge0\), \(F_T(t)=1-e^{-\lambda t}\) is between zero and one and increases with \(t\). In contrast, \(S_T(t)=e^{-\lambda t}\) is between zero and one and decreases with \(t\). Their sum must equal one at every threshold. These checks catch a reversed sign, an incorrect branch, or an accidental interchange of the functions.

1
Check the support.
If \(t<0\), then \(F_T(t)=0\) and \(S_T(t)=1\); do not use the positive-time formulas.
2
Match the units.
Multiply the rate per time unit by a time in that same unit, so \(\lambda t\) is dimensionless.
3
Choose the function that matches the event.
Use \(F_T(t)\) for \(T\le t\) and \(S_T(t)\) for \(T>t\).
4
Check the result.
For nonnegative \(t\), verify that both probabilities lie between zero and one and sum to one.

Common Sign and Interpretation Traps

A frequent integration error is to write

$$ \int_0^t \lambda e^{-\lambda u}\,du=e^{-\lambda t}-1. $$

This expression is negative for \(t>0\), so it cannot be a probability. The antiderivative is \(-e^{-\lambda u}\), not \(e^{-\lambda u}\). Applying the limits in order gives

$$ \left[-e^{-\lambda u}\right]_0^t =-e^{-\lambda t}+1 =1-e^{-\lambda t}. $$

Another error is to call \(e^{-\lambda t}\) the CDF. It is the survival function for nonnegative \(t\); the CDF is its complement. A final trap is to read a density value as the probability of an exact event time. For a continuous distribution, \(P(T=t)=0\), even though \(f_T(t)\) can be positive.

Key takeaway: For an exponential waiting time with rate \(\lambda>0\), the CDF is \(0\) below zero and \(1-e^{-\lambda t}\) at nonnegative \(t\); the survival function is \(1\) below zero and \(e^{-\lambda t}\) at nonnegative \(t\). Check the support, units, sign, and complement relationship before finalizing a calculation.

Check Your Understanding

For each problem, identify the relevant branch or function before evaluating it. Probabilities are rounded to four decimal places where needed.

  1. A hypothetical inspection wait \(T\) is exponential with rate \(0.5\) per hour. Find \(P(T\le2)\). (Answer: \(0.6321\))
  2. A hypothetical repair time \(T\) is exponential with rate \(0.08\) per day. Find the probability the repair takes more than 10 days. (Answer: \(0.4493\))
  3. A hypothetical alert time \(T\) is exponential with rate 3 per week. Find \(F_T(1/3)\) and \(f_T(1/3)\), stating the density’s units. (Answer: \(F_T(1/3)=0.6321\); \(f_T(1/3)=1.1036\) per week)
  4. A hypothetical waiting time \(T\) is exponential with rate \(0.6\) per hour. Find \(F_T(-0.5)\), \(S_T(-0.5)\), and \(F_T(0)\). (Answer: \(0,\ 1,\ 0\), respectively)
  5. A hypothetical device’s operating time \(T\) is exponential with rate \(0.25\) per month. Find the CDF and survival probabilities at 4 months and verify that they sum to one. (Answer: \(F_T(4)=0.6321\), \(S_T(4)=0.3679\); sum \(=1.0000\))
Clinical Biostats