topic Re: Global p-value of categorial variable after multiple imptutation??? in Statistical Procedures

Global p-value of categorial variable after multiple imptutation???

Alib — Wed, 13 Sep 2017 09:19:16 GMT

proc mi data=table1 seed=1305417 nimpute=10 out=imputedtable1;
   class slevel sex period city death;
   fcs logistic(slevel/link=glogit) logistic(period/link=glogit) reg(rate/details);
   var slevel sex period city death rate time_log ;
run;


Proc phreg data=pop1_imput1;
	format slevel slevel.;
	class slevel (param=ref ref='Primary');
	model delriskc1*riskc1(0)=slevel/eventcode=1;
by _Imputation_; 
ods output parameterEstimates  = a_mvn;
	run;

proc mianalyze parms(classvar=classval)=a_mvn;
class slevel;
modeleffects slevel;
RUN;
run;

data u_slevel2; set u_slevel1;
HR_est=EXP(ESTIMATE);
LCI_HR=HR_est*EXP(-1.96*STDERR);
UCI_HR=HR_est*EXP(+1.96*STDERR);
RUN;

proc print data=u_slevel2;
var Parm slevel Probt HR_est LCI_HR UCI_HR;
run;

After used the syntax above, I've had p-value and hazards ratio by categories of slevel but i need a global p-value of slevel (p-value for all categories), how i can get it?

Thanks

Re: Global p-value of categorial variable after multiple imptutation???

SAS_Rob — Wed, 13 Sep 2017 12:11:21 GMT

In this case it is recommended that you combine the actual Chi-Square statistics from the Type3 tests in PHREG. This cannot be done with Proc MIANALYZE, but there is a MACRO on Dr. Paul Allison's web site called %COMBCHI that can do this directly.
http://www.ssc.upenn.edu/~allison/#Macros

Re: Global p-value of categorial variable after multiple imptutation???

Alib — Wed, 13 Sep 2017 17:06:47 GMT

Many Thanks Rob