Tôi không biết liệu một mô hình hỗn hợp có phù hợp hay không (sử dụng các gói tiêu chuẩn trong đó cấu trúc hiệu ứng ngẫu nhiên là một công cụ dự báo tuyến tính), trừ khi bạn nghĩ rằng dữ liệu tại mọi thời điểm nên được trao đổi với nhau theo một nghĩa nào đó (trong trường hợp nào đó các khoảng không đều là một vấn đề) - nó thực sự sẽ không mô hình hóa sự tự tương quan thời gian một cách hợp lý. Có thể bạn có thể lừa lmer () thực hiện một số việc tự động nhưng chính xác là bạn sẽ làm thế nào để thoát khỏi tôi ngay bây giờ (tôi có thể không nghĩ thẳng). Ngoài ra, tôi không chắc chắn "biến nhóm" sẽ gây ra hiện tượng tự tương quan trong kịch bản mô hình hỗn hợp.
Nếu tự động tương quan thời gian là một tham số phiền toái và bạn không mong đợi nó cũng vậylớn, sau đó bạn có thể chia dữ liệu thành các kỷ nguyên về cơ bản tách rời nhau về mặt tương quan (ví dụ: tách chuỗi thời gian tại các điểm không có dữ liệu) và xem chúng là các bản sao độc lập. Sau đó, bạn có thể thực hiện một cái gì đó giống như một GEE trên tập dữ liệu đã sửa đổi này trong đó "cụm" được xác định theo thời gian bạn đang ở và các mục của ma trận tương quan làm việc là một chức năng của các quan sát được thực hiện cách nhau bao xa. Nếu hàm hồi quy của bạn là chính xác, thì bạn vẫn sẽ nhận được các ước tính nhất quán về các hệ số hồi quy, ngay cả khi cấu trúc tương quan bị sai. Điều này cũng sẽ cho phép bạn mô hình hóa nó dưới dạng dữ liệu đếm bằng cách sử dụng, ví dụ, liên kết nhật ký (như thường thấy trong hồi quy poisson). Bạn cũng có thể xây dựng trong một số mối tương quan khác biệt giữa các loài, trong đó mỗi điểm thời gian được xem như một vectơ đa biến của các loài được tính với một số liên kết phân rã theo thời gian giữa các điểm thời gian. Điều này sẽ yêu cầu một số tiền xử lý để lừa các gói GEE tiêu chuẩn thực hiện việc này.
Nếu tự động tương quan thời gian không phải là một tham số phiền toái, tôi sẽ thử một cái gì đó giống như một mô hình hiệp phương sai có cấu trúc trong đó bạn xem toàn bộ dữ liệu như một quan sát của một vectơ đa biến lớn sao cho hiệp phương sai giữa các quan sát trên các loài làYs,Ytu,v
cov(Ys,Yt)=fθ(s,t,u,v)
Trong đó là một số hàm tham số được biết đến với số lượng tham số hữu hạn, , cùng với một số tham số để chi phối cấu trúc trung bình. Bạn có thể cần phải "tự xây dựng" cho một mô hình như thế này, nhưng tôi cũng không ngạc nhiên nếu có các gói MPLUS để thực hiện những việc như thế này để đếm dữ liệu.fθ