enero 18, 2013

DCA con replicas desiguales (usando SAS)


Análisis de varianza para Diseño Completamente al Azar (DCA) con replicas desiguales (Tabla 2.3 en Procedimientos Estadísticos para Investigaciones Agrícolas. 1984. Gómez, K.A.; Gómez, A.A. pp. 17-20)

Tratamientos

Rendimiento en grano kg/ha 


Total
Media
Tipo
Proporcióna
Días después
de siembra DDS
hA/hB
2.0/0.25
21
3187
4610
3562
3217
14576
3644
hA/hC
3.0/1.00
28
3390
2875
2775
9040
3013
hA/hB
2.0/0.25
14
2797
3001
2505
3490
11793
2948
hA/hD
2.0/0.50
14
2832
3103
3448
2255
11638
2910
hA/hE
3.0/1.50
21
2233
2743
2727
7703
2568
hF
1.5
14
2952
2272
2470
7694
2565
hA/hB
2.0/0.25
28
2858
2895
2458
1723
9934
2484
hA/hG
3.0/1.00
28
2308
2335
1975
6618
2206
hA/hD
2.0/0.50
28
2013
1788
2248
2115
8164
2041
Manual
15 y 35
3202
3060
2240
2690
11192
2798
Control
1192
1652
1075
1030
4949
1237
Grand total (G)
103301
Grand media







2583
a = ingrediente activo
DDA = días después de siembra


/* Inician códigos en SAS
   Análisis de varianza para DCA con replicas desiguales*/

data herb;

input tipo $ proporcion $ DDS $ rendimiento;
cards;

hA/hB21       2.0/0.25      21     3187
hA/hB21       2.0/0.25      21     4610
hA/hB21       2.0/0.25      21     3562
hA/hB21       2.0/0.25      21     3217
hA/hC28       3.0/1.00      28     3390
hA/hC28       3.0/1.00      28     2875
hA/hC28       3.0/1.00      28     2775
hA/hC28       3.0/1.00      28     .
hA/hB14       2.0/0.25      14     2797
hA/hB14       2.0/0.25      14     3001
hA/hB14       2.0/0.25      14     2505
hA/hB14       2.0/0.25      14     3490
hA/hD14       2.0/0.50      14     2832
hA/hD14       2.0/0.50      14     3103
hA/hD14       2.0/0.50      14     3448
hA/hD14       2.0/0.50      14     2255
hA/hE21       3.0/1.50      21     2233
hA/hE21       3.0/1.50      21     2743
hA/hE21       3.0/1.50      21     2727
hA/hE21       3.0/1.50      21     .
hF            1.5           14     2952
hF            1.5           14     2272
hF            1.5           14     2470
hF            1.5           14     .
hA/hB28       2.0/0.25      28     2858
hA/hB28       2.0/0.25      28     2895
hA/hB28       2.0/0.25      28     2458
hA/hB28       2.0/0.25      28     1723
hA/hG28       3.0/1.00      28     2308
hA/hG28       3.0/1.00      28     2335
hA/hG28       3.0/1.00      28     1975
hA/hG28       3.0/1.00      28     .
hA/hD28       2.0/0.50      28     2013
hA/hD28       2.0/0.50      28     1788
hA/hD28       2.0/0.50      28     2248
hA/hD28       2.0/0.50      28     2115
Manual        15-y-35       .      3202
Manual        15-y-36       .      3060
Manual        15-y-37       .      2240
Manual        15-y-38       .      2690
Control       -             .      1192
Control       -             .      1652
Control       -             .      1075
Control       -             .      1030
;

ods html;
proc print;
run;

*Medias con errores estándar;
proc means sum mean stderr;
   class tipo;
run;

*Analisis de varianza;
proc anova;
   class tipo;
   model rendimiento = tipo;
run;

*Diagramas de cajas o boxplots;
proc boxplot;
   plot rendimiento * tipo;
run;quit;
ods html close;

SAS Output:







No hay comentarios.: