0

我正在使用 MuMin 包中的疏浚函数进行具有随机效果的游戏:

下面是使用包 mgcv 中的 bam 的全局模型。随机效应是s(StationID, bs = 're')

gam_global<-bam(Numberpertow ~ s(interval,k=6)+Stratum+ClosArea+CruiseID+s(interval, by=CruiseID,k=6)+offset(offset)+s(StationID, bs = 're'),data=l.data,method = "REML",family=nb(),na.action = "na.fail")

当我在全局模型上运行疏浚函数时,随机效应会从模型的子集中删除,因为疏浚正在运行带有和不带有某些变量的不同模型选项。我对确定是否需要随机效应并需要将其作为随机效应保留在所有模型中不感兴趣。我还需要在所有模型中保留偏移项。

进行网络搜索后,我根据这篇文章https://stackoverflow.com/questions/61064080/fixed-terms-in-dredge-function-mumin-in-gam[enter link description here] 1尝试了以下操作

我尝试了以下两种方法均未成功:

#genertric random effect name as shown in model summary output

gam_dredge<-dredge(gam_global,evaluate=TRUE,fixed=c("offset(offset)","s(StationID)"),rank="AIC")

#random effect name as in model

gam_dredge<-dredge(gam_global,evaluate=TRUE,fixed=c("offset(offset)","s(StationID,bs="re)"),rank="AIC")

我还将 bam 模型中的方法从REML更改为GCV.Cp

我对挖泥机的任何迭代都没有成功,并且在挖泥机卡盘运行时得到以下结果:

Fixed terms are "offset(offset)" and "(Intercept)"
4

1 回答 1

0

空白很重要。模型术语名称是's(StationID, bs = "re")'

根据?dredge,“出现在fixed并且subset必须采用 getAllTerms(global.model) 返回的确切形式的术语名称”。出于这个原因,您引用的接受的 SO 答案是不正确的(这就是为什么手册应该始终是主要信息来源的原因)。

但是,您可以通过给出fixed公式来使事情变得更容易,在这种情况下,空格是无关紧要的(在字符串之外):

~ offset(    offset) + s  (StationID,      bs  ="re"  )
于 2022-01-11T11:10:50.083 回答