Xin chào các chuyên gia thống kê và trình hướng dẫn lập trình R,
Tôi quan tâm đến việc mô hình hóa việc bắt động vật như là một chức năng của điều kiện môi trường và ngày trong năm. Là một phần của một nghiên cứu khác, tôi có số lần bị bắt trong ~ 160 ngày trong ba năm. Vào mỗi ngày, tôi có nhiệt độ, lượng mưa, tốc độ gió, độ ẩm tương đối, v.v ... Vì dữ liệu được thu thập liên tục từ cùng 5 ô, tôi sử dụng âm mưu như một hiệu ứng ngẫu nhiên.
Hiểu biết của tôi là nlme có thể dễ dàng giải thích cho sự tự tương quan thời gian trong phần dư nhưng không xử lý các hàm liên kết không phải Gaussian như lme4 (không thể xử lý tự động tương quan?). Hiện tại, tôi nghĩ rằng nó có thể hoạt động để sử dụng gói nlme trong R trên log (tính). Vì vậy, giải pháp của tôi ngay bây giờ sẽ là chạy một cái gì đó như:
m1 <- lme(lcount ~ AirT + I(AirT^2) + RainAmt24 + I(RainAmt24^2) + RHpct + windspeed +
sin(2*pi/360*DOY) + cos(2*pi/360*DOY), random = ~1|plot, correlation =
corARMA(p = 1, q = 1, form = ~DOY|plot), data = Data)
trong đó DOY = Ngày trong năm. Có thể có nhiều tương tác hơn trong mô hình cuối cùng, nhưng đây là ý tưởng chung của tôi. Tôi cũng có khả năng cố gắng mô hình hóa cấu trúc phương sai hơn nữa với một cái gì đó như
weights = v1Pow
Tôi không chắc có cách nào tốt hơn để thực hiện với hồi quy mô hình hỗn hợp Poisson hay không? Tôi vừa tìm thấy thảo luận toán học trong Chương 4 của "Mô hình hồi quy cho phân tích chuỗi thời gian" của Kedem và Fokianos. Nó là một chút ngoài tôi vào lúc này, đặc biệt là trong ứng dụng (mã hóa bằng R). Tôi cũng thấy một giải pháp MCMC trong Zuur et al. Sách mô hình hiệu ứng hỗn hợp (Chp 23) bằng ngôn ngữ BUGS (sử dụng winBUGS hoặc JAG). Đó có phải là lựa chọn tốt nhất của tôi? Có một gói MCMC dễ dàng trong R sẽ xử lý việc này không? Tôi không thực sự quen thuộc với các kỹ thuật GAMM hoặc GEE nhưng sẽ sẵn sàng khám phá những khả năng này nếu mọi người nghĩ rằng họ sẽ cung cấp cái nhìn sâu sắc hơn.Mục tiêu chính của tôi là tạo ra một mô hình để dự đoán các động vật bị bắt trong điều kiện môi trường. Thứ hai, tôi muốn giải thích những gì động vật phản ứng về hoạt động của chúng.
Bất kỳ suy nghĩ nào về cách tốt nhất để tiến hành (về mặt triết học), cách viết mã này trong R hoặc trong BUGS sẽ được đánh giá cao. Tôi còn khá mới với R và BUGS (winBUGS) nhưng đang học. Đây cũng là lần đầu tiên tôi từng cố gắng giải quyết vấn đề tự kỷ theo thời gian.
Cảm ơn, Dan