[R-br] erro na krigagem
Luis Iván Ortiz Valencia
liov2067 em gmail.com
Sábado Maio 16 01:17:12 BRT 2015
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 em 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 em 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
>> 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] <[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.html
>> > 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.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.
>
--
Luis Iván Ortiz Valencia
Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ
Estatístico Msc.
Spatial Analyst Msc.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150516/9eed1f96/attachment.html>
Mais detalhes sobre a lista de discussão R-br