您如何报告 Beta 回归(R 输出)的结果?

机器算法验证 r 报告 贝塔回归
2022-03-26 23:29:54

我正在寻找有关如何报告 beta 回归输出结果的建议/输入。我的数据以比率形式限制在 0 和 1 之间,我正在研究响应变量 (D_Ratio) 和连续的预测变量 (体长或 BL) 变量之间的简单关系。我使用了 R 中 betareg 包中的 betareg 函数。

例如,这是我的 R 输出:

Call:
    betareg(formula = D_Ratio ~ BL, data = wild, link = c("cloglog"))

    Standardized weighted residuals 2:
        Min      1Q  Median      3Q     Max 
     -1.4137 -0.6463 -0.1782  0.3970  2.6160 

Coefficients (mean model with cloglog link):
            Estimate Std. Error z value Pr(>|z|)    
(Intercept) -2.14147    0.51930  -4.124 3.73e-05 ***
     BL      0.05252    0.01673   3.139  0.00169 ** 

Phi coefficients (precision model with identity link): 

    Estimate Std. Error z value Pr(>|z|)
(phi)1.9522     0.2969   6.576 4.82e-11 ***
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Type of estimator: ML (maximum likelihood)
Log-likelihood: 8.766 on 3 Df
Pseudo R-squared: 0.2058
Number of iterations: 13 (BFGS) + 1 (Fisher scoring) 

首先,我注意到有两个表格需要考虑;来自平均模型链接的系数和来自精度模型的系数。我要报告哪些系数?我在其他线程中找到了不同的答案......现在我认为它应该是来自平均模型的伪 R 平方、Z 值、P 值......或者“估计”系数项是否意味着一些重要的东西,比如斜率? 我问是因为我的印象是这种关系不是一条直线。

不幸的是,我是一个相对较新的 R 用户,所以如果这里有编码问题,请告诉我。

1个回答

beta回归模型可以有两个子模型:(1)均值回归模型——类似于线性回归模型或二元回归模型;(2) 精度参数的回归模型——类似于线性回归模型中方差的倒数或 GLM 中的离散度。

到目前为止,您刚刚在 (1) 中使用了回归量,但在 (2) 中只使用了一个常数。我鼓励您检查在两个部分中都D_Ratio ~ BL | BL带有回归量的模型是否BL会导致更好的拟合。

如果不是,那么您可能最好从平均方程中报告系数,就像对二元回归模型一样。然后您可以添加精度参数估计(就像在线性回归中一样)、伪 R 平方和/或对数似然和/或 AIC/BIC。

如果回归量在模型的两个部分都起作用,那么可能会报告两组系数。

您还可以使用包中的函数mtable(betareg_object,...)memisc生成这样的表。也可以导出到 LaTeX。此外,您可能会考虑D_RATIO ~ BL使用拟合平均回归线加上可能的一些分位数(例如,5% 和 95%)的散点图。vignette("betareg", package = "betareg")有一些这样的例子。