# 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
biplot(mf)
predict(mf)
unclass(predict(mf)) - unclass(juracomp) # predict recovers the original composition
# }
Run the code above in your browser using DataLab