TALLER N.

4
Clasificación de modelos de análisis de varianza

Presentado por:
Diego Alejandro Ayala Jiménez

1. Suponga el modelo Yijk = μ+b +ai j(i) + ∈k(ij), donde k = 1, . . ., t; j = 1, . . ., r;

i = 1, . . ., m. Suponga que el único e
́ fecto aleatorio es el error experimental.
a. Presente el diagrama de estructura para llevar a cabo el análisis de
este experimento.
Yijk=μ + bi+ aj(i) + ek(ij)

b. Lleve a cabo el ANOVA especificando los grados de libertad, la suma

de cuadrados, las esperanzas de los cuadrados medios y los
estadísticos de prueba para las hipótesis de interés.
F.d.v G.l. S.C. i j k E(cm)

b m-1 ((1/jk)∑Yi..2) - (Y…2/mrt) 1 r t σ2e + ((rt)/m-1)∑(b̂)2i

a:b (r-1)m ((1/k)∑∑Yij.2) - 0 1 t σ2e + ((t/(r-

((1/jk)∑Yi..2) 1)m)∑(â)2j

error (t-1)mr (∑∑∑Yijk2) - 1 1 1 σ2e

((1/k)∑∑Yij.2)

total mrt-1 (∑∑∑Yijk2) - (Y…2/mrt)

Hipótesis:
1- Ho= μ1 = μ2 = … = μm; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(b)/CME; Rechazar si Ho si F> Fcrit α, m-1, (t-1)mr
2- Ho= μ1 = μ2 = … = μm; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(a)/CME; Rechazar si Ho si F> Fcrit α, (r-1)m, (t-1)mr

2. Para los ejercicios 5 al 8 del capıtulo 1

a. Presente el diagrama de estructura para llevar a cabo el análisis de
cada experimento.
 b. Lleve a cabo el ANOVA especificando los grados de libertad, la suma
de cuadrados, las esperanzas de los cuadrados medios y probar la
hipótesis de interés.
➢ Se investigó la pérdida de peso en porcentaje del peso inicial de la carne de
res tipo milanesa, después de cinco días de empacada en cuatro tipos de
envolturas: Icopor, Biopak, Cry-0-Pak y Shopak. Para tal finalidad se
realizaron cinco replicaciones de cada una.
Yij=μ + bi+ ej(i)

Donde ej(i) es el error experimental, bi es el efecto de las envolturas y μ es la media

general del experimento
F.d.v G.l. S.C. i J E(cm)

b 3 ((1/5)∑Yi.2) - (Y..2/20) 1 r σ2e + (r/3)∑(b)2i

error 16 (∑∑Yij2) - ((1/5)∑Yi.2) 1 1 σ2 e

total 19 (∑∑Yij2) - (Y..2/20)

Hipótesis:
1- Ho= μ1 = μ2 = μ3 = μ4; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(α)/CME; Rechazar si Ho si F> Fcrit α, 3, 16

➢ Se investigó la morfometria de J. Lamprotaenia (sardina azul) en las islas de

rosario (Bolívar). Entre otras características, midieron la longitud estándar de
sardinas hembras en tres estaciones de muestreo: fondo arenoso, manglar y
pasto marino.
Yij=μ + bi + ej(i)
Donde ej(i) es el error experimental, bi es el efecto de las tres estaciones de
muestreo y μ es la media general del experimento.
 F.d.v G.l. S.C. I J E(cm)

b 2 ((1/n)∑Yi.2) - (Y..2/3n) 1 N σ2e + (n/2)∑(b)2i

error (n-1)3 (∑∑Yij2) - ((1/n)∑Yi.2) 1 1 σ2e

total 3n-1 (∑∑Yij2) - (Y..2/3n)

Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j

F= CM(α)/CME; Rechazar si Ho si F> Fcrit α, 2, 3n-1

➢ Se midió la concentración de calcio en hojas de cinco especies de árboles

para saber si el nitrógeno en tres niveles diferentes afectaba el porcentaje de
peso seco de la hoja (esta es una forma de medir la concentración de calcio).
Dentro de cada árbol, se tomaron aleatoriamente cuatro hojas y se hicieron
cuatro lecturas de cada hoja.
Yijkl=μ + αi + δj(i) + ek(ij) + ηl (ijk)
Donde ηl (ijk) es el error muestral, ek(ij) es el error experimental, δj(i) es el efecto de las
cinco especies de plantas, αi es el efecto de los tres niveles de nitrógeno y μ es la
media general del experimento.

F.d.v G.l. S.C. i J k l E(cm)

a 2 ((1/80)∑Yi...2) - (Y….2/240) 0 5 4 4 σ2em +4 σ2ee + (80/2)∑(a)2i

b 12 ((1/16)∑∑Yij..2) - 1 0 4 4 σ2em +4 σ2ee + ((16/12)∑∑(b)2j

((1/80)∑Yi...2)

e.e. 45 ((1/4)∑∑∑Yijk.2 - 1 1 1 4 σ2em +4 σ2ee

((1/16)∑∑Yij.2

e.m. 180 (∑∑∑∑Yijkl2) - 1 1 1 1 σ2em

((1/4)∑∑∑Yijk.2)

total 239 (∑∑∑∑Yijkl2) - (Y….2/240)

Hipótesis:
1- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
 F= CM(α)/CM(ee); Rechazar si Ho si F> Fcrit α, 2, 45
2- Ho= μ1 = ... = μ5; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(δ)/CM(ee); Rechazar si Ho si F> Fcrit α, 12, 45
3- Ho= σ2em = σ2ee; H1= Por lo menos un par de σ2em = σ2ee
F= CM(e.e.) /CM(e.m.); Rechazar si Ho si F> Fcrit α, 45,180

➢ Un experimentador está estudiando el efecto de diluyentes de semen en la

tasa de concepción en vacas lecheras, usando las variables: Reducidor
(Citrato o Fosfato), Sulfanilamida (Si o No) y la cantidad de Penicilina (para la
cual se seleccionaron tres niveles entre 0 y 500 miligramos).
● Yijkl=μ + αi + Ԏ j + δk + (α Ԏ)ij+ (αδ)ik + (Ԏ δ)jk +(α Ԏ δ)ijk + el(ijk)
Donde el(ijk) es el error experimental, (α Ԏ δ)ijk es la interacción entre el reducidor, la
sulfanilamida y la cantidad de penicilina, ( Ԏ δ)jk es la interacción entre la
sulfanilamida y la cantidad de penicilina, (αδ) ik es la interacción entre el reducidor y
la cantidad de penicilina, (α Ԏ)ij es la interacción entre el reducidor y la sulfanilamida,
δk es el efecto de la cantidad de penicilina, Ԏ j es el efecto de la sulfanilamida, αi es el efecto
del reducidor y μ es la media general del experimento.

F.d.v G.l. S.C. i J k l E(cm)

α 2 ((1/9n)∑Yi...2) - (Y….2/27n) 0 3 3 n σ2ee + (9n/2)∑(α̂)2i

Ԏ 2 ((1/9n)∑Y.j..2) - (Y….2/27n) 3 0 3 n σ2ee + ((9n/2)∑(ˇԎ)2j

δ 2 ((1/9n)∑Y..k.2) - (Y….2/27n) 3 3 0 n σ2ee + ((9n/1)∑(δ̂)2k

e.e. (n-1) * 27 (∑∑∑∑Yijkl2)- ((1/n)∑∑∑Yijk.2 1 1 1 1 σ2ee

total 27n -1 (∑∑∑∑Yijkl2) - (Y….2/27n)

Hipótesis:
1- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(α)/CM(ee); Rechazar si Ho si F> Fcrit α, 2, ((n-1) *27)
2- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(Ԏ)/CM(ee); Rechazar si Ho si F> Fcrit Ԏ, 2, ((n-1) *27)

3- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j

F= CM(δ)/CM(ee); Rechazar si Ho si F> Fcrit δ, 2, ((n-1) *27)
 3. Suponga que (V) variedades de una leguminosa fueron plantadas en
diferentes fechas (F). De cada combinación fecha-variedad, fueron
seleccionadas algunas plantas (P), y ciertas hojas (L) fueron seleccionadas
aleatoriamente de cada planta. Cada hoja fue analizada por un numero de
métodos estándar (M), por diferentes analistas (R).
a. Presente el diagrama de estructura y su modelo asociado para llevar a
cabo el análisis de este experimento.
Yijklmno=μ + Vi + Fj + Pk(ij) + Ll(ijk)+ Mm(ijkl) +Rn(ijklm) + eo(ijklmn)
Donde eo(ijklmn) es el error experimental, Rn(ijklm) es el efecto de los diferentes analistas,
Mm(ijkl) es el efecto de los métodos estándar, Ll(ijk) es el efecto de las hojas (error
muestral), Pk(ij) es el efecto de las plantas (error experimental), F j es el efecto de las
fechas, Vi es el efecto de las variedades y μ es la media general del experimento.

b. Exponga la tabla de análisis de varianza.

F.d.v G.l. S.C. i J k l m n o E(cm)

Vi:μ v-1 ((1/fphers)∑Yi…...2) - (Y…...2/vfphers) 0 F p h e r s σ2e+ersσ2L+hersσ2P+

((fphers)/v-1)∑(V)2i
 Fj:μ f-1 ((1/vphers)∑Y.j.....2) - (Y….…2/vfphers) o 0 p h e r s σ2e+ersσ2L+hersσ2P+
((ophers)/f-1)∑(F)2j

(VF)ij:μ (v-1)(f-1) ((1/phers)∑∑Yij.....2)- 0 0 p h e r s σ2e+ersσ2L+hersσ2P+((phers/

((1/fphers)∑Yi…...2)- (v-1)(f-1)∑∑(VF)2ij
((1/vphers)∑Y.j.....2)+(Y….…2/vfphers)

Pk(ij):μV vf(p-1) ((1/hers)∑∑∑Yijk....2)- 1 1 1 h e r s σ2e+ersσ2L+hersσ2P

F ((1/phers)∑∑Yij.....2

Ll(ijk):μV vfp(h-1) ((1/ers)∑∑∑∑Yijkl...2)- 1 1 1 1 e r s σ2e+ersσ2L

FP ((1/hers)∑∑∑Yijk....2)

Mm(ijkl):μ vfph(e-1) ((1/rs)∑∑∑∑∑Yijklm..2)- 1 1 1 1 0 r s σ2e+rsσ2M

VFPL ((1/ers)∑∑∑∑Yijkl...2)

Rn(ijklm): vfphe(r-1) ((1/s)∑∑∑∑∑∑Yijklmn.2)- 1 1 1 1 1 0 s σ2e+sσ2R

μVFPL ((1/rs)∑∑∑∑∑Yijklm..2)

E0(ijklmn): vfpher(s-1) (∑∑∑∑∑∑∑Yijklmo2)- 1 1 1 1 1 1 1 σ 2E

μVFPR ((1/s)∑∑∑∑∑∑Yijklmn.2)

Total vfphers-1 (∑∑∑∑∑∑∑Yijklmo2) - (Y…...2/vfphers)

Hipótesis:
1- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(V)/CM(P); Rechazar si Ho si F> Fcrit α, v-1, vf(p-1)
2- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(F)/CM(P); Rechazar si Ho si F> Fcrit α, f-1, vf(p-1)
3- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(VF)/CM(P); Rechazar si Ho si F> Fcrit α;(V-1) (f-1); vf(p-1)
4- Ho= σ2P = 0; H1= Por lo menos un σ2P > 0
F= CM(P)/CM(L); Rechazar si Ho si F> Fcrit α; vf(p-1); vfp(h-1);
5- Ho= σ2L = 0; H1= Por lo menos un σ2L > 0
F= CM(L)/CM(E); Rechazar si Ho si F> Fcrit α; vfp(h-1); vfpher(s-1)
6- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(M)/CM(E); Rechazar si Ho si F> Fcrit α; vfph(e-1); vfpher(s-1)
7- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(R)/CM(E); Rechazar si Ho si F> Fcrit α; vfphe(r-1); vfpher(s-1)

4. Un experimento de caña de azúcar fue diseñado de tal modo que se

ensayaron los nutrientes mayores N, P y K con las siguientes dosis 0 y 200; 0
y 100; 0 y 80 Kg/Ha, respectivamente (Martínez 1988). El área experimental
era bastante uniforme y lo suficientemente grande como para que el material
experimental se ensaya en tres oportunidades. La respuesta observada fue el
rendimiento (Kg/Área de 44m2):

a. Encuentre el diagrama de estructura apropiado.

 b. Obtenga el modelo apropiado señalando todos los supuestos
necesarios que permitan realizar pruebas de hipótesis de efectos
principales e interacciones.
Yijlm=μ + Ni + Pj + Kl + (NP)ij+ (NK)ik + (PK)jk +(NPK)ijl + em(ijl)
Donde em(ijl) es el error experimental, (NPK)ijl es la interacción entre dosis de
nitrógeno, fósforo y potasio, (PK)jk es la interacción entre dosis de fósforo y potasio,
(NK)ik es la interacción entre dosis de nitrógeno y potasio, (NP) ij es la interacción
entre dosis de nitrógeno y fósforo, Kl es el efecto del potasio, Pj es el efecto del
fósforo, Ni es el efecto del nitrógeno y μ es la media general del experimento.
Supuestos:
● Yijlm es una variable aleatoria independiente que se distribuye N (μ ijlm,
σ2)
● eijlm ~ N (0, σ2)
● La sumatoria de los efectos de cada factor y de las interacciones es 0.
c. Construya la tabla ANOVA completa.
F.d.v G.l. S.C. i j l m E(cm)

N 1 ((1/12)∑Yi...2) - (Y….2/24) 0 2 2 3 σ2ee + (12/1)∑(N^)2i

P 1 ((1/12)∑Y.j..2) - (Y….2/24) 2 0 2 3 σ2ee + ((12/1)∑(P^)2j

K 1 ((1/12)∑Y..l.2) - (Y….2/24) 2 2 0 3 σ2ee + ((12/1)∑(K^)2l

e.e. 16 (∑∑∑∑Yijlm2)- ((1/3)∑∑∑Yijl.2 1 1 1 1 σ2ee

total 23 (∑∑∑∑Yijlm2) - (Y….2/24)

d. Obtenga un estimador de varianza para el error experimental.

● σ2ee = CME
5. Suponga un estudio donde se tiene el factor estado, el municipio, la
localización (rural-urbana) y las casas, en donde se desea estudiar la
presencia de agentes patógenos en la familia.
a. Presente el diagrama de estructura asociado con este estudio.
 b. Encuentre el modelo apropiado para el diagrama obtenido en a.
Yijkmno=μ + Pi +Tj(i) + Mk(ij) + Lm(ijk) + Cn(ijkm)+ eo(ijkmn)
Donde eo(ijkmn) es el error experimental, Cn(ijkm) es el efecto de las diferentes casas,
Lm(ijk) es el efecto de las diferentes localidades, Mk(ij) es el efecto de los municipios, Tj(i) es el
efecto de los estados, Pi es el efecto de los agentes patógenos en la familia y μ es la
media general del experimento.
c. Exponga la tabla de análisis de varianza suponiendo que el factor
municipio es aleatorio, estructurando bien las sumas de cuadrados,
esperanzas de los cuadrados medios y las pruebas F.
F.d.v G.l. S.C. i j k m n o E(cm)

Pi:μ e-1 ((1/ulcrb)∑Yi.....2) - (Y.....2/eulcrb) 0 u ℓ c r b σ2e+crbσ2M+(ulcrb/(e-

1))∑(P)2i

Tj(i):μP e(u-1) ((1/lcrb)∑∑Yij....2) - 1 0 ℓ c r b σ2e+crbσ2M+(lcrb/e(u-

((1/eulcrb)∑Yi…...2) 1))∑(T)2j(i)

Mk(ij):μP eu(ℓ-1) ((1/crb)∑∑∑Yijk...2)-((1/ℓcrb)∑∑∑Yij....2) 1 1 1 c r b σ2e+crbσ2M

T

Lm(ijk):μ euℓ(c-1) ((1/rb)∑∑∑∑Yijkm..2)- 1 1 1 0 r b σ2e+(rb/eul(c-1))∑(L)2m(ijk)

PTM ((1/crb)∑∑∑Yijk...2)

Cn(ijkm):μ euℓc(r-1) ((1/b)∑∑∑∑∑Yijkmn.2)- 1 1 1 1 0 b σ2e+(b/eulc(r-1))∑(C)2n(ijkm)

PTML ((1/rb)∑∑∑∑Yijkm..2)

Eo(ikjmn): euℓcr(b-1) (∑∑∑∑∑∑Yijkmno2)- 1 1 1 1 1 1 σ2e

μPTML ((1/b)∑∑∑∑∑Yijkmn.2)
C

Total euℓcrb-1 (∑∑∑∑∑∑Yijkmno2) - (Y.....2/eulcrb)

Hipótesis:
1- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(P)/CM(M); Rechazar si Ho si F> Fcrit α, e-1, eu(l-1)
2- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(T)/CM(M); Rechazar si Ho si F> Fcrit α, e(u-1), eu(l-1)
3- Ho= σ2M = 0; H1= Por lo menos un σ2M > 0
F= CM(M)/CM(E); Rechazar si Ho si F> Fcrit α; eu(l-1); eulcr(b-1)
4- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(L)/CM(E); Rechazar si Ho si F> Fcrit α; eul(c-1); eulcr(b-1)
5- Ho= μ1 = μ2 = μ3; H1= Por lo menos un par de μi = μj ∀ i ≠ j
F= CM(C)/CM(E); Rechazar si Ho si F> Fcrit α; eulc(r-1); eulcr(b-1)
6. En un estudio sobre crecimiento, desarrollo y adaptación de 6 materiales de
sorgo a las condiciones de los llanos orientales, se usó un diseño en bloques
completos aleatorizados con 6 genotipos, 3 replicaciones y 5 plantas por
replicación.
a. Presente el diagrama de estructura apropiado.

b. Presente el modelo adecuado señalando todos los supuestos

necesarios que permitan realizar pruebas de hipótesis de efectos
principales e interacciones.
Yijkl=μ + αi + βi + eij + ηk (ij)
Donde ηk (ij) es el error muestral, eij es el error experimental, βi es el efecto de los
bloques, αi es el efecto de los seis genotipos y μ es la media general del
experimento.
Supuestos:
● Yijlm es una variable aleatoria independiente que se distribuye N (μ ijlm, σ2)
● ee ~ N (0, σ2)
● Los efectos de bloques y tratamientos deben ser aditivos.
● em ~ N (0, σ2)
● Cov (ee, em) = 0
● E(ee2) = σ2ee
● E(em2) = σ2em
c. Construya la tabla ANOVA completa.
F.d.v G.l. S.C. i j k E(cm)

α 5 ((1/15)∑Yi..2) - (Y…2/90) 0 3 5 σ2em +5 σ2ee + (15/5)∑(α̂)2i

Β 2 ((1/30)∑Y.j.2) - (Y...2/90) 6 0 5 σ2em +5 σ2ee + ((30/2)∑(β)2j

e.e. 10 ((1/5)∑∑Yij.2) - ((1/15)∑Yi..2) 1 1 5 σ2em +5 σ2ee

- ((1/30)∑Y.j.2) + (Y…2/90)

e.m. 72 (∑∑∑Yijk2) - ((1/5)∑∑Yij.2) 1 1 1 σ2em

total 89 (∑∑∑Yijk2) - (Y…2/90)

7. Considere la estructura balanceada, donde se tienen 4 factores A, B, C y D,

donde A y B son cruzados, C está anidado en AB y D está anidado en C.
 a. Presente el diagrama de estructuras.

b. Exponga el modelo en notación estándar.

Yijlm=μ + Ai + Bj + (AB)ij+ Ck(ij) + Dm(ijk) + En(ijlkm)
Donde en(ijkm) es el error experimental, Dm(ijk) es el efecto del factor D, Ck(ij) es el efecto
del factor C, Ai es el efecto del factor A, Bj es el efecto del factor B, (AB)ij es la
interacción entre los factores A y B, y μ es la media general del experimento.
c. Escriba la tabla ANOVA, incluyendo la esperanza de los cuadrados
medios suponiendo todos los factores fijos.
F.d.v G.l. S.C. i j k m n E(cm)

Ai:μ a-1 ((1/bcde)∑Yi....2) - (Y….2/abcde) 0 b c d e σ2e + (bcde/(a-

1))∑(A)2i

Bj:μ b-1 ((1/acde)∑Y.j...2) - (Y….2/abcde) a 0 c d e σ2e + (acde/(a-

1))∑(B)2j

(AB)ij:μ (a-1)(b-1) ((1/cde)∑∑Yij...2) - 0 0 c d e σ2e + (cde/(a-1)(b-

((1/bcde)∑Yi....2) - 1))∑(AB)2ij
((1/acde)∑Y.j...2) + (Y….2/abcde)

Ck(ij):μAB ab(c-1) ((1/de)∑∑∑Yijk..2) - 1 1 0 d e σ2e + (de/(ab(c-

((1/cde)∑∑Yij...2) 1)))∑(C)2k(ij)

Dm(ijk):μABC abc(d -1) ((1/e)∑∑∑∑Yijkm.2) - 1 1 1 0 e σ2e + (e/(abc(d-

((1/de)∑∑∑Yijk..2) 1)))∑(D)2m(ijk)

En(ijkm):μABCD abcd(e-1) (∑∑∑∑∑Yijkmn2) - 1 1 1 1 1 σ2e

((1/e)∑∑∑∑Yijkm.2)

Total abcde-1 (∑∑∑∑∑Yijkmn2) - (Y….2/abcde)