r - 带有 beta 家族的 glmmTMB

Question

我正在尝试对取值从 0 到 0.7 的光合响应 (fvfm) 进行建模，我正在使用 beta 分布，因为我的通货膨胀率为零，并且我考虑了一个分散模型，因为我的值随着时间推移远离平均值，但是在运行这个模型时：

beta7<-glmmTMB(fvfm~Temperature*dia+(1|ID), zi=~1, disp=~dia, data=ds, family=beta_family(link="logit"))

这出现了：

Error in eval(family$initialize) : y values must be 0 < y < 1

我添加了响应变量 0.001 的所有值，因此我不再有零通货膨胀并且模型有效。这是对的？？或者有什么方法可以在不添加 0.001 的情况下运行？谢谢！！问候弗兰

Answer 1

如果您的结果包括 0，我认为零膨胀 beta 可能会起作用 - CRAN (1.0.0) 上的最新 glmmTMB 版本能够适应零膨胀 beta 回归。

否则，您可能还会考虑使用normalize()-function并设置include_boundsto 来“压缩”您的响应变量FALSE。