Остаточный анализ на основе результатов критерия хи-квадрат независимости в R

Question

Как я могу провести остаточный анализ на основе результатов критерия хи-квадрат независимости в R, как на картинке?

## From Agresti(2007) p.39
M <- as.table(rbind(c(762, 327, 468), c(484, 239, 477)))
dimnames(M) <- list(gender = c("F", "M"),
                    party = c("Democrat","Independent", "Republican"))
(Xsq <- chisq.test(M))  # Prints test summary\

Answer 1

Я нашел ответ:

df <- as.data.frame(M)
df

library(gmodels)

x <- CrossTable(df$gender, 
           df$party, 

           chisq=T, 
           expected=T, 
           asresid=T, 
           format="SPSS", 
           prop.t=F, 
           prop.chisq=F, 
           prop.r=F, 
           prop.c=F, 
           digits=2)

Answer 2

Вы можете получить доступ к стандартизированным остаткам в chsq$stdres, где chsq - это возвращаемый объект chisq.test.

Для использования примера из ?chisq.test

M <- as.table(rbind(c(762, 327, 468), c(484, 239, 477)))
dimnames(M) <- list(
    gender = c("F", "M"),
    party = c("Democrat","Independent", "Republican"))

chsq <- chisq.test(M)

chsq$stdres
#      party
#gender   Democrat Independent Republican
#     F  4.5020535   0.6994517 -5.3159455
#     M -4.5020535  -0.6994517  5.3159455