Я пытаюсь предсказать следующие 24 месяца на основе следующих данных из набора данных по электроэнергии в пакете TSA. Я сделал модель, используя линейную регрессию, на основе которой и месяц, и время.
library(TSA)
library(forecast)
data(electricity)
logdata = log(electricity)
month = season(electricity)
time = time(electricity)
model = lm(logdata~month +time(electricity))
summary(model)
Однако мне трудно прогнозировать будущее на 24 месяца.
newelec= ts(rep(0,24), start=end(electricity), freq=frequency(electricity))
month = season(newelec)
mat = as.matrix(model.matrix(~month+time(newelec)+0)[,-1])
pred = forecast(model,mat)
Когда я делаю это, я получаю следующую ошибку:
Error in `[[<-.data.frame`(`*tmp*`, length(tmpdata) + 1, value = c(1973, : replacement has 396 rows, data has 24
Любая помощь, направившая меня в правильном направлении, была бы великолепна!