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 <[email protected]> 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 <[email protected]>[email protected]> > Dourados, MS - [22 16.5'S, 54 49'W] > > > > _______________________________________________ > 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. >
_______________________________________________ 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.
