Hiểu thông báo cảnh báo cho mô hình hỗn hợp trong r l44

Question

Tôi đã xây dựng và chạy một hiệu ứng hỗn hợp mô hình hồi quy logistic trong gói lme4 cho r để ước tính khả năng chiếm dụng của các loài cá ở các vị trí khác nhau (ô/môi trường sống). Khung dữ liệu bao gồm 1.207.140 quan sát của 68 cá thể. Đối với mỗi cá nhân (mỗi ngày trong ~ 1 năm), nó mô tả số lần xuất hiện tại mỗi vị trí duy nhất liên quan đến tổng số lần xuất hiện ở tất cả các địa điểm.

Dưới đây là mô hình cơ sở:

m.base = glmer(cbind(N,t.move-N) ~ jdate + snSurface.Area + Restoration..P.A. +  
    Release.Location+ Sex + (1|Station) + (0 + jdate|ID), data=allfishdat, family=binomial) 
where N=# unique positions, t.move=total positions, jdate=julian date, Station=locations, ID=fish ID 

Tôi nhận được thông báo cảnh báo sau đây:

Warning messages: 
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : 
       Model failed to converge with max|grad| = 3349.26 (tol = 0.001) 
2: In if (resHess$code != 0) { : 
the condition has length > 1 and only the first element will be used 
3: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : 
Model is nearly unidentifiable: very large eigenvalue 
- Rescale variables?;Model is nearly unidentifiable: large eigenvalue ratio 
- Rescale variables? 

tôi làm một số tìm kiếm để cố gắng hiểu những gì những thông điệp có ý nghĩa và tác động của chúng trên mô hình, nhưng chưa hiểu các cảnh báo.

Answer 1

Nếu vấn đề đầu tiên phải làm với R cần phải trải qua nhiều lần lặp lại để đạt được sự hội tụ, mã bên dưới có thể hữu ích. Thay thế '20000' bằng bất kỳ số lặp tối đa nào có ý nghĩa đối với mô hình cụ thể của bạn. (Chú ý rằng mã mô hình ban đầu của bạn đã được sửa đổi ở cuối để bao gồm 'kiểm soát = my.control'.)

my.control=lmerControl(optCtrl=list(maxfun=20000); my.control 

m.base = glmer(cbind(N,t.move-N) ~ jdate + snSurface.Area + Restoration..P.A. + Release.Location+ Sex + (1|Station) + (0 + jdate|ID), data=allfishdat, family=binomial, control = my.control) 

Nó cũng có thể hữu ích cho bạn để xem xét lmeControls hiện tại của bạn với lệnh sau đây:

str(lmerControl()) 

Hơn nữa, câu trả lời trước đây có thể hữu ích cho bạn: increase iterations for new version of lmer?