Oh my friend, I finally got the output results without error by using your code. However, I have three questions about the existing code. Firstly, although the code can run properly, I am not quite sure what does the number 1 after the variable name indicate in the estimate statement, and how does the content after the estimate statement inform SAS to perform a joint test? Secondly, where should the p-value indicating the significance of joint test be located in the results of SAS output? I originally guessed that it should be in the F test at the end of the estimation section. However, in that section, the numerator degrees of freedom is 1 and other parts including p-values are blank (as shown in the figure below)." Thirdly, should the interaction term be included as a covariate in the model statement? The reason why I have this doubt is that I found that the interaction effects obtained by the following two code sections are exactly the same. proc surveyphreg data = work.nhanes_data varmethod = taylor nomcar;
weight WEIGHT;
strata SDMVSTRA;
cluster SDMVPSU;
class demo_age demo_sex;
model mor_peryear*mor_mortstat(0) = homa_ir demo_age demo_sex;
estimate homa_ir 1, demo_age 1 /joint e;
domain overall_selected;
run; proc surveyphreg data = work.nhanes_data varmethod = taylor nomcar;
weight WEIGHT;
strata SDMVSTRA;
cluster SDMVPSU;
class demo_age demo_sex;
model mor_peryear*mor_mortstat(0) = homa_ir demo_age demo_sex;
estimate homa_ir 1, demo_sex 1 /joint e;
domain overall_selected;
run; The only difference between the above two code sections is the variables specified in the estimate statement. So I guess the interaction term may also should be in the model statement as a covariate. Thank you very much for taking the time to answer my questions despite your busy schedule. Thank you so much!
... View more