## Not run:
# library(embryogrowth)
# data(nest)
# formated <- FormatNests(nest)
# # The initial parameters value can be:
# # "T12H", "DHA", "DHH", "Rho25"
# # Or
# # "T12L", "T12H", "DHA", "DHH", "DHL", "Rho25"
# x <- structure(c(118.768297442004, 475.750095909406, 306.243694918151,
# 116.055824800264), .Names = c("DHA", "DHH", "T12H", "Rho25"))
# # pfixed <- c(K=82.33) or rK=82.33/39.33
# pfixed <- c(rK=2.093313)
# resultNest_4p <- searchR(parameters=x, fixed.parameters=pfixed,
# temperatures=formated, derivate=dydt.Gompertz, M0=1.7,
# test=c(Mean=39.33, SD=1.92))
# data(resultNest_4p)
# pMCMC <- TRN_MHmcmc_p(resultNest_4p, accept=TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab