powered by
evaluation of the likelihood
vecchia_likelihood(z, vecchia.approx, covparms, nuggets, covmodel = "matern")
(multivariate normal) loglikelihood implied by the Vecchia approximation
the observed data
a vecchia object as generated by vecchia_specify()
covariance parameters as a vector
either a single (constant) nugget or a vector of nugget terms for the observations
covariance model, 'matern' by default
z=rnorm(5); locs=matrix(1:5,ncol=1); vecchia.approx=vecchia_specify(locs,m=3) vecchia_likelihood(z,vecchia.approx,covparms=c(1,2,.5),nuggets=.2)
Run the code above in your browser using DataLab