Hello, all   I used PROC LOGISTICS to run ordinal logistic models on 5 multiply imputed data (My outcome is ordinal with 3 categories).   I also included "the unequalslopes" function in order to check the proportional odds assumption. I was able to get pooled results on the 5 ordinal logit analyses but not on the 5 assumption test results. Is there a way to obtain pooled results on the proportional assumption tests (indicated in SAS output as "Linear Hypothesis Test Results").    I would greatly appreciate your help.    Here is my code:    proc logistic data=name; by _imputation_; model DV (order=data) = var1 var2 / unequalslopes MAXITER =1000 ; var1: test var1_4 = var1_3; var2: test var2_4 = var2_3; run; data perm.name2print; set perm.nameprint; parameter = compress (parameter); run; proc mianalyze parms (classvar=classval)= perm.name2print; modeleffects intercept var1 var2 ; ods output parameterestimates = mianalyze_parms; run; data OR; set mianalyze_parms; OR=exp(estimate); LCL_OR=exp(LCLMean); UCL_OR=exp(UCLMean); proc print; var parm OR LCL_OR UCL_OR; run; ods rtf close; ... View more

Dear all,    I am trying to create a time-dependent variable in a cox proportional hazard model, using the counting process, which is guided by http://support.sas.com/resources/papers/proceedings12/168-2012.pdf.    Unfortunately, I keep seeing this error message below.  ERROR 161-185: No matching DO/SELECT statement.   579 data perm.habc_cox10; 580 set perm.habc_cox9; 581 582 array symptotal_(*) symptotal_1 - symptotal_12; 583 array chng (11); 584 t= 1; 585 do i = 2 TO 12; 586 if symptotal_(i) NE symptotal_(i-1) THEN DO; /*NE = not equal to */ 587 chng(t) = i-1; 588 t = t+1; 589 end; 590 end; 591 run; 592 data perm.habc_cox11; 593 set perm.habc_cox10; 594 595 array symptotal_(*) symptotal_1 - symptotal_12; 596 array chng (*) chng1 - chng11; 597 start = 0; 598 censor2 = 0; 599 t = 1; 600 601 do i = 1 to time_event; 602 603 if chng(t) >. and (chng(t) < time_event) or i=time_event then; 604 605 if chng(t)>. then burden = symptotal_(chng(t)); 606 else burden = symptotal_(time_event); 607 608 stop = min(chng(t), time_event); 609 610 if i = time_event then censor2 = censor; 611 612 if t>1 then start = chng(t-1); 613 614 t = t+1; 615 output; 616 end; 617 end; --- 161 ERROR 161-185: No matching DO/SELECT statement. 618 run; I would really appreciate your help.    Thank you,      ... View more

Hello, all!   Do anyone know how to get pulled odds ratios for the PROC PLM procedure on "multiply imputed data"?    1) I just ran PROC MI to do multiple imputation;  2) ran PROC SURVEYLOGISTIC to run the model;  3) ran PROC MIANALYZE to get pulled results (including odds ratios).  4) Since I have 2 interaction terms in the model, I ran PROC PLM to conduct simple slope analyses;  5) ran PROC MIANALYZE to get pulled results for Step 4.    After Step 5, I only got pulled estimates but NOT pulled odds ratios. How do I get pulled odds ratios for Step 5? I would really really appreciate if you could share your insight with me!          ... View more

Dear all,    I have two interaction terms in my model and would like to run a simple slope analysis using PROC PLM. I already got pooled results for the mail effect but am having trouble with the PROC PLM part. Is there a way to get pooled results for PROC PLM on imputed data?    Please see my code below. Please note that I used PROC SURVEYLOGISTIC because the data is weighted and the outcome is ordinal.    (proc mi...........truncated)   proc surveylogistic data=perm.19m2; strata RAESTRAT; cluster RAEHSAMP; weight R11WTRESP; by _imputation_; model DV_r (order=data) = R11AGEY_B | stress  RAEDYRS h11itot self_health R11WGIWK R11JCTEN computer pres_retire  young_promo less_demand physical_eff gender race2_b race2_his race2_o occupation_1 occupation_2 occupation_3 reduce_hrs part_full /link=clogit clodds expb rsq; ods output parameterestimates = perm.19m2print; store perm.test; run;    (proc mianalyze....truncated)  ods graphics on; proc plm restore = perm.test; estimate 'stress slope, R11AGEY_B = mean-sd' stress 1 stress*R11AGEY_B 2.71 RAEDYRS 1 h11itot 1 self_health 1 R11WGIWK 1 R11JCTEN 1 computer 1 pres_retire 1 young_promo 1 less_demand 1 physical_eff 1 gender 1 race2_b 1 race2_his 1 race2_o 1 occupation_1 1 occupation_2 1 occupation_3 1 reduce_hrs 1 part_full 1, 'stress slope, R11AGEY_B = mean' stress 1 stress*R11AGEY_B 2.73 RAEDYRS 1 h11itot 1 self_health 1 R11WGIWK 1 R11JCTEN 1 computer 1 pres_retire 1 young_promo 1 less_demand 1 physical_eff 1 gender 1 race2_b 1 race2_his 1 race2_o 1 occupation_1 1 occupation_2 1 occupation_3 1 reduce_hrs 1 part_full 1, 'stress slope, R11AGEY_B = mean+sd' RAEDYRS 1 stress 1 stress*R11AGEY_B 2.75 RAEDYRS 1 h11itot 1 self_health 1 R11WGIWK 1 R11JCTEN 1 computer 1 pres_retire 1 young_promo 1 less_demand 1 physical_eff 1 gender 1 race2_b 1 race2_his 1 race2_o 1 occupation_1 1 occupation_2 1 occupation_3 1 reduce_hrs 1 part_full 1 /e ; run; ods graphics off;   (and then how to pool these results?)    I would really appreciate your help!       ... View more