Существуют ли парные постхокские сравнения для теста хи-квадрат в R?

Question

Мне интересно, существует ли в R пакет / функция для выполнения: «Последовательных парных сравнений для теста однородности пропорций хи-квадрат» (или его эквивалента), который описан здесь: http://epm.sagepub.com/cgi/content/abstract/53/4/951

В моей ситуации я просто делаю тест чи на матрице 2 на Х. Я нашел разницу, но я хочу знать, какой из столбцов «отвечает» за разницу.

Спасибо, Tal

Answer 1

См. Пакет fifer для функции chisq.post.hoc ()

Answer 2

"Критерий хи-квадрат" обычно генерируется как сумма квадратов индивидуальных отклонений ячеек от "ожидаемых" = произведений сумм строк и столбцов, деленная на общую сумму. Таким образом, можно сравнить вклады отдельных ячеек в сумму с критическим значением хи-квадрата с 1 d.f. Это довольно простая задача изменить код chisq.test () так, чтобы он возвращал хи-квадраты ячейки. Я просто добавил:

cell.chisq = (x - E)^2/E,

на вызов структуры в конце. Они не получат print () - ed, но вы можете присвоить результат объекту и использовать:

 obj$cell.chisq

Answer 3

Не уверен, что это работает для вашей проблемы, но я прочитал статью , в которой для выполнения парного теста chi-sqr

используется следующее

rcompanion::pairwiseNominalIndependence()