Существуют ли парные постхокские сравнения для теста хи-квадрат в R?
Мне интересно, существует ли в R пакет / функция для выполнения: "Последовательных парных сравнений для критерия хи-квадрат однородности пропорций" (или его эквивалента), который описан здесь: http://epm.sagepub.com/cgi/content/abstract/53/4/951
Моя ситуация заключается в том, что я просто делаю тест Чи на матрице 2 на Х. Я нашел разницу, но я хочу знать, какой из столбцов "отвечает" за разницу.
Спасибо тал
3 ответа
"Критерий хи-квадрат" обычно генерируется как сумма квадратов отклонений отдельных ячеек от "ожидаемых" = произведений сумм строк и столбцов, деленная на общую сумму. Таким образом, можно сравнивать вклад отдельных ячеек в сумму с критическим значением хи-квадрата с 1 df. Довольно просто изменить код chisq.test() так, чтобы он возвращал хи-квадраты ячеек. Я только добавил:
cell.chisq = (x - E)^2/E,
на вызов структуры в конце. Они не получат print()-ed, но вы можете присвоить результат объекту и использовать:
obj$cell.chisq
Не уверен, что это работает для вашей проблемы, но я прочитал статью, которая использует следующее для выполнения попарного теста chi-sqr
rcompanion::pairwiseNominalIndependence()