Я использую набор данных flu
из пакета astsa
, и он содержит месячные значения измерений с 1968 по 1978.
Теперь я хотел бы иметь последовательность того же диапазона (12 * 11= 132), которая содержит именно эти моменты времени измерения, так что я смогу построить кривую лёсса с линиями, как в коде ниже.
К сожалению, красная кривая Лёсса абсолютно плоская, даже если она выглядит хорошо, если я строю ее отдельно (с x = measure_points).
Проблема, безусловно, заключается в точках измерения времени, я просто не могу их правильно использовать в качестве индекса для линий (предсказать Лесс).Как я могу решить это правильно с моим набором данных?
install.packages("astsa")
library(astsa)
data(flu)
tm <- seq(flu)
infpneum <- ts(flu,start=1968,freq=12)
res <- loess(infpneum ~ tm,span=0.1)
pred <- predict(res)
measurement_points <- seq(as.Date("1968-01-01"),as.Date("1978-12-31"),by="months")
plot(infpneum)
lines(pred,type="l",col="red",x=measurement_points)