erro na validação cruzada com pacote geoR
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
Tal vez, confira no seus pontos amostrais na VC se existem pares de pontos "muito" próximos e elimine um ponto deles. Em 16 de outubro de 2013 17:45, Hélio Gallo Rocha <heliogallorocha@gmail.com
escreveu:
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 R-br@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.
-- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc.
Não importa o programa, se a matriz é singular não inverte! Agora o lance é entender porque está acontecendo isso. Da próxima vez que postar coloque também um exemplo com código reprodutível para podermos ver o que é. Não é possível nem começar a responder se não tem um exemplo. O máximo que eu posso dizer é que pelo erro que você postou não há um problema no programa, mas sim com seus dados. Flavio Barros [image: Facebook]<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.facebook.com%2Fflavio.barros.1650%3Fref%3Dtn_tnmn&sn=> [image: LinkedIn]<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.linkedin.com%2Fprofile%2Fview%3Fid%3D61839390%26trk%3Dtab_pro&sn=> Contact me: [image: Google Talk] flaviomargarito@gmail.com “"Only those who will risk going too far can possibly find out how far one can go." - T.S. Eliot<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.quotesdaddy.com%2Fquote%2F307858%2Fts-eliot%2Fonly-those-who-will-risk-going-too-far-can-possibly&sn=> ” Get this email app! <http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.wisestamp.com%2Fapps%2Fquotes%3Futm_source%3Dextension%26utm_medium%3Demail%26utm_term%3Dquotes%26utm_campaign%3Dapps&sn=> Create your free signature: <http://s.wisestamp.com/links?url=http%3A%2F%2Fr1.wisestamp.com%2Fr%2Flanding%3Fpromo%3D33%26dest%3Dhttp%253A%252F%252Fwww.wisestamp.com%252Femail-install%253Futm_source%253Dextension%2526utm_medium%253Demail%2526utm_campaign%253Dpromo_33&sn=> CLICK HERE!<http://s.wisestamp.com/links?url=http%3A%2F%2Fr1.wisestamp.com%2Fr%2Flanding%3Fpromo%3D33%26amp%3Bdest%3Dhttp%253A%252F%252Fwww.wisestamp.com%252Femail-install%253Futm_source%253Dextension%2526utm_medium%253Demail%2526utm_campaign%253Dpromo_33&sn=> 2013/10/16 Luis Iván Ortiz Valencia <liov2067@gmail.com>
Tal vez, confira no seus pontos amostrais na VC se existem pares de pontos "muito" próximos e elimine um ponto deles.
Em 16 de outubro de 2013 17:45, Hélio Gallo Rocha < heliogallorocha@gmail.com> escreveu:
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 R-br@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.
-- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc.
_______________________________________________ R-br mailing list R-br@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.
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
Senhores, bom dia! Considerando o quarto item colocado pelo Prof. Paulo Justiniano, sobre a real necessidade da validação cruzada, talvez tenham interesse no texto de Davis (1987), intitulado "Uses and abuses of cross-validation in geostatistics". Obviamente há referências mais atuais, mas essa aí é bastante didática porque, além de tudo, esclarece diferenças entre "cross validation" e "jacknifing". Segue um excerto do texto original: "Cross-validation cannot confirm that a particular model is or is not the optimum. It is a method to better examine and understand the phenomenon under study using the available data." Davis, B. M., 1987. Uses and abuses of cross-validation in geostatistics. Math. Geol., v. 17, p. 563–586. Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
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@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-paco... 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
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@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 19 de outubro de 2013 09:26, Hélio Gallo Rocha <heliogallorocha@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@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-paco... 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@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.
Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores. Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver. Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: < http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.h...
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt> Se importaria em dispor esses dados novamente? Atte. Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Eder, valeu pela força. Amanhã cedo vou analisar sua sugestão, mas de antemão adianto que sem a validação cruzada não tem como ter o dado do erro, nem tampouco fazer a krigagem, a não ser que o prof. Paulo de uma dica de como gerar o erro sem a validação. Tentei fazer o seguinte, mas não deu certo: 1. Fiz o ajuste (C0,C1 e alance) em cima do variograma sem tendencia de 1St 2. Como o pacote geoR não valida, plotei o variograma com os dados com tendencia e usei os dados (C0,C1 e alance) sem tendencia. 3. Fiz a validação e gerei o mapa. 4. Claro, o mapa não ficou como deveria, usei o sufer e arcgis para comparar. Isso deve ter acontecido pois os erros foram diferentes, pois os semivariogramas são muito diferentes. Pelo menos tudo isso serviu pro meu aprendizado, me fez pensar muito mais para ver uma solução, o que não ocorreria se não tivesse o erro. Abraço Hélio Em 20 de outubro de 2013 16:03, Eder Comunello [via R-br] < ml-node+s2285057n4660663h21@n4.nabble.com> escreveu:
Prezado Helio, boa tarde!
Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: < http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.h...
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4660663&i=0>[hidden email] <http://user/SendEmail.jtp?type=node&node=4660663&i=1>> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4660663&i=2> 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-paco... 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
ver post anterior sobre patch a ser incorporado em breve On Sun, 20 Oct 2013, Hélio Gallo Rocha wrote:
Eder, valeu pela força. Amanhã cedo vou analisar sua sugestão, mas de antemão adianto que sem a validação cruzada não tem como ter o dado do erro, nem tampouco fazer a krigagem, a não ser que o prof. Paulo de uma dica de como gerar o erro sem a validação.
Tentei fazer o seguinte, mas não deu certo: 1. Fiz o ajuste (C0,C1 e alance) em cima do variograma sem tendencia de 1St 2. Como o pacote geoR não valida, plotei o variograma com os dados com tendencia e usei os dados (C0,C1 e alance) sem tendencia. 3. Fiz a validação e gerei o mapa. 4. Claro, o mapa não ficou como deveria, usei o sufer e arcgis para comparar. Isso deve ter acontecido pois os erros foram diferentes, pois os semivariogramas são muito diferentes.
Pelo menos tudo isso serviu pro meu aprendizado, me fez pensar muito mais para ver uma solução, o que não ocorreria se não tivesse o erro.
Abraço
Hélio
Em 20 de outubro de 2013 16:03, Eder Comunello [via R-br] <ml-node+s2285057n4660663h21@n4.nabble.com> escreveu: Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: <http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.html#a4659536>
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <[hidden email][hidden email]> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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-paco... To unsubscribe from R-br, click here. NAML
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
Eder, De nada vai adiantar nossa análise com o pacote geoR, sem a validação cruzada não são gerados 1.os dados preditivos, 2. o erro, 3. nem faz krigagem!!! Isso quando temos de retirar tendencia, e ponto. Ninguém da lista pode dar uma luz, estou num mato sem cachorro. As dicas do pro. Paulo não me ajudaram,. Analisei a opção 4 dele, mas o mapa gerado pela verosimilhança tá longe e analisar erros, como???? Hélio Em 20 de outubro de 2013 16:03, Eder Comunello [via R-br] < ml-node+s2285057n4660663h21@n4.nabble.com> escreveu:
Prezado Helio, boa tarde!
Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: < http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.h...
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4660663&i=0>[hidden email] <http://user/SendEmail.jtp?type=node&node=4660663&i=1>> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4660663&i=2> 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-paco... 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
Os comentários estao ininteligiveis. A verossimilhanca nao gera mapa, e sim estima parâmetros. nem faz a krigagem -- nao dá apra entender o que quis dizer os codigos dos pacotes geoestatísticos (qq um deles) para obeter as quantidades desejadas sugiro ser mais rpeciso e eventualmente enviar um CMR On Mon, 21 Oct 2013, Hélio Gallo Rocha wrote:
Eder, De nada vai adiantar nossa análise com o pacote geoR, sem a validação cruzada não são gerados 1.os dados preditivos, 2. o erro, 3. nem faz krigagem!!! Isso quando temos de retirar tendencia, e ponto.
Ninguém da lista pode dar uma luz, estou num mato sem cachorro.
As dicas do pro. Paulo não me ajudaram,. Analisei a opção 4 dele, mas o mapa gerado pela verosimilhança tá longe e analisar erros, como????
Hélio
Em 20 de outubro de 2013 16:03, Eder Comunello [via R-br] <ml-node+s2285057n4660663h21@n4.nabble.com> escreveu: Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: <http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.html#a4659536>
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <[hidden email][hidden email]> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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-paco... To unsubscribe from R-br, click here. NAML
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
Caros Éder, Luis e Paulo, Muitíssimo obrigado pelos retornos. Estou causando um certo desconforto com meu problema, sei que ninguém é obrigado a resolver o problema dos outros. Como havia comentado, comecei a utilizar o R, em especial o pacote geoR. Aprendi muito nesta lista com a colaboração de muitos. Com certeza não expus de forma clara minhas dúvidas, gerando confusão. Resolvi, como sugestão do Prof. Paulo, postar um CRM e mais dois arquivos txt, um de dados e outro de recorte da área, nem seria necessário este ultimo. O CRM roda todo, apenas sendo necessária a atenção no ponto 6. https://www.dropbox.com/sh/jbd2pmvzm41hs6b/Hk2tm7fr5v Conto mais uma vez com a colaboração. Não sei há interesse por parte de alguém da lista ou do Prof. Paulo para ministrar um curso de geoestatística, não apenas para comando do geoR, mas conceitual também. Não há tempo, que tal um curso a distância? Em 22 de outubro de 2013 11:39, Paulo Justiniano [via R-br] < ml-node+s2285057n4660678h86@n4.nabble.com> escreveu:
Os comentários estao ininteligiveis. A verossimilhanca nao gera mapa, e sim estima parâmetros.
nem faz a krigagem -- nao dá apra entender o que quis dizer
os codigos dos pacotes geoestatísticos (qq um deles) para obeter as quantidades desejadas
sugiro ser mais rpeciso e eventualmente enviar um CMR
On Mon, 21 Oct 2013, Hélio Gallo Rocha wrote:
Eder, De nada vai adiantar nossa análise com o pacote geoR, sem a validação cruzada não são gerados 1.os dados preditivos, 2. o erro, 3. nem faz krigagem!!! Isso quando temos de retirar tendencia, e ponto.
Ninguém da lista pode dar uma luz, estou num mato sem cachorro.
As dicas do pro. Paulo não me ajudaram,. Analisei a opção 4 dele, mas o mapa gerado pela verosimilhança tá longe e analisar erros, como????
Hélio
Em 20 de outubro de 2013 16:03, Eder Comunello [via R-br] <[hidden email] <http://user/SendEmail.jtp?type=node&node=4660678&i=0>> escreveu: Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: < http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.html#a4659536>
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <[hidden email][hidden email]> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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-paco...
To unsubscribe from R-br, click here. NAML
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4660678&i=1> 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-paco... 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
Comentário abaixo
Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores.
... não necessáriamente. O interpolador é exato se: 1. nugget = 0 (tausq = 0 na geoR) 2. o nugget for interpretado com variaçaõ de micro escala Na geoR o nuget pode ser considerado erro de medida, variacao de micro escala ou uma istura de ambos em qualquer proporcao Os argumentos "signal" e "micro.scale" definem, o comportamento desejado krige.control(type.krige = "ok", trend.d = "cte", trend.l = "cte", obj.model = NULL, beta, cov.model, cov.pars, kappa, nugget, micro.scale = 0, dist.epsilon = 1e-10, aniso.pars, lambda)
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: <http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.html#a4659536>
Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
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@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Em 19 de outubro de 2013 09:26, Hélio Gallo Rocha <heliogallorocha@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@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-paco... To unsubscribe from R-br, click here. NAML
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list R-br@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.
participantes (6)
-
Elias T Krainski -
Flavio Barros -
Hélio Gallo Rocha -
Luis Iván Ortiz Valencia -
Paulo Justiniano -
Éder Comunello