model{ ## binomial likelihood for(i in 1:N){ count[i] ~ dbinom(p[i],n[i]) logit(p[i]) <- betaL[line[i]] + betaR[line[i]] * revert[i] + betaH[line[i]] * host[i] + betaI[line[i]] * revert[i] * host[i] } ## priors on main parameters for(j in 1:3){ ## line parameters betaL[j] ~ dnorm(gmn[1],tau[1]) betaR[j] ~ dnorm(gmn[2],tau[2]) betaH[j] ~ dnorm(gmn[3],tau[3]) betaI[j] ~ dnorm(gmn[4],tau[4]) } ## hyperpriors for(k in 1:4){ gmn[k] ~ dnorm(0,0.0001) tau[k] ~ dgamma(1,1) } } ## plot, prob. of acceptance pdf("binProbs.pdf",width=6,height=12) par(mar=c(3,5,3,1)) par(mfrow=c(2,1)) vioplot(c(outbin[[1]][,17],outbin[[2]][,17],outbin[[3]][,17]),c(outbin[[1]][,18],outbin[[2]][,18],outbin[[3]][,18]),c(outbin[[1]][,19],outbin[[2]][,19],outbin[[3]][,19]),c(outbin[[1]][,20],outbin[[2]][,20],outbin[[3]][,20]),c(outbin[[1]][,21],outbin[[2]][,21],outbin[[3]][,21]),c(outbin[[1]][,22],outbin[[2]][,22],outbin[[3]][,22]),col="gray",names=c("L1","L1R","L2","L2R","L3","L3R")) mtext("prob. of acceptance",2,cex=1.7,line=3) mtext("(a) lentil",3,cex=1.7,line=1) vioplot(c(outbin[[1]][,23],outbin[[2]][,23],outbin[[3]][,23]),c(outbin[[1]][,24],outbin[[2]][,24],outbin[[3]][,24]),c(outbin[[1]][,25],outbin[[2]][,25],outbin[[3]][,25]),c(outbin[[1]][,26],outbin[[2]][,26],outbin[[3]][,26]),c(outbin[[1]][,27],outbin[[2]][,27],outbin[[3]][,27]),c(outbin[[1]][,28],outbin[[2]][,28],outbin[[3]][,28]),col="gray",names=c("L1","L1R","L2","L2R","L3","L3R")) mtext("prob. of acceptance",2,cex=1.7,line=3) mtext("(b) mung bean",3,cex=1.7,line=1) dev.off()