lazer <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TAB15-1.DAT"),
# "C:\\E-drive\\Books\\ANREG2\\newdata\\tab5-6.dat",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\TAB15-1.dat",
		  sep="",col.names=c("O","Rel","Occ"))
attach(lazer)
lazer
laz <- xtabs(O~Rel+Occ)
laz

burn_in=1000
iterates=10000

or=O[Rel==1]
op=O[Rel==2]
oj=O[Rel==3]
Nr=sum(or)
Np=sum(op)
Nj=sum(oj)
ar=c(0.5,0.5,1.5,2.5)
ap=c(1,1,1.5,1.5)
aj=c(0.64,0.66,0.5,0.2)

occ.dat.jags<-list("or","op","oj","ar","ap","aj","Nr","Np","Nj")
occ.params<-c("pr","pp","pj")
inits <- function() {
  list( pr = ar/sum(ar), pp = ap/sum(ap), pj = aj/sum(aj))
}

occ.fit <- jags(data = occ.dat.jags, inits = inits,
  parameters.to.save = occ.params, n.chains = 1, n.thin=1,
n.iter = iterates+burn_in, n.burnin = burn_in,
model.file = "occupation-m.txt")



 x[] ~ dmulti(p[], N)
  p[] ~ ddirich(alpha[])
  ddirch
for (i in 1:I) {
x[i,1:K] ~ dmnorm(mu[], tau[,])
}

yr[] ~ dmulti(pr[],



model{
 for(k in 1:38){
   kills[k] ~ dpois(lambda[k]*trials[k])
   log(lambda[k]) <- beta[1]+beta[2]*(age[k]-mean(age[]))+
   beta[3]*(age[k]-mean(age[]))**2+eta[k] }
 for(i in 1:38){ eta[i] ~ dnorm(0,tau) }
 for(i in 1:11){
    log(ExKill[i]) <-  beta[1] + beta[2]*(pages[i]-mean(age[])) +
    beta[3]*(pages[i]-mean(age[]))**2 + 1/(2*tau)
   d[i] ~ dnorm(0,tau)
   log(r[i]) <- beta[1] + beta[2]*(pages[i]-mean(age[])) +
   beta[3]*(pages[i]-mean(age[]))**2  + d[i]
 }
 beta[1] ~ dnorm(0,0.0001)
 beta[2] ~ dnorm(0,0.0001)
 beta[3] ~ dnorm(0,0.0001)
 tau ~ dgamma(0.0001,0.0001)
}