Lmer और p-values ​​पर भ्रम: MCMC वालों की तुलना में मेमिस पैकेज से p-मान कैसे करते हैं?

Question

मैं इस धारणा के तहत था कि पैकेज lmer()में फ़ंक्शन lme4पी-वैल्यू (देखें lmer, पी-वैल्यू और सब कुछ ) का उत्पादन नहीं करता है ।

: मैं बजाय के रूप में इस सवाल का प्रति एमसीएमसी उत्पन्न पी मानों का उपयोग किया गया है में महत्वपूर्ण प्रभाव lme4मिश्रित मॉडल और इस सवाल: से उत्पादन में पी-मूल्यों नहीं खोजा जा सका lmer()में lm4में पैकेजR ।

हाल ही में मैंने एक पैकेज की कोशिश की जिसे मेमस्क कहा जाता है और यह getSummary.mer()अपने मॉडल के निश्चित प्रभावों को सीएसवी फ़ाइल में प्राप्त करने के लिए है। जैसे कि जादू से, ऐसा कॉलम pदिखाई देता है, जो मेरे MCMC पी-वैल्यू से बहुत निकटता से मेल खाता है (और उपयोग करने के साथ आने वाले प्रसंस्करण समय को नुकसान नहीं पहुँचाता है pvals.fnc())।

मैंने अस्थायी रूप से कोड में एक नज़र getSummary.merडाली है और पी-मान उत्पन्न करने वाली रेखा को देखा है:

p <- (1 - pnorm(abs(smry@coefs[, 3]))) * 2

क्या इसका मतलब यह है कि पी मानों को lmerचलाने के बजाय सीधे आउटपुट से उत्पन्न किया जा सकता है pvals.fnc? मुझे लगता है कि यह कोई संदेह नहीं है कि 'पी-वैल्यू बुतवाद' की बहस शुरू होगी, लेकिन मुझे जानने की दिलचस्पी है। मैंने पहले नहीं सुना है memiscजब यह आता है lmer।

अधिक सफल होने के लिए: एमसीएमसी पी-वैल्यू का उपयोग करने से क्या लाभ (यदि कोई है) उत्पन्न होने वालों पर getSummary.mer()?

Comments

मैं getSummary.merसमारोह का मूल लेखक हूं । सूचना -values केवल एक त्वरित जांच के रूप में इस्तेमाल किया जाना चाहिए। यदि मुझे याद है, तो मैंने वास्तव में इसे प्रदान किए गए ढांचे के भीतर काम करने के लिए केवल- शामिल किया है । लेकिन यह वास्तव में उपयोगकर्ता को एक उपयुक्त चेतावनी के साथ प्रदान किया जाना चाहिए, और मैं इसे बनाए रखने के बारे में देखने के लिए पैकेज अनुरक्षक से संपर्क करूंगा। मेरी सलाह का पालन करना है कि डग बेट्स द्वारा प्रदान किया गया: एमसीएमसी सुरक्षित शर्त है (यह मानते हुए कि दूसरों के पास बेहतर विकल्प नहीं हैं)। $p$$p$memisc

— जेसन मॉर्गन

@JasonMorgan मुझे इस सवाल का एक बहुत ही उचित जवाब लगता है।

— Glen_b -Reinstate Monica

@JasonMorgan मैं आपके कहे अनुसार सहमत हूं लेकिन वर्तमान mcmcsamp()में कई मुद्दों के कारण उपलब्ध नहीं है (एक Status of mcmcsampअधिक विवरण के लिए glmm.wikidot.com/faq में अनुभाग की जांच कर सकते हैं )। मुझे लगता है कि इस समय, शायद (पैरामीट्रिक?) बूटस्ट्रैपिंग एक व्यवहार्य-और लागू करने के लिए बहुत कठिन नहीं है- विकल्प; bootMer()functiom सेवा की जा सकती है।

— us --r11852

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

— Glen_b -Reinstate Monica

ओपी को एक अन्य तकनीकी / शब्दावली नोट। विधि का उपयोग करके प्राप्त memiscकिए गए पी-वैल्यू पीएल मान हैं जिन्हें वॉल्ड आंकड़ों के रूप में मनाया गया परीक्षण आँकड़ों से माना जाता है ( इस मामले में टी को वाल्ड जेड के रूप में माना जाता है )। ऐसा परीक्षण "बड़े नमूने" की धारणा पर निर्भर करता है और इसलिए अधिक से अधिक भरोसेमंद होता है क्योंकि आपके नमूना आकार बड़े होते हैं। मेरी जानकारी के लिए MCMC- आधारित मूल्य, इस तरह की धारणा पर निर्भर नहीं करता है। तो वैसे भी, Wald परीक्षणों और उनके विकल्पों के बारे में थोड़ा पढ़ना आपके प्रश्न पर और प्रकाश डालने में मदद कर सकता है।

— जेक वेस्टफॉल

Answer 1

मैं getSummary.merसमारोह का मूल लेखक हूं । सूचना -values केवल एक त्वरित जांच के रूप में इस्तेमाल किया जाना चाहिए। यदि मुझे याद है, तो मैंने वास्तव में इसे प्रदान किए गए ढांचे के भीतर काम करने के लिए केवल- शामिल किया है । लेकिन यह वास्तव में उपयोगकर्ता को एक उपयुक्त चेतावनी के साथ प्रदान किया जाना चाहिए, और मैं इसे बनाए रखने के बारे में देखने के लिए पैकेज अनुरक्षक से संपर्क करूंगा। मेरी सलाह का पालन करना है कि डग बेट्स द्वारा प्रदान किया गया: एमसीएमसी सुरक्षित शर्त है (यह मानते हुए कि दूसरों के पास बेहतर विकल्प नहीं हैं)।$p$$p$memisc