R Ошибка модели OpenBugs

Я хочу запустить иерархическую модель Пуассона для данных о хоккейных воротах. вот модель, настроенная в ошибках:

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

Любые идеи о том, где я ошибся...?


person jmy    schedule 04.05.2017    source источник


Ответы (1)


Я хотел бы просто прокомментировать, но у меня недостаточно репутации. В любом случае, я собрал некоторые данные и запустил упрощенную версию вашей модели (то есть только с возрастом и сыгранными минутами). Я запустил его в версии OpenBUGS с графическим интерфейсом, и это не сработало. Затем я изменил априорные коэффициенты для вас на dnorm(0,0.01), и они обновились. Так что я бы посоветовал сменить приору. Те, что у вас есть, очень, очень расплывчаты. Изменение их, как предлагается, не повлияет на ваш вывод, и вы можете заставить модель работать.

person amg    schedule 09.05.2017