मैं आर पैकेज glmerमें एक मॉडल फिट कर रहा हूं lme4। मैं पी-मान के साथ एक एनोवा तालिका की तलाश में हूं, लेकिन मुझे ऐसा कोई भी पैकेज नहीं मिल रहा है जो इसे फिट करे। क्या आर में ऐसा करना संभव है?
मैं जिस मॉडल की फिटिंग कर रहा हूं वह फॉर्म का है:
model1<-glmer(dmn~period*teethTreated+(1|fullName),
family="poisson",
data=subset(dataset,
group=='Four times a year'),
control=glmerControl(optimizer="bobyqa"))जवाबों:
यदि आप वाल्ड परीक्षणों के लिए समझौता करना चाहते हैं तो यह काम करना चाहिए:
library(lme4)
library(car)
gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp, family = binomial)
Anova(gm1,type="III")हालाँकि, ध्यान दें ?Anovaकि (से ):
पदनाम "टाइप- II" और "टाइप- III" एसएएस से उधार लिए गए हैं, लेकिन यहां उपयोग की जाने वाली परिभाषाएं एसएएस द्वारा नियोजित लोगों के लिए सटीक रूप से मेल नहीं खाती हैं। टाइप- II परीक्षण की गणना सीमांतता के सिद्धांत के अनुसार की जाती है, शब्द के उच्च-क्रम वाले रिश्तेदारों की अनदेखी के अलावा, अन्य सभी के बाद प्रत्येक शब्द का परीक्षण; तथाकथित प्रकार- III परीक्षण सीमांतता का उल्लंघन करते हैं, अन्य सभी के बाद मॉडल में प्रत्येक शब्द का परीक्षण करते हैं। टाइप- II परीक्षणों की यह परिभाषा एसएएस द्वारा विश्लेषण-के-विचरण मॉडल के लिए निर्मित परीक्षणों से मेल खाती है, जहां सभी भविष्यवक्ता कारक हैं, लेकिन आम तौर पर अधिक नहीं (यानी, जब मात्रात्मक भविष्यवक्ता होते हैं)। टाइप- III परीक्षण के लिए मॉडल तैयार करने में बहुत सावधानी बरतें, या परीक्षण की गई परिकल्पना का कोई मतलब नहीं होगा।
मैं यह सुनिश्चित करने के लिए आपके परिणामों की बहुत सावधानी से जांच करूँगा कि वे क्या समझ रहे हैं!
वैकल्पिक रूप से, आप afex::mixedसंभावना अनुपात परीक्षण या पैरामीट्रिक बूटस्ट्रैप के माध्यम से अनुरूप तालिकाओं को प्राप्त करने के लिए उपयोग कर सकते हैं ; उत्तरार्द्ध सबसे सटीक है, लेकिन अब तक सबसे धीमा है।
जीएलएमएम के संदर्भ में पी-मूल्य संगणना की अधिक सामान्य चर्चा के लिए पैकेज ?pvaluesमें देखें lme4।