Variograma ajustado

Para o ajuste do variograma anterior é necessário fazer um chute do efeito pepita, patamar e alcance para o modelo esférico. Os valores escolhidos foram 2,5, 14 e 1,1, respectivamente.

Programa:

data(jura)

vgm1=variogram(Co~1,locations=~Xloc+Yloc,data=prediction.dat)

f<-fit.variogram(vgm1,vgm(14,"Sph",1.1,2.5))

ff<-variogramLine(f,maxdist=2,n=500,min=0.000001)

plot(ff,col="blue",main="Variograma ajustado (Cobalto)",xlab="distância")

points(vgm1[,2],vgm1[,3],col="red")

Resultado:

Valores ajustados dos parâmetros (comando fit.variogram):

model psill range

1 Nug 1.304963 0.000000

2 Sph 12.524976 1.183541

Sendo Sph referente ao modelo esférico onde psill é o patamar e range é o alcance. Nug é o efeito de pepita.