r - 如何从 stan_glm 中的系数中提取标准误差

Question

我有一个看起来像这样的模型：

fit <- stan_glm(switch ~ arsenic + dist100, 
                family=binomial(link="logit"), 
                data=wells, 
                refresh=0)

输出如下所示：

stan_glm
 family:       binomial [logit]
 formula:      switch ~ arsenic + dist100
 observations: 3020
 predictors:   3
------
            Median  MAD_SD 
(Intercept)  0.0054  0.0810
arsenic      0.4618  0.0415
dist100     -0.9025  0.1040

------

我想要得到的拟合的 MAD_SD 系数。我通常使用 se(fit) 这给了我系数，但由于某种原因我得到了错误

Error in se(fit) : argument "predicted" is missing, with no default

我该如何纠正这个错误？

标签： rglmstan

Answer 1

如果有人想知道这个问题的答案，stan_glm 对象基本上就像将名称映射到元素的大标题。在这种情况下，元素被命名为“ses”并且它是一个列表。所以你可以做

fit$ses[0]

索引它

编辑：事实证明，我完全需要这样做的原因是因为我在 stan_glm 之后导入了一个包，该包也有一个名为“se”的函数并且它被屏蔽了。