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) 

}" 

내가 다음 컴파일하고 데이터로드 : handleRes (고해상도)에서 오류 : 나는 오류가이 시점에서

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) 

을 여기에 버그에서 설정 한 모델이다 NA

내가 어디로 잘못 갔는지 ...에 대한 아이디어가 있습니까?

Answer 1

나는 단지 의견을 말할 수는 있었지만 나는 충분한 평판이 없다고 생각한다. 어쨌든 일부 데이터를 구성하고 모델의 단순화 된 버전을 실행했습니다 (예 : 연령 및 분만 사용). OpenBUGS의 GUI 버전에서 실행했는데 작동하지 않았습니다. 나는 당신의 계수를 dnorm(0,0.01으로 바꾸었다. 그래서 전임자를 변경하는 것이 좋습니다. 당신이 가진 것들은 매우 모호합니다. 제안 된대로 변경해도 추측에는 영향을 미치지 않으므로 모델을 실행할 수 있습니다.