embryogrowth (version 6.4)

predict.tsd: Estimate sex ratio according to constant incubation temperature

Description

Estimate sex ratio according to constant incubation temperature The data.frame has the following components: temperatures, SE, sexratio, CI.minus.sexratio, CI.plus.sexratio, range.CI Note that incubation duration is a very bad proxy for sex ratio. See Georges, A., Limpus, C. J. & Stoutjesdijk, R. 1994. Hatchling sex in the marine turtle Caretta caretta is determined by proportion of development at a temperature, not daily duration of exposure. J. Exp. Zool., 270, 432-444.

Usage

"predict"(object, temperatures = NULL, durations = NULL, SE = NULL, quantiles = c(0.025, 0.5, 0.975), replicates = 1000, progressbar = FALSE, ...)

Arguments

object

A result file generated by tsd

temperatures

A vector of temperatures

durations

A vector of durations

SE

The standard error for temperatures or durations

quantiles

The quantiles to be returned, default=c(0.025, 0.5, 0.975)

replicates

Number of replicates to estimate CI

progressbar

Logical. Does a progression bar must be shown?

...

Not used

Value

A data.frame with informations about sex-ratio

Details

predict.tsd Estimate sex ratio according to constant incubation temperature

See Also

Other Functions.for.temperature.dependent.sex.determination: DatabaseTSD, TSP.list, stages, tsd_MHmcmc_p, tsd_MHmcmc, tsd

Examples

## Not run: 
# library(embryogrowth)
# m <- c(10, 14, 7, 4, 3, 0, 0)
# f <- c(0, 1, 2, 4, 15, 10, 13)
# t <- c(25, 26, 27, 28, 29, 30, 31)
# result <- tsd(males=m, females=f, temperatures=t)
# plot(result)
# predict(result, temperatures=c(25, 31))
# predict(result, temperatures=c(25, 31), SE=c(1, 2))
# d <- c(72, 70, 65, 63, 62, 60, 59)
# result <- tsd(males=m, females=f, durations=d)
# predict(result, durations=c(67, 68))
# ## End(Not run)