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

[Paulo Justiniano]

Comentário acerca de umponto abaixo

On Sat, 19 Oct 2013, Éder Comunello wrote:

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

Isto mesmo
os residuos podem ser analisados da mesma forma que a xvalid() faz
Uma laternativa é montar os residuos na forma de objeto da classe 
xvalid(0 e se aproveitar das funcoes prontas


Alems disto - o Elias me enviou um "patch" que devo colocar na xvalid para 
resolver isto

em breve disponível  em uma nova versao da geoR



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