Caros listeiros. Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem: KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0)) deu esse erro: Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated"). Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
Hélio, boa tarde! Poderia disponibilizar o CMR e os dados que geram o erro? Também seria bom o sessionInfo() com geoR carregada, pra ver a versão. Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma instrução de uma versão antiga da geoR. Mas não deve causar maiores problemas. Até mais, Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Caros listeiros. Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem: KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0)) deu esse erro: Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated"). Agradeço desde já -- 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.
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que está usnado modelo com covariáveis (trend) e que estes possuem valores muito grandes ou pequenos Se suas coordenadas forem em UTM considere dividi-las por 100 quanto ao warning é inofensivo mas isto foi resolvido na ultima versão disponível no CRAN On Fri, 15 May 2015, Éder Comunello wrote:
Hélio, boa tarde!
Poderia disponibilizar o CMR e os dados que geram o erro?
Também seria bom o sessionInfo() com geoR carregada, pra ver a versão. Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma instrução de uma versão antiga da geoR. Mas não deve causar maiores problemas.
Até mais,
Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Caros listeiros.
Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem:
KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0))
deu esse erro:
Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated").
Agradeço desde já
-- 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.
Caros Paulo e Eder, Realmente estou utilizando modelos com trend. Dividi as coordenadas por 1000. Já utilizei muito o CRM em questão e nunca havia dado este erro. Os dados são referentes a 400 estações meteorológicas no Paraná que serão comparadas a um modelo ECMWF. Quanto a sessionInfo() : R - 3.1.3 geoR - 1.7-4.1 Em 15 de maio de 2015 17:02, Paulo Justiniano [via R-br] < ml-node+s2285057n4664494h82@n4.nabble.com> escreveu:
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que está usnado modelo com covariáveis (trend) e que estes possuem valores muito grandes ou pequenos
Se suas coordenadas forem em UTM considere dividi-las por 100
quanto ao warning é inofensivo mas isto foi resolvido na ultima versão disponível no CRAN
On Fri, 15 May 2015, Éder Comunello wrote:
Hélio, boa tarde!
Poderia disponibilizar o CMR e os dados que geram o erro?
Também seria bom o sessionInfo() com geoR carregada, pra ver a versão.
Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma instrução de uma
versão antiga da geoR. Mas não deve causar maiores problemas.
Até mais,
Éder Comunello <[hidden email] <http:///user/SendEmail.jtp?type=node&node=4664494&i=0>> Dourados, MS - [22 16.5'S, 54 49'W]
Caros listeiros.
Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem:
KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0))
deu esse erro:
Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated").
Agradeço desde já
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664494&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.
_______________________________________________ R-br mailing list [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664494&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/R-br-erro-na-krigagem-tp4664489p4664494.ht... 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
ceritfique-se que dividiu porm100 tabto as cos dados qto as do grid de predicao se persistir, experimente subtrair os minimos (globais) das coordenadas (isto nao altera as distnacias) gr[,1] <- gr[,1] - min(gr[,1]) gr[,2] <- gr[,2] - min(gr[,2]) gd$coords[,1] <- gd$coords[,1] - min(gr[,1]) gd$coords[,2] <- gd$coords[,2] - min(gr[,2]) On Fri, 15 May 2015, Hélio Gallo Rocha wrote:
Caros Paulo e Eder, Realmente estou utilizando modelos com trend. Dividi as coordenadas por 1000.
Já utilizei muito o CRM em questão e nunca havia dado este erro.
Os dados são referentes a 400 estações meteorológicas no Paraná que serão comparadas a um modelo ECMWF.
Quanto a sessionInfo() : R - 3.1.3 geoR - 1.7-4.1
Em 15 de maio de 2015 17:02, Paulo Justiniano [via R-br] <ml-node+s2285057n4664494h82@n4.nabble.com> escreveu: Quanto ao erro pode ser vevido a diversas coisas mas suspeito que está usnado modelo com covariáveis (trend) e que estes possuem valores muito grandes ou pequenos
Se suas coordenadas forem em UTM considere dividi-las por 100
quanto ao warning é inofensivo mas isto foi resolvido na ultima versão disponível no CRAN
On Fri, 15 May 2015, Éder Comunello wrote:
> > Hélio, boa tarde! > > Poderia disponibilizar o CMR e os dados que geram o erro? > > Também seria bom o sessionInfo() com geoR carregada, pra ver a versão. Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma instrução de uma > versão antiga da geoR. Mas não deve causar maiores problemas. > > Até mais, > > Éder Comunello <[hidden email]> > Dourados, MS - [22 16.5'S, 54 49'W] > > Caros listeiros. > > > Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem: > > KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0)) > > deu esse erro: > > Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : > sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 > Além disso: Mensagens de aviso perdidas: > 'overlayPointsWithPolygons' is deprecated. > Use 'over' instead. > See help("Deprecated") and help("sp-deprecated"). > > > Agradeço desde já > > > -- > 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. > > _______________________________________________ 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/R-br-erro-na-krigagem-tp4664489p4664494.ht... To unsubscribe from R-br, click here. NAML
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
Paulo, Atualizei o geoR, quanto ao grid de predição, está limitado pela borda( do estado), veja: gr <- pred_grid(gd$borders, by=5) que ao ser plotado apresenta coordenadas div. por 1000. Após a atualização do geoR, apresentou a parte superior do erro anterior: Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 após rodar as linhas sugeridas, não houve diferença no valor das coordenadas dos dados(gd) Em 15 de maio de 2015 17:20, Paulo Justiniano [via R-br] < ml-node+s2285057n4664496h44@n4.nabble.com> escreveu:
ceritfique-se que dividiu porm100 tabto as cos dados qto as do grid de predicao
se persistir, experimente subtrair os minimos (globais) das coordenadas (isto nao altera as distnacias)
gr[,1] <- gr[,1] - min(gr[,1]) gr[,2] <- gr[,2] - min(gr[,2])
gd$coords[,1] <- gd$coords[,1] - min(gr[,1]) gd$coords[,2] <- gd$coords[,2] - min(gr[,2])
On Fri, 15 May 2015, Hélio Gallo Rocha wrote:
Caros Paulo e Eder, Realmente estou utilizando modelos com trend. Dividi as coordenadas por
Já utilizei muito o CRM em questão e nunca havia dado este erro.
Os dados são referentes a 400 estações meteorológicas no Paraná que
serão comparadas a um modelo ECMWF.
Quanto a sessionInfo() : R - 3.1.3 geoR - 1.7-4.1
Em 15 de maio de 2015 17:02, Paulo Justiniano [via R-br] <[hidden email]
<http:///user/SendEmail.jtp?type=node&node=4664496&i=0>> escreveu:
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que
está
usnado modelo com covariáveis (trend) e que estes possuem valores
muito
grandes ou pequenos
Se suas coordenadas forem em UTM considere dividi-las por 100
quanto ao warning é inofensivo mas isto foi resolvido na ultima
versão
disponível no CRAN
On Fri, 15 May 2015, Éder Comunello wrote:
> > Hélio, boa tarde! > > Poderia disponibilizar o CMR e os dados que geram o erro? > > Também seria bom o sessionInfo() com geoR carregada, pra ver a
versão. Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma
instrução de uma > versão antiga da geoR. Mas não deve causar maiores problemas. > > Até mais, > > Éder Comunello <[hidden email]> > Dourados, MS - [22 16.5'S, 54 49'W] > > Caros listeiros. > > > Rodando krigagem de clima, com retirada de tendência polinomial
de segunda ordem:
> > KC.2nd <-
krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0))
> > deu esse erro: > > Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data)))
:
> sistema é computacionalmente singular: condição recíproca
número = 4.2826e-23
> Além disso: Mensagens de aviso perdidas: > 'overlayPointsWithPolygons' is deprecated. > Use 'over' instead. > See help("Deprecated") and help("sp-deprecated"). > > > Agradeço desde já > > > -- > 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.
> > _______________________________________________ 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/R-br-erro-na-krigagem-tp4664489p4664494.ht...
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=4664496&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/R-br-erro-na-krigagem-tp4664489p4664496.ht... 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
Veja se tem alguma estacões "muito próximas" ou "repetidas" , aparentemente 2 colunas na sua matriz de semivariogramas sao "muito" parecidas. Se tiver retire uma. Em 15 de maio de 2015 17:55, Hélio Gallo Rocha <heliogallorocha@gmail.com> escreveu:
Paulo, Atualizei o geoR, quanto ao grid de predição, está limitado pela borda( do estado), veja: gr <- pred_grid(gd$borders, by=5)
que ao ser plotado apresenta coordenadas div. por 1000.
Após a atualização do geoR, apresentou a parte superior do erro anterior:
Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23
após rodar as linhas sugeridas, não houve diferença no valor das coordenadas dos dados(gd)
Em 15 de maio de 2015 17:20, Paulo Justiniano [via R-br] < ml-node+s2285057n4664496h44@n4.nabble.com> escreveu:
ceritfique-se que dividiu porm100 tabto as cos dados qto as do grid de predicao
se persistir, experimente subtrair os minimos (globais) das coordenadas (isto nao altera as distnacias)
gr[,1] <- gr[,1] - min(gr[,1]) gr[,2] <- gr[,2] - min(gr[,2])
gd$coords[,1] <- gd$coords[,1] - min(gr[,1]) gd$coords[,2] <- gd$coords[,2] - min(gr[,2])
On Fri, 15 May 2015, Hélio Gallo Rocha wrote:
Caros Paulo e Eder, Realmente estou utilizando modelos com trend. Dividi as coordenadas por
Já utilizei muito o CRM em questão e nunca havia dado este erro.
Os dados são referentes a 400 estações meteorológicas no Paraná que
serão comparadas a um modelo ECMWF.
Quanto a sessionInfo() : R - 3.1.3 geoR - 1.7-4.1
Em 15 de maio de 2015 17:02, Paulo Justiniano [via R-br] <[hidden
email] <http:///user/SendEmail.jtp?type=node&node=4664496&i=0>> escreveu:
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que
está
usnado modelo com covariáveis (trend) e que estes possuem valores
muito
grandes ou pequenos
Se suas coordenadas forem em UTM considere dividi-las por 100
quanto ao warning é inofensivo mas isto foi resolvido na ultima
versão
disponível no CRAN
On Fri, 15 May 2015, Éder Comunello wrote:
> > Hélio, boa tarde! > > Poderia disponibilizar o CMR e os dados que geram o erro? > > Também seria bom o sessionInfo() com geoR carregada, pra ver a
versão. Esse warning é pra avisar que o comando da {sp} mudou. Poderia ser uma
instrução de uma > versão antiga da geoR. Mas não deve causar maiores problemas. > > Até mais, > > Éder Comunello <[hidden email]> > Dourados, MS - [22 16.5'S, 54 49'W] > > Caros listeiros. > > > Rodando krigagem de clima, com retirada de tendência polinomial
de segunda ordem:
> > KC.2nd <-
krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0))
> > deu esse erro: > > Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data)))
:
> sistema é computacionalmente singular: condição recíproca
número = 4.2826e-23
> Além disso: Mensagens de aviso perdidas: > 'overlayPointsWithPolygons' is deprecated. > Use 'over' instead. > See help("Deprecated") and help("sp-deprecated"). > > > Agradeço desde já > > > -- > 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.
> > _______________________________________________ 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/R-br-erro-na-krigagem-tp4664489p4664494.ht...
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=4664496&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/R-br-erro-na-krigagem-tp4664489p4664496.ht... 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.
-- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc.
participantes (4)
-
Hélio Gallo Rocha -
Luis Iván Ortiz Valencia -
Paulo Justiniano -
Éder Comunello