# load Quatre Montagnes dataset
data(quatre_montagnes)
# build ABA model for basal area, with three metrics as predictors
model_aba <- aba_build_model(quatre_montagnes$G_m2_ha,
quatre_montagnes[, c("zpcum8", "ipcumzq70", "p_hmin")],
transform = "log", nmax = 2
)
# plot field values VS predictions in cross-validation
aba_plot(model_aba, main = "Basal area")
Run the code above in your browser using DataLab