## Not run:
# library(embryogrowth)
# eo <- subset(STSRE_TSD, Species=="Emys orbicularis", c("Males", "Females",
# "Incubation.temperature"))
# eo_logistic <- with(eo, tsd(males=Males, females=Females,
# temperatures=Incubation.temperature))
# pMCMC <- tsd_MHmcmc_p(eo_logistic, accept=TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab