1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
70
71 /* mvtubes.sas */
72
73 title 'Multivariate tests on Little Fungus Tube data';
74
75 data mould;
76 infile '/home/brunner0/441s20/LittleTubes2.data.txt' firstobs=2;
77 input line tube mcg length10 weight;
78 /* Make indicator dummy variables using arrays */
79 /* v1-v6 will have these values for every case. */
80 v1=198; v2=205; v3=213; v4=221; v5=223; v6=225;
81 array dummy{6} mcg198 mcg205 mcg213 mcg221 mcg223 mcg225;
82 array value{6} v1-v6;
83 do j = 1 to 6;
84 if mcg = . then dummy{j} = .;
85 else if mcg = value{j} then dummy{j}=1;
86 else dummy{j}=0;
87 end;
88 /* Effect coding */
89 array fungus{5} fungus1-fungus5;
90 do i = 1 to 5;
91 if mcg=225 then fungus{i} = -1;
92 else fungus{i}=dummy{i};
93 end;
94
NOTE: The infile '/home/brunner0/441s20/LittleTubes2.data.txt' is:
Filename=/home/brunner0/441s20/LittleTubes2.data.txt,
Owner Name=brunner0,Group Name=oda,
Access Permission=-rw-r--r--,
Last Modified=19Feb2020:11:30:30,
File Size (bytes)=1775
NOTE: 24 records were read from the infile '/home/brunner0/441s20/LittleTubes2.data.txt'.
The minimum record length was 69.
The maximum record length was 69.
NOTE: The data set WORK.MOULD has 24 observations and 24 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 808.96k
OS Memory 31912.00k
Timestamp 02/19/2020 04:34:38 PM
Step Count 34 Switch Count 2
Page Faults 0
Page Reclaims 169
Page Swaps 0
Voluntary Context Switches 19
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 264
95 proc freq;
96 title2 'Check dummy variable creation';
97 tables (mcg198--mcg225 fungus1-fungus5)*mcg / norow nocol nopercent missing;
98
NOTE: There were 24 observations read from the data set WORK.MOULD.
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.26 seconds
user cpu time 0.26 seconds
system cpu time 0.00 seconds
memory 3793.46k
OS Memory 33712.00k
Timestamp 02/19/2020 04:34:39 PM
Step Count 35 Switch Count 5
Page Faults 0
Page Reclaims 468
Page Swaps 0
Voluntary Context Switches 25
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 608
99 proc glm;
100 title2 'Manova with proc glm';
101 class mcg;
102 model length10 weight = mcg;
103 /* Test equality of mcgs excluding 198: a COLLECTION of contrasts */
104 contrast 'AllBut198' mcg 0 1 -1 0 0 0,
105 mcg 0 0 1 -1 0 0,
106 mcg 0 0 0 1 -1 0,
107 mcg 0 0 0 0 1 -1;
108 manova H = _all_;
109
NOTE: PROCEDURE GLM used (Total process time):
real time 0.56 seconds
user cpu time 0.28 seconds
system cpu time 0.03 seconds
memory 16517.65k
OS Memory 47168.00k
Timestamp 02/19/2020 04:34:39 PM
Step Count 36 Switch Count 3
Page Faults 0
Page Reclaims 3974
Page Swaps 0
Voluntary Context Switches 1132
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 1304
110 proc reg plots=none;
111 title2 'Manova with proc reg and effect coding;';
112 model length10 weight = fungus1-fungus5;
113 overall: mtest fungus1=fungus2=fungus3=fungus4=fungus5=0;
114
NOTE: PROCEDURE REG used (Total process time):
real time 0.09 seconds
user cpu time 0.10 seconds
system cpu time 0.00 seconds
memory 2399.06k
OS Memory 47552.00k
Timestamp 02/19/2020 04:34:39 PM
Step Count 37 Switch Count 2
Page Faults 0
Page Reclaims 343
Page Swaps 0
Voluntary Context Switches 22
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 96
115 proc reg plots=none;
116 title2 'Manova with proc reg and cell means coding';
117 model length10 weight = mcg198 mcg205 mcg213 mcg221 mcg223 mcg225 / noint;
118 Overall: test mcg198=mcg205=mcg213=mcg221=mcg223=mcg225;
119 mOverall: mtest mcg198=mcg205=mcg213=mcg221=mcg223=mcg225;
120 AllBut198: test mcg205=mcg213=mcg221=mcg223=mcg225;
121 mAllBut198: mtest mcg205=mcg213=mcg221=mcg223=mcg225;
122
123
124
125
126
127 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
138