I have a longitudinal dataset with 1000 observations. Each participant has 2 to 10 timepoints repeated measurement of BMI. /*id is the unique number for each participant*/ /*here I want to exam the different decline among people who survive and died separately*/ proc mixed data=bmidata PLOTS(MAXPOINTS= 20000) NOCLPRINT; class t id death_status (ref='0'); model bmi = death_status back_timescale death_status*back_timescale age_at_time0 back_timescale_square sex educ /outp=prediction solution ddfm =bw; REPEATED t/ type=vc subject=id r rcorr; random intercept back_timescale /type=VC subject=projid ; run;   I got the results.   So now, I could write the formulation about BMI and these variables. I am wondering how to calculate the BMI and 95% confidence interval if I valued death_status = 1, back_timescale=-5, age_at_time0 = 75 (around the mean of that timepoint), sex = 1 (female), education = 18 .   Furthermore, I want to compare two BMI values with different death_status but same values for other variables.         ... View more

I used proc reg and want to output parameterestimates.  The first time, it works well and output the table I want. 74 ods output ParameterEstimates=pest2; 75 ods listing close; 76 proc reg data=rush.slope_cha_merge_mri; 77 model total_GMV_n =bmi_slope apoe_e4/clb; 78 run; NOTE: The data set WORK.PEST2 has 3 observations and 10 variables. However,  when I rerun the code, there is a warning and parameterestimates could not be output. And the warning is still here even I changed the name of output table.   Does anyone know how to resolve this?   80 ods output ParameterEstimates=pest2; 81 ods listing close; WARNING: Output 'ParameterEstimates' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used. WARNING: The current ODS SELECT/EXCLUDE/OUTPUT statement was cleared because the end of a procedure step was detected. Probable causes for this include the non-termination of an interactive procedure (type quit; to end the procedure) and a run group with no output. NOTE: PROCEDURE REG used (Total process time): real time 14.61 seconds cpu time 1.01 seconds 82 proc reg data=rush.slope_cha_merge_mri; 83 model total_GMV_n =bmi_slope apoe_e4/clb; 84 run; 85 86 ods output ParameterEstimates=pest3; 87 ods listing close; WARNING: Output 'ParameterEstimates' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used. WARNING: The current ODS SELECT/EXCLUDE/OUTPUT statement was cleared because the end of a procedure step was detected. Probable causes for this include the non-termination of an interactive procedure (type quit; to end the procedure) and a run group with no output. NOTE: PROCEDURE REG used (Total process time): real time 21.98 seconds cpu time 0.71 seconds 88 proc reg data=rush.slope_cha_merge_mri; 89 model total_GMV_n =bmi_slope apoe_e4/clb; 90 run; ... View more

I am analysing the slope of one continuous_variable during follow-up across groups. In the following code, t=fu_year. I used PROC MIXED. In the procedure, I used ESTIMATE to calculate the difference between slopes of the corresponding group and the reference group (group 4 here).   proc mixed data=dataset order=formatted PLOTS(MAXPOINTS= 100000); class t group (ref='4'); model continuous_var=fu_year group group*fu_year/ solution ddfm=kr; repeated t /subject=projid_mixed type=un; estimate "slope difference for group 1, ref=group4" fu_year 0 group*fu_year 1 0 0 -1 0 0/cl; estimate "slope difference for group 2, ref=group4" fu_year 0 group*fu_year 0 1 0 -1 0 0/cl; estimate "slope difference for group 3, ref=group4" fu_year 0 group*fu_year 0 0 1 -1 0 0/cl; estimate "slope difference for group 5, ref=group4" fu_year 0 group*fu_year 0 0 0 -1 1 0/cl; estimate "slope difference for group 6, ref=group4" fu_year 0 group*fu_year 0 0 0 -1 0 1/cl; run;   The point estimate is supposed to be same from the table named "solution for fixed effects"  and the table named "estimate". however, my results varied in these two tables.   In the "solution for fixed effects" Group 1,2,3 have a significantly faster decline (fu_year*group) over time than group 4 .     However, in the table for Estimate, both the estimate and the P-value changed. Is anyone know why the estimates were different between those two tables? Is this normal?   ... View more

I want to calculate the HR and 95%CI with one unit increase in variable A using proc surveyphreg. However, there is no output for the confidence interval.   In some files, they suggested to use HAZARDRATIO to output CI, however, there is no this function in my SAS version.   Is there any other solutions? Thanks in advance. ... View more