Olá Éder, entendi até aí eu me dou bem, o problema é que queria fazer uma calibração nesse preenchimento e usar essa calibração para os dados que realmente serão preenchidos, tentei criar uma função de aproximação approxfun(sim,obs), mas não obtive resultados melhores, atualmente estou estudando o hydroPSO para ver se consigo algo. Obrigado pela ajuda Abraço Em 4 de maio de 2015 18:13, Éder Comunello <comunello.eder@gmail.com> escreveu:
Caro Wagner, boa tarde!
Não tenho experiência com esse pacote, mas tentei aplicar o procedimento que utilizo pra avaliar o preenchimento em outras situações. O dataset é restrito, mas deve servir pra testar...
### <code r> require(mtsdi) data(miss) notNA <- which(!is.na(miss), arr.ind=T) ### valores presentes set.seed(333); sel.pos <- notNA[sample(nrow(notNA), 20),] ### reservando 20 valores sel.obs <- miss[sel.pos] ### valores reservados
miss2 <- miss miss2[sel.pos] <- NA ### "exclui" valores reservados f <- ~c31+c32+c33+c34+c35 i <- mnimput(f,miss2,eps=1e-3,ts=TRUE, method="spline",sp.control=list(df=c(7,7,7,7,7))) summary(i) imput <- data.frame(predict(i)) sel.pre <- imput[sel.pos] ### predição referente aos 20 valores reservados
cbind(sel.pre, sel.obs) ### comparação library(hydroGOF) ggof(sel.pre,sel.obs) ### </code>
[image: Imagem inline 1]
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
