I am trying to extract the semivariance values associated with a given semivariogram model developed in gstat, the end goal being to compare modelled semivariance with observed semivariance at defined distances. Using the function variogramLine to get semivariance values from my semivariogram model, I find that the values don't match the model when plotted over the semivariogram itself.
A sample of the data is available here: https://www.dropbox.com/sh/mjqhika3stqjv6m/AAAIN_89Hpaiye0G7f_J-Hy_a?dl=0. My code looks as follows:
library(raster)library(sp)library(gstat)library(DiceEval)## Read in sample datasample_region_sub=read.csv('sample_region_sub_zeroed.csv')coordinates(sample_region_sub) = ~x+yxx=coordinates(sample_region_sub)[,1]yy=coordinates(sample_region_sub)[,2]## Create anisotropic variogramanisotropy_major=60anisotropy_minor=0.4aniso
A sample of the data is available here: https://www.dropbox.com/sh/mjqhika3stqjv6m/AAAIN_89Hpaiye0G7f_J-Hy_a?dl=0. My code looks as follows:
library(raster)library(sp)library(gstat)library(DiceEval)## Read in sample datasample_region_sub=read.csv('sample_region_sub_zeroed.csv')coordinates(sample_region_sub) = ~x+yxx=coordinates(sample_region_sub)[,1]yy=coordinates(sample_region_sub)[,2]## Create anisotropic variogramanisotropy_major=60anisotropy_minor=0.4aniso