该模型的代码是,
data {
int<lower=0> N;
int<lower=0> K;
matrix[N,K] y;
matrix[N,K] X;
}
parameters {
matrix[K,K] lambda_hat;
matrix[K,K] sigma_0;
matrix[K,K] L;
}
transformed parameters {
matrix[N,K] mu;
vector[4] v;
for (i in 1:K)
v[i]=1;
for (i in 1:N)
mu[i,]'=lambda_hat*x[i,]';
}
model {
L~ wishart(K,diag_matrix(v));
sigma_0~ wishart(K,L);
lambda_hat~ wishart(K,L);
for (t in 1:N)
y[t,]'~ multi_normal(y[(t-1),]'+mu[t,]',sigma_0);
}
运行此代码时出现以下错误,
rt = stanc("mystancode.stan")
得到的错误是,*SYNTAX ERROR, MESSAGE(S) FROM PARSER: Illegal statement started with non-void expression parsed as transpose(mu[i, :]) 不是合法的赋值、采样或函数语句。注意
- 赋值语句只允许左侧的变量(带有可选索引);
- 采样语句允许在左侧使用任意值表示表达式。
- 用作语句的函数必须声明为具有 void 返回
第 20 行第 7 列的“model2e850c17b4_mystancode”错误
18:
19: for (i in 1:N)
20: mu[i,]'=lambda_hat*x[i,]';
^
21:
解析器预期:stanc 中的错误(“mystancode.stan”):由于上述错误,无法解析 Stan 模型“mystancode”。*