require(R2jags)

tadpole_BUGS_fit <- function(data) {
tfit_jags <- jags(model="tadpole_bugs.txt",
                  data=c(N=nrow(data),as.list(data)),
                  parameters.to.save=c("c","d","g"),
                  n.chains=1,
                  inits=list(list(c=0.45,d=13,g=1)))
}

## dropdev <- function(x) {
##   if (length(dim(x))==2) {
##     if ("deviance" %in% rownames(x)) {
##       x[rownames(x)!="deviance",]
##     } else if ("deviance" %in% colnames(x)) {
##       x[,colnames(x)!="deviance"]
##     } else x
##   } else {
##     if (is.null(names(x))) x else x[names(x)!="deviance"]
##   }
## }