β族glmmTMB

Question

我试图模拟光合响应(fvfm)，它的值从0到0.7，我使用的是β分布，因为我的通货膨胀为零，我考虑了一个色散模型，因为我的值随着时间的推移而偏离平均值，但是在运行这个模型时：

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-最新的glmmTMB版本的CRAN (1.0.0)能够适应零膨胀的β-回归，那么零膨胀的β可能会起作用。

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