Чи можна використовувати тест Манна-Вітні для постсоціальних порівнянь після Крускала-Уолліса?


9

У мене є симуляція, коли тварина поміщена у вороже середовище і приурочена до того, щоб побачити, як довго вона може вижити, використовуючи певний підхід до виживання. Існує три підходи, якими він може скористатися, щоб вижити. Я провів 300 моделювання тварини, використовуючи кожен підхід до виживання. Усі симуляції проводяться в одному середовищі, але є певна випадковість, тому вона щоразу відрізняється. Я час, скільки секунд виживає тварина в кожному моделюванні. Жити довше - краще. Мої дані виглядають так:

Approach 1, Approach 2, Approach 2
45,79,38
48,32,24
85,108,44
... 300 rows of these

Я не впевнений у всьому, що роблю після цього моменту, тому дайте мені знати, чи роблю щось дурне і неправильне. Я намагаюся з’ясувати, чи є статистична різниця в тривалості життя, використовуючи певний підхід.

Я провів тест Шапіро на кожному з зразків, і вони повернулися з крихітними значеннями р, тому я вважаю, що дані не нормалізуються.

Дані про рядки не мають стосунків один до одного. Випадкове насіння, використовуване для кожного моделювання, було різним. Як результат, я вважаю, що дані не спарені.

Оскільки дані не нормалізуються, не спарюються і було більше двох зразків, я провів тест Крускала Уолліса, який повернувся з р-значенням 0,048. Потім я перейшов до посади, обравши Манна Вітні. Насправді не впевнений, чи варто тут використовувати Манна Вітні.

Я порівнював кожен підхід до виживання один з одним, виконуючи тест Манна Вітні, тобто {(підхід 1, підхід 2), (підхід 1, підхід 3), (підхід 2, підхід 3)}. Не було виявлено статистичної значущості між парами (підхід 2, підхід 3) за допомогою двосхилого тесту, але була різниця значущості, виявлена ​​за допомогою односхилого тесту.

Проблеми:

  1. Я не знаю, чи має таке використання Манна Вітні.
  2. Я не знаю, чи варто мені скористатися одним або двома хвостиками Манна Вітні.

Чи є у вас апріорні гіпотези щодо відносної сили різних підходів (наприклад, підхід1> підхід2> підхід3)? Це дуже важливо для відповіді на ваші запитання.
амеба

У мене є середнє, середнє та стандартне відхилення, і, схоже, підхід 3 є кращим, оскільки він має вищу медіану та середню, але він також має значно більш високе стандартне відхилення, тому я не впевнений. Але я не міг цього знати.
Phlox Midas

Або він також відомий як корекція Бонферроні?
Флокс Мідас

Флокс: якщо раніше "не було відомо про це", вам абсолютно не слід використовувати однохвостий тест, а лише двосторонній (як @Alexis згадував у своїй відповіді).
амеба

6
@amoeba "її";)
Олексій

Відповіді:


15

Ні, не слід використовувати Манн-Вітні U випробування на цю обставину.

Ось чому: тест Данна є підходящим постфактум тест * після відхилення тесту Круськала-Уолліса. Якщо перейти від відмови від Крускала-Уолліса до виконання звичайних парних тестів (наприклад, Вілкоксона або Манна-Уітні), то виходять дві проблеми: (1) ранги, які використовуються для парних тестів на рангову суму не чини, використані тестом Крускала-Уолліса; та (2) у тестах рангової суми не використовується об'єднана дисперсія, що передбачається нульовою гіпотезою Крускала-Уолліса. У тесту Данна цих проблем немає

Пост- тести після відхилення тесту Крускала-Уолліса, які були скориговані для декількох порівнянь, можуть не відхилити всі парні тести для заданої сімейної помилки або заданої помилкової частоти виявлення, що відповідає заданійαдля тесту омнібус, як і будь-який інший сценарій багаторазового порівняння (omnibus / post-hoc) .

Якщо у вас немає підстав вважати, що час виживання однієї групи довший або коротший, ніж апріорний , вам слід скористатися двосторонніми тестами.

Тест Данна можна виконати в Stata, використовуючи dunntest (тип net describe dunntest, from(https://www.alexisdinno.com/stata)), а в R за допомогою пакету dunn.test .

Також мені цікаво, чи можете ви скористатися підходом до аналізу виживання, щоб оцінити, чи і коли тварина гине на основі різних умов?


* Кілька менш відомих постспеціальних парних тестів для відхилення Крускала-Уолліса включають Коновер-Іман (як Данн, але на основі розподілу t , а не розподілу z , реалізованого для Stata в пакеті conovertest , і для R у пакеті conover.test ), і для тестів Dwass-Steel-Citchlow-Fligner.


Дякую за вашу відповідь. Тест Данна також відомий як тест Неменій-Даміко-Вулф-Данн чи це окремий тест?
Флокс Мідас

Я запитую, бо не можу знайти жодної реалізації тесту Данна.
Phlox Midas

@PhloxMidas Я не знаю про "тест Неменія-Даміко-Вулфа-Данна", але Вікіпедія передбачає, що це відповідний постспеціальний тест після відмови від тесту омнібуса в проекті повторних заходів - наприклад, після тесту Фрідмана. Також дивіться мій коментар щодо Stata.
Олексій

7

Уніфікуючим узагальненням Крускал-Уолліс / Вілкоксона є модель пропорційних шансів, яка допускає загальні контрасти з точковими або одночасними довірчими інтервалами для коефіцієнтів шансів. Це реалізовано в моїх rmsпакетах ormі contrast.rmsфункціях R.


1

Ви також можете використовувати критичну різницю після Коновера або критичну різницю після Шайха і Хамерле. Перший є більш ліберальним, тоді як останній є точним, але йому не вистачає сил. Обидва способи проілюстровані на моєму веб-сайті svetstat.com, а webapp яскравих даних також дозволяє обчислити ці критичні відмінності та виконати пост-спеціальні тести відразу. Крускал-Уолліс на svetstat.com


-1

Якщо ви використовуєте SPSS, виконайте пост-хан Манна-Вітні з корекцією Бонферроні (значення p, поділене на кількість груп).


"Ман-Вітні" страждає від двох проблем, які я виявляю у своїй відповіді, і є невідповідним пост-спеціальним тестом для Крускала-Уолліса.
Олексій
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.