library(bbmle)

load("tadpole_sims.RData")
load("../R/fit_fits.RData")
parnames <- c("c","d","g")
simdat <- tadpole_sims

truevals <- melt(coef(fits[[1]]))
truevals <- data.frame(value=truevals,variable=rownames(truevals))


source("../../TOOLS/sim_plotfuns.R")
S <- simplots(simdat,parnames,truevals)

lapply(S,print)