FitVariogram: Fit Varigoram using autofitVariogram modified
In AybarCL/Dorado: Easy Rainfall and Temperature Interpolation

Description Usage Arguments Details Value Author(s) See Also

This function is a very slight variation about the function autofitVariogram of packages automap. Only added the param boundaries so the user can define it.

FitVariogram(formula, input_data, model = c("Sph", "Exp", "Gau", "Ste"),
  kappa = c(0.05, seq(0.2, 2, 0.1), 5, 10), fix.values = c(NA, NA, NA),
  verbose = FALSE, GLS.model = NA, start_vals = c(NA, NA, NA),
  miscFitOptions = list(), boundaries, ...)

`formula`	that defines the dependent variable as a linear model of independent variables; suppose the dependent variable has name 'z', for ordinary and simple kriging use the formula 'z~1'; for simple kriging also define 'beta' (see below); for universal kriging, suppose 'z' is linearly dependent on 'x' and 'y', use the formula 'z~x+y'.
`input_data`	An object of SpatialPointsDataFrame-class.
`model`	The list of variogrammodels that will be tested.
`kappa`	Smoothing parameter of the Matern model. Provide a list if you want to check more than one value.
`fix.values`	Can be used to fix a variogram parameter to a certain value. It consists of a list with a length of three. The items describe the fixed value for the nugget, range and sill respectively. They need to be given in that order. Setting the value to NA means that the value is not fixed.
`verbose`	logical, if TRUE the function will give extra feedback on the fitting process
`GLS.model`	If a variogram model is passed on through this parameter a Generalized Least Squares sample variogram is calculated.
`start_vals`	Can be used to give the starting values for the variogram fitting. The items describe the fixed value for the nugget, range and sill respectively. They need to be given in that order. Setting the value to NA means that the value will be automatically chosen.
`miscFitOptions`	A list with named arguments that provide additional control over the fitting process. For example: `list(merge.small.bins = TRUE)`. If the list is empty, autofitVariogram uses default values. The following parameters can be set: `merge.small.bins`: logical, when TRUE, the function checks if there are bins with less than 5 points. If so, the first two bins are merged and the check is repeated. This is done until all bins have more than `min.np.bin` points. `min.np.bin`: integer, the minimum number of points allowed in a bin before we start merging bins. See also `merge.small.bins`.
`boundaries`	Are the Boundaries for the bins in km.
`...`	parameters that are passed on to variogram variogram when calculating the sample variogram

Geostatistical routines are used from package gstat. A few simple choices are made when estimating the inital guess for fit.variogram. The initial sill is estimated as the mean of the max and the median of the semi-variance. The inital range is defined as 0.10 times the diagonal of the bounding box of the data. The initial nugget is defined as the min of the the semi-variance. There are five different types of models that are often used:

Sph: A shperical model.
Exp: An exponential model.
Gau: A gaussian model.
Mat: A model of the Matern familiy
Ste: Matern, M. Stein's parameterization

A list of all permitted variogram models is available by typing vgm() into the R console. autofitVariogram iterates over the variogram models listed in model and picks the model that has the smallest residual sum of squares with the sample variogram. For the Matern model, all the kappa values in kappa are tested.

Note that when using the power model, and not specifying starting values yourself, the sill is set to 1, the range to 1 and the nugget to 0. This is because the normal initial values for those paramters don't work well with the power model. I consider this a temporary solution, any suggestions are appreciated.

It is possible to pass anisotropy parameters to autofitVariogram. However, autofitVariogram does not fit anisotropic variogram models. The function sees the anisotropic sample variogram as one big sample variogram. So it fits an average isotropic variogram model from the anisotropic sample variogram. A warning is issued when a users passes alpha to autofitVariogram.

An object of type autofitVariogram is returned. This object contains the experimental variogram, the fitted variogram model and the sums of squares (sserr) between the sample variogram and the fitted variogram model.

Paul Hiemstra

fit.variogram, autoKrige, posPredictionInterval

AybarCL/Dorado documentation built on Jan. 13, 2020, 1:13 p.m.