library(bbmle)

load("../../../R/fit_fits.RData")
Rfit <- fits[[1]]
d <- as.data.frame(Rfit@data)

simdata <- function(seed) {
  if (!missing(seed)) set.seed(seed)
  transform(d,Kill=simulate(Rfit))
}

dat<-simdata()

# Write data in the same filename and format as the real data 

write.table(dat,file="tadpole.dat",row.names=FALSE)