R OpenBugsモデルエラー

Question

私は、ホッケー目標のデータを扱うために階層ポアソンモデルを実行することを検討しています。バグに設定してこちらのモデルは次のとおりです。

modelString <- "model { 

    for(i in 1:n_obs){ 

    hockey_goals[i] ~ dpois(p[i]) 


    log(p[i]) <- p_inter + p_age * age[i] + p_sv_pct * sv_pct[i] + p_team * team[i] + p_win_pct * win_pct[i] + log(n_mins[i]) 

    } 
    p_inter ~ dnorm(0,0.00001) 
    p_age ~ dnorm(0, 0.00001) 
    p_sv_pct ~ dnorm(0, 0.00001) 
    p_team ~ dnorm(0, 0.00001) 
    p_win_pct ~ dnorm(0, 0.00001) 

}" 

私は、データをコンパイルしてロード：

season_goals <- data$GA 
n_mins <- data$MIN 
age <- data$Age 
sv_pct <- data$SV. 
team <- data$Tm 
win_pct <- data$W/data$GP 
data <- list(n_obs=length(season_goals),n_mins=n_mins,hockey_goals=season_goals,age=age, 
sv_pct=sv_pct,team=team,win_pct=win_pct) 

# Get the data into BUGS: 
modelData(bugsData(data)) 
#------------------------------------------------------------------------------ 
# INTIALIZE THE CHAINS. 

nchain = 1 
modelCompile(numChains=nchain) 
modelGenInits() 

#------------------------------------------------------------------------------ 
# RUN THE CHAINS. 
samplesSet(c("p_age","p_sv_pct","p_team","p_win_pct")) 
# R command defines a new variable that specifies an arbitrary chain length: 
chainLength = 10000 
# BRugs tells BUGS to generate a MCMC chain: 
modelUpdate(chainLength) 

この時点で、私はエラーを取得する：handleResで エラー（RES）：NA

私がどこに間違っていたかに関するアイデア...？

Answer 1

私はちょうどコメントすることができますが、私は十分な評判を持っていません。とにかく、いくつかのデータを作り、モデルの簡略化されたバージョンを実行しました（つまり、再生された年齢と分のみ）。私はOpenBUGSのGUIバージョンでそれを実行し、それは動作しませんでした。私はあなたの係数をdnorm(0,0.01に変更して更新しました。 それで、私は前任者の変更を提案します。あなたが持っているものは非常に曖昧です。あなたの推論に影響することはありません示唆されたそれらを変更し、実行するモデルを得ることがあります。