我想用' glmulti'进行模型选择。封装
我使用以下代码:
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)
我收到以下错误:
Error in fn(x, ...) : Downdated VtV is not positive definite
什么是VtV?