Пакет DLM в R

Я по сути скопировал пример Nile из пакета DLM в R, но не совсем понимаю, как он работает / как получаются определенные входные данные?

Кто-нибудь использовал указанный пакет, и если да, то могут ли они, возможно, подробнее рассказать о том, что делает каждый элемент из нижеприведенного, и какие соответствующие цифры будут получены?

 nileBuild <- function(par) {
  dlmModPoly(1, dV = exp(par[1]), dW = exp(par[2]))
}
nileMLE <- dlmMLE(Nile, rep(0,2), nileBuild); #nileMLE$conv
nileMod <- nileBuild(nileMLE$par)
V(nileMod)
W(nileMod)
nileFilt <- dlmFilter(Nile, nileMod)
nileSmooth <- dlmSmooth(nileFilt)
plot(cbind(Nile, nileFilt$m[-1], nileSmooth$s[-1]), plot.type="s",
     col=c("black","red","blue"), ylab="Level", main="Nile river", lwd=c(1,2,2))

s <- dlmSmooth(Nile, dlmModPoly(1, dV = 15100, dW = 1470))
plot(Nile, type ="s")
lines(dropFirst(s$s), col = "red")

Я особенно не понимаю строку nileBuild и то, что она делает... вместе с тем, почему 0 и 2 были выбраны в 'rep(0,2)'? Я прочитал документы, которые идут в комплекте, но я изо всех сил пытаюсь понять это.

Ура

0 ответов

Другие вопросы по тегам