Displaying 1 result from an estimated 1 matches for "veg_cover".
2012 Aug 30
0
storage mode error question (R2winBUGS)
...alpha[j]~dnorm(mu.alpha,tau.alpha)
}
mu.alpha~dnorm(0,0.01)
tau.alpha<-1/(sd.alpha*sd.alpha)
sd.alpha~dunif(0,5)
#Likelihood
for (i in 1:n){
for(j in 1:nsite)
C[i,j]~dpois(lambda[i,j])
log(lambda[i,j])<-mu+beta1*HABITAT[i,j]+alpha[j]
} #J
}#i
}
",fill=TRUE)
sink()
win.data<-list(C=Veg_cover$COLONY.SIZE,n=length(Veg_cover$COLONY.SIZE),
nsite=Veg_cover$SITE,HABITAT=Veg_cover$HABITAT)
inits<-function()list(mu=runif(1,0,4),beta1=runif(1,-1,1),alpha=runif(5,-2,2))
params<-c("mu","beta1","alpha")
ni<-6000
nt<-5
nb<-1000
nc<-3...