具有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的运行方式？ 谢谢！！ 关于Fran

Answer 1

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

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