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>
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.