Tôi đang chạy một glmm với một biến phản ứng nhị thức và một công cụ dự đoán phân loại. Hiệu ứng ngẫu nhiên được đưa ra bởi thiết kế lồng nhau được sử dụng để thu thập dữ liệu. Dữ liệu trông như thế này:
m.gen1$treatment
[1] sucrose control protein control no_injection .....
Levels: no_injection control sucrose protein
m.gen1$emergence
[1] 1 0 0 1 0 1 1 1 1 1 1 0 0....
> m.gen1$nest
[1] 1 1 1 2 2 3 3 3 3 4 4 4 .....
Levels: 1 2 3 4 5 6 8 10 11 13 15 16 17 18 20 22 24
Mô hình đầu tiên tôi chạy trông như thế này
m.glmm.em.<-glmer(emergence~treatment + (1|nest),family=binomial,data=m.gen1)
Tôi nhận được hai cảnh báo giống như thế này:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.0240654 (tol = 0.001, component 4)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?
Tóm tắt mô hình cho thấy rằng một trong các phương pháp điều trị có lỗi tiêu chuẩn lớn bất thường, bạn có thể thấy ở đây:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 2.565 1.038 2.472 0.0134 *
treatmentcontrol -1.718 1.246 -1.378 0.1681
treatmentsucrose 16.863 2048.000 0.008 0.9934
treatmentprotein -1.718 1.246 -1.378 0.1681
Tôi đã thử các trình tối ưu hóa khác nhau từ điều khiển ánh sáng và các chức năng từ các gói khác và tôi nhận được một đầu ra tương tự. Tôi đã chạy mô hình bằng cách sử dụng glm bỏ qua hiệu ứng ngẫu nhiên và vấn đề vẫn còn. Trong khi khám phá dữ liệu tôi nhận ra rằng việc điều trị với Std cao. lỗi chỉ thành công trong biến phản ứng. Chỉ để kiểm tra xem điều đó có thể gây ra sự cố hay không, tôi đã thêm một điểm dữ liệu giả với "lỗi" cho việc xử lý đó và mô hình chạy trơn tru và đưa ra lỗi tiêu chuẩn hợp lý. Bạn có thể nhìn thấy cái đó ở kia:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.4090 1.6712 2.040 0.0414 *
treatmentcontrol -1.8405 1.4290 -1.288 0.1978
treatmentsucrose -0.2582 1.6263 -0.159 0.8738
treatmentprotein -2.6530 1.5904 -1.668 0.0953 .
Tôi đã tự hỏi nếu trực giác của tôi là đúng về việc thiếu thất bại cho điều trị đó ngăn ngừa một ước tính tốt, và làm thế nào tôi có thể làm việc xung quanh vấn đề này.
Cảm ơn trước!