Bazı eksik gözlemleri olan bir dizi iklim verisinde anormal değerleri tespit etmeye çalışıyorum. İnternette arama yaparken birçok mevcut yaklaşım buldum. Bunlardan stl ayrışması, eğilim ve mevsimsel bileşenleri ortadan kaldırmak ve geri kalanını incelemek anlamında çekici görünüyor. STL Okuma : Loess'e Dayalı Mevsimsel Trend Ayrıştırma Prosedürü , stl
değişkenlik atama ayarlarını belirleme konusunda, aykırı değerlerden etkilenmeyen ve eksik değerlere rağmen uygulanabilecek esnek görünmektedir. Ancak, R
dört yıllık gözlemlerle uygulamaya çalışarak ve tüm parametreleri http://stat.ethz.ch/R-manual/R-patched/library/stats/html/stl.html'ye göre tanımlayarak , karşılaşıyorum hata:
"time series contains internal NAs"
(ne zaman na.action=na.omit
) ve
"series is not periodic or has less than two periods"
(ne zaman na.action=na.exclude
).
Frekansın doğru tanımlandığını iki kez kontrol ettim. Bloglarda alakalı sorular gördüm, ancak bunu çözebilecek herhangi bir öneri bulamadım. stl
Eksik değerleri olan bir seride başvurmak mümkün değil mi? Artefaktları tanıtmak istemediğimden (ve sonuç olarak ... tespit etmek istemediğimden) bunları enterpolasyon konusunda çok isteksizim. Aynı nedenden ötürü, bunun yerine ARIMA yaklaşımlarını kullanmanın ne kadar tavsiye edilebilir olacağını bilmiyorum (ve eksik değerlerin hala bir sorun olup olmadığı).
stl
Eksik değerleri olan bir seriye başvurmanın bir yolunu biliyorsanız veya seçimlerimin metodolojik olarak sağlam olmadığını düşünüyorsanız veya daha iyi bir öneriniz varsa lütfen paylaşın . Alanda oldukça yeniyim ve ilgili bilginin yığınları tarafından boğulmuşum.