ऑफसेट के साथ GLM पॉज़िशन की भविष्यवाणी करें


9

मुझे पता है कि यह शायद एक बुनियादी सवाल है ... लेकिन मुझे इसका जवाब नहीं मिल रहा है।

मैं एक Poisson परिवार के साथ GLM फिट कर रहा हूं, और फिर भविष्यवाणियों पर एक नज़र डालने की कोशिश की, हालांकि ऑफसेट को ध्यान में रखा गया है:

model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003),
offset=(log(population)), data=data, subset=28:36, family=poisson())

predict (model_glm, type="response")

मुझे मामले मिलते हैं दर नहीं ...

मैंने भी कोशिश की है

model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+
offset(log(population)), data=data, subset=28:36, family=poisson())

उसी परिणाम के साथ। हालांकि जब मैं GAM से अनुमान लगाता हूं, mgcv का उपयोग करके, भविष्यवाणियां ऑफसेट (मुझे दर मिलती हैं) पर विचार करती हैं।

मुझे कुछ याद आ रहा है?


1
कृपया यहाँ और आर-हेल्प लिस्ट पर क्रॉस-पोस्ट न करें ... और यदि आप स्टैकओवरफ़्लो / स्टैकएक्सचेंज फोरम पर पोस्ट करने जा रहे हैं, तो मुझे लगता है कि एसओ बेहतर होगा (यह एक तकनीकी आर सवाल है, न कि एक आँकड़े सवाल ...)
बेन बोल्कर

जवाबों:


12

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


1
मुझे जवाब देने के लिए धन्यवाद। मुझे यह अनुमान लगाने में अजीब नहीं लगता कि मामलों की भविष्यवाणी करना, मुझे लगा कि दरों (मामलों / आबादी) के लिए भविष्यवाणी निर्धारित करने के लिए मैं कुछ याद कर रहा हूं। चूंकि GAM में मुझे इसके (मामलों / आबादी) की भविष्यवाणी करने के लिए कुछ और जोड़ना नहीं था।
सैंड्रा
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.