OPTIONS NOCENTER ; TITLE1 'RIESBY DATA - FIT OF MODELS'; PROC IML; /* random intercepts model */ xt = { 1 1 1 1 1 1 , 0 1 2 3 4 5 }; x = T(xt); zt = { 1 1 1 1 1 1 }; z = T(zt); beta = { 23.5518 , -2.3757 }; var_ups = { 16.1544 }; var_err = { 19.0378 }; means = T(x*beta); varcov = z*var_ups*t(z) + var_err*I(6); sd = sqrt(vecdiag(varcov)); sdt = T(sd); s = 1/sd; corr = diag(s)*varcov*diag(s); PRINT 'random intercepts model'; PRINT 'Estimated Means across time', means [FORMAT=10.4]; PRINT 'Estimated SDs across time', sdt [FORMAT=10.4]; PRINT 'Estimated Correlations', corr [FORMAT=10.4]; /* random intercepts, week, week^2, & week^3 (orthogonalized) */ xt = { 1 1 1 1 1 1 , 0 1 2 3 4 5 , 0 1 4 9 16 25 , 0 1 8 27 64 125 }; x = T(xt)*INV(ROOT(xt*T(xt))); zt = { 1 1 1 1 1 1 , 0 1 2 3 4 5 , 0 1 4 9 16 25 , 0 1 8 27 64 125 }; z = T(zt)*INV(ROOT(zt*T(zt))); beta = { 43.2403 , -9.8813 , 0.3286 , 0.6162 }; var_ups = {110.9063 34.9184 -11.1415 -7.2492 , 34.9184 36.2505 0.2580 -4.3317 , -11.1415 0.2580 9.2360 4.0265 , -7.2492 -4.3317 4.0265 5.0425 }; var_err = { 8.9202 }; means = T(x*beta); varcov = z*var_ups*t(z) + var_err*I(6); sd = sqrt(vecdiag(varcov)); sdt = T(sd); s = 1/sd; corr = diag(s)*varcov*diag(s); PRINT 'random intercepts, week, week^2, & week^3 (orthogonalized)'; PRINT 'Estimated Means across time', means [FORMAT=10.4]; PRINT 'Estimated SDs across time', sdt [FORMAT=10.4]; PRINT 'Estimated Correlations', corr [FORMAT=10.4]; ----------- OUTPUT ----------- RIESBY DATA - FIT OF MODELS random intercepts model Estimated Means across time MEANS 23.5518 21.1761 18.8004 16.4247 14.0490 11.6733 Estimated SDs across time SDT 5.9323 5.9323 5.9323 5.9323 5.9323 5.9323 Estimated Correlations CORR 1.0000 0.4590 0.4590 0.4590 0.4590 0.4590 0.4590 1.0000 0.4590 0.4590 0.4590 0.4590 0.4590 0.4590 1.0000 0.4590 0.4590 0.4590 0.4590 0.4590 0.4590 1.0000 0.4590 0.4590 0.4590 0.4590 0.4590 0.4590 1.0000 0.4590 0.4590 0.4590 0.4590 0.4590 0.4590 1.0000 random intercepts, week, week^2, & week^3 (orthogonalized) Estimated Means across time MEANS 23.5076 21.4816 18.8741 16.1446 13.7523 12.1565 Estimated SDs across time SDT 4.5022 4.7341 5.3118 6.3720 7.0648 7.3175 Estimated Correlations CORR 1.0000 0.5026 0.4079 0.2940 0.2111 0.1382 0.5026 1.0000 0.5659 0.4481 0.3568 0.2906 0.4079 0.5659 1.0000 0.6777 0.6170 0.4612 0.2940 0.4481 0.6777 1.0000 0.7713 0.5799 0.2111 0.3568 0.6170 0.7713 1.0000 0.7023 0.1382 0.2906 0.4612 0.5799 0.7023 1.0000