[R-br] erro na validação cruzada com pacote geoR

[Éder Comunello]

Caro Hélio, boa tarde!

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 <http://cran-r.c3sl.ufpr.br/web/views/Spatial.html>. 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...

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...

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?

Atte.,

Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]


Em 19 de outubro de 2013 09:26, Hélio Gallo Rocha <heliogallorocha em gmail.com
> escreveu:

> Caros da lista,
>
> Obrigado pelas dicas, se puderem me elucidar em dois pontos:
>
> 1. Quais outros pacotes trabalham com geoestatistica.
> 2. Com analisar o residuo se ele é gerado na validação?
>
> grato
>
> Hélio
>
>
> Em 17 de outubro de 2013 09:20, Paulo Justiniano [via R-br] <
> ml-node+s2285057n4660630h98 em n4.nabble.com> escreveu:
>
>> Há um problema na geoR para fazer a validacao cruzada com trend (model de
>> media nao constante.
>>
>> Algumas opcoes:
>> 1. fazer um  pequeno programa que remova uma observacao por vez e chame a
>> krige.conv()
>>
>> 2. Usar outros pacotes do R
>>
>> 3. Fazer a valizacao nos "residuos"
>>
>> 4. Pensar se a validacao cruzada é de fato nbecessária. Embora seja
>> usada é bem menos informativa e util que outras medidadas de avaliacao do
>> ajuste baseadas na verossimilhanã e análise dos residuos
>>
>>
>> On Wed, 16 Oct 2013, Hélio Gallo Rocha wrote:
>>
>> > Caros da lista.
>> > Estou uma situação bastante complicada.
>> >
>> > Estou no fim do mestrado, onde utilizei o R todo o tempo.
>> >
>> > Mas estou numa situação que talvez tenha de deixar de usar o R e passar
>> para um "pago".
>> >
>> > Já postei sobre  este problema e vi outros com o mesmo.
>> >
>> > Após a retirada de tendencia do semivariograma e fazer a validação
>> cruzada para escolha do melhor
>> > modelo o R dá esse erro:
>> > xErro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) :
>> >   sistema é computacionalmente singular: condição recíproca número =
>> 1.18641e-024
>> >
>> > Alguém poderia me dar uma luz?
>> >
>> > Não tem solução? é isso?
>> >
>> > grato a todos
>> >
>> > --
>> > Hélio Gallo Rocha
>> > IFSULDEMINAS - Câmpus Muzambinho
>> >
>> >
>> _______________________________________________
>> R-br mailing list
>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4660630&i=0>
>> 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.
>>
>> ------------------------------
>>  If you reply to this email, your message will be added to the
>> discussion below:
>>
>> http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pacote-geoR-tp4660619p4660630.html
>>  To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
>> .
>> NAML<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>
>>
>
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131019/02697661/attachment.html>