Trying to run a model with all 6 terms does not work as the model is rank-deficient: #library(MuMIn) # not runįorm<-formula(paste0('Y~',paste0(vars,collapse='+'),'+(1|group)'))įmod<-glmer(form,data=dfdat,family='poisson')įixed-effect model matrix is rank deficient so dropping 1 column / coefficientÄ®rror: pwrssUpdate did not converge in (maxit) iterations To give a reproducible example, let's generate a random data set: dfdat<-factor(sample(1:2,nrow(dfdat),replace=T))Äfdat$Y<-rpois(nrow(dfdat),10)+rpois(nrow(dfdat),as.numeric(dfdat$group)) However, to use dredge() effectively one needs to have a full model including all terms - which can lead to full model to be rank-deficient. Since several variables are collinear I use the subsetting function of dredge to avoid correlated variables. I am trying to run variable selection on Poisson mixed-effect models using glmer() and dredge().
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |