5.5.3 Weibull Fit

Summary

Weibull fit is a kind of parameter method to analyze the relationship between the survival function and failure time. After analysis, we can get parameter estimates, which can determine survival function and hazard function of Weibull distribution .

Weibull distribution:

f(x)=\frac{c}{\sigma }({\frac{x-\theta }{\sigma }})^{c-1}\exp (-\left({\frac{x-\theta }{\sigma }}\right)^{c}) where x >\theta, for c,\sigma > 0

Survival function:

S(x)=\exp(-\left({\frac{x-\theta }{\sigma }}\right)^{c})

Hazard function:

h(x)=\frac{c}{\sigma }\left({\frac{x-\theta }{\sigma }}\right)^{c-1}

where c is the shape parameter, \sigma is the scale parameter, and \theta is the location parameter. In Origin, Weibull fit only discusses c and \sigma, and assumes \theta = 0.

If c > 1 the hazard increases , if c = 1 the hazard is constant(exponential model), if c < 1 the hazard decreases.

Minimum Origin Version Required: OriginPro 9.1 SR0

What you will learn

This tutorial will show you how to:

  • Perform Weibull Fit
  • How to explain analysis report

Perform Weibull Fit

  1. Start with an empty worksheet. Select Help: Open Folder: Sample Folder... to open the "Samples" folder. In this folder, open the Statistics subfolder and find the file weibull fit.dat. Drag-and-drop this file into the empty worksheet to import it.
  2. Select Statistics: Survival Analysis: weibull fit to open the dialog.
  3. Put the A(X) column into Time Range. Similarly, put the B(Y) column into Censor Range.
  4. Select 0 as the censoring value form censoring value(s) drop-down list.
  5. Expand Plots branch and check Survival Plot and Hazard Plot.
    W01.png
  6. Click the OK button to perform the Weibull fit analysis.
W04.png

Interpreting the Result

Go to worksheet WeibullFit1 for the analysis report.

  • From the "Summary of Event and Censored Values" table , we can see that censored =19 and percent Censored =0.2111.
W02.png
  • From the "Analysis of parameter estimates", we can get all parameter estimates for Weibull distribution.
    Intercept=\theta=4.1959, (\theta is the intercept of the small extreme distribution, \theta = ln(Weibull Scale))
    Weibull Scale= \sigma =66.4153,
    Weibull Shape=c=2.0204,
    Scale =0.495 (scale = 1/c).
W03.png
c > 1,so we can conclude that the hazard increases as time.
  • Furthermore, we can get the survival function and hazard function:
W05.png W06.png