HelpersMG (version 5.1)

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

Description

Plot the estimates that best describe lethality of doses.

Usage

# S3 method for LD50
plot(
  x,
  ...,
  las.x = 1,
  las.y = 1,
  lab.PT = "LD50",
  at = NULL,
  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 dose

at

Position of ticks in x-axis

lab.TRD

Label to describe transitional range of dose

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

See Also

Other LD50 functions: LD50_MHmcmc_p(), LD50_MHmcmc(), 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, xlim=c(0, 300))
# }