Tôi đang chạy mô hình hồi quy logistic ở dạng:
lmer(response~1+(1|site), family=binomial, REML = FALSE)
Thông thường tôi sẽ tính toán ICC từ phương sai đánh chặn và dư, nhưng tóm tắt của mô hình không bao gồm phương sai dư. Làm thế nào để tôi tính toán này?
1
Tại sao bạn tính toán ICC?
—
AdamO
Để kiểm tra giả định rằng hồi quy logistic thông thường không hợp lệ đối với các dữ liệu này, là bằng chứng cho thấy tôi nên sử dụng GLMM. Tôi đã tìm thấy một phương trình: ICClogit = phương sai đánh chặn ^ 2 / (phương sai đánh chặn ^ 2 + pi ^ 2/3). Điều này có vẻ hợp lý?
—
Megan
Bạn đang sử dụng phương pháp khả năng tối đa đầy đủ. Bạn có thể làm một bài kiểm tra tỷ lệ khả năng với 1 bậc tự do so với mô hình hiệu ứng cố định không?
—
AdamO
@Megan: Bạn nói đúng. Trong thực tế, Zeger et al. (1988) gợi ý công trình tốt hơn so với π 2 / 3 là đúng còn lại cho các mô hình hồi quy logistic, mặc dù hai đang rất gần. Xem SL Zeger, KY Liang và PS Albert. Các mô hình cho dữ liệu theo chiều dọc: một cách tiếp cận phương trình ước lượng tổng quát. Sinh trắc học, 44: 1049-1060 1988.
—
Randel
@Megan: Đó là
—
Wolfgang
intercept_variance / (intercept_variance + pi^2/3)
- vì vậy đừng bình phương phương sai.