load.model <- function () {
    list(name="lmer.LOGLIN.MAT.MAP.intra.r.set.species",
         var.BLUP = 'set.id',
         lmer.formula.tree.id=formula("logG~1+logD+sumBn.inter+MAT+MAP+sumBn.intra+(1+logD||species.id)+(1|plot.id)+(1+sumBn.inter+sumBn.intra||set.id)"))
}