Rstan：如何读入从 R 到 Stan 的多个向量

Question

假设我在 Stan 中有一些向量

vector[100] mu[20]

这意味着mu是 20 个向量，每个向量的长度为 100。

我是从 R 读为mu = martix(NA,nrow=100,ncol=20)还是读为mu = array(NA,dim=c(100,20)) ？ （我希望我没有混淆尺寸）

或者是不正确的，你这样做有什么不同吗？

Answer 1

因此，您有两种选择来提供matrix / array或向量list 。 在您的示例中，您必须更改行和列的顺序。

rstan错误消息很清楚，因此您将知道尺寸是否有问题。

library(rstan)
# dummy parameters
scode <- "
data {
  vector[100] mu[20];
}
parameters {
  real y[2]; 
}
"
# Both fit pass without errors
fit1 <- stan(
  model_code = scode,
  data = list(mu = matrix(rnorm(20 * 100), nrow = 20, ncol = 100))
) 
fit2 <- stan(
  model_code = scode, 
  data = list(mu = lapply(1:20, function(x) rnorm(100)))
)