1          OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 72         
 73         /* monkey1.sas (2018 version) */
 74         title 'Primate hippocampal function: Zola-Morgan and Squire (1990)';
 75         /* Science, Vol. 250 (12 Oct. 1990) , Pages 288-290 */
 76         title2 'Multivariate approach to repeated measures (within-cases)';
 77         
 78         data memory;
 79              infile '/folders/myfolders/441s18/Lecture/monkey.data.txt' firstobs=2;
 80              input monkey $ treatmnt $ week2 week4 week8 week12 week16;
 81         
 
 NOTE: The infile '/folders/myfolders/441s18/Lecture/monkey.data.txt' is:
       Filename=/folders/myfolders/441s18/Lecture/monkey.data.txt,
       Owner Name=root,Group Name=vboxsf,
       Access Permission=-rwxrwx---,
       Last Modified=08Mar2018:12:50:29,
       File Size (bytes)=854
 
 NOTE: 18 records were read from the infile '/folders/myfolders/441s18/Lecture/monkey.data.txt'.
       The minimum record length was 42.
       The maximum record length was 42.
 NOTE: The data set WORK.MEMORY has 18 observations and 7 variables.
 NOTE: DATA statement used (Total process time):
       real time           0.00 seconds
       cpu time            0.00 seconds
       
 82         proc means data=memory mean;
 
 83              class treatmnt;
 84              var week2 -- week16;
 85         run;
 
 NOTE: There were 18 observations read from the data set WORK.MEMORY.
 NOTE: PROCEDURE MEANS used (Total process time):
       real time           0.05 seconds
       cpu time            0.05 seconds
       
 
 86         
 87         proc glm data=memory;
 88              class treatmnt;
 89              model week2 -- week16 = treatmnt;
 90              repeated time profile /  short summary mean;
 91         run;
 
 NOTE: The Huynh-Feldt epsilon and the corresponding adjusted p-value have been enhanced to include a correction based on Lecoutre 
       (1991). Use the UEPSDEF=HF option on the REPEATED statement to revert to the previous definition.
 92         
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           4.78 seconds
       cpu time            0.72 seconds
       
 93         proc glm data=memory noprint;
 
 94              title3 'Replicate test for main effect of treatment: F=8.08, p=0.0118';
 95              class treatmnt;
 96              model week2 -- week16 = treatmnt;
 97              manova H = treatmnt
 98                     M = week2+week4+week8+week12+week16 / short;
 99              /* M is a matrix of coefficients for transforming the DVs */
 100        run;
 
 101        
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.00 seconds
       cpu time            0.00 seconds
       
 102        proc glm data=memory noprint;
 
 103             title3 'Replicate tests for main effect of time: Lambda=0.84009249';
 104             title4 'And time by treatment interaction: Lambda=0.44106117';
 105             class treatmnt;
 106             model week2 -- week16 = treatmnt;
 107             manova H = intercept treatmnt
 108                    M = week2-week4, week4-week8, week8-week12, week12-week16
 109                       / short;
 110        run;
 
 111        
 112        /* But the real point is that the treatment only affects recent memories,
 113        not older ones. A basic MANOVA is really more to the point. Follow up
 114        with Bonferroni-corrected univariate tests. 0.05/5 = 0.01. */
 115        
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.00 seconds
       cpu time            0.01 seconds
       
 116        proc glm data=memory;
 
 117             title3 'MANOVA, no repeated measures';
 118             class treatmnt;
 119             model week2 -- week16 = treatmnt;
 120             manova h = treatmnt;
 121        run;
 
 122        
 123        
 124        OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 137