1          OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 72         
 73         /* monkey1.sas */
 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 '/home/u1407221/441s24/data/monkey.data.txt' firstobs=2;
 80              input monkey $ treatmnt $ week2 week4 week8 week12 week16;
 81              label week2  = 'Trained 2 weeks before injection'
 82                    week4  = 'Trained 4 weeks before injection'
 83                    week8  = 'Trained 8 weeks before injection'
 84                    week12 = 'Trained 12 weeks before injection'
 85                    week16 = 'Trained 16 weeks before injection';
 86         
 
 NOTE: The infile '/home/u1407221/441s24/data/monkey.data.txt' is:
       Filename=/home/u1407221/441s24/data/monkey.data.txt,
       Owner Name=u1407221,Group Name=oda,
       Access Permission=-rw-r--r--,
       Last Modified=25Mar2024:13:47:17,
       File Size (bytes)=854
 
 NOTE: 18 records were read from the infile '/home/u1407221/441s24/data/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
       user cpu time       0.00 seconds
       system cpu time     0.00 seconds
       memory              761.40k
       OS Memory           30120.00k
       Timestamp           03/25/2024 06:33:18 PM
       Step Count                        40  Switch Count  2
       Page Faults                       0
       Page Reclaims                     158
       Page Swaps                        0
       Voluntary Context Switches        21
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           264
       
 
 87         proc means data=memory mean;
 88              class treatmnt;
 89              var week2 -- week16;
 90         run;
 
 NOTE: There were 18 observations read from the data set WORK.MEMORY.
 NOTE: PROCEDURE MEANS used (Total process time):
       real time           0.04 seconds
       user cpu time       0.04 seconds
       system cpu time     0.01 seconds
       memory              9553.96k
       OS Memory           37052.00k
       Timestamp           03/25/2024 06:33:18 PM
       Step Count                        41  Switch Count  1
       Page Faults                       0
       Page Reclaims                     1843
       Page Swaps                        0
       Voluntary Context Switches        12
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           8
       
 
 91         
 92         proc glm data=memory;
 93              class treatmnt;
 94              model week2 -- week16 = treatmnt;
 95              repeated time profile /  short summary mean;
 96         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.
 97         
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.93 seconds
       user cpu time       0.52 seconds
       system cpu time     0.04 seconds
       memory              25080.75k
       OS Memory           53948.00k
       Timestamp           03/25/2024 06:33:19 PM
       Step Count                        42  Switch Count  2
       Page Faults                       0
       Page Reclaims                     8048
       Page Swaps                        0
       Voluntary Context Switches        1818
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           2912
       
 
 98         proc glm data=memory noprint;
 99              title3 'Replicate test for main effect of treatment: F=8.08, p=0.0118';
 100             class treatmnt;
 101             model week2 -- week16 = treatmnt;
 102             manova H = treatmnt
 103                    M = week2+week4+week8+week12+week16 / short;
 104             /* M is a matrix of coefficients for transforming the Response vars */
 105        run;
 
 106        
 
 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              1566.40k
       OS Memory           52408.00k
       Timestamp           03/25/2024 06:33:19 PM
       Step Count                        43  Switch Count  3
       Page Faults                       0
       Page Reclaims                     218
       Page Swaps                        0
       Voluntary Context Switches        22
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           288
       
 
 107        proc glm data=memory noprint;
 108             title3 'Replicate tests for main effect of time: Lambda=0.84009249';
 109             title4 'And time by treatment interaction: Lambda=0.44106117';
 110             class treatmnt;
 111             model week2 -- week16 = treatmnt;
 112             manova H = intercept treatmnt
 113                    M = week2-week4, week4-week8, week8-week12, week12-week16
 114                       / short;
 115        run;
 
 116        
 117        /* But the real point is that the treatment only affects recent memories,
 118        not older ones. A basic MANOVA is really more to the point. Follow up
 119        with Bonferroni-corrected univariate tests. 0.05/5 = 0.01. */
 120        
 
 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              1567.12k
       OS Memory           52408.00k
       Timestamp           03/25/2024 06:33:19 PM
       Step Count                        44  Switch Count  4
       Page Faults                       0
       Page Reclaims                     218
       Page Swaps                        0
       Voluntary Context Switches        32
       Involuntary Context Switches      0
       Block Input Operations            0
       Block Output Operations           288
       
 
 121        proc glm data=memory;
 122             title3 'MANOVA, no repeated measures';
 123             class treatmnt;
 124             model week2 -- week16 = treatmnt;
 125             manova h = treatmnt;
 126        run;
 
 127        quit;
 
 NOTE: PROCEDURE GLM used (Total process time):
       real time           0.77 seconds
       user cpu time       0.43 seconds
       system cpu time     0.03 seconds
       memory              4940.21k
       OS Memory           54204.00k
       Timestamp           03/25/2024 06:33:20 PM
       Step Count                        45  Switch Count  3
       Page Faults                       0
       Page Reclaims                     2879
       Page Swaps                        0
       Voluntary Context Switches        1813
       Involuntary Context Switches      4
       Block Input Operations            0
       Block Output Operations           2560
       
 
 128        
 129        OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 141