BookmarkSubscribeRSS Feed
Alib
Calcite | Level 5
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

2 REPLIES 2
SAS_Rob
SAS Employee

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

hackathon24-white-horiz.png

The 2025 SAS Hackathon has begun!

It's finally time to hack! Remember to visit the SAS Hacker's Hub regularly for news and updates.

Latest Updates

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.

Discussion stats
  • 2 replies
  • 1798 views
  • 0 likes
  • 2 in conversation