Sorry about that. My code and more of the data is provided below.
proc mixed data=be.time04 method=reml covtest;
class acrostic;
model weight = spmrkt_grcy psi time county spmrkt_grcy*time spmrkt_grcy*time*psi /solution ddfm=sat;
random intercept / subject = acrostic type=un;
lsmeans spsis / e diff cl adjust=bon;
run;
Psi timepoints | Psi |
Baseline_arm_1 | 125 |
Month_4_arm_1 | 88 |
Baseline_arm_1 | 110 |
Month_4_arm_1 | 110 |
Baseline_arm_1 | 59 |
Month_4_arm_1 | 76 |
Baseline_arm_1 | 103 |
Month_4_arm_1 | 113 |
Baseline_arm_1 | 84 |
Month_4_arm_1 | 92 |
Baseline_arm_1 | 117 |
Month_4_arm_1 | 92 |
Baseline_arm_1 | 84 |
Month_4_arm_1 | 92 |
Baseline_arm_1 | 123 |
Month_4_arm_1 | 118 |
Baseline_arm_1 | 123 |
Month_4_arm_1 | 131 |
Baseline_arm_1 | 90 |
Month_4_arm_1 | 87 |
Month_4_arm_1 | 125 |
Baseline_arm_1 | 123 |
Month_4_arm_1 | 131 |
Baseline_arm_1 | 90 |
Month_4_arm_1 | 87 |
Baseline_arm_1 | 125 |