library(rjags) # read data islands=read.csv("island-data.csv") #make a sequence to get predicted values x=seq(1,60) #create a list of data data=list(PA=islands[,1], y=islands[,2], n=nrow(islands), x=x, n.pred=length(x)) # specify initial values inits=list(list(a=0,b=0)) # call jags jm=jags.model("Island_JAGS.R",data=data,inits,n.chains=length(inits), n.adapt = 10000) zm=coda.samples(jm,variable.names=c("a","b","p.hat"), n.iter=50000, n.thin=1) zj=jags.samples(jm,variable.names=c("a","b","p.hat", "p.hat[10]"), n.iter=25000, n.thin=1) #get the quantiles for each island 1-60 s<-summary(zj$p.hat,quantile,c(.025,.5,.975))$stat #plot median and credible intervals plot(x,s[2,], xlab="Perimeter to Area Ratio", ylab= "Probability Occupied", typ='l') lines(x,s[1,], lty="dashed") lines(x,s[3,], lty="dashed")