Anova Anidado
Anova Anidado
Anova Anidado
redes
□ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ R4 A Oligo 35,7 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
□ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ R5 A Oligo 22,6 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
□ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ R6 B Eutr 61,3 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
1
ANOVA anidado – ejemplo 1 ANOVA anidado – ejemplo 2
Nos interesa saber si la cantidad de algas azules cambia en los
lagos eutrofizados. Se escogieron aleatoriamente 10 lagos Una bióloga de la UMSA está estudiando si la
eutróficos y 10 lagos oligotróficos en el dept. Beni. Se usaron frecuencia del ataque del murciélago vampiro
los redes especiales para cuantificar la cantidad de las algas (Desmodus rotundus) varia dependiendo del tipo
azules.
Unidad experimental (subgrupos): lagos (Ο) de vegetación alrededor de la comunidad en el
Unidad de muestreo: redes para algas (□) depto. Potosí.
La cantidad de algas azules observada en cada red Ella escogió 7 comunidades donde se han
– submuestra. reportados las ataques de D. rotundus.
1)Usar ANOVA anidado Estas 7 comunidades tiene diferente vegetación
2)Usar la media de algas por lago como var. dep. alrededor:
Red Lago Tipo Algas oligotróficos eutróficos Bosque nativo – 2 comunidades
R1 A Oligo 23,5 Ο Ο ΟΟ Ο Ο ΟΟ Ο Ο Ο ΟΟ Ο Ο Ο ΟΟ Ο Ο Bosque de Eucalyptus – 3 comunidades
R2 A Oligo 31,8
□ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ Sin vegetación – 2 comunidades
R3 A Oligo 13,7 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
¿Hay diferencia de las ataques de D. rotundus
redes
R4 A Oligo 35,7
□ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
R5 A Oligo 22,6 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ entre 3 tipos de vegetación?
R6 B Eutr 61,3 □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □ □
Com. 6
Com. 3 ***
Com. 6 ****
Com. 1
Com. 2
Com. 7
Com. 4
Com. 5
Com. 1 ***
Com. 2 ***
Com. 7 ***
Com. 4 ***
Com. 5 ****
C.5 C.6 C.7 C.5 C.6 C.7
BE Com. 6 ****
BN BE SV
S Com. 5 ****
V Com. 4 ***
Com. 7 ***
2
ANOVA anidado ANOVA anidado
ejemplo - tucanes ejemplo - tucanes
e.g., Quisiéramos saber si hay diferencia en el e.g., Quisiéramos saber si hay diferencia en el
número de especies de tucanes según la altura número de especies de tucanes según la altura
de cerros en Sudamérica (Alto, Medio y Bajo) de cerros en Sudamérica (Alto, Medio y Bajo)
Se escogieron 8 cerros (2-3 cerros para cada Se escogieron 8 cerros (2-3 cerros para cada
altura). altura).
En cada cerro, se establecieron 4 transectos de 1 En cada cerro, se establecieron 4 transectos de 1
km para realizar el censo. km para realizar el censo.
Diseño de este estudio: Modelo: yijk = μ + αi + β j(i) + εk(ij)
Variable dep.: # de spp. de tucanes Hipótesis:
Factor (Grupo): Altura H0: αi = 0 – efecto de grupo (altura)
3
ANOVA anidado ANOVA anidado
- Cálculo de SS - Cálculo de SS
e.g., número de especies de tucanes entre 3 alturas. e.g., número de especies de tucanes entre 3 alturas.
3 alturas (grupos) -> 8 cerros (subgrupos) -> 32 transectos 3 alturas (grupos) -> 8 cerros (subgrupos) -> 32 transectos
yijk = μ + αi + β j(i) + εk(ij) yijk = μ + αi + β j(i) + εk(ij)
SST = SSA + SSB(A) + SSE SST = SSA + SSB(A) + SSE
SS total = SS entre grupos + SS entre subgrupos SST = ∑(yijk - y...)2 – SS total
(dentro de grupos) + SS dentro de subgrupos
4
ANOVA anidado ANOVA anidado
- Cálculo de SS y gl - Cálculo de SS y gl
e.g., número de especies de tucanes entre 3 alturas. e.g., número de especies de tucanes entre 3 alturas.
3 alturas (grupos) -> 8 cerros (subgrupos) -> 32 transectos 3 alturas (grupos) -> 8 cerros (subgrupos) -> 32 transectos
yijk = μ + αi + β j(i) + εk(ij) yijk = μ + αi + β j(i) + εk(ij)
SST = ∑(yijk - y...)2 glT = N - 1 SST = ∑(yijk - y...)2 glT = N - 1
SSA = ∑(yi.. - y...)2 glA = a - 1 = 3 - 1 SSA = ∑(yi.. - y...)2 glA = a - 1
SSB(A) = ∑(yij. - yi..)2 SSB(A) = ∑(yij. - yi..)2 glB(A) = ∑(bi-1) = ∑bi-a
SSE = ∑(yijk - yij.)2 SSE = ∑(yijk - yij.)2 =8–3
Grupo ALTO MEDIO BAJO Grupo ALTO MEDIO BAJO
Subgrupo S.Mal Agua Saja Apolo S.Lui Chul. Malta Boqu Subgrupo S.Mal Agua Saja Apolo S.Lui Chul. Malta Boqu
5 5 3 5 3 7 2 2 5 5 3 5 3 7 2 2
7 10 7 7 6 5 4 5 7 10 7 7 6 5 4 5
8 9 10 10 9 9 5 6 8 9 10 10 9 9 5 6
11 13 14 13 12 3 10 8 11 13 14 13 12 3 10 8
Subgrupo y11. y12. y13. y21. y22. y31. y32. y33. Subgrupo y11. y12. y13. y21. y22. y31. y32. y33.
Grupo y1.. y2.. y3.. y... Grupo y1.. y2.. y3.. y...
5
ANOVA anidado
- Cálculo de F Tabla de ANOVA anidado
e.g., número de especies de tucanes entre 3 alturas.
Fuente SS gl MS F P
3 alturas (grupos) -> 8 cerros (subgrupos) -> 32 transectos
yijk = μ + αi + β j(i) + εk(ij) Grupo SSA a-1 MSA MSA/MSB(A) ***
(Factor)
SST = SSA + SSB(A) + SSE
Subgrupo SSB(A) ∑bi - a MSB(A) MSB(A)/MSE ***
SS total = SS entre grupos + SS entre subgrupos
(dentro de grupos) + SS dentro de subgrupos Error SSE N - ∑bi MSE
Total SST N-1
{H0: β j(i) = 0} – efecto de subgrupo
¿# de tucanes varia entre subgrupos (cerros)? Para probar el efecto de Grupo (Factor),
F = MS entre subgrupos / MS dentro de subgrupos utilizar MSB(A) (∑bi - a) en lugar de MSE
F = MSB(A) / MSE (N - ∑bi) como denominador.
P = ***, glB(A) y glE {Para probar el efecto de Subgrupo,
MS denominador es MSE (MS dentro de
utilizar MSE.} SPSS – ANOVA anidado
subgrupos).
F y P no
sirven!
Hay diferencia
/DESIGN = Grupo Subgrupo(Grupo) Origen SS gl MS F P
significativa en la
riqueza de spp. de
/TEST = Grupo VS Subgrupo(Grupo). Altura 61,6 2 30.8 16,9 0,006 tucanes entre diferentes
Cerro(Altura) 9,1 5 1,8 0,7 0,97 alturas (F2,5 = 17, P <
Error 271 24 11,3 0,05)
ANOVA anidado post hoc (SPSS) ANOVA anidado post hoc (SPSS)
Editor de sintaxis
Hipótesis:
/POSTHOC=Grupo(TUKEY) VS Subgrupo(Grupo)
H0: μALTO = μBAJO
H0: μALTO = μMEDIO
H0: μMEDIO = μBAJO
6
ANOVA anidado como ANOVA de
una vía con las medias
ANOVA anidado con Tucán
ANOVA con diseños
más complejos
ANOVA de 1 vía con las medias de Tucán
ANOVA anidado con más de una ANOVA anidado con más de una
jerarquía (subgrupo, subsubgrupo) jerarquía (subgrupo, subsubgrupo)
Nos interesa saber si la cantidad de algas azules Nos interesa saber si la cantidad de algas azules
cambia en los lagos eutrofizados. Se escogieron cambia en los lagos eutrofizados. Se escogieron
aleatoriamente 10 lagos eutróficos y 10 lagos aleatoriamente 10 lagos eutróficos y 10 lagos
oligotróficos en el dept. Beni. oligotróficos en el dept. Beni.
Para muestrear todo el lago, cada lago fue Para muestrear todo el lago, cada lago fue
dividido a 2-6 zonas según su tamaño. dividido a 2-6 zonas según su tamaño.
En cada zona, se usaron 3 redes de plánctones. En cada zona, se usaron 3 redes de plánctones.
Var.dep., Unid.Exp., Unid.Muestreo, Grupo, Var.dep., Unid.Exp., Unid.Muestreo, Grupo,
Subgrupo, etc.? Subgrupo, etc.?
Var.dep.: Cantidad de algas azules
Grupo (Factor): Tipo de lago
Subgrupo (Unid.Exp.): Lago
Subsubgrupo: Zona
Unid.Muestreo: Red de plánctones
7
ANOVA anidado con más de una
jerarquía (subgrupo, subsubgrupo) ANOVA anidados donde hay
Origen SS gl MS F P factores en ambos niveles de
Tipo de SSA a-1 MSA MSA ***
lago
Lago SSB(A) ∑(b-1) MSB(A)
/MSB(A)
MSB(A) ***
grupos y subgrupos
/MSC(B(A))
Zona SSC(B(A)) ∑(c-1) MSC(B(A)) MSC(B(A)) ***
/MSE
Error SSE MSE
yijkl = μ + αi + β j(i) + γ k(j(i)) + εijkl
αi: efecto de grupo (tipo de lago)
β j(i): efecto de subgrupo (lago)
γk(j(i)): efecto de subsubgrupo (zona)
SST = SSA + SSB(A) + SSC(B(A)) + SSE
ANOVA anidados donde hay factores en ANOVA anidados donde hay factores en
ambos niveles de grupos y subgrupos ambos niveles de grupos y subgrupos
Una bióloga de la UMSA está estudiando si la La incidencia del murciélago vampiro varia según las
frecuencia del ataque del murciélago vampiro características de comunidad y de corral?
(Desmodus rotundus) varia según las Incidencia fue medida en corral.
características de comunidad y de corral en el Las 12 comunidades:
depto. Potosí. 3 tipos de vegetación (B.nativo, B.Euc. y SV)
Ella escogió 12 comunidades y la abundancia 2 niveles de distancia al bosque (lejos y cerca).
(incidencia) de D. rotundus fue medido en 3-4 3-4 corrales/comunidad:
corrales en cada comunidad. 2 tipos de techo (abierto, cerrado).
Las comunidades tiene uno de 3 tipos de
vegetación (bosque nativo, bosque de Eucaliptos
y sin vegetación) y 2 niveles de distancia al
bosque (lejos y cerca).
Hay diferentes tipos de corrales: techo (abierto,
cerrado).
ANOVA anidados donde hay factores en ANOVA anidados donde hay factores en
ambos niveles de grupos y subgrupos ambos niveles de grupos y subgrupos
La incidencia del murciélago vampiro varia según las Origen SS gl MS F P
características de comunidad y de corral? Vegetación SSA a-1 MSA MSA/MSB **
Incidencia fue medida en corral.
Distancia SSC c-1 MSC MSC/MSB **
Las 12 comunidades:
3 tipos de vegetación (B.nativo, B.Euc. y SV) Int.V*D SSAC (a-1)(c-1) MSAC MSAC/MSB **
2 niveles de distancia al bosque (lejos y cerca). Comunidad SSB(AC) ∑∑(b-1) MSB MSB(AC)/MSE **
3-4 corrales/comunidad: Techo SSF f-1 MSF MSF/MSE **
2 tipos de techo (abierto, cerrado). Error SSE MSE
Var.dep., Grupo, Subgrupo, Unidades?
Var. dep.: Incidencia de vampiro yijklm = μ + αi + γ j + (αγ)ij + β k(ij) + δl + εijklm
Grupo (factor en el nivel de unid.exp.): vegetación, αi, γj: efecto de factores principales
distancia al bosque (αγ)ij: efecto de interacciones entre A y C
Factor en el nivel de unid.muestreo: techo de corral β k(ij): efecto de subgrupo
Subgrupo (Unid.Exp.): Comunidad δl: efecto de factor en el nivel de unid.muest.
Unid.Muestreo: Corral SST = SSA + SSC + SSAC + SSB(AC) + SSF + SSE