Learn R Programming

embryogrowth (version 6.1.1)

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

## S3 method for class 'tsd':
predict(object, temperatures = NULL, durations = NULL,
  SE = NULL, range.CI = 0.95, 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
range.CI
The range of confidence interval for estimation, default=0.95
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: STSRE_TSD; TSP.list; tsd_MHmcmc_p; tsd_MHmcmc; tsd

Examples

Run this code
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))

Run the code above in your browser using DataLab