library(rjags) ## mcmc settings burn<-2000 nc<-3 stps<-40000 thn<-20 dat<-read.csv("SINGLESEEDASSAY.csv") y<-which(dat[,4]=='Y') count<-as.numeric(dat[y,5]) n<-count+as.numeric(dat[-y,5]) line<-as.numeric(dat[y,1]) revert<-as.numeric(dat[y,2]=='R') host<-as.numeric(dat[y,3]=='LENTIL') mod<-jags.model("binmod",data=list(count=count,line=line,revert=revert,host=host,n=n,N=length(count)),n.chains=nc) update(mod,n.iter=burn) outbin<-coda.samples(model=mod,variable.names=c("p","betaL","betaR","betaH","betaI","gmn"),n.iter=stps,thin=thn)