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