Se utilizarán los datos de contenido de proteína (expresados en porcentaje), obtenidos en el experimento sobre amaranto realizado por Alfaro Villatoro (1985).
Supongamos que al realizar los análisis de proteína, la muestra correspondiente al corte a 25 días en el 4º bloque fue extraviada.
Bloque
|
Edad al corte (días)
| ||
25
|
40
|
60
| |
1
|
30.6
|
22.4
|
13.6
|
2
|
28.5
|
24.1
|
13.1
|
3
|
29.1
|
26.3
|
15.8
|
4
|
24.3
|
20.2
| |
5
|
30.6
|
18.6
|
13.8
|
6
|
30.6
|
20.9
|
12
|
7
|
28.3
|
21.1
|
12.7
|
8
|
28.7
|
23.7
|
13.7
|
Procedimiento en SAS:
Title1 'Diseño en Bloques Completos al Azar Desbalanceado';
Title2 'Contenidos de proteina para 3 edades de corte';
DATA proteina;
INPUT Edadcorte $ Bloque Proteina; /* El simbolo $ despues de una
variable indica que es de caracter alphanumerico*/
datalines;
25d 1 30.6
25d 2 28.5
25d 3 29.1
25d 5 30.6
25d 6 30.6
25d 7 28.3
25d 8 28.7
40d 1 22.4
40d 2 24.1
40d 3 26.3
40d 4 24.3
40d 5 18.6
40d 6 20.9
40d 7 21.1
40d 8 23.7
60d 1 13.6
60d 2 13.1
60d 3 15.8
60d 4 20.2
60d 5 13.8
60d 6 12.0
60d 7 12.7
60d 8 13.7
;
proc glm data=proteina;
class Edadcorte Bloque;/*Variables clasificatorias*/
model Proteina= Edadcorte Bloque / ss3;/*modelo / suma de cuadrados TIPO3*/
means Edadcorte / Tukey lines;/* Tukey=comparacion entre los tratamientos,
lines=evita comparaciones 1 a 1 para Tukey cuando hay datos faltantes */
run;
quit;
No hay comentarios:
Publicar un comentario