data(KNMIdata)
locations <- KNMIdata$coord
pairIndices <- selectPairIndices(locations, maxDistance = 0.5)
Mestimator(KNMIdata$data, locations, pairIndices, k = 50, model="BR",
iterate=FALSE, covMat = FALSE)$thetaRun the code above in your browser using DataLab