[R-br] erro na krigagem

Sexta Maio 15 17:55:59 BRT 2015

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150515/ead8d289/attachment.html>