embryogrowth (version 8.0)

plot.tsd: Plot results of tsd() that best describe temperature-dependent sex determination

Description

Plot the estimates that best describe temperature-dependent sex determination.

Usage

# S3 method for tsd
plot(
  x,
  ...,
  show.observations = TRUE,
  show.model = TRUE,
  males.freq = TRUE,
  show.PTRT = TRUE,
  las.x = 1,
  las.y = 1,
  lab.PT = paste("Pivotal ", x$type),
  resultmcmc = NULL,
  chain = 1,
  l = 0.05,
  replicate.CI = 10000,
  range.CI = 0.95,
  mar = c(4, 4, 4, 1) + 0.4,
  temperatures.plot = seq(from = 25, to = 35, by = 0.1),
  durations.plot = seq(from = 40, to = 70, by = 0.1),
  lab.TRT = paste0("Transitional range of ", x$type, "s l=", x$l * 100, "%"),
  col.TRT = "gray",
  col.TRT.CI = rgb(0.8, 0.8, 0.8, 0.5),
  col.PT.CI = rgb(0.8, 0.8, 0.8, 0.5),
  show.CI = TRUE,
  warn = TRUE
)

Arguments

x

A result file generated by tsd()

...

Parameters for plot()

show.observations

Should the observations be shown

show.model

Should the model be shown

males.freq

Should the graph uses males relative frequency [TRUE] or females [FALSE]

show.PTRT

Should the P and TRT information be shown

las.x

las parameter for x axis

las.y

las parameter for y axis

lab.PT

Label to describe pivotal temperature

resultmcmc

A result of tsd_MHmcmc()

chain

What chain to be used is resultmcmc is provided

l

Sex ratio limits to define TRT are l and 1-l

replicate.CI

replicate.CI replicates from the hessian matrix to estimate CI

range.CI

The range of confidence interval for estimation, default=0.95

mar

The par("mar") parameter

temperatures.plot

Temperatures used for showing curves of sex ratio

durations.plot

Durations used for showing curves of sex ratio

lab.TRT

Label to describe transitional range of temperature

col.TRT

The color of TRT

col.TRT.CI

The color of CI of TRT based on range.CI

col.PT.CI

The color of CI of PT based on range.CI

show.CI

Do the CI for the curve should be shown

warn

Do the warnings must be shown ? TRUE or FALSE

Value

Nothing

Details

plot.tsd plot result of tsd() that best describe temperature-dependent sex determination

References

Girondot, M. 1999. Statistical description of temperature-dependent sex determination using maximum likelihood. Evolutionary Ecology Research, 1, 479-486.

Godfrey, M.H., Delmas, V., Girondot, M., 2003. Assessment of patterns of temperature-dependent sex determination using maximum likelihood model selection. Ecoscience 10, 265-272.

Hulin, V., Delmas, V., Girondot, M., Godfrey, M.H., Guillon, J.-M., 2009. Temperature-dependent sex determination and global change: are some species at greater risk? Oecologia 160, 493-506.

Girondot M., Submited. On the concept of embryological thermosensitive period for sex determination in reptiles.

See Also

Other Functions for temperature-dependent sex determination: DatabaseTSD.version(), DatabaseTSD, P_TRT(), TSP.list, predict.tsd(), stages, tsd_MHmcmc_p(), tsd_MHmcmc(), tsd()

Examples

# NOT RUN {
CC_AtlanticSW <- subset(DatabaseTSD, RMU=="Atlantic, SW" & 
                          Species=="Caretta caretta" & (!is.na(Sexed) & Sexed!=0))
tsdL <- with (CC_AtlanticSW, tsd(males=Males, females=Females, 
                                 temperatures=Incubation.temperature-Correction.factor, 
                                 equation="logistic"))
plot(tsdL)
# }