Programming the statistical procedures from SAS

P-value for Stratified log-rank test

Accepted Solution Solved
Reply
Occasional Contributor
Posts: 6
Accepted Solution

P-value for Stratified log-rank test

 

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
Grand Advisor
Posts: 16,842

Re: P-value for Stratified log-rank test

[ Edited ]

You move from a one sided p-value to a 2 sided p-value by mulitplying 2....so that might be what's happening here

View solution in original post


All Replies
Grand Advisor
Posts: 10,026

Re: P-value for Stratified log-rank test

I think you want as part of the proc code:

 

ods output homtests= youroutputset;

Occasional Contributor
Posts: 6

Re: P-value for Stratified log-rank test

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!
Solution
‎11-23-2015 03:54 PM
Grand Advisor
Posts: 16,842

Re: P-value for Stratified log-rank test

[ Edited ]

You move from a one sided p-value to a 2 sided p-value by mulitplying 2....so that might be what's happening here

Grand Advisor
Posts: 10,026

Re: P-value for Stratified log-rank test

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.

Grand Advisor
Posts: 16,842

Re: P-value for Stratified log-rank test

Do you see the value you want in the output from the code above?
☑ This topic is SOLVED.

Need further help from the community? Please ask a new question.

Discussion stats
  • 5 replies
  • 1043 views
  • 0 likes
  • 3 in conversation