[R-br] erro na validação cruzada

[Hélio Gallo Rocha]

Fernando, ai vai o CRM.
Não consegui disponibilizar o arquivo via R, os comandos sugeridos na lista
não deram certo, então terá de salvar num txt, para não atachar arquivo.

Dados:
 https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt

library(geoR)
#direcionando para dir onde está o arquivo dados250.txt
setwd("k:/2013/r/")

#carregando arquivo:
dados=read.geodata("dados250.txt ",coords.col=1:2,data.col=3,sep =
"",header=T);summary(dados);n=summary(dados)[[1]];n

#distancia máxima
hmax=summary(dados)[[3]][[2]]*.8 ;hmax

#variograma com direção de 135 graus
vario135= variog(dados,uvec=12,pairs.min=100,max.dist=hmax, direction=(pi/4
+ pi/2),estimator="classical")
plot(vario135, xlab="cento e trinta e cinco",main="numero de
pontos");text(vario135$u,vario135$v,vario135$n,pos=1);lines(vario135$u,vario135$v)

#Ajustar variograma experimental ao teórico
as135=eyefit(vario135)

# valores do Patamar(sill), alcance e nugget , só para comparar
sillp135=as135[[1]]$cov.pars[1];sillp135   ##C1
alcance135=as135[[1]]$cov.pars[2];alcance135 ##phi
pepita135=as135[[1]]$nugget;pepita135    ##C0

#validação cruzada
xv.as135=xvalid(dados,model=as135)

ai dá o erro:
xvalid: performing cross-validation at location ... 1, Erro em
solve.default(Vcov, trend.d) :
  sistema é computacionalmente singular: condição recíproca número =
2.24178e-20

Grato

Hélio



Em 3 de junho de 2013 23:05, Fernando Mayer [via R-br] <
ml-node+s2285057n4659534h22 em n4.nabble.com> escreveu:

> Como não houveram respostas ainda, eu sugiro que vc forneça mais
> detalhes como comandos e pacotes que vc utilizou, e um CMR para que as
> pessoas possam reproduzir o erro. Veja também que nem todos sabem o
> que eh exatamente "validação cruzada na direção 45 graus", então uma
> breve expilcação também vai ajudar.
>
> Sobre o erro, está claro pela mensagem que a matriz Vcov é singular
> (não possui inversa), então o solve() não vai resolver. Sem fornecer
> mais detalhes é o que pode ser concluído...
>
> []s,
>
>
>
> ---
> Fernando de Pol Mayer
> Doutorando em Estatística e Experimentação Agronômica
> Escola Superior de Agricultura "Luiz de Queiroz" - ESALQ
> Universidade de São Paulo - USP
> URL: http://fernandomayer.github.com
> e-mail: fernando.mayer [@] {gmail.com, usp.br}
>
>
> 2013/6/3 Hélio Gallo Rocha <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659534&i=0>>:
>
>
> > Caros da Lista,
> >
> > Está dando este erro ao tentar rodar a validação cruzada na direção de
> 45
> > graus:
> >
> > xvalid: performing cross-validation at location ... 1, Erro em
> > solve.default(Vcov, trend.d) :
> >   sistema é computacionalmente singular: condição recíproca número =
> > 8.6266e-21
> >
> >
> >
> > Grato a todos
> >
> > --
> > Hélio Gallo Rocha
> > IFSULDEMINAS - Câmpus Muzambinho
> >
> > _______________________________________________
> > R-br mailing list
> > [hidden email] <http://user/SendEmail.jtp?type=node&node=4659534&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=4659534&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-validacao-cruzada-tp4659526p4659534.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/20130604/451ff985/attachment.html>