Manera apropiada de lidiar con una tabla de contingencia de 3 niveles


12

Tengo una tabla de contingencia de tres niveles, con datos de conteo para varias especies, la planta huésped de la cual fueron recolectadas y si esa recolección ocurrió en un día lluvioso (¡esto realmente importa!). Usando R, los datos falsos podrían ser algo como esto:

count    <- rpois(8, 10)
species  <- rep(c("a", "b"), 4)
host     <- rep(c("c","c", "d", "d"), 2)
rain     <- c(rep(0,4), rep(1,4))
my.table <- xtabs(count ~ host + species + rain)


, , rain = 0

    species
host  a  b
   c 12 15
   d 10 13

, , rain = 1

    species
host  a  b
   c 11 12
   d 12  7

Ahora, quiero saber dos cosas: ¿Están las especies asociadas con las plantas hospederas? ¿La lluvia o no afecta esta asociación? Solía loglm()de MASSpara esto:

 # Are species independent to host plants, given the effect of rain?
loglm(~species + host + rain + species*rain + host*rain, data=my.table)

 # Given any relationship between host plants and species, does rain change it?
loglm(~species + host + rain + species*host)

Esto está un poco fuera de mi nivel de comodidad, y quería comprobar que había configurado correctamente los modelos y que esta era la mejor manera de abordar estas preguntas.

Respuestas:


10

Hay dos formas de interpretar su primera pregunta, que se reflejan en las dos formas en que la hizo: "¿Las especies están asociadas con las plantas hospederas?" y, "¿Son las especies independientes de las plantas hospederas, dado el efecto de la lluvia?"

La primera interpretación corresponde a un modelo de independencia conjunta , que establece que las especies y los huéspedes son dependientes, pero conjuntamente independientes de si llovió:

pshr=pshpr

pshr(s,h,r)shrpsh(s,h,)pr

La segunda interpretación corresponde a un modelo de independencia condicional , que establece que las especies y los huéspedes son independientes dado que llovió:

psh|r=ps|rph|rpshr=psrphr/pr

psh|r(s,h,r)r

Puede probar estos modelos en R ( loglintambién funcionaría bien, pero estoy más familiarizado con esto glm):

count <- c(12,15,10,13,11,12,12,7)
species <- rep(c("a", "b"), 4)
host <- rep(c("c","c", "d", "d"), 2)
rain <- c(rep(0,4), rep(1,4))
my.table <- xtabs(count ~ host + species + rain)
my.data <- as.data.frame.table(my.table)
mod0 <- glm(Freq ~ species + host + rain, data=my.data, family=poisson())
mod1 <- glm(Freq ~ species * host + rain, data=my.data, family=poisson())
mod2 <- glm(Freq ~ (species + host) * rain, data=my.data, family=poisson())
anova(mod0, mod1, test="Chi") #Test of joint independence
anova(mod0, mod2, test="Chi") #Test of conditional independence

mod1mod2mod0pshr=psphprsummary(mod2)

Una forma diferente de abordar su primera pregunta sería realizar la prueba exacta de Fischer ( fisher.test(xtabs(count ~ host + species))) en la tabla colapsada 2x2 (primera interpretación) o la prueba Mantel-Haenszel ( mantelhaen.test(xtabs(count ~ host + species + rain))) para tablas 2x2 estratificadas 2 o escribir una prueba de permutación que respete la estratificación (segunda interpretación).

Parafraseando su segunda pregunta, ¿la relación entre las especies y el huésped depende de si llovió?

mod3 <- glm(Freq ~ species*host*rain - species:host:rain, data=my.data, family=poisson())
mod4 <- glm(Freq ~ species*host*rain, data=my.data, family=poisson())
anova(mod3, mod4, test=”Chi”)
pchisq(deviance(mod3), df.residual(mod3), lower=F)

El modelo completo mod4está saturado, pero puede probar el efecto en cuestión observando la desviación de mod3lo que he hecho anteriormente.


Gracias Lockedoff, especialmente por ayudarme a resolver mi propio pensamiento con respecto a la diferencia entre el modelo de independencia condicional y conjunta
David w

1

hostrainhostrain

El comando R sería:

glm (fórmula = especie ~ huésped + lluvia, familia = binomio (logit), pesos = recuentos)

p


1
La regresión logística parece estar bien, pero tiene la restricción adicional de que el total de filas y columnas es fijo. Este puede no ser el caso con los datos de Poisson. Creo que las respuestas no diferirán mucho.
suncoolsu

1

Inicialmente sugerí probar una de las técnicas de ordenación restringida del veganpaquete, pero pensándolo bien, dudo que esto sea útil, ya que en realidad tienes 2 tablas de contingencia. Espero que la segunda parte de este ejemplo [PDF: Demostración R - Análisis categórico] sea útil.


Piensa que el enlace está roto, ¿te refieres a este categórico aquí ? Eso fue útil, gracias!
David w

Sí, parece que el espacio en la URL lo rompe.
ILS
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.