0


R で R2OpenBUGS を実行する際に特定の問題がありますが、OpenBUGS では問題なく動作します。
問題が何であるかを理解したい。

これが私のコードです:

モデルのボラティリティ;
定数 n=180;
{

# 尤度:
(t in 1:n) の ys の同時分布
{
ysigma2[t] <- 1/exp(theta[t]);
y[t] ~ dnorm(0,ysigma2[t]);
}

以前の分布

ムー ~ dnorm(0,0.1);
ファイスター ~ dbeta(20,1.5);
itau2 ~ dgamma(2.5,0.025);
ベータ <- exp(mu/2);
ファイ <- 2*フィスター-1;
タウ <- sqrt(1/イタウ 2);
theta0 ~ dnorm(mu, itau2);
thmean[1] <- mu + phi*(theta0-mu);
theta[1] ~ dnorm(thmean[1],itau2);
for (t in 2:n)
{
thmean[t] <- mu + phi*(theta[t-1]-mu);
theta[t] ~ dnorm(thmean[t],itau2);

}
}




R コード:
svm.sim <- bugs(data, inits, model.file = "C:/Documents and Settings/code.txt", parameters = c("mu", "phi", "tau"), n. chain = 1, n.iter = 1000, codaPkg = TRUE,debug=TRUE)

そして、R から実行したときに OpenBUGS のログに表示されるエラーは次のとおりです。



モデルは構文的に正しい
データをロード した
モデルをコンパイル
しました ノードのこのコンポーネントは確率的ファイ エラー pos 25 では ありません
タイプ UpdaterNormal.StdUpdater のノードの初期値を生成できません
モデルを初期化する
必要があります モニターを使用する前にモデルを初期化する必要があります モニターを使用する前に
モデルを初期化する必要があります
モニターの前に使用 モデル を初期化する必要があります モニターの前に使用モデルを初期化する必要があります DIC をモニターする前に使用モデルを初期化する必要があります モニター
を 更新する前に モデルを初期化する必要があります モニターを使用する前にモデルを初期化する必要があります 使用する DIC モニターが設定され て いませんよろしく








ディナーカル

4

1 に答える 1

0

コードを少し変更して実行したところ (以下を参照)、正常に動作するように見えました。値 n をデータとして指定し、OpenBUGS 自体に初期値を生成させました。init の指定方法に問題がありますか?

    sink("stack_model.txt")
    cat("
    model{
    # likelihood: joint distribution of ys 
    for (t in 1:n) 
    { 
    yisigma2[t] <- 1/exp(theta[t])
    y[t] ~ dnorm(0,yisigma2[t])
  }
  mu ~ dnorm(0,0.1)
  phistar ~ dbeta(20,1.5) 
  itau2 ~ dgamma(2.5,0.025)
  beta <- exp(mu/2)
  phi <- 2*phistar-1
  tau <- sqrt(1/itau2)
  theta0 ~ dnorm(mu,itau2) 
  thmean[1] <- mu + phi*(theta0-mu) 
  theta[1] ~ dnorm(thmean[1],itau2) 
  for (t in 2:n) 
  { 
    thmean[t] <- mu + phi*(theta[t-1]-mu) 
    theta[t] ~ dnorm(thmean[t],itau2) 

  } 
}",fill=TRUE)
sink()
win.data<-list(n=180)
svm.sim <- bugs(win.data, inits=NULL, model.file = "stack_model.txt", parameters = c("mu", "phi", "tau"), n.chains = 2, n.iter = 1000, DIC=F, debug=TRUE,  codaPkg =TRUE)
于 2015-07-28T01:00:30.817 に答える