tempConst | R Documentation |
Timeseries of temperatures for nests
tempConst
A dataframe with raw data.
Timeseries of constant temperatures for nests
Marc Girondot marc.girondot@universite-paris-saclay.fr
## Not run:
library(embryogrowth)
# Same as:
# GenerateConstInc(durations = rep(104*60*24, 11),
# temperatures = 25:35,
# names = paste0("T",25:35))
data(tempConst)
tempConst_f <- FormatNests(tempConst)
data(nest)
formated <- FormatNests(nest)
x <- structure(c(109.683413821537, 614.969219372661, 306.386903812694,
229.003478775323), .Names = c("DHA", "DHH", "T12H", "Rho25"))
# See the stages dataset examples for justification of M0 and rK
pfixed <- c(rK=1.208968)
resultNest_4p_SSM <- searchR(parameters=x, fixed.parameters=pfixed,
temperatures=formated, integral=integral.Gompertz, M0=0.3470893,
hatchling.metric=c(Mean=39.33, SD=1.92))
plotR(result=resultNest_4p_SSM, show.hist = TRUE,
ylim=c(0, 8), curve="ML quantiles")
# Now use the fited parameters from resultNest_4p_SSM with
# the constant incubation temperatures:
plot(resultNest_4p_SSM, temperatures=tempConst_f,
stop.at.hatchling.metric=TRUE, series="T30", xlim=c(0,50),
ylimT=c(22, 32), hatchling.metric=c(Mean=39.33, SD=1.92),
embryo.stages="Caretta caretta.SCL")
plot(resultNest_4p_SSM, temperatures=tempConst_f,
stop.at.hatchling.metric=TRUE, series="T25", xlim=c(0,120),
ylimT=c(22, 32), hatchling.metric=c(Mean=39.33, SD=1.92),
embryo.stages="Caretta caretta.SCL")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.