JAGS编译错误

Question

我看似直截了当的模型一直给我带来奇怪的编译错误。

运行时错误： 第17行编译错误。 未知变量Nx 要么用数据提供这个变量的值，要么在关系的左边定义它。

这是我的代码：

dataList = list(
  x = x ,
  y = y ,
  xP = xProbe
)

Here's the JAGS model:  
1    data{  
2        Ntotal <- length(y)  
3        ym <- mean(y)  
4        ysd <- sd(y)  
5        for ( i in 1:Ntotal ) {  
6          zy[i] <- ( y[i] - ym ) / ysd  
7        }  
8        Nx <- length(xName)  
9        for ( j in 1:Nx ) {  
10         xm[j]  <- mean(x[,j])  
11         xsd[j] <-   sd(x[,j])  
12         for ( i in 1:Ntotal ) {  
13           zx[i,j] <- ( x[i,j] - xm[j] ) / xsd[j]  
14         }  
15       }  
16       Nprobe <- Ntotal  
17       for ( j in 1:Nx ) {  
18         xPm[j]  <- mean(xP[,j])  
19         xPsd[j] <-   sd(xP[,j])  
20         for (i in 1:Nprobe ) {  
21           zxP[i,j] <- ( xP[i,j] - xPm[j] ) / xPsd[j]  
22         }  
23       }  
24   }  
25   model{  
26       for ( i in 1:Ntotal ) {  
27         zy[i] ~ dt( zbeta0 + sum(zbeta[1:Nx]*zx[i,1:Nx]), 1/zsigma^2, nu)  
28       }  
29      # Priors vague on standardized scale:  
30       zbeta0 ~ dnorm( 0 , 1/2^2 )     
31       for ( j in 1:Nx ) {  
32         zbeta[j] ~ dnorm( 0 , 1/2^2 )  
33       }  
34       zsigma ~ dunif( 1.0E-5 , 1.0E+1 )  
35       nu ~ dexp(1/30.0)  
36      # Transform to original scale:  
37       beta[1:Nx] <- ( zbeta[1:Nx] / xsd[1:Nx] )*ysd  
38       beta0 <- zbeta0* ysd + ym - sum(zbeta[1:Nx]*xm[1:Nx]/xsd[1:Nx])*ysd  
39       sigma <- zsigma*ysd  
40      # Predicted y values as xProbe:  
41       for ( i in 1:Nprobe ) {  
42         zyP[i] ~ dt(zbeta0+sum(zbeta[1:Nx]*zxP[i,1:Nx]),1/zsigma^2, nu)  
43         yP[i] <- zyP[i] * ysd + ym  
44       }  
45   }  

Answer 1

似乎您还没有在第8行中定义xName是什么，这是Nx所依赖的。在dataList中定义它，您应该做得很好(关于修复这个单一错误)。