使用JAGS调整线性模型时遇到问题

时间:2019-11-30 21:07:34

标签: r bayesian jags

我正在尝试使用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 ,:运行时错误:第8行的编译错误。   y子集表达式中的维不匹配。

我不知道这是什么错误:/请有人帮我。

0 个答案:

没有答案