随机效应的MuMin疏浚gam

Question

我使用MuMin包中的疏浚函数对具有随机效果的gam：

下面是使用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](https://stackoverflow.com/questions/61064080/fixed-terms-in-dredge-function-mumin-in-gam%5Benter) link description here的方法

我尝试了以下两种方法，但都没有成功：

#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)"

Answer 1

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

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

但是，您可以通过将fixed作为一个公式来简化，在这种情况下，空格是不相关的(字符串之外)：

~ offset(    offset) + s  (StationID,      bs  ="re"  )