11-18-2015 10:55 AM
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;
STRATA stratum1 stratum2 stratum3 / GROUP=trt;
/* 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 ;
STRATA STRVAL1 STRVAL2 STRVAL3 / GROUP=TRT01PN ;
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?
11-18-2015 02:13 PM
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.
Need further help from the community? Please ask a new question.