小弟初学R2WinBUGS
按照课本的指点，依葫芦画瓢做的一下程序

library(R2WinBUGS)
library(lattice)
library(coda)
sink("lme.model1.txt")
cat("
model {
# Priors
for (i in 1:ngroups){
alpha ~ dnorm(mu.int, tau.int) # Random intercepts
}
mu.int ~ dnorm(0, 0.001)# Mean hyperparameter for random intercepts
tau.int <- 1 / (sigma.int * sigma.int)
sigma.int ~ dunif(0, 100)# SD hyperparameter for random intercepts
beta ~ dnorm(0, 0.001) # Common slope
tau <- 1 / ( sigma * sigma) # Residual precision
sigma ~ dunif(0, 100) # Residual standard deviation
# Likelihood
for (i in 1:n) {
mass ~ dnorm(mu, tau)# The random variable
mu <- alpha[pop] + beta* length # Expectation
}
}
",fill=TRUE)
sink()

# Bundle data
win.data <- list(mass = as.numeric(mass), pop = as.numeric(pop),
length = length, ngroups = max(as.numeric(pop)), n = n)
到了这里时怎么把自己的数据读成win.data，请大侠指点
自己的数据data在附件中