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 <[email protected]>
Dourados, MS - [22 16.5'S, 54 49'W]
Em 19 de outubro de 2013 09:26, Hélio Gallo Rocha <[email protected]>
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]
<[email protected]> 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
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.
