Greetings, I am trying to estimate the means of a continuous outcome called 'sbpchangemed' at different levels of a categorical variable (called acescore with 5 levels) (exposure) at three different values of a continuous variable (called hsiscore). I found the mean within each tertile of hsiscore. Tertile 1, low, has a mean of 2.0. Tertile 2, medium, has a mean of 10.3, and Tertile 3, high, has a mean of 29.5. I would like to estimate the change in sbpchangemed by acescore at a hsiscore of 2, a hsiscore of 10.3, and a hsiscore of 29.5. I have added the table shell below for context. I obtain estimates when I run the code below. However, I have not interpreted estimate statements in quite some time. Question: Could someone please review the code below and assess whether this outputs the table shell data correctly? proc surveyreg data=Analyze; cluster psu2; stratum strata2; weight weight2; domain keep2('1'); class acescore female ethnic incgt30k1 hsedu married insurance1 engpref1 nativity ; model sbpchangemed = acescore age1 followuptime2 female ethnic incgt30k1 hsedu married insurance1 engpref1 nativity bmi1 sbp1med hsiscore acescore*hsiscore /solution clparm anova; estimate "1 vs 0 hsiscore =2" acescore -1 1 0 0 0 hsiscore 2 ; estimate "2 vs 0 hsiscore =2" acescore -1 0 1 0 0 hsiscore 2 ; estimate "3 vs 0 hsiscore =2" acescore -1 0 0 1 0 hsiscore 2 ; estimate "4+ vs 0 hsiscore =2" acescore -1 0 0 0 1 hsiscore 2 ; estimate "1 vs 0 hsiscore =10.3" acescore -1 1 0 0 0 hsiscore 10.3 ; estimate "2 vs 0 hsiscore =10.3" acescore -1 0 1 0 0 hsiscore 10.3 ; estimate "3 vs 0 hsiscore =10.3" acescore -1 0 0 1 0 hsiscore 10.3 ; estimate "4+ vs 0 hsiscore =10.3" acescore -1 0 0 0 1 hsiscore 10.3 ; estimate "1 vs 0 hsiscore = 29.5" acescore -1 1 0 0 0 hsiscore 29.5 ; estimate "2 vs 0 hsiscore = 29.5" acescore -1 0 1 0 0 hsiscore 29.5 ; estimate "3 vs 0 hsiscore = 29.5" acescore -1 0 0 1 0 hsiscore 29.5 ; estimate "4+ vs 0 hsiscore = 29.5" acescore -1 0 0 0 1 hsiscore 29.5 ; run;
... View more