Merhaba istatistiksel gurular ve R programlama sihirbazları,
Hayvan yakalamalarını çevresel koşulların ve yılın gününün bir fonksiyonu olarak modellemekle ilgileniyorum. Başka bir çalışmanın parçası olarak, üç yıl içinde ~ 160 gün boyunca yakalama sayım var. Bu günlerin her birinde sıcaklık, yağış, rüzgar hızı, bağıl nem vb. Var. Veriler aynı 5 parselden tekrar tekrar toplandığından rastgele bir etki olarak arsa kullanıyorum.
Anladığım kadarıyla, nlme artıklarda geçici otokorelasyonu kolayca hesaplayabilir, ancak lme4 (otokorelasyonu işleyemeyen?) Gibi Gauss dışı bağlantı işlevlerini işlemez. Şu anda, günlükte (sayım) R'de nlme paketini kullanmanın işe yarayabileceğini düşünüyorum. Yani benim çözüm şu anda şöyle bir şey çalıştırmak olacaktır:
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)
burada DOY = Yılın Günü. Son modelde daha fazla etkileşim olabilir, ancak bu benim genel fikrim. Ayrıca, varyans yapısını benzer bir şeyle daha fazla modellemeye çalışabilirim
weights = v1Pow
Poisson karma model regresyonu veya başka bir şeyle daha iyi bir yol olup olmadığından emin değilim? Kedem ve Fokianos'un "Zaman Serisi Analizi için Regresyon Modelleri" Bölüm 4'te matematiksel bir tartışma buldum. Şu anda, özellikle uygulamada (R'de kodlama) biraz ötesindeydi. Ayrıca Zuur ve ark. BUGS dilinde Karışık Efekt Modelleri kitabı (Bölüm 23) (winBUGS veya JAG kullanarak). Bu benim en iyi seçeneğim mi? R'de bunun üstesinden gelebilecek kolay bir MCMC paketi var mı? GAMM veya GEE tekniklerine gerçekten aşina değilim, ancak insanlar daha iyi kavrayış sağlayacaklarını düşünürlerse bu olasılıkları araştırmaya istekli olacaklardı.Temel hedefim, çevre koşulları göz önüne alındığında hayvan yakalamalarını tahmin etmek için bir model oluşturmak. İkinci olarak, hayvanların aktiviteleri açısından neye yanıt verdiğini açıklamak istiyorum.
İlerlemenin en iyi yolu hakkındaki düşünceler (felsefi olarak), bunu R veya HATA'da nasıl kodlayacağınız takdir edilecektir. R ve BUGS (winBUGS) için oldukça yeniyim ama öğreniyorum. Bu aynı zamanda geçici otokorelasyonu ilk kez ele almaya çalıştım.
Teşekkürler, Dan