在 GAM mgcv 中运行 Betar distriubtion 时出错

问题描述

我正在尝试使用 mgcv 包运行 gam，其中响应变量是比例数据。数据过度分散，所以最初我使用了拟二项分布。但是，因为我使用的模型选择不是特别有用，因为它不会产生 AIC 分数。

相反，我正在尝试使用 beta 发行版，因为我读到它可能是合适的。

mRI_br <- bam(ri ~ SE_score + s(deg,k=7) + s(gs,k=7) + TL + species + sex + season + year + s(code,bs = 're') + s(station,bs = 're'),family=betar(),data=node_dat,na.action = "na.fail")

但是，我在运行模型时收到此警告。

Warning messages:
1: In estimate.theta(theta,family,y,mu,scale = scale1,... :
  step failure in theta estimation

当我尝试检查模型摘要时，我收到此错误。

> summary(mRI_br)
Error in chol.default(diag(p) + crossprod(S2 %*% t(R2))) : 
  the leading minor of order 62 is not positive definite

我想知道：

1. 导致这些错误和警告的原因是什么，如何解决？
2. 如果没有，还有任何其他分布可以与比例数据一起使用，使我能够随后使用模型选择技术（例如 MuMIn 包中的疏浚功能。

可以找到数据集的副本here

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）