Я пытаюсь создать график для данных временных рядов. Красная линия изображает мою модель сегментированной регрессии. Чтобы изменение уровня на графике выглядело более отчетливо, я создал новый фрейм данных с 0,1 единицами времени, чтобы предсказать синюю линию. Вместо 24 временных точек у меня теперь 240.
model <- lm(y ~ time + cons + tsi)
datanew <- data.frame(cons=rep(c(0, 1), c(120, 120)), time=1:240/10, tsi=0)
datanew$tsi[121:240] <- 1:120/10
pred1 <- predict(model, type="response", datanew)
lines((1:240/10), pred1, col="4")
Красная / синяя линия изображает модель без контроля автокорреляции. Зеленая линия - это модель с контролем автокорреляции:
Для зеленой линии я хотел бы сделать то же самое, что и длякрасная линия. Однако, поскольку я контролировал автокорреляцию с использованием функции prais_winsten
prais <- prais_winsten(y ~ time + cons + tsi,
data=ARtsi)
с использованием функции прогнозирования, по-видимому, это невозможно, поскольку я получаю следующее сообщение об ошибке:
«Ошибка в UseMethod (« предикат »): нет применимого метода для« предиката », примененного к объекту класса« похвалы »» *
Кто-нибудь знает лучший способ улучшить мой график при использовании prais_winsten функция?