আমি আর প্যাকেজে একটি 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) যে:
"টাইপ -2" এবং "টাইপ-III" উপাধিগুলি এসএএস থেকে ধার করা হয়েছে, তবে এখানে ব্যবহৃত সংজ্ঞাগুলি এসএএস দ্বারা নিযুক্তদের সাথে যথাযথভাবে মেলেনি। টাইপ -২ পরীক্ষাগুলি মার্জিনালটির নীতি অনুসারে গণ্য করা হয়, প্রতিটি পদটির উচ্চতর অর্ডার আত্মীয়দের অগ্রাহ্য করা ছাড়া অন্য সকলের পরে পরীক্ষা করে; তথাকথিত ধরণের III পরীক্ষাগুলি প্রান্তিকতা লঙ্ঘন করে, প্রতিটি ক্ষেত্রেই অন্যদের সমস্তের পরে মডেলটিতে পরীক্ষা করে। প্রকার -২ পরীক্ষার এই সংজ্ঞাটি এসএএস-এর-বিশ্লেষণ-বৈকল্পিক মডেলগুলির জন্য উত্পাদিত পরীক্ষার সাথে মিলে যায়, যেখানে ভবিষ্যদ্বাণীকারীরা সমস্ত কারণ হয়, তবে সাধারণত হয় না (অর্থাত্ যখন পরিমাণগত ভবিষ্যদ্বাণী থাকে)। প্রকার III পরীক্ষার জন্য মডেল গঠনে খুব সতর্কতা অবলম্বন করুন বা পরীক্ষিত হাইপোথিসগুলি অর্থবোধ করবে না।
আপনার ফলাফলগুলি বুদ্ধিমান হওয়ার জন্য আমি খুব যত্ন সহকারে আপনার ফলাফলগুলি পরীক্ষা করব!
বিকল্পভাবে, আপনি afex::mixedসম্ভাবনা অনুপাত পরীক্ষা বা প্যারামেট্রিক বুটস্ট্র্যাপের মাধ্যমে অ্যালবাম টেবিলগুলি পেতে ব্যবহার করতে পারেন ; পরবর্তীটি সর্বাধিক নির্ভুল তবে দূর থেকেও ধীর।
দেখুন ?pvaluesএ lme4GLMMs প্রেক্ষাপটে P-মান গুনতি আরও সাধারণ আলোচনার জন্য প্যাকেজ।