error de sintaxis misteriosa [...] cerca de "datos" en jags

votos
0

Cuando se ejecuta mi modelo jags me sale el siguiente mensaje de error

Error en el archivo modelo de análisis: error de sintaxis en la línea 5, cerca de datos

En resumen, lo que tengo es de dos variables de respuesta ( cmaxy cmd) y una variable predictora ( dbh). Mi idea es estimar la correlación entre cmaxy cmdque no se explica por dbh, desde los simples correlaciones entre estas variables son triviales.

A continuación se muestra el código:

# Model input
dat = list(N = nrow(data), # number of observations
       data = as.matrix(log(data[,c(dbh,cmax,cmd)])), # log-transform variables
       T = diag(2)/1000, # var-covar matrix for non-informative priors
       r = 2, # number of variables
       m = c(0,0)) # means for non-informative priors
inits = list(P=diag(2)/1000, A=c(0,0), B=c(0,0))

# JAGS model
cat(model{

### Likelihood
for(i in 1:N){
  M[i,1:r] <- A[1:r] + B[1:r]*data[i,1]
  data[i,2:3] ~ dmnorm(M[i,1:r],P[1:r,1:r])
  }

### Priors
P[1:r,1:r] ~ dwish(T,r)
A[1:r] ~ dmnorm(m,T)
B[1:r] ~ dmnorm(m,T)

### Statistics
V <- inverse(P)
sigmaH <- sqrt(V[1,1])
sigmaW <- sqrt(V[2,2])
covHW <- V[1,2]
corHW <- covHW/(sigmaH*sigmaW)

},
  file=Ch1/BM2.txt)

# Run JAGS
res = jags.model(file=Ch1/BM2.txt, data=dat, inits=inits, n.chains=1, 
             n.adapt=500)
Publicado el 07/11/2018 a las 22:32
fuente por usuario
En otros idiomas...                            


1 respuestas

votos
0

Por lo tanto, cuando vine a saber que había otro módulo ( "datos"), aparte de "modelo" en un archivo de modelo jags, empecé a sospechar que llamar a mis "datos" de entrada trama de datos no era la idea más inteligente. Después de cambiar el nombre de esta trama de datos todo parece haber funcionado bien.

Respondida el 09/11/2018 a las 03:12
fuente por usuario

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more