Usage
gp_mle(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
Value
negative log-likelihood value, or maxval if the computation fails