एक मिश्रित मॉडल से प्रतिगमन आउटपुट की व्याख्या करते समय जब श्रेणीबद्ध चर के बीच बातचीत शामिल होती है


14

मेरे पास मिश्रित मॉडल / लैमर के मेरे उपयोग के बारे में एक प्रश्न है। मूल मॉडल यह है:

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

समूह और स्थिति दोनों कारक हैं: समूह के दो स्तर हैं (groupA, groupB) और स्थिति के तीन स्तर हैं (condition1, condition2, condition3)। यह मानव विषयों का डेटा है, इसलिए प्रत्येक व्यक्ति के लिए 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

अब, मुझे पता है कि सूचीबद्ध पंक्तियाँ कारकों के प्रत्येक स्तर को संदर्भ स्तर से तुलना करती हैं। समूह के लिए, संदर्भ GroupA है और स्थिति के लिए, संदर्भ कंडिशन 1 है।

क्या मैं निम्नलिखित तरीके से इस आउटपुट की व्याख्या करने में सही होगा:

  • समूहों के बीच कोई समग्र अंतर नहीं है (इसलिए groupB> ap05 का एपी)
  • हालत 1 और हालत 2 के बीच और हालत 1 और हालत 3 के बीच कुल अंतर।
  • GroupA, कंडिशन 1 बनाम ग्रुपबी, कंडीशन 2 और ग्रुपए, कंडिशन 1 बनाम ग्रुप बी, कंडीशन 3 के बीच अंतर।

क्या वो सही है? मुझे लगता है कि मैं इस बारे में थोड़ा भ्रमित हूं कि दो अलग-अलग कारकों के स्तरों के बीच बातचीत के संबंध में इसकी व्याख्या कैसे करें।

मैंने यहाँ विभिन्न प्रश्नों को पढ़ा है और कुछ वेब खोजों को किया है, और glht के साथ विरोधाभास प्राप्त करने में कामयाब रहे: क्या यह समूहों और स्थितियों के बीच के अंतर को देखने का एक बेहतर तरीका होगा? मुझे लगा कि इस मामले में बातचीत के संकेत यहां मौजूद होंगे।


लेकिन, अगर हम Group = B को संदर्भ स्तर Group = A से तुलना करना चाहते हैं तो Condition = 2 (या 3)? हो सकता? और, मुझे लगता है कि अगर "कंडीशन 1 और कंडीशन 2 के बीच का अंतर अलग है जब समूह = ए बनाम समूह = बी" की तुलना समान है अगर "समूह = ए और समूह = बी के बीच का अंतर तब अलग है जब कंडीशन 1 बनाम कंडीशन 2 "। क्या वो सही है? अन्यथा, उनके पी-मान क्या हैं?

यह प्रश्न के उत्तर की तरह नहीं दिखता है। बल्कि, आपके पास एक नया प्रश्न है। ऐसे पोस्ट करने के लिए सबसे अच्छा है।
निक कॉक्स

जवाबों:


21

दिए गए प्रतिगमन तालिका का उपयोग करते हुए, हम DVदो कारकों के प्रत्येक संयोजन के लिए , आश्रित चर के अपेक्षित मूल्य की तालिका की गणना कर सकते हैं, जो इसे और अधिक स्पष्ट कर सकता है (नोट मैंने साधारण अनुमानों का उपयोग किया है, एमसीएमसी अनुमानों का नहीं):

जीआरयूपीजीआरयूपीबीसीnमैंटीमैंn16.13726.0758सीnमैंटीमैंn26.25226.0853सीnमैंटीमैंn36.23726.1149

मैं इस तालिका का संदर्भ देते हुए आपकी व्याख्याओं का जवाब देकर आपके प्रश्न का उत्तर दूंगा।

समूहों के बीच कोई समग्र अंतर नहीं है (इसलिए groupB> ap05 का एपी)

पीConditionCondition=16.13726.0758

यह परीक्षण नहीं कर रहा है कि क्या समूहों के बीच एक समग्र अंतर है। उस परीक्षण को करने के लिए, आपको Conditionपूरी तरह से मॉडल से बाहर निकलना होगा और इसके महत्व का परीक्षण करना होगा Group

हालत 1 और हालत 2 के बीच और हालत 1 और हालत 3 के बीच कुल अंतर।

Condition2Condition3Condition1Group=A6.1372 । स्थिति चर में समग्र अंतर के लिए परीक्षण करने के लिए, आपको Groupमॉडल से बाहर निकलने और conditionअकेले परीक्षण करने की आवश्यकता होगी ।

GroupA, कंडिशन 1 बनाम ग्रुपबी, कंडीशन 2 और ग्रुपए, कंडिशन 1 बनाम ग्रुप बी, कंडीशन 3 के बीच अंतर।

बातचीत की शर्तें यह परखती हैं कि क्या एक चर का प्रभाव दूसरे चर के स्तर पर निर्भर करता है।

groupB:condition2Condition1Condition2Group=AGroup=B

6.2522-6.1372=0.115
6.0853-6.0758=0.0095
Condition2Condition1GroupA GroupBCondition3

1
यह एक शानदार प्रतिक्रिया है: बहुत से, इसे एक साथ रखने के लिए समय निकालने के लिए बहुत धन्यवाद! आपकी राय में, इस तरह से कुछ के लिए अनुवर्ती विरोधाभासों को चलाने में बहुत कम बिंदु होगा?
विज्जरो

2
आपका बहुत बहुत स्वागत है @vizzero! इस मामले में, ऐसा लगता है कि सभी हितों की तुलना मॉडल में है, इसलिए मुझे यकीन नहीं है कि पोस्ट-हॉक परीक्षणों का उद्देश्य क्या होगा। इसके अलावा, चूंकि हम महत्वपूर्ण बातचीत करते हैं, इसलिए समूह के साधनों की तुलना का महत्व (जैसे समूह ए बनाम समूह बी, स्थिति की अनदेखी) मेरे लिए स्पष्ट नहीं है।
मैक्रो

शानदार प्रतिक्रिया, @ मार्को। क्या आप किसी फ़ंक्शन को जानते हैं जो एक मॉडल में निर्दिष्ट सभी भविष्यवाणियों के समग्र प्रभाव का परीक्षण करता है और प्रत्येक सबमॉडल को हाथ से निर्दिष्ट किए बिना परीक्षण करता है?
crsh

(x1,...,xp)y
E(y|x1,...,xp)=f(x1,...,xp)
f(x1,...,xp)cसी
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.