I am using gstat package for deriving a variogram using exponential model. My R code goes like this:

g <- gstat(id= "HV",formula=HV~X+Y,data=vari)
vari.var <- variogram(g)
plot(vari.var,col="black",pch=16,cex=1.3,xlab="Distance",ylab="Semivariance",main="Omnidirectional variogram NDVI")
model.vari <- fit.variogram(vari.var,vgm("Exp",range=range,nugget=0),fit.method=2)
plot(vari.var,model=model.vari,col="black",pch=16,cex=1.3,main="Omnidirectional variogram NDVI").

I am not sure whether it is correct or not.

Also I am not sure how to calculate sill values.

1 Answer 1


Yes, this is correct. When you print the model by typing


you'll see sill values, split up in a nugget component (the offset) and the exponential component. The sum of these two is usually indicated by "the sill value" (i.e., around 25).

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.