このモデルの後部からサンプリングしようとしています: ここに画像の説明を入力してください
x は 10x1 ベクトル、mu は 10x1 ベクトル、sigma は 10x10 行列、psi_0 は 10x10 行列、太字の 1 は 5x1 ユニティ ベクトル、残りはスカラーです。F と E は、分析における 2 つの株式グループの添え字です (ポートフォリオの最適化を行っています)。次のコードを実行しようとしています。
stancode <- "data {
int<lower=1> T;
matrix[T,10] data_total;
real<lower=0> sigma_squared;
real k;
real<lower=0> r_squared;
matrix[10,10] psi;
}
parameters {
vector[10] mu_var;
matrix[10,10] Sigma_var;
real alpha_energy_var;
real alpha_fin_var;
real beta;
}
transformed parameters {
vector[10] mu_mean;
mu_mean = [alpha_fin_var,alpha_fin_var,alpha_fin_var,alpha_fin_var,alpha_fin_var,alpha_energy_var,alpha_energy_var,alpha_energy_var,alpha_energy_var,alpha_energy_var]';
}
model {
target += normal_lpdf(beta | k, sqrt(r_squared));
target += inv_wishart_lpdf(Sigma_var | 12, psi);
target += normal_lpdf(alpha_energy_var | beta, sqrt(sigma_squared));
target += normal_lpdf(alpha_fin_var | beta, sqrt(sigma_squared));
target += multi_normal_lupdf(mu_var | mu_mean, Sigma_var);
target += multi_normal_lupdf(data_total | mu_var,Sigma_var);
}"
stanmodel <- stan_model(model_code = stancode, model_name="stanmodel")
(data_total はデータ マトリックスであり、"_var" サフィックスは、コードの残りの部分との一貫性のために追加されています)
最後の行でエラーが表示されます: PARSER EXPECTED: "(", モデル コードの最後から 2 番目の行に関して。私はスタックで、この種のエラーが発生してより多くの専門家の開発者を困惑させているのを見てきました。
誰でも私を助けることができますか?何が悪かったのかよくわかりません。