在GAMLSS的BEINF回归中获取newdata均值的置信区间

我使用gamlss包的BEINF（零-一膨胀beta分布）族，对健康结局**出院存活天数占比（%DAOH）**建模。该结局在0和1处存在点质量，且聚类于最大值附近，目前用marginaleffects工具开展组间比较。

当尝试在newdata上预测结局的期望均值时，BEINF参数化下的均值计算公式为(mu+tau)/(1+nu+tau)。我能手动算出均值，但无法顺利获取对应的置信区间——因为marginaleffects需要通过what =参数分别预测mu、sigma、nu、tau四个参数，无法直接给出均值的置信区间。

delta法理论上可以解决这个问题，但我没能成功实现。另外，STATA和SAS中获取零-一膨胀beta分布结局的预测均值置信区间操作很简便；bootstrap也是可选思路，但我没找到适配自身工作流的便捷实现方式。

有没有人尝试过从GAMLSS膨胀beta族模型中获取均值的置信区间？求可行的实现思路！

内容的提问来源于stack exchange，提问作者Kelvin