Tôi đang xem lại một bài báo về sự thụ phấn, trong đó dữ liệu được phân phối nhị thức (quả chín hoặc không). Vì vậy, tôi đã sử dụng glmer
với một hiệu ứng ngẫu nhiên (cây riêng lẻ) và một hiệu ứng cố định (điều trị). Một nhà phê bình muốn biết liệu thực vật có ảnh hưởng đến bộ quả hay không - nhưng tôi gặp khó khăn trong việc diễn giải glmer
kết quả.
Tôi đã đọc trên web và dường như có thể có vấn đề với việc so sánh trực tiếp glm
và glmer
mô hình, vì vậy tôi không làm điều đó. Tôi đã tìm ra cách đơn giản nhất để trả lời câu hỏi sẽ là so sánh phương sai hiệu ứng ngẫu nhiên (1.449, bên dưới) với tổng phương sai hoặc phương sai được giải thích bằng cách xử lý. Nhưng làm thế nào để tôi tính toán các phương sai khác? Chúng dường như không được bao gồm trong đầu ra dưới đây. Tôi đã đọc một cái gì đó về phương sai dư không được đưa vào nhị thức glmer
- làm thế nào để tôi giải thích tầm quan trọng tương đối của hiệu ứng ngẫu nhiên?
> summary(exclusionM_stem)
Generalized linear mixed model fit by maximum likelihood (Laplace
Approximation) [glmerMod]
Family: binomial ( logit )
Formula: cbind(Fruit_1, Fruit_0) ~ Treatment + (1 | PlantID)
AIC BIC logLik deviance df.resid
125.9 131.5 -59.0 117.9 26
Scaled residuals:
Min 1Q Median 3Q Max
-2.0793 -0.8021 -0.0603 0.6544 1.9216
Random effects:
Groups Name Variance Std.Dev.
PlantID (Intercept) 1.449 1.204
Number of obs: 30, groups: PlantID, 10
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.5480 0.4623 -1.185 0.2359
TreatmentD -1.1838 0.3811 -3.106 0.0019 **
TreatmentN -0.3555 0.3313 -1.073 0.2832
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
(Intr) TrtmnD
TreatmentD -0.338
TreatmentN -0.399 0.509