I tried to compare the result obtained by fitting the spherical variogram model using fit.variogram and nls functions. The large difference in the results is a great concern for me knowing that the two functions use Weighted Least Squares (WLS) approach for estimating parameters. The steps taken are as follows:

## Fit empirical semivariogram using gstat

empvar