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

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

Моя ситуация заключается в том, чтобы просто провести тест хи на матрице 2 на X. Разницу нашел, но хочу узнать, какая из колонок "отвечает" за разницу.

Спасибо, Тал


r statistics posthoc chi-squared
person Tal Galili    schedule 17.03.2010    source источник
comment
Таль, у меня ощущение, что это здесь не по теме. Вы знаете, где найти r-help, так как в любом случае последние несколько дней вы занимаетесь кросспостингом.   -  person Dirk Eddelbuettel    schedule 17.03.2010
comment
Привет, Дирк, я ценю твое предложение и подозреваю, что твои чувства окажутся правильными. Причина, по которой я разместил его здесь, заключается в том, что я уже задавал этот вопрос в R-help несколько месяцев назад и не получил ответа. Поэтому я подумал проверить, не нашел ли кто-нибудь здесь решение. Бест, Таль   -  person Tal Galili    schedule 17.03.2010
comment
Привет, Шейн, кажется, и ты, и Дирк правы. Кроме того, я не знал о mathoverflow - спасибо за ссылку!   -  person Tal Galili    schedule 17.03.2010
comment
@Tal: возможно, вы можете улучшить этот вопрос теперь, когда у вас есть ученая степень в области статистики? Я вижу, что цитируемая статья подверглась критике как статистически некорректная: semanticscholar.org/paper/ (это не редкость, когда речь идет о множественных сравнениях.)   -  person IRTFM    schedule 08.03.2021

Ответы (3)


arrow_upward
3
arrow_downward

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

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

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

 obj$cell.chisq
person IRTFM    schedule 20.03.2010
comment
Спасибо Двин. Тем не менее, мой вопрос был связан со сравнением двух строк (или столбцов), чтобы выяснить, какие пары являются значимыми. Можно было бы попробовать все попарные таблицы столбцов и вычислить их значения P, но как тогда скорректировать множественность? - person Tal Galili; 21.03.2010
comment
Суммируйте строки obj$cell.chisq и примените корректировку Бонферрони к p-значениям, полученным из критических значений chisq. Ключевым моментом является то, что тесты chisq для таблиц можно разложить по ячейкам, по строкам, по столбцам или по их комбинациям. Люди вводятся в заблуждение своим курсом внутренней статистики, заставляя думать, что критерий хи-квадрат означает только одно, хотя на самом деле он означает много вещей. Аналитик по-прежнему должен четко держать в голове, что именно тестируется и сколько степеней свободы потребовал весь его анализ. - person IRTFM; 04.04.2010

arrow_upward
1
arrow_downward

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

person user2521252    schedule 16.01.2015
comment
Этот пакет больше не поддерживается, начиная с версии R 3.5.0. - person Nakx; 02.02.2019

arrow_upward
0
arrow_downward

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

rcompanion::pairwiseNominalIndependence()
person Tianyi    schedule 08.10.2017