r - 将 Robust SE 用于 gee 包中的 plot_models

Question

我注意到使用包 sjPlot 中的 plot_models 会根据朴素标准错误给出置信区间。我希望它使用 Robust SE。有简单的解决方法吗？

Answer 1

目前，sjPlot 不支持此选项，但是，它计划在即将到来的更新中。sjPlot使用parameters包来计算模型参数 - 如果你不介意从 GitHub更新参数包（并安装see包），你已经可以使用这个功能：

library(parameters)
library(gee)

data(warpbreaks)
model <- gee(breaks ~ tension, id = wool, data = warpbreaks)
#> Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27
#> running glm to get initial regression estimate
#> (Intercept)    tensionM    tensionH 
#>    36.38889   -10.00000   -14.72222

mp <- model_parameters(model)

mp
#> Parameter   | Coefficient |   SE |          95% CI |     z | df |      p
#> ------------------------------------------------------------------------
#> (Intercept) |       36.39 | 2.80 | [ 30.90, 41.88] | 12.99 | 51 | < .001
#> tension [M] |      -10.00 | 3.96 | [-17.76, -2.24] | -2.53 | 51 | 0.015 
#> tension [H] |      -14.72 | 3.96 | [-22.48, -6.96] | -3.72 | 51 | < .001

plot(mp)

mp <- model_parameters(model, robust = TRUE)

mp
#> Parameter   | Coefficient |   SE |          95% CI |     z | df |      p
#> ------------------------------------------------------------------------
#> (Intercept) |       36.39 | 5.77 | [ 25.07, 47.71] |  6.30 | 51 | < .001
#> tension [M] |      -10.00 | 7.46 | [-24.63,  4.63] | -3.94 | 51 | 0.186 
#> tension [H] |      -14.72 | 3.73 | [-22.04, -7.41] | -1.34 | 51 | < .001

plot(mp)

^{由reprex 包（v0.3.0）于 2019 年 12 月 23 日创建}