I understand we are migrating to gstlearn, but I have an RGeostats question on fitting madograms or rodograms. It appears that model.auto does not seem to work with madograms, so I am using model.fit. It also doesn't seem to work. I get an error message, "Error in model.input(ndim = ndim, nvar = nvar, flag.sill = flag.sill, :
could not find function "model.input". Here is my code:
madfit <- # madfit is the name I want to assign for the modeled madogram
model.fit(
mad_Omni, # the name of the experimental variogram model I wish to fit - created in earlier code
struct = struct, # struct = the exponential variogram - assigned in earlier code
flag.fit = TRUE,
draw = TRUE,
verbose = 1,
title = paste(property, "Madogram Model Omnidirectional"),
# pos.legend = 1,
xlab = "Lag distance",
ylab = expression(paste("Variance (", gamma, "(h))", sep = ""))
)
Where am I going wrong?