如何解决在 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
我想知道:
可以找到数据集的副本here
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。