<div dir="ltr">Caro Hélio, boa tarde!<div><br></div><div>Acerca dos outros pacotes que trabalham com geoestatística, um bom ponto de partida é a "CRAN Task View: Analysis of Spatial Data" que você pode acessar em <<a href="http://cran-r.c3sl.ufpr.br/web/views/Spatial.html">http://cran-r.c3sl.ufpr.br/web/views/Spatial.html</a>>. Tem um tópico "Geostatistics" onde há links para os principais pacotes. Contudo, na minha opinião particular, não sendo possível usar o 'geoR', eu partiria para o 'gstat'. Acredito que seja o pacote geoestatístico mais frequente na lista R-sig-geo, outro local onde dá pra conseguir muita informação...</div>
<div><br></div><div>Acredito que o resíduo referido pelo Prof. Paulo seja aquele dado pela diferença entre os dados observados e a predição da krigagem. Quando vc roda o 'krige.conv' você utiliza o modelo ajustado para predizer pontos em um grade definida (p. ex. loc=grid), mas parece ser possível trocar o grid pelas coordenadas dos pontos originais (p. ex. loc=pts). Feito isso, poderia obter a diferença observados-preditos (resíduo). Se for isso mesmo, obter os resíduos não seria tão difícil, mas também não tenho muita clareza em como prosseguir daí pra frente...</div>
<div><br></div><div>Eu tenho uns estudos geoestatísticos em 'standby' que pretendo retomar em breve. Neles eu me deparei com o mesmo problema que você vivencia agora. Sendo assim, tenho grande interesse nessa discussão, já que terei que lidar com isso também. Das estratégias oferecidas pelo Prof. Paulo, eu particularmente pretendo utilizar a primeira, que julgo ser a mais fácil de implementar a curto prazo. Considerou essa possibilidade?</div>
<div><br></div><div>Atte.,</div></div><div class="gmail_extra"><br clear="all"><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>
<br><br><div class="gmail_quote">Em 19 de outubro de 2013 09:26, Hélio Gallo Rocha <span dir="ltr"><<a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@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">Caros da lista,<div><br></div><div>Obrigado pelas dicas, se puderem me elucidar em dois pontos:</div><div><br></div><div>1. Quais outros pacotes trabalham com geoestatistica.</div><div>2. Com analisar o residuo se ele é gerado na validação?</div>
<div><br></div><div>grato</div><div><br></div><div>Hélio</div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 17 de outubro de 2013 09:20, Paulo Justiniano [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4660630h98@n4.nabble.com" target="_blank">ml-node+s2285057n4660630h98@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Há um problema na geoR para fazer a validacao cruzada com trend (model de
<br>media nao constante.
<br><br>Algumas opcoes:
<br>1. fazer um pequeno programa que remova uma observacao por vez e chame a
<br>krige.conv()
<br><br>2. Usar outros pacotes do R
<br><br>3. Fazer a valizacao nos "residuos"
<br><br>4. Pensar se a validacao cruzada é de fato nbecessária. Embora seja
<br>usada é bem menos informativa e util que outras medidadas de avaliacao do
<br>ajuste baseadas na verossimilhanã e análise dos residuos
<br></div></div><div><div><div class="h5"><br><br>On Wed, 16 Oct 2013, Hélio Gallo Rocha wrote:
<br><div><div><br>> Caros da lista.
<br>> Estou uma situação bastante complicada.
<br>>
<br>> Estou no fim do mestrado, onde utilizei o R todo o tempo.
<br>>
<br>> Mas estou numa situação que talvez tenha de deixar de usar o R e passar para um "pago".
<br>>
<br>> Já postei sobre este problema e vi outros com o mesmo.
<br>>
<br>> Após a retirada de tendencia do semivariograma e fazer a validação cruzada para escolha do melhor
<br>> modelo o R dá esse erro:
<br>> xErro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) :
<br>> sistema é computacionalmente singular: condição recíproca número = 1.18641e-024
<br>>
<br>> Alguém poderia me dar uma luz?
<br>>
<br>> Não tem solução? é isso?
<br>>
<br>> grato a todos
<br>>
<br>> --
<br>> Hélio Gallo Rocha
<br>> IFSULDEMINAS - Câmpus Muzambinho
<br>>
<br>></div></div></div></div>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4660630&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><div class="im"><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
</div></div><hr noshade size="1" color="#cccccc">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div>
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pacote-geoR-tp4660619p4660630.html" target="_blank">http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pacote-geoR-tp4660619p4660630.html</a>
</div><div><div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><div class="im"><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</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>