BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
NancyForget
Calcite | Level 5

Hello all,

Usually for Kaplan-Meier estimates, the client is interested of having the 25%, 50% and 75% quartiles for survival estimates.  Now the client would like to have also the 10% percentile, since we don’t have enough events.

In proc lifetest, we can easily retrieve the 25%, 50% and 75% quartiles with the associated Confidence Interval (CI).

    

Does anyone know how to get the CI for another percentile, let’s say 10%?  From outsurv dataset, we can easily pull out the 10th percentile but not the associated CI to that value.  I do not want the CI for the curve, but the CI associated to the 10th percentile value.

OBSDUR_CENSOR_SURVIVALSDF_LCLSDF_UCL
10.1.0000001.000001.00000
2800.9200000.716390.97937
31300.880000

0.67256

0.95964

41800.8400000.628060.93673
52300.8000000.584450.91146
65200.7600000.542050.88428
76300.6800000.460930.82527
87000.6400000.422150.79378

So here's the 10th percentile is DUR=13 (using loglog).    So how do we get the CI associated to DUR=13?

Usually proc lifestest gives the 25%, 50% and 75% percentile automatically.

Thanks for your help

1 ACCEPTED SOLUTION

Accepted Solutions
JacobSimonsen
Barite | Level 11

I dont think you can get a CI of the percentile that is estimated by a non-parametric method like the KM-curve. That method is aiming on giving estimates and CI in terms of an interval around the curve, not an interval on the time-axis.

A way to get an interval on the time axis is to use an accelarated failure model instead. Here you model the time to event by assuming some parametric distribution on the time-to-event. See the documentation for PROC LIFEREG. For example, if you assume a weibull distribution on the time-to event, you can then find an analytic expression for the 10th percent percentile. Then you can plug in the confidence limits of the parameters into that expression.

View solution in original post

7 REPLIES 7
Reeza
Super User

You only have estimates for time points that you have. Usually estimates are reported for specific time points, not percentiles.

That being said the usual protocol is to pick the time closes to the time of interest and use that as the estimate which is in line with the step curve nature of KM estimates.  Extending that to percentiles, pick the time point that's closest to 10th percentile and use that. 

NancyForget
Calcite | Level 5

Yes I am ok with getting the time point closest to 10th percentile.  In the above example, it will be DUR=13.

Now, I need to calculate the 95% CI of that DUR=13. 

Another example from SAS (http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_lifetest_sec...).  On this link, you will find the Quartile Estimates,  The quartiles 25%, 50% and 75% are given in proc lifestest (see below on print screen).   But if I want the 10th percentile, how can I get the associated 95% CI?

KM quartiles estimates.jpg

NancyForget
Calcite | Level 5

The oursurv gives the CI of the curve, not the CI of the percentile.  For example, on this print screen below, the 10th percentile is AVAL=4.  The SDF_LCL and SDF_UCL is the CI of the survival estimate (here 0.8947368421) and not the CI of AVAL=4.

I'm looking for the CI of the percentile, not the curve.  Thanks

outsur v.jpg

Reeza
Super User

Ok, you want a CI on the survival time, the 10th percentile of the time of survival.

Your first sentence says survival estimates but those are time points not survival function estimates therefore my confusion.

I think this is the description of what you're looking for and I *think* it has to be calculated manually Smiley Sad

SAS/STAT(R) 9.2 User's Guide, Second Edition

JacobSimonsen
Barite | Level 11

I dont think you can get a CI of the percentile that is estimated by a non-parametric method like the KM-curve. That method is aiming on giving estimates and CI in terms of an interval around the curve, not an interval on the time-axis.

A way to get an interval on the time axis is to use an accelarated failure model instead. Here you model the time to event by assuming some parametric distribution on the time-to-event. See the documentation for PROC LIFEREG. For example, if you assume a weibull distribution on the time-to event, you can then find an analytic expression for the 10th percent percentile. Then you can plug in the confidence limits of the parameters into that expression.

NancyForget
Calcite | Level 5

Correct.  Does anyone know how to calculate this manually?

THanks in advance

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 7 replies
  • 9989 views
  • 0 likes
  • 3 in conversation