Karışık modellerde zamanın bir yordayıcı olarak yer almasına izin veriliyor mu?


10

Her zaman zamanın regresyonlarda (gamlar dahil) bir yordayıcı olarak kullanılmaması gerektiğine inanıyordum, çünkü o zaman, kişi sadece trendin kendisini "tanımlayacaktır". Bir çalışmanın amacı, bir hayvanın aktivitesindeki varyansı açıklayan sıcaklık vb. Çevresel parametreleri bulmaksa, o zaman merak ediyorum, zaman nasıl faydalı olabilir? ölçülemeyen parametreler için bir proxy olarak?

Liman porpoises aktivite verileri üzerinde zaman içinde bazı eğilimler burada görülebilir: -> GAMM yaparken bir zaman serisi boşlukları nasıl ele alınır?

Benim sorunum: Modelime zaman eklediğimde (julian günlerde ölçülür), diğer tüm parametrelerin% 90'ı önemsiz hale gelir (mgcv'den daha pürüzsüz ts-shrinkle onları dışarı atın). Eğer zaman aşımı bırakırsam, bazıları önemlidir.

Soru şudur: Bir öngörücü olarak zamana izin verilir mi (hatta ihtiyaç duyulabilir mi?) Veya analizimi bozuyor mu?

şimdiden çok teşekkürler

Yanıtlar:


12

Zamana izin verilir; gerekip gerekmediği, modellemeye çalıştığınız şeye bağlı olacaktır? Sahip olduğunuz sorun, zamanın da yapabileceği, ancak daha az serbestlik derecesini kullanarak verilerdeki eğilime uygun görünen ortak değişkenlere sahip olmanızdır - bu nedenle Zaman yerine bırakılırlar.

Eğer ilgi sistemi modellemekse, tepkinin zaman içinde nasıl değiştiğini modellemek yerine, tepki ile ortak değişkenler arasındaki ilişki zaman içinde değişmez. Amaç, cevabın ortalama seviyesindeki değişikliği modellemekse, Zamanı dahil edin, ancak eşdeğeri dahil etmeyin. Ne dediğine itibaren, sen, eski değil ikincisi istediğiniz görünür ve gerektiği değil Modelinizdeki Süresi dahil. (Ancak aşağıdaki ekstra bilgileri göz önünde bulundurun.)

Gerçi birkaç uyarı var. Teorinin geçerli olması için artıklar iid olmalıdır (veya bir korelasyon yapısı kullanarak bağımsızlık varsayımını gevşetirseniz id). Cevabı ortak değişkenlerin bir fonksiyonu olarak modelliyorsanız ve verilerdeki herhangi bir eğilimi yeterince modellemiyorlarsa, kalan korelasyon yapısı bu eğilim ile başa çıkmadıkça, teorinin varsayımlarını ihlal eden bir eğilime sahip olacaktır.

Tersine, sadece yanıttaki eğilimi modelleyiyorsanız (sadece Zaman dahil), artıklarda (uygun eğilim hakkında) trendle (Zaman) açıklanmayan sistematik varyasyon olabilir ve bu da varsayımları ihlal edebilir kalanlar için. Bu gibi durumlarda, kalıntıları iid hale getirmek için başka ortak değişkenler eklemeniz gerekebilir

neden bu bir sorundur? Örneğin, trend bileşeninin önemli olup olmadığını veya ortak değişkenlerin etkilerinin önemli olup olmadığını test ederken, kullanılan teori kalıntıların iid olduğunu varsayarsa, eğer değilse, varsayımlar karşılanmaz ve p değerleri önyargılı olacaktır.

Tüm bunların amacı, verilerin tüm çeşitli bileşenlerini, artıkların kullandığınız teori için geçerli olacak, takılan bileşenlerin önemli olup olmadığını test edecek, geçerli olacak şekilde modellemeniz gerektiğidir.

Örnek olarak, mevsimsel verileri ele alalım ve verilerdeki uzun vadeli varyasyonu, eğilimi tanımlayan bir modele uymak istiyoruz. Eğilimi sadece mevsimsel döngüsel varyasyonu değil, modellersek, mevcut eğilimin önemli olup olmadığını test edemeyiz çünkü artıklar işe yaramayacaktır Bu tür veriler için hem mevsimsel bileşen hem de trend içeren bir modele uymamız gerekir bileşeni ve yalnızca sezonluk bileşeni içeren boş bir model. Daha sonra, takılan eğilimin önemini değerlendirmek için genelleştirilmiş bir olasılık oranı testi kullanarak iki modeli karşılaştırırız. Bu, kullanılan iki modelin bileşenleri anova()üzerinde yapılır .$lmegamm()


Sevgili Gavin, çok yararlı yorumlarınız için çok teşekkürler. Umarım yakında size de yardımcı olabilirim;) GLRT'yi anova ile denediğimde bana "nesne" sabit "bulunamadı ':(
Jens

1
@Jens çağrısı olmalı anova(mod1$lme, mod2$lme). Gauss olmayan bir model takıyorsanız, PQL yöntemlerinde gerçek bir günlük olasılığı olmadığı için bu işe yaramayabilir, PQL adında yarı olabilirlik. Bu, gamm4'ü kullanmanın bir nedenidir , ancak lme4'ün izin vermediği için korelasyon yapısı hakkında bir şeyler yapmanız gerekir .
Gavin Simpson
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.