HelpersMG (version 3.2)

plot.LD50: Plot results of LD50() that best describe LD50

Description

Plot the estimates that best describe lethality of exposures.

Usage

# S3 method for LD50
plot(x, ..., las.x = 1, las.y = 1, lab.PT = "LD50",
  lab.TRD = paste0("Transitional range of doses l=", l * 100, "%"),
  col.TRD = "gray", col.TRD.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)

Arguments

x

A result file generated by IC50()

...

Parameters for plot()

las.x

las parameter for x axis

las.y

las parameter for y axis

lab.PT

Label to describe pivotal temperature

lab.TRD

Label to describe transitional range of temperature

col.TRD

The color of TRD

col.TRD.CI

The color of CI of TRD 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

Value

Nothing

Details

plot.LD50 plot result of IC50() that best describe IC50

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 LD50 functions: LD50, logLik.LD50, predict.LD50

Examples

# NOT RUN {
#' data <- data.frame(Doses=c(80, 120, 150, 150, 180, 200),
Alive=c(10, 12, 8, 6, 2, 1),
Dead=c(0, 1, 5, 6, 9, 15))
LD50_logistic <- LD50(data, equation="logistic")
predict(LD50_logistic, doses=c(140, 170))
plot(LD50_logistic)
# }