[R-es] Crear funcion para identificar residuales en tabla de contingencia

[Jorge I Velez]

Hola Manuel,
Interesante.
El inconveniente es qué haces con la celda que tiene residuales
significativos. La remueves?  La haces cero?  Igualas su frecuencia al
valor esperado?
Adjunto un ejemplo donde se utiliza la prueba de independencia chi^2 en
lugar del modelo loglineal y la celda "problema" se reemplaza por su valor
esperado bajo independencia. La función la construí con ayuda de Gemini.
Espero sea de ayuda.
Saludos,
Jorge Vélez.-


On Sun, Apr 12, 2026 at 2:17 PM MANUEL MORALES ORTIZ <morales using us.es> wrote:

> Hola, quiero construir una función en la que me identifique los residuales
> significativos de una tabla de contingencia y me vaya analizando la tabla
> hasta que se acepte el modelo de independencia y no queden residuales
> significativos. Un ejemplo tentativo de lo que se pretende sería el
> siguiente:
>
> library(MASS)
> x<- c(13,6,7,9,9,25,5,6,6,7,6,14,7,8,14,46)
> datos<- matrix(x,4,4,byrow=T)
> inicio<- matrix(rep(1,16),4,4)
> res<- loglm(~1+2,data= datos,fitted=T,start = inicio)
> resid<- residuals(res)
> resid<-abs(as.vector(resid))
> max<- which.max(resid)
>
> # Asigno el valor 0 a la caslla de la tabla que tiene el mayor residual
> datos<- as.vector(datos)
> datos[6]<- 0
> datos<- matrix(datos,4,4)
>
> # Vuelvo a analizar la tabla
> inicio[2,2] <- 0
> res<- loglm(~1+2,data= datos,fitted=T,start = inicio)
> chi2<-res$pearson
> aux<- chi2 > qchisq(.95,9)
>
> # A continuación, se calculan los residuales,
>
> resid<- residuals(res)
>
> # Repetir el procedimiento hasta que no queden residuales significativos.
>
> ¿Alguna sugerencia?
>
> Saludos
>
>
>
>
>
>
>
> --
>
>
>
>
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20260413/61189791/attachment-0001.html>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: manuel-residuales.pdf
Type: application/pdf
Size: 153974 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20260413/61189791/attachment-0001.pdf>