Displaying 1 result from an estimated 1 matches for "geodata1".
Did you mean:
geodata
2009 Sep 17
0
geoR, variofit
....
Thanks, Sascha Bellaire!
Parts of code, just multiplying wl1$PSI with 1000 changed it.
#Calculate the extent of the current "sampling" design
locations <- cbind(wl1$X,wl1$Y)
extent <- max(dist(locations))
max.dist <- extent/2
data1 <- cbind(wl1$X,wl1$Y,wl1$PSI*1000)
geodata1 <- as.geodata(data1, coords.col = 1:2, data.col = 3)
# Calculate the sample variogram
vario_var1 <- variog(geodata1,
uvec = nlags,
trend = trend,
option = "bin",
estimator.type = "modulus",
max.dist = max.dist,...