lapply(c("rjags","R2jags","arm","coda","superdiag","R2WinBUGS"),library,
    character.only=TRUE)
source("http://jeffgill.org/files/jeffgill/files/abortion.panel-data.list_.txt")
names(attitudes.list) <- c("N","P","r","n","x1")

attitudes.mod <- function() {
    for (i in 1:N)  {
        for (j in 1:P)  {
            logit(p[i,j]) <-  b0[j] + b1[i]*x1[i]
            r[i,j] ~ dbin(p[i,j], n[i])
        }
        b1[i] ~ dnorm(mu1, nu1)
    }
    for (j in 1:P)  {
        b0[j] ~ dnorm(mu0, nu0)
    }
    mu0  ~ dnorm(0.0,0.1)
    mu1  ~ dnorm(0.0,0.1)
    nu0 ~ dgamma(1,0.1)
    nu1 ~ dgamma(1,0.1)
}

attitudes.params <- c("b0","b1")
attitudes.out <- jags(data=attitudes.list, parameters.to.save=attitudes.params,
    n.iter=10000, model=attitudes.mod,n.burnin=2500, n.thin=1,n.chains=3)
update(attitudes.out, n.iter=50000, n.burnin=0, n.thin=1,n.chains=3)

attitudes.mcmc <- as.mcmc(attitudes.out)
attitudes.mat <- attitudes.mcmc[[1]]
means.vec <- apply(attitudes.mat,2,mean)
sd.vec <- apply(attitudes.mat,2,sd)