<div dir="ltr">Wagner, boa tarde!<div><br></div><div>Eu não recomendo usar <font face="monospace, monospace">splinefun()</font> porque a calibração do modelo feita na execução de  <span style="font-family:monospace,monospace;font-size:12.8000001907349px">mnimput() </span><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">é bem mais sofisticada</font></span><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">. O objeto </font><font face="monospace, monospace">imput </font><font face="arial, helvetica, sans-serif">é  o próprio objeto </font><font face="monospace, monospace">miss2 </font><font face="arial, helvetica, sans-serif">com as falhas preenchidas com base na covariância e comportamento temporal das séries. Já é o produto final, não precisa ajustar por </font><font face="monospace, monospace">splinefun()</font><font face="arial, helvetica, sans-serif">.</font></span></div><div><br></div><div><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">Se você fizer </font><font face="monospace, monospace">imput-miss2 </font><font face="arial, helvetica, sans-serif">verá que a diferença é 0 para todos valores iniciais e um novo valor será imputado para todos os casos onde havia NA. </font></span></div><div><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif"><br></font></span></div><div><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">O procedimento que postei faz uma avaliação geral da qualidade do ajuste pras 5 séries porque exige menor quantidade de pontos de teste. Daria pra fazer uma avaliação de cada série, mas teria que reservar uma quantidade significativa de pontos de cada série, o que pode impactar no modelo.</font></span></div><div><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif"><br></font></span></div><div>Além do mais, se você utilizar <font face="monospace, monospace">splinefun</font>() pra predizer valores utilizados na calibração, vai retornar os próprios valores utilizados...</div><div><span style="font-family:monospace,monospace;font-size:12.8000001907349px">fun <- splinefun(sel.obs, sel.pre)</span><br></div><div><div style="font-family:monospace,monospace;font-size:12.8000001907349px">all.equal(fun(sel.obs), sel.pre) ### [1] TRUE</div><div style="font-family:monospace,monospace;font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">Notar que:</font></div><div style="font-family:monospace,monospace"><div><span style="font-size:12.8000001907349px">args(splinefun)</span></div><div><span style="font-size:12.8000001907349px">### function (x, y = NULL, method = c("fmm", "periodic", "natural", </span></div><div><span style="font-size:12.8000001907349px">###    "monoH.FC", "hyman"), ties = mean) </span></div><div><br></div></div></div><div><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">A forma que segue é coerente, mas os resultados não são tão animadores quanto os retornados por </font><font face="monospace, monospace">mtsdi::</font></span><span style="font-size:12.8000001907349px"><font face="monospace, monospace">mnimput()</font></span></div><div><div><font face="monospace, monospace"><span style="font-size:12.8000001907349px">fun2 <- splinefun(miss[notNA]) ### ajusta uma spline em função dos dados disponíveis</span></font></div><div><font face="monospace, monospace"><span style="font-size:12.8000001907349px">sel.pre2 <- fun2(sel.obs)      ### usa a splinefun pra estimar os dados reservados</span></font></div><div><font face="monospace, monospace"><span style="font-size:12.8000001907349px">ggof(sel.pre2,sel.obs)         ### avalia predições comparando com dados reservados ("verdade")</span></font></div></div><div><span style="font-family:monospace,monospace;font-size:12.8000001907349px"><br></span></div><div><span style="font-family:monospace,monospace;font-size:12.8000001907349px"><img src="https://mail.google.com/mail/u/0/?ui=2&ik=c48285784a&view=fimg&th=14d25a2d239939ab&attid=0.1&disp=emb&realattid=ii_14d25a04f2d2a87e&attbid=ANGjdJ-mzt08d4nd91iahQNqUudeV3NvZjV1XBYpkKj6SVx4dSA7O3cqdiySp7GgwPzJLCYJTysQ75OPgZcds2y6qv2ELRJsKZt0NvADvNp2HT5giJKjBcBHnW_jMU8&sz=s0-l75&ats=1430855537815&rm=14d25a2d239939ab&zw" alt="Imagem inline 1" width="519" height="361"><br></span></div><div><span style="font-family:monospace,monospace;font-size:12.8000001907349px"><br></span></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><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">Em 5 de maio de 2015 13:57, Wagner Wolff <span dir="ltr"><<a href="mailto:wwolff@usp.br" target="_blank">wwolff@usp.br</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><font face="monospace, monospace">Olá Éder vou tentar esclarecer com seu exemplo, minhas modificações estão em negrito.<br></font></div><span class=""><div><font face="monospace, monospace"><br>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)<br><b><br></b></font></div></span><div><font face="monospace, monospace"><b>fun.cal <- splinefun(sel.pre,sel.obs)<br>     ## função para calibrar<br>sel.pre.cal <- fun.cal(sel.obs)</b>           <b>## calcular valores preditos após calibração<br><br></b></font></div><div><font face="monospace, monospace"><b>## Assim usaria a função fun.cal para calibrar outros valores imputados do mesmo banco de dados.<br></b></font></div><div><font face="monospace, monospace"><b>Entretanto, não foi útil, pois quando uso para outros valores imputados houve aumento no erro<br></b></font></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">Em 5 de maio de 2015 10:40, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Wagner, bom dia!<div><br></div><div>Não entendi exatamente o que você está buscando. Teria dados pra um exemplo mais prático?</div><div><br></div><div>Por exemplo, você quer completar dados de uma série de chuva  em particular e dispõe dados de quatro estações próximas, assumindo algum grau de correlação/covariância entre elas? Sua preocupação é saber quão "boa" foi a predição ou algo assim?</div><span><div><br></div><div><br><div class="gmail_extra"><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></span></div>
<br></div></div><span class="">_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></span></blockquote></div><br></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></div>