febrero 14, 2013

Análisis de varianza para un diseño balanceado lattice de 4 x 4 (en SAS)

Tabla 2.10 Número de cultivadores por metro cuadrado de 16 tratamientos de fertilizante examinados en un diseño balanceado lattice de 4 X 4. (Fuente: paginas 39–52 en Statistical Procedures for Agricultural Research. 1984. Gómez, K.A.; Gómez, A.A. Segunda Edición; IRRI—John Wiley & Sons.)

Bloque Número
Cultivador, no./m2
Bloque Total (B)
Bloque Número
Cultivador, no./m2
Bloque Total (B)
Rep. I
Rep. II
(1)
(2)
(3)
(4)
(1)
(5)
(9)
(13)
1
147
152
167
150
616
5
140
165
182
152
639
(5)
(6)
(7)
(8)
(10)
(2)
(14)
(6)
2
127
155
162
172
616
6
97
155
192
142
586
(9)
(10)
(11)
(12)
(7)
(15)
(3)
(11)
3
147
100
192
177
616
7
155
182
192
192
721
(13)
(14)
(15)
(16)
(16)
(8)
(12)
(4)
4
155
195
192
205
747
8
182
207
232
162
783
Rep. total R1
2595
Rep. total R2
2729
Rep. III
Rep. IV
(1)
(6)
(11)
(16)
(1)
(14)
(7)
(12)
9
155
162
177
152
646
13
220
202
175
205
802
(5)
(2)
(15)
(12)
(13)
(2)
(11)
(8)
10
182
130
177
165
654
14
205
152
180
187
724
(9)
(14)
(3)
(8)
(5)
(10)
(3)
(16)
11
137
185
152
152
626
15
165
150
200
160
675
(13)
(10)
(7)
(4)
(9)
(6)
(15)
(4)
12
185
122
182
192
681
16
155
177
185
172
689
Rep. total R3
2607
Rep. total R4
2890
Rep. V
(1)
(10)
(15)
(8)
17
147
112
177
147
583
(9)
(2)
(7)
(16)
18
180
205
190
167
742
(13)
(6)
(3)
(12)
19
172
212
197
192
773
(5)
(14)
(11)
(4)
20
177
220
205
225
827
Rep. total R5
2925
Nota: Los valores en paréntesis corresponden a los números de tratamientos.

* Tabla 2.10 Número de cultivadores por metro cuadrado de 16 tratamientos
  de fertilizante examinados en un diseño balanceado lattice de 4 X 4;

data latt;
input bloque rep trt yield;
cards;
1      1      1      147
1      1      2      152
1      1      3      167
1      1      4      150
2      1      5      127
2      1      6      155
2      1      7      162
2      1      8      172
3      1      9      147
3      1      10     100
3      1      11     192
3      1      12     177
4      1      13     155
4      1      14     195
4      1      15     192
4      1      16     205
5      2      1      140
5      2      5      165
5      2      9      182
5      2      13     152
6      2      10     97
6      2      2      155
6      2      14     192
6      2      6      142
7      2      7      155
7      2      15     182
7      2      3      192
7      2      11     192
8      2      16     182
8      2      8      207
8      2      12     232
8      2      4      162
9      3      1      155
9      3      6      162
9      3      11     177
9      3      16     152
10     3      5      182
10     3      2      130
10     3      15     177
10     3      12     165
11     3      9      137
11     3      14     185
11     3      3      152
11     3      8      152
12     3      13     185
12     3      10     122
12     3      7      182
12     3      4      192
13     4      1      220
13     4      14     202
13     4      7      175
13     4      12     205
14     4      13     205
14     4      2      152
14     4      11     180
14     4      8      187
15     4      5      165
15     4      10     150
15     4      3      200
15     4      16     160
16     4      9      155
16     4      6      177
16     4      15     185
16     4      4      172
17     5      1      147
17     5      10     112
17     5      15     177
17     5      8      147
18     5      9      180
18     5      2      205
18     5      7      190
18     5      16     167
19     5      13     172
19     5      6      212
19     5      3      197
19     5      12     192
20     5      5      177
20     5      14     220
20     5      11     205
20     5      4      225
;

ods html;
proc means sum maxdec = 2; *summation;
   class bloque; var yield;
proc means sum maxdec = 2;
   class trt; var yield;
proc means sum maxdec = 2;
   var yield;
run;

proc mixed method = type3;
   class rep bloque trt;
   model yield = trt;
   random rep bloque(rep);
   lsmeans trt/pdiff;
run; quit;
ods html close;

SAS Output:

Class Level Information
ClassLevelsValues
rep51 2 3 4 5
bloque201 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
trt161 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16


Dimensions
Covariance Parameters3
Columns in X17
Columns in Z25
Subjects1
Max Obs Per Subject80


Number of Observations
Number of Observations Read80
Number of Observations Used80
Number of Observations Not Used0


Type 3 Analysis of Variance
SourceDFSum of SquaresMean SquareExpected Mean SquareError TermError DFF ValuePr > F
trt15183961226.412500Var(Residual) + Q(trt)MS(Residual)453.800.0003
rep45946.0500001486.512500Var(Residual) + 4 Var(bloque(rep)) + 16 Var(rep)1.25 MS(bloque(rep)) - 0.25 MS(Residual)12.5251.710.2091
bloque(rep)1511382758.789167Var(Residual) + 3.2 Var(bloque(rep))MS(Residual)452.350.0138
Residual4514533322.962500Var(Residual)....


Covariance Parameter Estimates
Cov ParmEstimate
rep38.6729
bloque(rep)136.20
Residual322.96


Fit Statistics
-2 Res Log Likelihood596.1
AIC (smaller is better)602.1
AICC (smaller is better)602.5
BIC (smaller is better)600.9


Type 3 Tests of Fixed Effects
EffectNum DFDen DFF ValuePr > F
trt15454.33<.0001

 Resultado de F-calculada mostrado en página 51.




No hay comentarios.: