r - Is there any way to fit a `glm()` so that all levels are included (i.e. no reference level)? -
कोड पर विचार करें:
x & lt; - read.table ("http: / [/ Data.princeton.edu/wws509/datasets/cuse.dat ", हैडर = TRUE] [, 1: 2] फिट & lt; - glm (शिक्षा ~ उम्र, परिवार =" binomial ", डेटा = x) सारांश (फिट)
जहां आयु के चार स्तर हैं: "& lt; 25" "25-29" "30-39" "40-49"
परिणाम निम्न हैं:
तो डिफ़ॉल्ट रूप से, स्तरों में से एक एक संदर्भ स्तर के रूप में प्रयोग किया जाता है क्या सभी 4 स्तर + अवरोधन (अर्थात कोई संदर्भ स्तर नहीं है) के लिए glm आउटपुट गुणांक है? एसएएस जैसे सॉफ्टवेयर पैकेज डिफ़ॉल्ट रूप से करते हैं, इसलिए मैं सोच रहा था कि इसके लिए कोई विकल्प था।
धन्यवाद!
सूत्र
, विशेष रूप से, + 0
को अपने मॉडल विनिर्देशन में शामिल करने का अर्थ देखें ...
< कोड> # नमूना डेटा - व्याख्यात्मक चर (सतत) x & lt; - runif (100) # स्पष्टीकरणीय डेटा, 3 स्तरों के साथ कारक एफ & lt; - as.factor (नमूना (3, 100, TRUE)) # परिणाम डेटा y & lt; # इंटरफ़ेस सारांश (जीएलएम (y ~ x + f + 0)) के साथ # मॉडल: (# 100) + rnorm (100) + rnorm (100, माध्य = सी (1,3,6)) # कॉल: # = Y ~ x + f + 0) # डायवेंविस अवशिष्ट: # मिन 1 क्यू मेडियान 3 क्यू मैक्स # -5.7316 -1.8 923 0.0195 1.8918 5.9520 # क्यूफीएफ़ेंटर: # अनुमान स्टैट त्रुटि टी मूल्य पीआर (& gt; | टी |) #x 0.3216 0.9772 0.329 0.743 # एफ 1 3.4493 0.6823 5.055 2.06e-06 *** # एफ 2 3.634 9 0.6959 5.223 1.02e-06 *** # एफ 3 3.1 9 62 0.6598 4.844 4.87 ई-06 ***
Comments
Post a Comment