form_glmulti = as.formula(paste("rrF ~ management + mNdep + mMAT + mMAP"))
lmer.glmulti<-function(formula,data,random="",...){
   lmer(paste(deparse(formula),random),data=regdata,REML=F,...)
}
glmulti_lmm<-glmulti(form_glmulti,random="+(1|PlotNR)",data=regdata,method="h",
                     fitfunc=lmer.glmulti, intercept=TRUE,marginality=FALSE,level=2)