Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Brookmeyer and Crowley Confidence Intervals?

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 11-02-2019 11:00 PM
(20102 views)

Does anybody know how to generate 95% Confidence Intervals "based on the Brookmeyer and Crowley method" ? I'm programming a table that contains stats that come from PROC LIFETEST, so I was assuming (or maybe just hoping) that there's an ODS data set containing these. I've never heard of Brookmeyer and Crowley CI's before, unfortunately.

Thanks!!!

...dave

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

According to: (page 33 of the PDF) https://support.sas.com/documentation/onlinedoc/stat/142/lifetest.pdf,

"Brookmeyer and Crowley constructed the confidence interval for the median survival time based on the confidence interval for S(t)..."

This method is then extended to estimate any pth percentile. There are numerous options in the LIFETEST procedure but it looks like the confidence intervals for the quartile estimates use the extended Brookmeyer and Crowley method. The ODS table is called Quartiles. Below is a sample from the sas documentation:

```
ods output Quartiles;
proc lifetest
data=sashelp.BMT
conftype=linear;
time T * Status(0);
strata Group /
test=logrank
adjust=sidak;
run;
```

I just added the option CONFTYPE=LINEAR to see if this code would correspond to what was shown in the linked PDF above. If you run this code and compare observation 3 in the output to the table at the bottom of page 34 of the PDF, you'll see that the LIFETEST procedure uses the described methodology by default. It's up to you to determine which g-transformation (i.e. linear, logit, etc.) to use in generating these estimates.

best,

-unison

-unison

5 REPLIES 5

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

According to: (page 33 of the PDF) https://support.sas.com/documentation/onlinedoc/stat/142/lifetest.pdf,

"Brookmeyer and Crowley constructed the confidence interval for the median survival time based on the confidence interval for S(t)..."

This method is then extended to estimate any pth percentile. There are numerous options in the LIFETEST procedure but it looks like the confidence intervals for the quartile estimates use the extended Brookmeyer and Crowley method. The ODS table is called Quartiles. Below is a sample from the sas documentation:

```
ods output Quartiles;
proc lifetest
data=sashelp.BMT
conftype=linear;
time T * Status(0);
strata Group /
test=logrank
adjust=sidak;
run;
```

I just added the option CONFTYPE=LINEAR to see if this code would correspond to what was shown in the linked PDF above. If you run this code and compare observation 3 in the output to the table at the bottom of page 34 of the PDF, you'll see that the LIFETEST procedure uses the described methodology by default. It's up to you to determine which g-transformation (i.e. linear, logit, etc.) to use in generating these estimates.

best,

-unison

-unison

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Thanks, Unison! I will look into that as well!

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I found a macro online which actually calculates these things when passed a standard TTE data set.

Written by Xiao Chen, it ran flawlessly. I passed the results to the statistician on my project and will report his comments back here.

Brookmeyer-Crowley 95% Confidence Interval

Lower Upper

newgroup total Estimate Limit Limit bc_lower bc_upper

1 88 155.000 154.000 156.000 154

2 175 120.000 115.000 121.000 115

3 70 99.000 75.000 155.000 75 130

4 194 121.000 117.000 121.000 117

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

That's great, Dave. Looking forward to it.

Best of luck!

-unison

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hello @daveconifer1 and welcome to the SAS Support Communities!

I've just seen that @unison was faster in answering your question (while I was writing the text below), so this is rather a confirmation of what he has written ...

As far as I know, this refers to the 1982 article by R. Brookmeyer and J. Crowley "A Confidence Interval for the Median Survival Time" in *Biometrics* (see the reference in the PROC LIFETEST documentation). This confidence interval (and the corresponding CIs for the other two quartiles) are available in ODS table Quartiles of PROC LIFETEST.

I think for the original CI as per Brookmeyer and Crowley you need to specify the option CONFTYPE=**LINEAR** in the PROC LIFETEST statement because the default has been changed to CONFTYPE=LOGLOG in SAS version 9.1, but this is rather based on a generalization of the original CI. The formula and more details can be found in section "Breslow, Fleming-Harrington, and Kaplan-Meier Methods" of the PROC LIFETEST documentation.

Note (in the procedure output "Quartile Estimates") that the upper confidence limit is not included in the interval, which can be easily overlooked when only the ODS output dataset is used.

**Available on demand!**

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

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.