1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;7071 /* mvtubes.sas */7273 title 'Multivariate tests on Little Fungus Tube data';7475 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;94NOTE: 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)=1775NOTE: 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 secondsuser cpu time 0.00 secondssystem cpu time 0.00 secondsmemory 808.96kOS Memory 31912.00kTimestamp 02/19/2020 04:34:38 PMStep Count 34 Switch Count 2Page Faults 0Page Reclaims 169Page Swaps 0Voluntary Context Switches 19Involuntary Context Switches 0Block Input Operations 0Block Output Operations 26495 proc freq;96 title2 'Check dummy variable creation';97 tables (mcg198--mcg225 fungus1-fungus5)*mcg / norow nocol nopercent missing;98NOTE: There were 24 observations read from the data set WORK.MOULD.NOTE: PROCEDURE FREQ used (Total process time):real time 0.26 secondsuser cpu time 0.26 secondssystem cpu time 0.00 secondsmemory 3793.46kOS Memory 33712.00kTimestamp 02/19/2020 04:34:39 PMStep Count 35 Switch Count 5Page Faults 0Page Reclaims 468Page Swaps 0Voluntary Context Switches 25Involuntary Context Switches 0Block Input Operations 0Block Output Operations 60899 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_;109NOTE: PROCEDURE GLM used (Total process time):real time 0.56 secondsuser cpu time 0.28 secondssystem cpu time 0.03 secondsmemory 16517.65kOS Memory 47168.00kTimestamp 02/19/2020 04:34:39 PMStep Count 36 Switch Count 3Page Faults 0Page Reclaims 3974Page Swaps 0Voluntary Context Switches 1132Involuntary Context Switches 0Block Input Operations 0Block Output Operations 1304110 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;114NOTE: PROCEDURE REG used (Total process time):real time 0.09 secondsuser cpu time 0.10 secondssystem cpu time 0.00 secondsmemory 2399.06kOS Memory 47552.00kTimestamp 02/19/2020 04:34:39 PMStep Count 37 Switch Count 2Page Faults 0Page Reclaims 343Page Swaps 0Voluntary Context Switches 22Involuntary Context Switches 0Block Input Operations 0Block Output Operations 96115 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;122123124125126127 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;138