« View all new features in
Mathematica
9
◄
previous
|
next
►
New in
Mathematica
9
›
Survival Analysis
Estimate Confidence Limits and Bands
Use
SurvivalModelFit
to compute pointwise, equal-precision, and Hall-Wellner 95% confidence bands for the survival curve for some censored data.
In[1]:=
X
\[ScriptCapitalD] = EventData[{{2081, \[Infinity]}, {1602, \[Infinity]}, {1496, \ \[Infinity]}, {1462, \[Infinity]}, {1433, \[Infinity]}, {1377, \ \[Infinity]}, {1330, \[Infinity]}, {996, \[Infinity]}, {226, \ \[Infinity]}, {1199, \[Infinity]}, {1111, \[Infinity]}, {530, \ \[Infinity]}, {1182, \[Infinity]}, {1167, \[Infinity]}, 418, 383, 276, 104, 609, 172, 487, 662, 194, 230, 526, 122, 129, 74, 122, 86, 466, 192, 109, 55, 1, 107, 110, 332}];
Fit a model to compute a 95% confidence interval about the survival at 300.
In[2]:=
X
model = SurvivalModelFit[\[ScriptCapitalD], ConfidenceRange -> {100, 600}]
Out[2]=
In[3]:=
X
model["PointwiseBands", "SF"][300]
Out[3]=