Analysis of variance

Variate: DI1

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 277.8 138.9 1.00

REP.*Units* stratum
TRT 5 590.3 118.1 0.85 0.545
Residual 10 1388.9 138.9

Total 17 2256.9

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 1 *units* 2 19.4 s.e. 8.8

Tables of means

Variate: DI1

Grand mean 6.9

TRT 1 2 3 4 5 6
8.3 8.3 0.0 8.3 16.7 0.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 9.62

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 21.44

66 IF _ibalance.eq.0 .OR. _ibalance.eq.1

67 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
68 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
69 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
70 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
71 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 8.333 -16.68 33.35
2 8.333 -16.68 33.35
3 0.000 -25.02 25.02
4 8.333 -16.68 33.35
5 16.667 -8.35 41.68
6 0.000 -25.02 25.02

72 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
73 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 0.000 a
6 0.000 a
1 8.333 a
4 8.333 a
2 8.333 a
5 16.667 a

74 ENDIF
75 SET [IN=*]
81 "One-way design in randomized blocks"
82 DELETE [REDEFINE=yes] _ibalance
83 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
84 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI2; SAVE=_a2save
Analysis of variance

Variate: DI2

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 625.0 312.5 1.00

REP.*Units* stratum
TRT 5 1562.5 312.5 1.00 0.465
Residual 10 3125.0 312.5

Total 17 5312.5

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DI2

Grand mean 29.2

TRT 1 2 3 4 5 6
33.3 41.7 16.7 33.3 33.3 16.7

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 14.43

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 32.16

85 IF _ibalance.eq.0 .OR. _ibalance.eq.1

86 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
87 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
88 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
89 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
90 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 33.33 -4.193 70.86
2 41.67 4.140 79.19
3 16.67 -20.860 54.19
 4 33.33 -4.193 70.86
5 33.33 -4.193 70.86
6 16.67 -20.860 54.19

91 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
92 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 16.67 a
6 16.67 a
1 33.33 a
4 33.33 a
5 33.33 a
2 41.67 a

93 ENDIF
94 SET [IN=*]
100 "One-way design in randomized blocks"
101 DELETE [REDEFINE=yes] _ibalance
102 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
103 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI3; SAVE=_a2save
Analysis of variance

Variate: DI3

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 1458.3 729.2 1.52

REP.*Units* stratum
TRT 5 312.5 62.5 0.13 0.982
Residual 10 4791.7 479.2

Total 17 6562.5

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 3 *units* 3 33.3 s.e. 16.3

Tables of means

Variate: DI3

Grand mean 54.2

TRT 1 2 3 4 5 6
58.3 50.0 50.0 58.3 50.0 58.3

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 17.87

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 39.82

104 IF _ibalance.eq.0 .OR. _ibalance.eq.1

105 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
106 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
107 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
108 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
109 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 58.33 11.865 104.8
2 50.00 3.532 96.5
3 50.00 3.532 96.5
4 58.33 11.865 104.8
5 50.00 3.532 96.5
6 58.33 11.865 104.8

110 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
111 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
2 50.00 a
3 50.00 a
5 50.00 a
1 58.33 a
4 58.33 a
6 58.33 a

112 ENDIF
113 SET [IN=*]
119 "One-way design in randomized blocks"
120 DELETE [REDEFINE=yes] _ibalance
121 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
122 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI4; SAVE=_a2save
Analysis of variance

Variate: DI4

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 208.3 104.2 0.20

REP.*Units* stratum
TRT 5 833.3 166.7 0.32 0.890
Residual 10 5208.3 520.8

Total 17 6250.0

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DI4

Grand mean 75.0

TRT 1 2 3 4 5 6
83.3 66.7 66.7 83.3 75.0 75.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 18.63

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 41.52

123 IF _ibalance.eq.0 .OR. _ibalance.eq.1

124 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
125 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
126 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
127 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
128 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 83.33 34.89 131.8
2 66.67 18.22 115.1
3 66.67 18.22 115.1
 4 83.33 34.89 131.8
5 75.00 26.55 123.4
6 75.00 26.55 123.4

129 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
130 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
2 66.67 a
3 66.67 a
5 75.00 a
6 75.00 a
1 83.33 a
4 83.33 a

131 ENDIF
132 SET [IN=*]
138 "One-way design in randomized blocks"
139 DELETE [REDEFINE=yes] _ibalance
140 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
141 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI5; SAVE=_a2save
Analysis of variance

Variate: DI5

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 486.1 243.1 0.66

REP.*Units* stratum
TRT 5 590.3 118.1 0.32 0.889
Residual 10 3680.6 368.1

Total 17 4756.9

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DI5

Grand mean 81.9

TRT 1 2 3 4 5 6
83.3 83.3 75.0 83.3 91.7 75.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 15.66

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 34.90

142 IF _ibalance.eq.0 .OR. _ibalance.eq.1

143 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
144 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
145 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
146 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
147 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 83.33 42.61 124.1
2 83.33 42.61 124.1
3 75.00 34.27 115.7
 4 83.33 42.61 124.1
5 91.67 50.94 132.4
6 75.00 34.27 115.7

148 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
149 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 75.00 a
6 75.00 a
1 83.33 a
2 83.33 a
4 83.33 a
5 91.67 a

150 ENDIF
151 SET [IN=*]
157 "One-way design in randomized blocks"
158 DELETE [REDEFINE=yes] _ibalance
159 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
160 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI6; SAVE=_a2save
Analysis of variance

Variate: DI6

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 69.4 34.7 0.22

REP.*Units* stratum
TRT 5 590.3 118.1 0.74 0.611
Residual 10 1597.2 159.7

Total 17 2256.9

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 2 *units* 6 -19.4 s.e. 9.4

Tables of means

Variate: DI6

Grand mean 93.1

TRT 1 2 3 4 5 6
91.7 83.3 100.0 100.0 91.7 91.7

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 10.32

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 22.99

161 IF _ibalance.eq.0 .OR. _ibalance.eq.1

162 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
163 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
164 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
165 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
166 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 91.67 64.84 118.5
2 83.33 56.50 110.2
3 100.00 73.17 126.8
4 100.00 73.17 126.8
5 91.67 64.84 118.5
6 91.67 64.84 118.5

167 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
168 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
2 83.33 a
1 91.67 a
5 91.67 a
6 91.67 a
3 100.00 a
4 100.00 a

169 ENDIF
170 SET [IN=*]
176 "One-way design in randomized blocks"
177 DELETE [REDEFINE=yes] _ibalance
178 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
179 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DI7; SAVE=_a2save
Analysis of variance

Variate: DI7

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 0. 0.

REP.*Units* stratum
TRT 5 0. 0.
Residual 10 0. 0.

Total 17 0.

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DI7

Grand mean 100.0

TRT 1 2 3 4 5 6
100.0 100.0 100.0 100.0 100.0 100.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. *
s.e.d. 0.00

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. *
l.s.d. *

180 IF _ibalance.eq.0 .OR. _ibalance.eq.1

181 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
182 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
183 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
184 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
185 DF=_rdf
Fault 1, code UF 1, statement 12 in procedure CONFIDENCE

The estimate of VARIANCE is zero or negative.

186 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
187 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Warning 2, code UF 2, statement 159 in procedure AMCOMPARISON

Multiple comparisons cannot be calculated for TRT as its standard errors are zero.

188 ENDIF
Fault 3, code PC 9, statement 1 on line 188

Command: ENDIF
Invalid control structure: unmatched ENDIF statement.
ENDIF does not have preceding IF

189 SET [IN=*]

195 "One-way design in randomized blocks"
196 DELETE [REDEFINE=yes] _ibalance
197 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
198 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS1; SAVE=_a2save
Analysis of variance

Variate: DS1

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 44.44 22.22 1.00

REP.*Units* stratum
TRT 5 94.44 18.89 0.85 0.545
Residual 10 222.22 22.22

Total 17 361.11

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 1 *units* 2 7.8 s.e. 3.5

Tables of means

Variate: DS1

Grand mean 2.8

TRT 1 2 3 4 5 6
3.3 3.3 0.0 3.3 6.7 0.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 3.85

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 8.58

199 IF _ibalance.eq.0 .OR. _ibalance.eq.1

200 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
201 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
202 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
203 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
204 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 3.333 -6.674 13.34
2 3.333 -6.674 13.34
3 0.000 -10.007 10.01
4 3.333 -6.674 13.34
5 6.667 -3.340 16.67
6 0.000 -10.007 10.01

205 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
206 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 0.000 a
6 0.000 a
1 3.333 a
4 3.333 a
2 3.333 a
5 6.667 a

207 ENDIF
208 SET [IN=*]
214 "One-way design in randomized blocks"
215 DELETE [REDEFINE=yes] _ibalance
216 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
217 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS2; SAVE=_a2save
Analysis of variance

Variate: DS2

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 100.00 50.00 1.00

REP.*Units* stratum
TRT 5 250.00 50.00 1.00 0.465
Residual 10 500.00 50.00

Total 17 850.00

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DS2

Grand mean 11.7

TRT 1 2 3 4 5 6
13.3 16.7 6.7 13.3 13.3 6.7

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 5.77

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 12.86

218 IF _ibalance.eq.0 .OR. _ibalance.eq.1

219 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
220 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
221 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
222 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
223 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 13.33 -1.677 28.34
2 16.67 1.656 31.68
3 6.67 -8.344 21.68
 4 13.33 -1.677 28.34
5 13.33 -1.677 28.34
6 6.67 -8.344 21.68

224 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
225 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 6.67 a
6 6.67 a
1 13.33 a
4 13.33 a
5 13.33 a
2 16.67 a

226 ENDIF
227 SET [IN=*]
233 "One-way design in randomized blocks"
234 DELETE [REDEFINE=yes] _ibalance
235 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
236 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS3; SAVE=_a2save
Analysis of variance

Variate: DS3

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 486.1 243.1 2.36

REP.*Units* stratum
TRT 5 77.8 15.6 0.15 0.975
Residual 10 1030.6 103.1

Total 17 1594.4

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 3 *units* 3 17.2 s.e. 7.6

Tables of means

Variate: DS3

Grand mean 25.6

TRT 1 2 3 4 5 6
28.3 25.0 23.3 28.3 23.3 25.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 8.29

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 18.47

237 IF _ibalance.eq.0 .OR. _ibalance.eq.1

238 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
239 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
240 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
241 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
242 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 28.33 6.783 49.88
2 25.00 3.450 46.55
3 23.33 1.783 44.88
4 28.33 6.783 49.88
5 23.33 1.783 44.88
6 25.00 3.450 46.55

243 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
244 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 23.33 a
5 23.33 a
2 25.00 a
6 25.00 a
1 28.33 a
4 28.33 a

245 ENDIF
246 SET [IN=*]
252 "One-way design in randomized blocks"
253 DELETE [REDEFINE=yes] _ibalance
254 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
255 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS4; SAVE=_a2save
Analysis of variance

Variate: DS4

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 36.1 18.1 0.07

REP.*Units* stratum
TRT 5 456.9 91.4 0.36 0.865
Residual 10 2547.2 254.7

Total 17 3040.3

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DS4

Grand mean 41.4

TRT 1 2 3 4 5 6
46.7 40.0 31.7 45.0 40.0 45.0

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 13.03

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 29.04

256 IF _ibalance.eq.0 .OR. _ibalance.eq.1

257 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
258 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
259 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
260 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
261 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 46.67 12.786 80.55
2 40.00 6.120 73.88
3 31.67 -2.214 65.55
 4 45.00 11.120 78.88
5 40.00 6.120 73.88
6 45.00 11.120 78.88

262 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
263 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 31.67 a
2 40.00 a
5 40.00 a
4 45.00 a
6 45.00 a
1 46.67 a

264 ENDIF
265 SET [IN=*]
271 "One-way design in randomized blocks"
272 DELETE [REDEFINE=yes] _ibalance
273 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
274 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS5; SAVE=_a2save
Analysis of variance

Variate: DS5

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 136.1 68.1 0.23

REP.*Units* stratum
TRT 5 290.3 58.1 0.20 0.957
Residual 10 2947.2 294.7

Total 17 3373.6

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DS5

Grand mean 49.7

TRT 1 2 3 4 5 6
50.0 48.3 43.3 48.3 56.7 51.7

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 14.02

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 31.23

275 IF _ibalance.eq.0 .OR. _ibalance.eq.1

276 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
277 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
278 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
279 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
280 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 50.00 13.56 86.44
2 48.33 11.89 84.78
3 43.33 6.89 79.78
 4 48.33 11.89 84.78
5 56.67 20.22 93.11
6 51.67 15.22 88.11

281 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
282 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
3 43.33 a
2 48.33 a
4 48.33 a
1 50.00 a
6 51.67 a
5 56.67 a

283 ENDIF
284 SET [IN=*]
290 "One-way design in randomized blocks"
291 DELETE [REDEFINE=yes] _ibalance
292 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
293 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS6; SAVE=_a2save
Analysis of variance

Variate: DS6

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 322.3 161.2 0.63

REP.*Units* stratum
TRT 5 185.8 37.2 0.15 0.977
Residual 10 2558.3 255.8

Total 17 3066.5

Information summary

All terms orthogonal, none aliased.

Message: the following units have large residuals.

REP 2 *units* 6 -31.2 s.e. 11.9

Tables of means

Variate: DS6

Grand mean 63.8

TRT 1 2 3 4 5 6
66.7 58.3 63.3 65.0 68.0 61.7

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 13.06

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 29.10

294 IF _ibalance.eq.0 .OR. _ibalance.eq.1

295 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
296 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
297 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
298 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
299 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

 Mean Lower Upper
TRT
1 66.67 32.71 100.62
2 58.33 24.38 92.29
3 63.33 29.38 97.29
4 65.00 31.05 98.95
5 68.00 34.05 101.95
6 61.67 27.71 95.62

300 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
301 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
2 58.33 a
6 61.67 a
3 63.33 a
4 65.00 a
1 66.67 a
5 68.00 a

302 ENDIF
303 SET [IN=*]
309 "One-way design in randomized blocks"
310 DELETE [REDEFINE=yes] _ibalance
311 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
312 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS7; SAVE=_a2save
Analysis of variance

Variate: DS7

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 11.11 5.56 0.06

REP.*Units* stratum
TRT 5 177.78 35.56 0.38 0.852
Residual 10 938.89 93.89

Total 17 1127.78

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DS7

Grand mean 76.1

TRT 1 2 3 4 5 6
76.7 75.0 80.0 71.7 80.0 73.3

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 7.91

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 17.63

313 IF _ibalance.eq.0 .OR. _ibalance.eq.1

314 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
315 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
316 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
317 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
318 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 76.67 56.10 97.24
2 75.00 54.43 95.57
3 80.00 59.43 100.57
 4 71.67 51.10 92.24
5 80.00 59.43 100.57
6 73.33 52.76 93.90

319 AMCOMPARISON [PRINT=critical,description,letter; METHOD=bonferroni;

DIRECTION=ascending;\
320 PROB=0.05; SAVE=_a2save['save']] TRT

Bonferroni test

Experiment-wise error rate = 0.05

TRT

Comparison-wise error rate = 0.0033

Critical values for t

5% 3.827
1% 4.855
0.1% 6.526

Mean
4 71.67 a
6 73.33 a
2 75.00 a
1 76.67 a
3 80.00 a
5 80.00 a

321 ENDIF
322 SET [IN=*]
328 "One-way design in randomized blocks"
329 DELETE [REDEFINE=yes] _ibalance
330 A2WAY [PRINT=aovtable,information,means; TREATMENTS=TRT; BLOCKS=REP; FPROB=yes;
PSE=diff,\
331 lsd; LSDLEVEL=5; PLOT=*; COMBINATIONS=present; EXIT=_ibalance] DS7; SAVE=_a2save
Analysis of variance

Variate: DS7

Source of variation d.f. s.s. m.s. v.r. F pr.

REP stratum 2 11.11 5.56 0.06

REP.*Units* stratum
TRT 5 177.78 35.56 0.38 0.852
Residual 10 938.89 93.89

Total 17 1127.78

Information summary

All terms orthogonal, none aliased.

Tables of means

Variate: DS7

Grand mean 76.1

TRT 1 2 3 4 5 6
76.7 75.0 80.0 71.7 80.0 73.3

Standard errors of differences of means

Table TRT
rep. 3
d.f. 10
s.e.d. 7.91

Least significant differences of means (5% level)

Table TRT
rep. 3
d.f. 10
l.s.d. 17.63

332 IF _ibalance.eq.0 .OR. _ibalance.eq.1

333 DELETE [REDEFINE=yes] _mean, _rep, _var, _rdf
334 AKEEP [SAVE=_a2save['save']] TRT; MEAN=_mean; REP=_rep; VARIANCE=_var;
RTERM=_resid
335 AKEEP [SAVE=_a2save['save']] #_resid; DF=_rdf
336 CONFIDENCE [METHOD=smm; PROB=0.05] MEANS=_mean; REPLICATION=_rep;
VARIANCE=_var;\
337 DF=_rdf

Studentized Maximum Modulus 95.0% confidence intervals

Equal number of observations per mean. (Input as scalar.)

MEAN, LOWER, UPPER are tables.

Mean Lower Upper

TRT
1 76.67 56.10 97.24
2 75.00 54.43 95.57
3 80.00 59.43 100.57
 4 71.67 51.10 92.24
5 80.00 59.43 100.57
6 73.33 52.76 93.90

338 AMCOMPARISON [PRINT=critical,description,letter; METHOD=fplsd;

DIRECTION=ascending;\
339 PROB=0.05; STUDENTIZE=no; SAVE=_a2save['save']] TRT

Fisher's protected least significant difference test

Experiment-wise error rate = 0.05

Comparison-wise error rate = 0.05

TRT

Warning 4, code UF 2, statement 194 in procedure AMCOMPARISON

Fisher's protected LSD is not calculated as variance ratio for TRT is not significant.

340 ENDIF
341 SET [IN=*]