turn on suggestions

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

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- P-value for Stratified log-rank test

Topic Options

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

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 10:55 AM

Hello,

I'm trying to derive the p-value for a stratified log-rank test using lifetest. The specs are such:

A stratified log-rank test will be used to test the equivalence

hypothesis of survival distributions. The stratification factors are

the randomization stratification factors. The LIFETEST procedure in

SAS with the TIME statement including a variable with survival times

and a (right) censoring variable, with STRATA statement for the

stratifying variable and the GROUP option for identifying treatment

arms will be used to obtain one sided p-value from log-rank score

test. The implementation is as follows.

PROC LIFETEST data=dataset METHOD=KM CONFTYPE=LOGLOG;

TIME survtime*censor(1);

STRATA stratum1 stratum2 stratum3 / GROUP=trt;

RUN;

/* stratum1 stratum2 stratum3 represents the 3 stratum variables (to

be included for stratified analysis only);

survtime represents variable containing event/censor times;

censor represents censoring variable (1=censored, 0=event);

trt represents treatment arm variable; */

I need to capture the ods output into a dataset for reporting. This is the current code:

ODS TRACE ON;

ODS OUTPUT ;

PROC LIFETEST data=pop3_sall METHOD=KM CONFTYPE=LOGLOG ;

Time months*event(0);

STRATA STRVAL1 STRVAL2 STRVAL3 / GROUP=TRT01PN ;

RUN;

ODS OUTPUT CLOSE;

ODS TRACE OFF;

Which ODS output dataset needs to be captured and does it require extra code to derive the p-value?

Thank you!

Accepted Solutions

Solution

11-23-2015
03:54 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 02:17 PM - edited 11-23-2015 04:37 PM

All Replies

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 01:20 PM

I think you want as part of the proc code:

ods output homtests= youroutputset;

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 02:13 PM

Thank you for responding. I originally used "homtests" for the output dataset, but the p-value doesn't seem to make sense based on the table I'm checking. The discription in the footnote says:

"P-value is obtained from the one-sided stratified logrank test". The p-value based on homtests is: 0.0898 and the table reports: 0.045. Any ideas?

Thank you for your help!

"P-value is obtained from the one-sided stratified logrank test". The p-value based on homtests is: 0.0898 and the table reports: 0.045. Any ideas?

Thank you for your help!

Solution

11-23-2015
03:54 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 02:17 PM - edited 11-23-2015 04:37 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-19-2015 10:54 AM

Without your data it's a bit difficult to diagnose any specific values.

For the example procedure I run the values as appearing in the output were in that table, unless you're looking at a different table.

You would have to post the table with the values you're referencing and preferably the data that generated the output. The data need not be real or complete but sufficient to generate example output and should be in the form of datastep code.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-18-2015 01:21 PM

Do you see the value you want in the output from the code above?