问题标签 [rstanarm]

我写信是为了找出对象"linear.predictors"返回的 内容stan_glm()。

显然，"linear.predictors"与用户提供的预测变量不同（文档没有帮助）。

无论如何，有没有办法从stan_glm()对象中获取预测值？

这是一个单一的预测器（即 mom_iq）示例：

我想知道是否有可能（并且可能推荐）从包中获得标准化系数？（在文档中没有找到任何具体内容）stan_glm()rstanarm

我可以像正常回归一样标准化所有变量吗？（见下文）

例子：

标准化：

stan_glm()我想知道如何根据包中的分组变量获得后验预测rstanarm？

例如，如果我的数据中有一个二进制(0, 1)编码的分组变量"vs"（基本 R 数据：）mtcars，我如何获得何时vs == 0和何时的预测vs == 1？

这是我的 R 代码：

背景： McEearth（2016 年）在他的重新思考书籍第 158-159 页中，使用索引变量而不是虚拟编码来对称为“clade”的 3 类变量进行预测“kcal.per.g”（线性回归）。

问题：我想知道我们是否可以在 中应用相同的方法"rstanarm"？我为下面的可能演示提供了数据和 R 代码。

我正在使用出色的绘图库bayesplot来可视化我正在估计的模型的后验概率区间rstanarm。我想通过将系数的后验间隔放到同一个图上来以图形方式比较来自不同模型的绘图。

例如，想象一下，beta1, beta2, beta3对于两个不同的模型，我有 1000 次后验图，用于三个参数：

Bayesplot 为单个模型绘制进行了出色的可视化，它是 ggplot2 '在引擎盖下'，所以我可以随意定制：

但是我想要实现的是将这两个模型一起绘制在同一个图上，这样对于每个系数我都有两个区间，并且可以通过将颜色映射到模型来区分哪个区间是哪个区间。但是我无法弄清楚如何做到这一点。一些不起作用的东西：

关于我如何做到这一点的任何想法？或者如何在给定后绘制矩阵的情况下手动完成？

^{由reprex 包（v0.2.1）于 2018 年 10 月 18 日创建}

我正在拟合一个模型（为了重现性而非常简化），如下所示：

库（数据集） 库（rstanarm）

我想我不应该期望两个相等的“模拟”数据集（例如 0.2216209 和 0.2177802）有相同的中位数？然而，上面的骨架代码是否是模拟不同场景的正确方法（例如 Petal.Length = 1.4 与 Petal.Length = 2.4）？

是否有可能在 rstanarm 中提取用于 MCMC 采样的 stan 代码？

我想将我自己的模型参数化和先前的选择与 rstanarm 中使用的进行比较。

我将不胜感激使用修改后的 brms 生成的 stan 模型更新我的 brmsfit 对象，因为我想以 brms 可能尚不支持的方式将各种权重列传递给可能性。

我的目标是在对此处讨论的权重分布进行边缘化后获得效果的后验分布。

我需要在 brms 或 stanarm 中而不是直接在 stan 中执行此操作，因为我想使用stanfit 对象当前不支持的https://github.com/mjskay/tidybayes的功能。

这个想法可能类似于使用修改后的 brms 生成的 Stan 模型创建 brmsfit 对象中讨论的想法，但我不确定在实践中如何做到这一点，因此非常感谢任何帮助。

此外，我不明白如何在实践中实施@paul.buerkner解释的第二行和第三行：

通过 brm(...,chains = 0) 创建模型的“空” brmsfit 对象。生成并修改 Stan 代码并将 slot fit 替换为 rstan::stan_model 的输出（您可能还想更改 slot model 中的 Stan 代码）。运行更新（，重新编译 = FALSE）

#要修改和更新的brms模型

#修改后的 brms 生成的 stan 模型modified_brms.stan

#我能够成功地做到这一点，但不确定如何从这里开始：

这个问题也已经发布在这里。提前感谢您的帮助。

我将不胜感激在下面指定我的 brms 模型以便能够将多列权重传递给模型，stan code如下所示。

我需要在 brms 或 stanarm 中执行此操作，而不是直接在 stan 中执行此操作，因为我想使用stanfit 对象当前不支持的https://github.com/mjskay/tidybayes的功能。

＃样本数据：

#我的尝试：

#斯坦代码：

提前感谢您的帮助。

我需要将数据作为具有不同长度的向量列表和具有相同行数但具有不同列数的矩阵列表传递。有没有办法在 Rstan 中传递数据？