问题标签 [rstan]

Stan 提供的功能

用于获取对称矩阵 A 的特征值和特征向量，但是如果您的矩阵是非对称的怎么办？你还能得到特征值和特征向量吗？另外，有没有办法测试我是否得到正确的特征向量？如何测试 Stan 程序的各个部分？

我正在尝试在 Rstan 中运行三级嵌套线性模型，但一直收到错误消息。

我从这个三层嵌套模型中获得了灵感：

http://rstudio-pubs-static.s3.amazonaws.com/64315_bc3a395edd104095a8384db8d9952f43.html

本质上，我重复测量了四年的平均宽带速度数据（总共 696 次观测），嵌套在 174 个地方当局中，最后嵌套在 11 个区域中。

这是数据：

代码：

错误：

我一直在研究一个三级 RStan 模型，其中重复的宽带测量（年份 ID = yrid）嵌套在地方当局（LA ID = 铺设）内，最终嵌套在区域内（区域 ID = rnid）。（记录的）因变量是速度，（记录的）预测变量是人口密度（pd）和超高速宽带渗透（sfbb）。目前在地方当局和地区一级（2 级和 3 级）有随机拦截。

如何扩展模型以在 1 级或 2 级具有随机斜率？

这是数据的一个子集，RStan 模型和整个 R 代码。任何帮助将非常感激。

拉出 stan 代码以提高可读性：

我正在使用 RStan 来拟合两个法线的混合。

我想添加一个条件，即 mu[1] 小于 mu[2]。我怎样才能做到这一点？

谢谢你的帮助

我在让 Rstan 在我的 OSX 设置上工作时遇到了很多麻烦。对于背景信息，我的系统如下：

• OSX 优胜美地 v. 10.10.5
• R 3.2.1 GUI 1.66 小牛

每次通过 R stan 在 R 中运行 stan 模型时，都会出现以下错误（无论我正在运行的模型如何，因此已经检查了 stan 代码的语法错误）

其他相关信息：
我可以在 PyStan 中成功运行测试程序。这似乎主要是 Rstan 的一个问题。

谁能给我任何关于如何解决这个问题的线索？

我有一个字符串，say"ti[123]"和一个对应的值，say1.2 我应该在 R 中使用什么命令来分配1.2给 vector 的第 123 个位置ti？我试过assign("ti[123]",1.2)了，但它创建了一个带有名称的新变量ti[123]，这不是我所期望的。

我需要这样做的原因是我使用 rstan 的optimizing函数来优化模型，它返回一个命名向量，具有这样的名称（如果您碰巧在模型中使用参数向量）。

I'm trying the following model with the lme4 package:

I now try to fit the same model using rstan through the glmer2stan package:

this is the result

I don't think that's the same model. Is my glmer2stan specification wrong? I know that glmer2stan is not actively developed any more but it should handle this simple model, shouldn't it?

UPDATE:

thanks to the tip by Roland I changed the Machine factor levels to dummies and it now works:

我下载了包 'rstan' 用于使用 R 进行贝叶斯工作。我收到了这个错误

我如何使 rstan 工作？

我已经使用 devtools 创建了我的包骨架，然后添加了一堆 R 代码、元数据、文档等。我想在这个包中使用 rstan。我知道 rstan::rstan.package.skeleton 创建了一个包骨架来促进这一点。那么，使用该包中使用 rstan 所需的结构来增强现有包的最佳实践是什么？谢谢你。

根据parallel::detectCores(). 我的防火墙已关闭。我正在stan_lmer从包中加载执行以下调用rstanarm。

确切的错误是：