получить метки соответствующей категории для значимых остатков Пирсона - PullRequest
Новая
       9

получить метки соответствующей категории для значимых остатков Пирсона

0 голосов
/ 27 ноября 2018

Я рассчитываю Chi2-тест, чтобы проверить, есть ли корреляция между двумя категориальными переменными.После этого я также проверяю, какой из стандартизированных остатков Пирсона (который можно извлечь из chisq.test ()) становится значимым, чтобы найти локальные различия.Моя цель - получить вектор с соответствующими ярлыками категорий для каждого значащего значения ячейки.Вот пример кода: 

set.seed(43)
A <- sample(c("small", "med", "large"), 200, replace = T)
B <- sample(c("x", "y", "z"), 200, replace = T, prob = c(0.7, 0.15, 0.15))

chisq_table <- table(A,B)

chisq_result <- chisq.test(chisq_table)

stdres <- chisq_result$stdre

alpha <- 0.05

# get critical upper and lower z-value
upper_z <-  qnorm(1-alpha/2)
lower_z <- qnorm(alpha/2)

testResiduals <- function(residual){

  if (residual < lower_z){
    return(TRUE)
  } else if(residual > upper_z){
    return(TRUE)
  } else {
    return(FALSE)
  }
}

posthoc_result <- sapply(stdres,FUN=testResiduals)
posthoc_result

posthoc_result вернет логический вектор, и, как вы можете видеть, одна комбинация двух переменных стала значимой.Я хочу знать, какая комбинация это.Поэтому оптимально я хотел бы получить вектор, такой как, например,> essential_residuals: "x & medium"

PS: я знаю, что chisq.test () не получил значительный и что он не будетв этом случае целесообразно интерпретировать стандартизированные остатки Пирсона.Я также знаю, что было бы лучше также исправить альфа-уровень при проверке остатков на значимость.Я упустил эти «статистические проблемы» для лучшей читабельности.

...