我試圖實現Chebyshev濾波器來平滑時間序列,但不幸的是,在數據序列中有NAs。R濾波器()處理NA
例如,
t <- seq(0, 1, len = 100)
x <- c(sin(2*pi*t*2.3) + 0.25*rnorm(length(t)),NA, cos(2*pi*t*2.3) + 0.25*rnorm(length(t)))
我用的切比雪夫濾波器:cf1 = cheby1(5, 3, 1/44, type = "low")
我試圖過濾器系列排除NAS上的時間,但不會弄亂訂單/位置。所以,我已經試過na.rm=T
,但似乎沒有這樣的說法。 Then
z <- filter(cf1, x) # apply filter
謝謝你們。
我只是想知道complete.case和na.omit是否一樣。另外,由於我使用的是觀測到的SST時間序列,我不確定輸入缺失值是否是一個好主意。 – 2012-07-18 13:15:40
希望這個更新解決了這個問題。 – chandler 2012-07-19 11:37:15