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)))
)