Saya memasang glmermodel dalam lme4paket R. Saya mencari tabel anova dengan nilai p yang ditunjukkan di sini, tetapi saya tidak dapat menemukan paket yang cocok. Apakah mungkin melakukannya di R?
Model yang saya pas adalah dalam bentuk:
model1<-glmer(dmn~period*teethTreated+(1|fullName),
family="poisson",
data=subset(dataset,
group=='Four times a year'),
control=glmerControl(optimizer="bobyqa"))
Jawaban:
Jika Anda bersedia menerima tes Wald ini harusnya berhasil:
library(lme4)
library(car)
gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp, family = binomial)
Anova(gm1,type="III")Namun, perhatikan (dari ?Anova) bahwa:
Sebutan "tipe-II" dan "tipe-III" dipinjam dari SAS, tetapi definisi yang digunakan di sini tidak sesuai persis dengan yang digunakan oleh SAS. Tes Tipe-II dihitung sesuai dengan prinsip marginalitas, menguji setiap istilah setelah yang lainnya, kecuali mengabaikan kerabat tingkat tinggi istilah; yang disebut tes tipe-III melanggar marginalitas, menguji setiap istilah dalam model setelah semua yang lain. Definisi tes Tipe-II ini sesuai dengan tes yang diproduksi oleh SAS untuk model analisis varians, di mana semua prediktor adalah faktor, tetapi tidak lebih umum (yaitu, ketika ada prediktor kuantitatif). Berhati-hatilah dalam merumuskan model untuk tes tipe-III, atau hipotesis yang diuji tidak masuk akal.
Saya akan memeriksa hasil Anda dengan sangat hati-hati untuk memastikan hasilnya masuk akal!
Atau, Anda dapat menggunakan afex::mixeduntuk mendapatkan tabel analog melalui uji rasio kemungkinan atau bootstrap parametrik; yang terakhir adalah yang paling akurat, tetapi juga yang paling lambat sejauh ini.
Lihat ?pvaluesdalam lme4paket untuk diskusi yang lebih umum tentang perhitungan nilai-p dalam konteks GLMM.