# NOT RUN {
library(embryogrowth)
eo <- subset(Database_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)
# }
Run the code above in your browser using DataLab