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:
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