Description Usage Arguments Value References Examples
Calculates lethal time (LT) and its fiducial confidence limits (CL) using a probit analysis according to Finney 1971, Wheeler et al. 2006, and Robertson et al. 2007.
1 2 3 4 5 6 7 8 9 10 11 12 13 |
formula |
an object of class |
data |
an optional data frame, list or environment (or object coercible by as.data.frame to a data frame) containing the variables in the model. If not found in data, the variables are taken from environment(formula), typically the environment from which |
p |
Lethal time (LT) values for given p, example will return a LT50 value if p equals 50. If more than one LT value wanted specify by creating a vector. LT values can be calculated down to the 1e-16 of a percentage (e.g. LT99.99). However, the tibble produced can and will round to nearest whole number. |
weights |
vector of 'prior weights' to be used in the fitting process. Only needs to be supplied if you are taking the response / total for your response variable within the formula call of |
subset |
allows for the data to be subseted if desired. Default set to |
log_base |
default is |
log_x |
default is |
het_sig |
significance level from person's chi square goodness-of-fit test that is used to decide if a heterogeneity factor is used. |
conf_level |
Adjust confidence level as necessary or |
conf_type |
default is |
long_output |
default is |
Returns a tibble with predicted LT for given p level, lower CL (LCL), upper CL (UCL), LCL, Pearson's chi square goodness-of-fit test (pgof), slope, intercept, slope and intercept p values and standard error, and LT variance.
Finney, D.J., 1971. Probit Analysis, Cambridge University Press, Cambridge, England, ISBN: 052108041X
Wheeler, M.W., Park, R.M., and Bailey, A.J., 2006. Comparing median lethal concentration values using confidence interval overlap or ratio tests, Environ. Toxic. Chem. 25(5), 1441-1444.10.1897/05-320R.1
Robertson, J.L., Savin, N.E., Russell, R.M. and Preisler, H.K., 2007. Bioassays with arthropods. CRC press. ISBN: 9780849323317
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | head(lamprey_time)
results <- LT_probit((response / total) ~ log10(hour),
p = c(50, 99),
weights = total,
data = lamprey_time,
subset = c(month == "May"))
# view calculated LT50 and LT99 for seasonal
# toxicity of a piscicide, 3-trifluoromethyl-4-nitrophenol, to lamprey in 2011
results
# dose-response curve can be plotted using 'ggplot2'
|
# A tibble: 6 x 7
nominal.dose dose month hour response survive total
<dbl> <dbl> <chr> <int> <int> <int> <int>
1 1.1 1.25 May 1 0 58 58
2 1.1 1.25 May 2 0 58 58
3 1.1 1.25 May 3 0 58 58
4 1.1 1.25 May 4 1 57 58
5 1.1 1.25 May 5 5 53 58
6 1.1 1.25 May 6 14 44 58
# A tibble: 2 x 19
p n time LCL UCL LCL_dis UCL_dis chi_square df pgof_sig h
<dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int> <dbl> <dbl>
1 50 11 9.99 8.93 11.7 1.06 1.75 13.5 9 0.141 1.50
2 99 11 38.1 26.3 73.7 11.8 35.6 13.5 9 0.141 1.50
# … with 8 more variables: slope <dbl>, slope_se <dbl>, slope_sig <dbl>,
# intercept <dbl>, intercept_se <dbl>, intercept_sig <dbl>, z <dbl>,
# var_m <dbl>
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.