在時間序列上應用黃土平滑

Question

我想通過應用黃土函數來平滑時間曲線，但我無法使其工作。 一個例子：

mydat <- runif(50) 
day1 <- as.POSIXct("2012-07-13", tz = "UTC") 
day2 <- day1 + 49*3600*24 
pdays <- seq(day1, day2, by = "days") 
lo <- loess(mydat ~ pdays) 

我得到以下信息：

Error: NA/NaN/Inf in foreign function call (arg 2) 

是否有可能應用黃土平滑的時間序列

任何幫助或指導，不勝感激！

Answer 1

我認爲這裏的想法是以數字形式（使用as.numeric）轉換您的時間序列，以便您可以執行操作。

mydat <- runif(50) 
day1 <- as.POSIXct("2012-07-13", tz = "UTC") 
day2 <- day1 + 49*3600*24 
pdays <- seq(day1, day2, by = "days") 
lo <- loess(mydat ~ as.numeric(pdays)) 

# And then if you want to plot the result: 
plot(pdays,mydat) 
lines(pdays, lo$fitted)