1          OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 70         
 71         /* monkey1.sas (2020 version) */
 72         title 'Primate hippocampal function: Zola-Morgan and Squire (1990)';
 73         /* Science, Vol. 250 (12 Oct. 1990) , Pages 288-290 */
 74         title2 'Multivariate approach to repeated measures (within-cases)';
 75         
 76         data memory;
 77              infile '/home/brunner0/441s20/monkey.data.txt' firstobs=2;
 78              input monkey $ treatmnt $ week2 week4 week8 week12 week16;
 79         
 
 NOTE: The infile '/home/brunner0/441s20/monkey.data.txt' is:
       Filename=/home/brunner0/441s20/monkey.data.txt,
       Owner Name=brunner0,Group Name=oda,
       Access Permission=-rw-r--r--,
       Last Modified=01Mar2020:17:38:09,
       File Size (bytes)=854
 
 NOTE: 18 records were read from the infile '/home/brunner0/441s20/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.02 seconds
       user cpu time       0.00 seconds
       system cpu time     0.00 seconds
       memory              764.21k
       OS Memory           30888.00k
       Timestamp           03/01/2020 10:41:02 PM
       Step Count                        122  Switch Count  2
       Page Faults                       0
       Page Reclaims                     161
       Page Swaps                        0
       Voluntary Context Switches        21
       Involuntary Context Switches      0
       Block Input Operations            8
       Block Output Operations           264
       
 
 80         proc means data=memory mean;
 81              class treatmnt;
 82              var week2 -- week16;
 83         run;
 
 NOTE: There were 18 observations read from the data set WORK.MEMORY.
 NOTE: PROCEDURE MEANS used (Total process time):
       real time           0.03 seconds
       user cpu time       0.04 seconds
       system cpu time     0.01 seconds
       memory              9280.03k
       OS Memory           37820.00k
       Timestamp           03/01/2020 10:41:02 PM
       Step Count                        123  Switch Count  1
       Page Faults                       0
       Page Reclaims                     1804
       Page Swaps                        0
       Voluntary Context Switches        10
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           8
       
 
 84         
 85         proc glm data=memory;
 86              class treatmnt;
 87              model week2 -- week16 = treatmnt;
 88              repeated time profile /  short summary mean;
 89         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.
 90         
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           3.56 seconds
       user cpu time       0.55 seconds
       system cpu time     0.03 seconds
       memory              17366.75k
       OS Memory           47680.00k
       Timestamp           03/01/2020 10:41:06 PM
       Step Count                        124  Switch Count  2
       Page Faults                       0
       Page Reclaims                     5563
       Page Swaps                        0
       Voluntary Context Switches        1684
       Involuntary Context Switches      32
       Block Input Operations            0
       Block Output Operations           1888
       
 
 91         proc glm data=memory noprint;
 92              title3 'Replicate test for main effect of treatment: F=8.08, p=0.0118';
 93              class treatmnt;
 94              model week2 -- week16 = treatmnt;
 95              manova H = treatmnt
 96                     M = week2+week4+week8+week12+week16 / short;
 97              /* M is a matrix of coefficients for transforming the Response vars */
 98         run;
 
 99         
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.00 seconds
       user cpu time       0.00 seconds
       system cpu time     0.00 seconds
       memory              1572.75k
       OS Memory           46008.00k
       Timestamp           03/01/2020 10:41:06 PM
       Step Count                        125  Switch Count  3
       Page Faults                       0
       Page Reclaims                     220
       Page Swaps                        0
       Voluntary Context Switches        24
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           288
       
 
 100        proc glm data=memory noprint;
 101             title3 'Replicate tests for main effect of time: Lambda=0.84009249';
 102             title4 'And time by treatment interaction: Lambda=0.44106117';
 103             class treatmnt;
 104             model week2 -- week16 = treatmnt;
 105             manova H = intercept treatmnt
 106                    M = week2-week4, week4-week8, week8-week12, week12-week16
 107                       / short;
 108        run;
 
 109        
 110        /* But the real point is that the treatment only affects recent memories,
 111        not older ones. A basic MANOVA is really more to the point. Follow up
 112        with Bonferroni-corrected univariate tests. 0.05/5 = 0.01. */
 113        
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.00 seconds
       user cpu time       0.01 seconds
       system cpu time     0.01 seconds
       memory              1588.87k
       OS Memory           46008.00k
       Timestamp           03/01/2020 10:41:06 PM
       Step Count                        126  Switch Count  4
       Page Faults                       0
       Page Reclaims                     216
       Page Swaps                        0
       Voluntary Context Switches        34
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           288
       
 
 114        proc glm data=memory;
 115             title3 'MANOVA, no repeated measures';
 116             class treatmnt;
 117             model week2 -- week16 = treatmnt;
 118             manova h = treatmnt;
 119        run;
 
 120        
 121        
 122        OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 133