Вот образец моих данных, которые я использую.
Я сначала помещаю эти данные в линейную модель, а затем в dlmModReg модель (после расчета dV, dW, m0, C0).
model<-lm(Y1~X1[,1]+X1[,2])
model<-dlmModReg(X1,dV=a,dW=b,m0=c,C0=matrix(d,3,3,byrow=T))
где Y1 - «Разряд», X1 содержит «rain_102» и «rain_111». Затем я перехожу к фильтрации Калмана на этих данных
filtermod <- dlmFilter(Y1, model)
и перехожу к построению графика
plot(Y1, type = "o", col = "seagreen")
lines(dropFirst(filtermod$m), type = "o", pch = 20, col = "brown")
Как ожидается, прогноз не содержит информации о таймсерии. У меня вопрос: где данные временной серии «dt.tm» go, чтобы гарантировать, что я получаю правильные графики и прогнозы фильтрации Калмана.
Я попытался включить переменную «dt.tm» в Y1 и X1; однако я получаю ошибки.