Інтерпретація виходу регресії зі змішаної моделі, коли включаються взаємодії між категоричними змінними


14

У мене є питання щодо мого використання змішаної моделі / лмера. Основна модель така:

lmer(DV ~ group * condition + (1|pptid), data= df)

Група і стан - це обидва фактори: група має два рівні (група A, група B), а стан має три рівні (умова1, умова2, умова3). Це дані від людських предметів, тому pptid є випадковим ефектом для кожної людини.

Модель виявила наступне із значенням p вихід:

                                   Estimate MCMCmean HPD95lower HPD95upper  pMCMC Pr(>|t|)
(Intercept)                          6.1372   6.1367     6.0418     6.2299 0.0005   0.0000
groupB                              -0.0614  -0.0602    -0.1941     0.0706 0.3820   0.3880
condition2                           0.1150   0.1151     0.0800     0.1497 0.0005   0.0000
condition3                           0.1000   0.1004     0.0633     0.1337 0.0005   0.0000
groupB:condition2                   -0.1055  -0.1058    -0.1583    -0.0610 0.0005   0.0000
groupB:condition3                   -0.0609  -0.0612    -0.1134    -0.0150 0.0170   0.0148

Тепер я знаю, що перелічені рядки порівнюють кожен рівень факторів із еталонним рівнем. Для групи посиланням є група А, а для умови - посиланням умова1.

Чи правильно я трактую цей вихід таким чином:

  • Немає загальних відмінностей між групами (отже, група B має ap> 0,05)
  • Загальні відмінності між умовою 1 та умовою 2 та між умовою 1 та умовою 3.
  • Відмінності між групою A, умовою 1 проти групи B, умовою 2, а також між групою A, умовою 1 проти групи B, умовою 3.

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

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


Але, якщо ми хочемо порівняти групу = B з еталонним рівнем Group = A, коли Умова = 2 (або 3)? Можливо? І я вважаю, що порівняйте, якщо "різниця між умовою1 та умовою2 є різною, коли група = А проти групою = В" є однаковою для порівняння, якщо "різниця між групою = А та групою = В є різною, коли Умова1 проти Умовою2 ". Це правильно? в іншому випадку, які їх p-значення?

Це не схоже на відповідь на запитання. Швидше у вас нове запитання. Найкраще публікувати як таке.
Нік Кокс

Відповіді:


21

Використовуючи дану таблицю регресії, ми можемо обчислити таблицю очікуваного значення залежної змінної DVдля кожної комбінації двох факторів, що може зробити це більш зрозумілим (Примітка. Я використовував звичайні оцінки, а не оцінки MCMC):

GroupAGroupBCondition16.13726.0758Conditiон26.25226.0853Сонгiтiон36.23726.1149

Я відповім на ваше запитання, відповідаючи на ваші тлумачення, посилаючись на цю таблицю.

Немає загальних відмінностей між групами (отже, група B має ap> 0,05)

pConditionCondition=16.13726.0758

Це не тестування, чи є загальна різниця між групами. Щоб зробити цей тест, вам доведеться повністю вийти Conditionз моделі і перевірити важливість Group.

Загальні відмінності між умовою 1 та умовою 2 та між умовою 1 та умовою 3.

Condition2Condition3Condition1Group=A6.1372Groupcondition

Відмінності між групою A, умовою 1 проти групи B, умовою 2, а також між групою A, умовою 1 проти групи B, умовою 3.

Умови взаємодії перевіряють, чи залежить вплив однієї змінної від рівня іншої змінної.

groupB:condition2Condition1Condition2Group=AGroup=B

6.25226.1372=.115
6.08536.0758=.0095
Condition2Condition1GroupA GroupB , і це, як я інтерпретувати це. Схоже, подібна динаміка відбувається, в меншій мірі, і стосовно Condition3.

1
Це фантастична відповідь: багато, велике спасибі за те, що знайшли час, щоб скласти це разом! На вашу думку, чи не буде мало сенсу застосовувати наступні контрасти щодо подібного?
vizzero

2
Вас дуже вітає @vizzero! У такому випадку схоже, що всі порівняння інтересів є у моделі, тому я не впевнений, якою була б мета спеціальних тестів. Крім того, оскільки ми бачимо значну взаємодію, важливість порівняння групових засобів (наприклад, група А проти групи В, ігнорування стану) для мене незрозуміла.
Макрос

Чудова відповідь, @Marco. Чи знаєте ви функцію, яка перевіряє загальний ефект усіх прогнозів, визначених у моделі, автоматично, не вказуючи та перевіряючи кожен підмодель вручну?
crsh

(x1,...,xp)y
E(y|x1,...,xp)=f(x1,...,xp)
f(x1,...,xp)cc
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.