# NOT RUN {
data("jura", package="gstat")
X = jura.pred[,1:2]
Zc = compositions::acomp(jura.pred[,7:9])
vg = logratioVariogram(data=Zc, loc=X, azimuth=c(0:18)*10,
azimuth.tol=22.5)
image(vg)
image(vg, jointColor=TRUE)
# }
Run the code above in your browser using DataLab