Cubic correlation function.
If length(theta) = 1, then the model is isotropic:
all theta_j = theta.
corrcubic(theta, d, ret = "all")parameters in the correlation function
m*n matrix with differences between given data points
A string. If set to "all" or "dr", the derivative of r (dr) will be returned, else dr is NA.
returns a list with two elements:
rcorrelation
drm*n matrix with the Jacobian of r at x. It is
assumed that x is given implicitly by d[i,] = x - S[i,],
where S[i,] is the i'th design site.