Проблема настройки линейной модели с помощью JAGS - PullRequest
Новая
       55

Проблема настройки линейной модели с помощью JAGS

0 голосов
/ 01 декабря 2019

Я пытаюсь настроить линейную модель с помощью JAGS, но у меня проблемы с кодом. Я пишу: 

library(R2jags)
library(BEST)

base<-data.table::data.table(read.csv("/Users/franco/Documents/Todo/UNAM/Facultad\ de\ Ciencias/Asignaturas\ Actuaría/Análisis\ Bayesiano\ de\ Datos/Tareas/Tarea-Examen\ 2/FootballLeague.csv"))
X <- cbind(1,as.matrix(base[,-c(1,2,12)]))
y <-as.matrix(base[,2])
n <- length(y)
m <- ncol(X)

model.jags <- function(){
  tau ~ dgamma(0.01, 0.01)
  for(i in 1:m){
    beta[i] ~ dnorm(0,0.001)
  }

  for (i in 1:n){

    y[i] ~ dnorm(x[i,]%*%beta,tau)

  } 

  sigma <- pow(tau,-1)
}
jags.params <- c("beta","sigma")
jags.modelo <- jags(model.file=model.jags,parameters.to.save=jags.params,
                    data = list('n' = n,
                                'y' = y,
                                'x' = X,
                                'm'=m),
                    n.chains = 2,
                    n.thin=1,
                    DIC=FALSE,
                    n.burnin = 10000,
                    n.iter = 20000)

И R выдает эту ошибку:

Ошибка в jags.model (model.file, data = data, inits = init.values, n.chains = n.chains,: RUNTIME ERROR: Ошибка компиляции в строке 8. Несоответствие размеров в подмножестве выражения y.

Я не знаю, в чем заключается ошибка: / Может кто-нибудь мне помочь,пожалуйста.

...