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

[Proyecto R-UCA]

Hola, Manuel:

Si haces repetidos contrastes a un nivel de confianza 1-alpha el resultado final no tendrá ese nivel de confianza. Hay que ajustar los
niveles usando la corrección de Bonferroni.

Por otra parte, el test que estás haciendo no es sobre una casilla es sobre la tabla entera. Si quieres analizar las casillas tendrás que
utilizar otro tipo de test. Tal vez uno basado en errores estandarizados.

No recuerdo ni encuentro ninguno en este momento. Un hilo del que tirar podría ser:
https://pubmed.ncbi.nlm.nih.gov/24788323/

Un saludo.


El dom, 12-04-2026 a las 19:17 +0000, MANUEL MORALES ORTIZ escribió:
> Hola, quiero construir una funcin 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
> sera 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 continuacin, 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://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-help-es__;!!D9dNQwwGXtA!W0VL6NufuD4gWEkqY5ZCuu6w5RV4JvHth-kVWEWTvR9OyPbMbVMVItbUdLmpq24bnGLzc3YPGzI$
>