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
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
Many Thanks Rob
April 27 – 30 | Gaylord Texan | Grapevine, Texas
Walk in ready to learn. Walk out ready to deliver. This is the data and AI conference you can't afford to miss.
Register now and lock in 2025 pricing—just $495!
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.