5

レス関数を適用して、プロットした時間曲線を滑らかにしたいのですが、動作させることができません。例:

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)

時系列にレス平滑化を適用することは可能ですか?

どんな助けやガイダンスも大歓迎です!

4

1 に答える 1

15

ここでの考え方は、時系列を(を使用して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)
于 2012-07-30T09:09:19.300 に答える