Options ps=54 ls=76; Data vasey_multi; Input patno lbrow1 lbrow2 lbrow3 lbrow4; Cards; 1 143 368 345 772 2 142 155 161 178 3 109 167 356 956 4 123 135 137 187 5 276 216 232 307 6 235 386 398 425 7 208 175 207 293 8 267 358 698 771 9 183 193 631 403 10 245 268 572 1383 11 324 507 556 504 12 148 378 342 796 13 130 142 150 173 14 119 171 333 1062 15 102 94 93 69 16 279 204 229 299 17 244 365 392 406 18 196 168 199 287 19 279 358 822 671 20 167 183 731 203 21 345 238 572 1652 22 524 507 520 504 ; Run; Data vasey_uni; set vasey_multi; Keep patno time lbrow; Time = 10; lbrow=lbrow1; output; Time = 2; lbrow=lbrow2; output; Time = 3; lbrow=lbrow3; output; Time = 4; lbrow=lbrow4; output; Run; ods PDF file="psychophysiology_rm.pdf"; ods graphics on; Title "GLM ANOVA and MANOVA"; Proc glm data = vasey_multi; Model lbrow1 lbrow2 lbrow3 lbrow4 = / nouni; Repeated time (1 2 3 4) contrast(1) /short summary printe; Manova/printe; Run; Title "Repeated - Compound Symmetry"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=cs subject=patno r; Run; Title "Repeated - Type H"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=hf subject=patno r; Run; Title "Repeated - AR1"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=AR(1) subject=patno r; Run; Title "Repeated - Unstructured"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=un subject=patno r; Run; Title "Repeated - Heterogeneous Compound Symmetry"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=csh subject=patno r; Run; Title "Repeated - Heterogeneous AR1"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=ARH(1) subject=patno r; Run; Title "Repeated - Heterogeneous Toeplitz"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated / type=toeph subject=patno r; Run; Title "Random Intercepts"; Proc Mixed data = vasey_uni; Class time patno; Model lbrow = time/solution; Repeated intercept/ subject=patno; Run; ods graphics off; ods PDF close;