This function calculates the liklelihood of the geostatistical model in presence of positional error using Quasi Monte Carlo Integration. It is used internally in the optimisation routine.
qmci(data, sequence, sigma2, phi, nugget, mu, delta, kappa)