Olá pessoal da lista, utilizei um pacote chamado mtsdi para preenchimento de dados em séries temporais. Gostaria de saber como posso proceder para fazer a calibração destes preenchimentos, por exemplo, eu tiro da minha série um período amostrado e com dados vizinhos correlacionados eu faço o preenchimento deste período assim eu comparo o observado com o preenchido, até aí tudo bem, mas eu gostaria de fazer uma calibração neste preenchimento e assim utilizá-la para os dados que realmente serão preenchidos. Segue um exemplo hipotético do uso do pacote.
library(mtsdi) library(hydroGOF) x0 <- rnorm(500) x1 <- rnorm(500) x2 <- rnorm(500) x3 <- rnorm(500) x4 <- c(x3[1:299],rep(NA,tim=400-299),x3[401:500]) x <- data.frame(x0,x1,x2,x4) f <- ~x0+x1+x2+x4 i <- mnimput(f,x,eps=1e-3,ts=T,log=F, method="spline",sp.control=list(df=c(7,7,7,7))) imput <- data.frame(predict(i)) (pred <- imput$x4[300:400]) (obs <- x3[300:400]) ggof(pred,obs)
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
