Skip to contents

This function plots the probit-transformed mortalities (probit.trans() function) as a function of the log10 of the dose, the regressions predicted by the resist.ratio() function, with or without confidence levels, if the dose-mortality responses are linear (option).

Usage

mort.plot(
  data,
  strains = NULL,
  plot.conf = TRUE,
  conf.level = 0.95,
  LD.value = c(25, 50, 95),
  test.validity = TRUE,
  legend.par = c("bottomleft"),
  ...
)

Arguments

data

a data frame of probit transformed mortality data using the function probit.trans()

strains

character. list of test strains to be plotted. If not provided, the function will plot all the strains in the data set.

plot.conf

logical. Whether to plot the confidence intervals for each strain, default TRUE

conf.level

numerical. The confidence interval to be plotted

LD.value

numerical. Level of lethal dose to be tested. default=c(25,50,95)

test.validity

logical. When TRUE (default), if a strain mortality-dose response fails the chi-square test for linearity in the resist.ratio() function, no regression will be plotted, only the observed data.

legend.par

multi-type. Arguments to be passed to the legend as in legend. default position bottomleft. If no legend desired use FALSE. Note: if pch, lty, and col are passed to the plot, they don't need to be passed to legend()

...

parameters to be passed on to graphics for the plot (e.g. col, pch)

Value

A plot of dose-mortality responses for bioassays

Author

Piyal Karunarathne, Pascal Milesi, Pierrick Labbé

Examples

data(bioassay)
transd<-probit.trans(bioassay$assay2)
data<-transd$tr.data
strains<-levels(data$strain)
mort.plot(data,strains)
#> Warning: longer object length is not a multiple of shorter object length
#> Warning: longer object length is not a multiple of shorter object length