## 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