# NOT RUN {
data("jura", package="gstat")
juracomp = compositions::acomp(jura.pred[, -(1:6)])
lrvg = logratioVariogram(data=juracomp, loc=jura.pred[,1:2])
mf = Maf(juracomp, vg=lrvg)
mf
compositions::coloredBiplot(mf, xlabs.col=as.integer(jura.pred$Rock)+2)
# }
Run the code above in your browser using DataCamp Workspace