Question

我试图向我们打包＆＃34; AICcmodavg＆＃34;使用函数＆＃34; glmer＆＃34;在一组候选混合模型中进行选择。在包＆＃34; lme4＆＃34;下使用二项式链接功能。但是，当我尝试运行＆＃34; aictab＆＃34;在这个模型列表上的函数，我收到此错误消息＆＃34; formatCands中的错误（cand.set）：函数不支持模型类和＃34;的混合。我的代码发布在下面。

Table<- read.delim("C:/Analysis/Data Files/Text Files/QuadSEFull.txt",  header=TRUE)
Table<- subset(Table,Species!="None")

library("lme4", lib.loc="C:/Users/mcs7gb/Documents/R/win-library/3.0")
library("AICcmodavg",lib.loc="C:/Users/mcs7gb/Documents/R/win-library/3.0")

y<-cbind(Table$Visual,Table$Excavation)

mod<-list()

mod[[1]]<-glmer(y~(1|Species)+Discharge,family=binomial,data=Table)
mod[[2]]<-glmer(y~(1|Species)+Visibility,family=binomial,data=Table)
mod[[3]]<-glmer(y~(1|Species)+Depth,family=binomial,data=Table)
mod[[4]]<-glmer(y~(1|Species)+AvgDepth,family=binomial,data=Table)
mod[[5]]<-glmer(y~(1|Species)+MaxDepth,family=binomial,data=Table)
mod[[6]]<-glmer(y~(1|Species)+DomSub,family=binomial,data=Table)
mod[[7]]<-glmer(y~(1|Species)+SubDomSub,family=binomial,data=Table)
mod[[8]]<-glmer(y~(1|Species)+Embed,family=binomial,data=Table)
mod[[9]]<-glmer(y~(1|Species)+CoreArea,family=binomial,data=Table)
mod[[10]]<-glmer(y~(1|Species)+BedArea,family=binomial,data=Table)
mod[[12]]<-glmer(y~(1|Species)+Temp,family=binomial,data=Table)
mod[[14]]<-glmer(y~(1|Species)+Bottom_Velocity,family=binomial,data=Table)
mod[[15]]<-glmer(y~(1|Species)+AvgVelocity,family=binomial,data=Table)
mod[[17]]<-glmer(y~(1|Species)+Macrophytes,family=binomial,data=Table)
mod[[18]]<-glmer(y~(1|Species)+WoodyDebris,family=binomial,data=Table)
mod[[19]]<-glmer(y~(1|Species)+Stream,family=binomial,data=Table)
mod[[20]]<-glmer(y~(1|Species)+JulDate,family=binomial,data=Table)
mod[[21]]<-glmer(y~(1|Species)+Site,family=binomial,data=Table)
mod[[22]]<-glmer(y~(1|Species)+Discharge+Visibility+Depth+AvgDepth
      +MaxDepth+Table$DomSub+SubDomSub+Embed+CoreArea+BedArea
      +Temp+Bottom_Velocity+AvgVelocity+Table$Macrophytes+WoodyDebris
      +JulDate+Site+Stream,family=binomial,data=Table)
Modnames<- paste("mod", 1:length(mod), sep = " ")
z<-aictab(cand.set=mod,modnames=Modnames,second.ord=TRUE,nobs=NULL,sort=TRUE)

带有glmer的R AIC型号选择表（包lme4）

0 个答案: