जब परिणाम भिन्नात्मक (दो संख्याओं का अनुपात) हो तो R में लॉजिस्टिक रिग्रेशन कैसे करें?


24

मैं एक कागज की समीक्षा कर रहा हूं जिसमें निम्नलिखित जैविक प्रयोग हैं। एक उपकरण का उपयोग द्रव कतरनी तनाव की अलग-अलग मात्रा में कोशिकाओं को उजागर करने के लिए किया जाता है। जैसा कि अधिक कतरनी तनाव कोशिकाओं पर लागू होता है, उनमें से अधिक सब्सट्रेट से अलग होने लगते हैं। कतरनी तनाव के प्रत्येक स्तर पर, वे उन कोशिकाओं को गिनते हैं जो संलग्न रहती हैं, और चूंकि उन्हें पता है कि कोशिकाओं की कुल संख्या जो शुरुआत में जुड़ी हुई थी, वे एक भिन्नात्मक लगाव (या टुकड़ी) की गणना कर सकते हैं।

यदि आप पक्षपाती अंश बनाम कतरनी तनाव की साजिश करते हैं, तो परिणाम एक तार्किक वक्र है। सिद्धांत रूप में, प्रत्येक व्यक्तिगत सेल एक एकल अवलोकन है, लेकिन जाहिर है कि हजारों या दसियों हजार कोशिकाएं होती हैं, इसलिए डेटा सेट विशाल होगा, अगर इसे सामान्य तरीके से स्थापित किया गया था (प्रत्येक पंक्ति एक अवलोकन के साथ)।

इसलिए, स्वाभाविक रूप से, मेरा प्रश्न (जैसा कि शीर्षक में कहा गया है) अब समझ में आना चाहिए। हम DV के रूप में भिन्नात्मक परिणाम का उपयोग करके लॉजिस्टिक रिग्रेशन कैसे करते हैं? क्या कुछ स्वचालित परिवर्तन है जो glm में किए जा सकते हैं?

इसी तर्ज पर, यदि संभावित 3 या अधिक (भिन्नात्मक) माप थे, तो कोई बहुराष्ट्रीय लॉजिस्टिक रिग्रेशन के लिए ऐसा कैसे करेगा?


यहाँ एक बहुभिन्नरूपी परिणाम लॉजिस्टिक रिग्रेशन के बारे में कुछ उदाहरण दिए गए हैं:http://www.ats.ucla.edu/stat/r/dae/mlogit.htm
marbel

1
आप जो वर्णन करते हैं वह ध्वनि नहीं करता है जैसे आपके पास स्वतंत्र अवलोकन होंगे (चूंकि कोशिकाएं संभवतः नहीं बढ़ती हैं क्योंकि कतरनी बढ़ जाती है, फिर भी प्रत्येक उच्च तनाव सेटिंग में संलग्न संख्या पिछले संख्या से अधिक नहीं होनी चाहिए); इस निर्भरता को ध्यान में रखा जाना चाहिए। (यह विकास घटता के साथ स्थिति की याद दिलाता है।) --- आप बस संख्याओं को एक GLM में प्लग नहीं कर सकते हैं जैसे कि वे स्वतंत्र थे ... और कोई भी उत्तर इस मुद्दे को संबोधित नहीं करता है।
Glen_b -Reinstate मोनिका

2
@Glen_b प्रत्येक प्रयोग विभिन्न कोशिकाओं के साथ किया जाएगा, अर्थात "100%" लगाव और कतरनी तनाव के एक अलग मूल्य को लागू करने के साथ।
Thecity2

आह ठीक है। वह स्वतंत्र परिणाम देगा।
Glen_b -Reinstate Monica

जवाबों:


37

glmमें समारोह Rएक रसद प्रतिगमन मॉडल के लिए सूत्र निर्दिष्ट करने के लिए 3 तरीके अनुमति देता है।

सबसे आम है कि डेटा फ़्रेम की प्रत्येक पंक्ति एक एकल अवलोकन का प्रतिनिधित्व करती है और प्रतिक्रिया चर या तो 0 या 1 है (या 2 स्तरों वाला कारक, या केवल 2 अद्वितीय मानों के साथ अन्य चर)।

एक अन्य विकल्प 2 कॉलम मैट्रिक्स का उपयोग करने के लिए प्रतिक्रिया चर के रूप में है जिसमें पहला कॉलम 'सफलताओं' की गिनती और दूसरा स्तंभ 'विफलताओं' की गिनती है।

आप प्रतिक्रिया को 0 और 1 के बीच के अनुपात के रूप में भी निर्दिष्ट कर सकते हैं, फिर एक अन्य कॉलम को 'वजन' के रूप में निर्दिष्ट करें जो कुल संख्या देता है जो अनुपात से है (इसलिए 0.3 की प्रतिक्रिया और 10 का वजन 3 के समान है) सफलताओं 'और 7' विफलताओं ')।

या तो अंतिम 2 तरीके आप क्या करने की कोशिश कर रहे हैं फिट होगा, आखिरी आपके डेटा का वर्णन करने के लिए सबसे प्रत्यक्ष लगता है।


9

एक शुरुआत के रूप में, यदि आपके पास एक आश्रित चर है जो एक अनुपात है, तो आप बीटा प्रतिगमन का उपयोग कर सकते हैं। यह (मेरे सीमित ज्ञान के साथ) कई अनुपातों तक विस्तारित नहीं होता है।

बीटा रिग्रेशन अवलोकन और आर कार्यान्वयन के लिए बिटारेग की जाँच करें ।


धन्यवाद! ऐसा लगता है कि मुझे द्विपद मामले के लिए क्या चाहिए।
thecity2

2

मैं nnet::multinomएक समान उद्देश्य के लिए उपयोग कर रहा हूँ (पैकेज nnet MASS का हिस्सा है), यह [0, 1] में निरंतर इनपुट को स्वीकार करता है।

यदि आपको एक संदर्भ की आवश्यकता है: C. Beleites et.al . : एस्ट्रोसाइटोमा ऊतकों के रमन स्पेक्ट्रोस्कोपिक ग्रेडिंग: सॉफ्ट संदर्भ जानकारी का उपयोग करना। गुदा बायोएनल केम, 2011, वॉल्यूम। 400 (9), पीपी। 2801-2816


महान! मेरे पास वह पैकेज है, और यह महसूस नहीं किया कि उसमें यह क्षमता थी।
thecity2

@ कबीलेइट्स: क्या यह आश्रित को [0,1] की अनुमति देता है? हालांकि मैं एक मामूली आश्रित के लिए एक समारोह था (भविष्यवक्ताओं को [0,1] तक बढ़ाया जाना चाहिए ...
B_Miner

@B_Miner: हाँ, आश्रित [0, 1] में हो सकता है। फ़ंक्शन छिपी हुई परत के बिना और लाजिस्टिक सिग्मॉइड के साथ एक कृत्रिम तंत्रिका नेटवर्क को फिट करता है। और हां, बेहतर अनुमान के लिए भविष्यवक्ताओं को मोटे तौर पर [0, 1] के पैमाने पर रखने की सिफारिश की जाती है।
cbeleites
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.