En un experimento de campo, el investigador puede no tener tiempo para cosechar (totalmente) cada unidad experimental. De esta manera, podrá seleccionar al azar varios cuadros por parcela y cosechar el grano en cada cuadro seleccionado. Describiríamos estas observaciones como “muestras dentro de unidades experimentales”.
Tratamientos
|
Submuestra
|
Repeticiones
| |||
I
|
II
|
III
|
IV
| ||
A
|
1
|
5.6
|
5
|
5.5
|
5.3
|
2
|
5.7
|
5.1
|
5.4
|
5.5
| |
B
|
1
|
6.7
|
4.7
|
5.7
|
6.2
|
2
|
8.7
|
3.7
|
6.5
|
5.8
| |
C
|
1
|
7.6
|
7.4
|
7.5
|
5.7
|
2
|
7.8
|
7.2
|
7.6
|
6.7
| |
t = 3 ( A, B, C)
r = 4 (I, II, III, IV)
m= 2 (submuestras)
Procedimiento en SAS:
Title1 'Diseño Completamente al Azar con Submuestras';
Title2 'Tratamientos, Repeticiones, Submuestra, Produccion';
data produccion;
input Tratamientos $ Repeticiones Submuestra Produccion;/*El simbolo $ despues de
una variable indica que es de caracter alphanumerico*/
datalines;
A 1 1 5.6
A 1 2 5.7
B 1 1 6.7
B 1 2 8.7
C 1 1 7.6
C 1 2 7.8
A 2 1 5.0
A 2 2 5.1
B 2 1 4.7
B 2 2 3.7
C 2 1 7.4
C 2 2 7.2
A 3 1 5.5
A 3 2 5.4
B 3 1 5.7
B 3 2 6.5
C 3 1 7.5
C 3 2 7.6
A 4 1 5.3
A 4 2 5.5
B 4 1 6.2
B 4 2 5.8
C 4 1 5.7
C 4 2 6.7
;
Proc sort data=produccion;
by tratamientos;
/*el procedimiento ordena los datos de acuerdo con lo que estipule el parametro
by, para este caso ordena los datos de A hacia C en la varable tratamiento*/
run;
proc glm data=produccion;
class Tratamientos Repeticiones; /*Variables clasificatorias*/
model Produccion= Tratamientos Tratamientos*Repeticiones / ss3;/*modelo / suma de cuadrados TIPO3*/
test h=Tratamientos e=Tratamientos*Repeticiones; /*con la instruccion test se realiza la prueba de
hipotesis 'h' a los tratamientos con el error 'e=Tratamientos*Repeticiones'*/
means Tratamientos / lsd lines e=Tratamientos*Repeticiones; /*lsd=comparacion entre los tratamientos,
lines=evita comparaciones 1 a 1 para lsd cuando hay datos faltantes; las cuales se presentan por
omision, e=Tratamientos*Repeticiones en este caso le indica a lsd usar el CME deel error experimental*/
run;
quit;
No hay comentarios:
Publicar un comentario