<div dir="ltr"><div><div><div>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.<br></div><br></div>Obrigado pela ajuda <br></div>Abraço<br></div><div class="gmail_extra"><br><div class="gmail_quote">Em 4 de maio de 2015 18:13, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Caro Wagner, boa tarde!</div><div><br></div><div>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...</div><div><br></div><div><div><font face="monospace, monospace">### <code r></font></div><div><font face="monospace, monospace">require(mtsdi)</font></div><div><font face="monospace, monospace">data(miss)</font></div><div><font face="monospace, monospace">notNA   <- which(!<a href="http://is.na" target="_blank">is.na</a>(miss), arr.ind=T)  ### valores presentes</font></div><div><font face="monospace, monospace">set.seed(333); sel.pos <- notNA[sample(nrow(notNA), 20),] ### reservando 20 valores</font></div><div><font face="monospace, monospace">sel.obs <- miss[sel.pos]                   ### valores reservados</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">miss2 <- miss</font></div><div><font face="monospace, monospace">miss2[sel.pos] <- NA                       ### "exclui" valores reservados</font></div><div><font face="monospace, monospace">f <- ~c31+c32+c33+c34+c35</font></div><div><font face="monospace, monospace">i <- mnimput(f,miss2,eps=1e-3,ts=TRUE, method="spline",sp.control=list(df=c(7,7,7,7,7)))</font></div><div><font face="monospace, monospace">summary(i)</font></div><div><font face="monospace, monospace">imput <- data.frame(predict(i))</font></div><div><font face="monospace, monospace">sel.pre <- imput[sel.pos]                  ### predição referente aos 20 valores reservados</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">cbind(sel.pre, sel.obs)                    ### comparação</font></div><div><font face="monospace, monospace">library(hydroGOF)</font></div><div><font face="monospace, monospace">ggof(sel.pre,sel.obs)</font></div><div><font face="monospace, monospace">### </code></font></div><div><br></div><div><img src="cid:ii_14d20c4c97bcbebc" alt="Imagem inline 1" height="412" width="519"><br></div><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div></div>
<br><div class="gmail_quote"><br></div></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>