[R-br] erro na krigagem

[Paulo Justiniano]

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 em 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.html
> To unsubscribe from R-br, click here.
> NAML
> 
> 
> 
> 
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
> 
>