zero mean version

gp_mle2(theta, d, nugget, myY, maxval = Inf)

Arguments

theta

proposed parameter vector: [mu, tauw, tauv, phi1...phiK]

d

spatial distance matrix

nugget

allows additional error in Y rather than fix interpolation to go through points

myY

vector of observed data

maxval

maximum value