Hi all,
I ran the following system of probit models using proc qlim (find syntax bellow). The journal is requesting I report confidence intervals instead of standard errors for the parameter estimates.
Is there a way to either a) get CI using SAS syntax; or b) formula to generate CI using parameter estimates and SE; or c) an explanation that reviewers would accept as to why my model doesn't produce CI?
Thank you!
Maria
Title "Qulim toddler";
proc qlim data=Parent_survey; where age ne 1;
class inc race;
model M = mon inc race ;
model TM = mon inc race ;
model IF = mon inc race ;
model NDM = mon inc race ;
endogenous IF TM M NDM ~ discrete;
test M.inc_1 = M.inc_2;
test M.race_1 = M.race_2;
test M.race_1 = M.race_3;
test M.race_1 = M.race_4;
test M.race_2 = M.race_3;
test M.race_2 = M.race_4;
test M.race_3 = M.race_4;
test TM.inc_1 = TM.inc_2;
test TM.race_1 = TM.race_2;
test TM.race_1 = TM.race_3;
test TM.race_1 = TM.race_4;
test TM.race_2 = TM.race_3;
test TM.race_2 = TM.race_4;
test TM.race_3 = TM.race_4;
test IF.inc_1 = IF.inc_2;
test IF.race_1 = IF.race_2;
test IF.race_1 = IF.race_3;
test IF.race_1 = IF.race_4;
test IF.race_2 = IF.race_3;
test IF.race_2 = IF.race_4;
test IF.race_3 = IF.race_4;
test NDM.inc_1 = NDM.inc_2;
test NDM.race_1 = NDM.race_2;
test NDM.race_1 = NDM.race_3;
test NDM.race_1 = NDM.race_4;
test NDM.race_2 = NDM.race_3;
test NDM.race_2 = NDM.race_4;
test NDM.race_3 = NDM.race_4;
output out=pred_tot_main predicted prob proball;
run ;
quit;
Title "predictions toddler";
proc tabulate data= pred_baby_main;
class inc race /missing;
var Prob2_IF Prob2_TM Prob2_M Prob2_NDM;
table inc,
Prob2_IF*(MEAN)
Prob2_TM*(MEAN)
Prob2_M*(MEAN)
Prob2_NDM*(MEAN);
table race,
Prob2_IF*(MEAN)
Prob2_TM*(MEAN)
Prob2_M*(MEAN)
Prob2_NDM*(MEAN) ;
run;
proc tabulate data= pred_tot_main;
class inc race /missing;
var
Prob2_IF
Prob2_TM
Prob2_M
Prob2_NDM
;
table
Prob2_M*(MEAN)
Prob2_TM*(MEAN)
Prob2_IF*(MEAN)
Prob2_NDM*(MEAN),inc
;
table
Prob2_M*(MEAN)
Prob2_TM*(MEAN)
Prob2_IF*(MEAN)
Prob2_NDM*(MEAN),race
;
run;
