From 13d2fa72cb1b31cc8d47d57e476f65c56942fe8f Mon Sep 17 00:00:00 2001 From: Sebastian Henz <sebastian.henz@ufz.de> Date: Mon, 18 May 2020 18:41:17 +0200 Subject: [PATCH] Last changes before manuscript submission. See NEWS for version 3.0.0 --- .Rbuildignore | 1 + .../merge_into_master.md | 2 +- DESCRIPTION | 22 +++--- NEWS.md | 7 +- R/data.R | 46 ++++++++++++ R/multi_tox.R | 68 +++++++++++------- R/plot_ecxsys.R | 2 +- R/stressaddition-package.R | 18 ++--- README.md | 2 +- data-raw/multiple_stress.R | 23 ++++++ data-raw/multiple_stress.csv | 59 +++++++++++++++ data/multiple_stress.rda | Bin 0 -> 524 bytes man/multi_tox.Rd | 68 +++++++++++------- man/multiple_stress.Rd | 34 +++++++++ man/plot_ecxsys.Rd | 2 +- man/stressaddition-package.Rd | 33 ++++----- 16 files changed, 297 insertions(+), 90 deletions(-) create mode 100644 R/data.R create mode 100644 data-raw/multiple_stress.R create mode 100644 data-raw/multiple_stress.csv create mode 100644 data/multiple_stress.rda create mode 100644 man/multiple_stress.Rd diff --git a/.Rbuildignore b/.Rbuildignore index 40f7f3d..99b26c0 100644 --- a/.Rbuildignore +++ b/.Rbuildignore @@ -7,3 +7,4 @@ ^LICENSE\.txt$ ^\.gitlab$ ^\.gitlab-ci\.yml$ +^data-raw$ diff --git a/.gitlab/merge_request_templates/merge_into_master.md b/.gitlab/merge_request_templates/merge_into_master.md index 7790e4c..b4b04c5 100644 --- a/.gitlab/merge_request_templates/merge_into_master.md +++ b/.gitlab/merge_request_templates/merge_into_master.md @@ -12,8 +12,8 @@ Comment on any notes that come up during devtools::check(). --> ### system details -* RStudio Version: * R Version: +* RStudio Version: * OS Version: ### checklist diff --git a/DESCRIPTION b/DESCRIPTION index a607fb5..0afe62b 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,8 +1,8 @@ Package: stressaddition Type: Package -Title: Modeling Tri-Phasic Concentration-Response Relationships +Title: Modelling Tri-Phasic Concentration-Response Relationships Version: 3.0.0 -Date: 2020-05-07 +Date: 2020-05-18 Authors@R: c(person("Sebastian", "Henz", role = c("aut", "cre"), @@ -12,20 +12,24 @@ Authors@R: c(person("Sebastian", "Liess", role = "aut", comment = c(ORCID = "0000-0002-3321-8909")), + person("Naeem", + "Shahid", + role = "ctb", + comment = c(ORCID = "0000-0001-6581-1654")), person("Helmholtz-Zentrum fuer Umweltforschung GmbH - UFZ", - role = c("cph", "fnd"))) -Description: The implementation of the tri-phasic concentration-response model - introduced in "Predicting low-concentration effects of pesticides" by - Liess, Henz and Knillmann (2019) <doi: 10.1038/s41598-019-51645-4>. It - allows modeling of ecotoxicological experiments where the response shows - signs of hormesis. + role = c("cph"))) +Description: The stress addition approach is an alternative to the traditional + concentration addition or effect addition models. It allows the modelling + of tri-phasic concentration-response relationships either as single toxicant + experiments, in combination with an environmental stressor or as mixtures of + two toxicants. License: GPL-3 Copyright: file inst/COPYRIGHTS URL: https://git.ufz.de/oekotox/stressaddition Encoding: UTF-8 LazyData: true Depends: - R (> 3.5) + R (>= 3.5) Imports: drc (>= 3.0), plotrix diff --git a/NEWS.md b/NEWS.md index 413ddf9..1422585 100644 --- a/NEWS.md +++ b/NEWS.md @@ -4,14 +4,17 @@ * Renamed all instances of "effect" to "survival". * Renamed all instances of "ec" to "lc". * Renamed `predict_mixture()`, which was a temporary development name, to `multi_tox()`. -* The argument `proportion_ca` in the mixture model `multi_tox()` was renamed and its value reversed. It is now called `sa_contribution` and specifies the proportion of stress addition in the calculation of toxicant stress. To convert your code from the old version use this equation: sa_contribution = 1 - proportion_ca. +* The argument `proportion_ca` in the mixture model `multi_tox()` was renamed and its value reversed. It is now called `sa_contribution` and specifies the proportion of stress addition in the calculation of toxicant stress. To convert your code from the old version use this equation: `sa_contribution = 1 - proportion_ca`. * Renamed `stress_tox_sam` to `stress_tox_sa` in the output of `multi_tox()`. ## Bugfixes * Fixed a bug where `plot_stress()` with argument `which = NULL` would result in an error. Now it correctly draws the axes without data. +* Fixed some errors in the documentation and improved the examples. ## New -* Exporte function `log10_ticks()` for calculating tick mark labels and positions on a base 10 logarithmic axis. +* Export function `log10_ticks()` for calculating tick mark labels and positions on a base 10 logarithmic axis. +* Add example data set `multiple_stress` for use with `multi_tox()`. +* Various minor changes to prepare for CRAN submission. # stressaddition 2.7.0 diff --git a/R/data.R b/R/data.R new file mode 100644 index 0000000..eb8c1d8 --- /dev/null +++ b/R/data.R @@ -0,0 +1,46 @@ +# Copyright (C) 2020 Helmholtz-Zentrum fuer Umweltforschung GmbH - UFZ +# See file inst/COPYRIGHTS for details. +# +# This file is part of the R package stressaddition. +# +# stressaddition is free software: you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation, either version 3 of the License, or +# (at your option) any later version. +# +# This program is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program. If not, see <https://www.gnu.org/licenses/>. + + +# Documentation for all included data sets. + + +#' Survival of Daphnia magna exposed to multiple stressors +#' +#' A data set of mixture toxicity experiments. Individuals of Daphnia magna were +#' exposed to combinations of food stress and multiple concentrations of the +#' pesticides esfenvalerate and prochloraz. The survival was recorded at 21 days +#' after contamination. +#' +#' @format A data frame with 58 rows and 4 variables: +#' \describe{ +#' \item{food}{The amount of food in percent.} +#' \item{esfenvalerate}{The concentration of esfenvalerate in µg/L.} +#' \item{prochloraz}{The concentration of prochloraz in µg/L.} +#' \item{survival}{The mean survival in percent.} +#' } +#' +#' @source +#' \href{https://doi.org/10.1021/acs.est.9b04293}{Shahid, N., Liess, M., +#' Knillmann, S., 2019. Environmental Stress Increases Synergistic Effects +#' of Pesticide Mixtures on Daphnia magna. Environ. Sci. Technol. 53, +#' 12586–12593.} +#' +#' Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects +#' of toxicant mixtures. Manuscript submitted for publication. +"multiple_stress" diff --git a/R/multi_tox.R b/R/multi_tox.R index 882c069..97ae971 100644 --- a/R/multi_tox.R +++ b/R/multi_tox.R @@ -24,11 +24,12 @@ #' details. #' #' The predictions are symmetric, i.e. it does not matter which of the toxicant -#' models is 1 or 2 as long as the concentration arguments are supplied in the -#' right order. See the example below. +#' models is a or b as long as the concentration arguments are supplied in the +#' same order. #' -#' This method is only suitable for experiments without environmental stress. -#' Any environmental stress in \code{model_a} or \code{model_b} is ignored. +#' This method is only suitable for experiments without or with low +#' environmental stress. Any environmental stress supplied as arguments to +#' \code{\link{ecxsys}} in \code{model_a} or \code{model_b} is ignored. #' #' @param model_a,model_b The ecxsys models of the toxicants. #' @param concentration_a,concentration_b The concentrations of the toxicants in @@ -45,31 +46,48 @@ #' @return A data frame with columns of the supplied concentrations and the #' corresponding mixture survival and stresses. #' -#' @examples toxicant_a <- ecxsys( -#' concentration = c(0, 0.05, 0.5, 5, 30), -#' hormesis_concentration = 0.5, -#' survival_tox_observed = c(90, 81, 92, 28, 0), +#' @examples +#' # Using a data set which is included in this package. See ?multiple_stress +#' ms <- multiple_stress +#' esfen <- ms[ms$food == 1 & ms$prochloraz == 0, ] +#' proch <- ms[ms$food == 1 & ms$esfenvalerate == 0, ] +#' +#' model_esfen <- ecxsys( +#' concentration = esfen$esfenvalerate, +#' survival_tox_observed = esfen$survival, +#' hormesis_concentration = 0.1 #' ) -#' toxicant_b <- ecxsys( -#' concentration = c(0, 0.1, 1, 10, 100, 1000), -#' hormesis_concentration = 10, -#' survival_tox_observed = c(26, 25, 24, 27, 5, 0), -#' survival_max = 30 +#' model_proch <- ecxsys( +#' concentration = proch$prochloraz, +#' survival_tox_observed = proch$survival, +#' hormesis_concentration = 100 #' ) -#' multi_tox( -#' toxicant_a , -#' toxicant_b , -#' c(0, 0.02, 0.2, 2, 20), -#' c(0, 0.03, 0.4, 5, 10) +#' +#' # Predict the survival at 8 different esfenvalerate concentrations +#' # but keep the prochloraz concentration constant at 32: +#' mt <- multi_tox( +#' model_esfen, +#' model_proch, +#' c(0, 0.0001, 0.001, 0.01, 0.1, 0.316, 1, 3.16), +#' 32, +#' sa_contribution = 0.8 #' ) +#' mt[1:3] # The concentrations and survival of the 8 mixtures. +#' +#' # Predict the survival at 4 different combinations +#' # of esfenvalerate and prochloraz: +#' mt <- multi_tox( +#' model_esfen, +#' model_proch, +#' c(0.1, 0.2, 0.3, 0.4), +#' c(0, 1, 32, 100), +#' sa_contribution = 0.8 +#' ) +#' mt[1:3] # The concentrations and survival of the 4 mixtures. +#' #' -#' # Example of symmetric prediction: -#' conc_a <- c(0, 0.03, 0.3, 3) -#' conc_b <- 5.5 -#' sa_contrib <- 0.75 -#' mix_a <- multi_tox(toxicant_a , toxicant_b , conc_a, conc_b, sa_contrib) -#' mix_b <- multi_tox(toxicant_b , toxicant_a , conc_b, conc_a, sa_contrib) -#' identical(mix_a$survival, mix_b$survival) +#' @references Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic +#' effects of toxicant mixtures. Manuscript submitted for publication. #' #' @export multi_tox <- function(model_a, diff --git a/R/plot_ecxsys.R b/R/plot_ecxsys.R index f0d9c0f..6d1ec9e 100644 --- a/R/plot_ecxsys.R +++ b/R/plot_ecxsys.R @@ -26,7 +26,7 @@ #' #' @name plot_ecxsys #' -#' @param model The list returned from \code{\link{ecxsys}}. +#' @param model The object returned from \code{\link{ecxsys}}. #' @param which A vector of names to plot. Allowed are the column names of the #' \code{model$curves} data frame. There is also #' \code{"survival_tox_observed"} and \code{"survival_tox_env_observed"} for diff --git a/R/stressaddition-package.R b/R/stressaddition-package.R index e36e1f8..c761fcb 100644 --- a/R/stressaddition-package.R +++ b/R/stressaddition-package.R @@ -27,21 +27,21 @@ # description. -#' @details See the publication linked below for more information including -#' equations. +#' @details This R package contains the definitions of the +#' \link[=ecxsys]{ECx-SyS} model and the \link[=multi_tox]{Multi-TOX} model. +#' See the publications linked below for more information including equations. #' -#' In the paper the model is introduced in the context of survival -#' experiments. However, we conjecture that it can also be used to model -#' other concentration dependent responses. -#' -#' This package is not on CRAN. It is hosted on the UFZ GitLab server. Visit -#' the repository linked below for the newest version. See the readme file -#' there for instructions on how to install and update the package. +#' Author contributions: M. Liess conceived the ECx-SyS and Multi-TOX models. +#' S. Henz developed this R package. N. Shahid contributed to the optimization +#' of the Multi-TOX model and provided the \link{multiple_stress} data set. #' #' @references \href{https://doi.org/10.1038/s41598-019-51645-4}{Liess, M., #' Henz, S. & Knillmann, S. Predicting low-concentration effects of #' pesticides. Sci Rep 9, 15248 (2019).} #' +#' Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects of +#' toxicant mixtures. Manuscript submitted for publication. +#' #' @import stats #' @import graphics #' @import grDevices diff --git a/README.md b/README.md index 0a1c02a..b92fa6f 100644 --- a/README.md +++ b/README.md @@ -5,7 +5,7 @@ This is the R implementation of the tri-phasic concentration-response model intr The EC<sub>x-SyS</sub> and Multi-TOX models from this package are also available as part of the [Indicate app](http://www.systemecology.eu/indicate) which offers a graphical user interface. ## Installation -Stressaddition is not yet on CRAN. You can install the most recent stable version from GitLab using the remotes package: +Stressaddition is not yet on CRAN. You can install the most recent development version from GitLab using the remotes package: ``` r install.packages("remotes") remotes::install_gitlab("oekotox/stressaddition", host = "git.ufz.de") diff --git a/data-raw/multiple_stress.R b/data-raw/multiple_stress.R new file mode 100644 index 0000000..856f0e8 --- /dev/null +++ b/data-raw/multiple_stress.R @@ -0,0 +1,23 @@ +# Copyright (C) 2020 Helmholtz-Zentrum fuer Umweltforschung GmbH - UFZ +# See file inst/COPYRIGHTS for details. +# +# This file is part of the R package stressaddition. +# +# stressaddition is free software: you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation, either version 3 of the License, or +# (at your option) any later version. +# +# This program is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program. If not, see <https://www.gnu.org/licenses/>. + + +# This script prepares the esfenvalerate-prochloraz-food data for packaging. + +multiple_stress <- read.csv("data-raw/multiple_stress.csv") +usethis::use_data(multiple_stress) diff --git a/data-raw/multiple_stress.csv b/data-raw/multiple_stress.csv new file mode 100644 index 0000000..1221185 --- /dev/null +++ b/data-raw/multiple_stress.csv @@ -0,0 +1,59 @@ +"food","esfenvalerate","prochloraz","survival" +1,0,0,84.44 +1,0,1,55.56 +1,0,32,53.33 +1,0,100,48.89 +1,0,1000,0 +1,1e-04,0,73.33 +1,1e-04,1,60 +1,1e-04,32,40 +1,1e-04,100,20 +1,0.001,0,66.67 +1,0.001,1,50 +1,0.001,32,26.67 +1,0.001,100,16.67 +1,0.01,0,53.33 +1,0.01,1,33.33 +1,0.01,32,23.33 +1,0.01,100,0 +1,0.1,0,46.67 +1,0.1,1,31.11 +1,0.1,32,16.67 +1,0.1,100,6.67 +1,0.316,0,20 +1,0.316,1,8.89 +1,0.316,32,0 +1,0.316,100,0 +1,1,0,0 +1,1,1,0 +1,1,32,0 +1,1,100,0 +1,3.16,0,0 +1,3.16,1,0 +1,3.16,32,0 +1,3.16,100,0 +100,0,0,93.33 +100,0,1,93.33 +100,0,32,88.89 +100,0,100,88.89 +100,0,3200,0 +100,0.001,0,93.33 +100,0.001,1,93.33 +100,0.001,32,84.44 +100,0.001,100,80 +100,0.1,0,80 +100,0.1,1,80 +100,0.1,32,66.67 +100,0.1,100,62.22 +100,0.316,0,55.56 +100,0.316,1,64.44 +100,0.316,32,28.89 +100,0.316,100,35.56 +100,1,0,31.11 +100,1,1,31.11 +100,1,32,0 +100,1,100,6.67 +100,3.16,0,0 +100,3.16,1,0 +100,3.16,32,0 +100,3.16,100,0 diff --git a/data/multiple_stress.rda b/data/multiple_stress.rda new file mode 100644 index 0000000000000000000000000000000000000000..5c2867452cbaabcf96d3fbdf0dd23a9d98f36087 GIT binary patch literal 524 zcmV+n0`vVsT4*^jL0KkKSy?>c_5cC#|NsAY%oSu~Rrc)v_r$+v|CmMqfG`X>0ALY- z9Dq1XV+Ftgy`f4D00E!@pbY>38ZiI^KmY&$13(6dnl&aP)cq-@nFpz$02%{9qaXkP z$N)6R0Knh?8UPvq&;S9W5CAj)0000q04ho<W~NOV0BB?YGyrL!05kv%8UPIq7$Y?_ z1|}~Xmu@`H?7dVmYpu8GZ@!nf%GU#5F)MAc4aFfosHp<ji;fI*up2>8ff0)5%+Yv6 z95W&yqzV-xwj>wVEO@*2650tcYDqQ}VBj*Ta^;tjFCkGPTuDjr7M4hAch*FhwF}xu zKv#G8b8RJ}^A01VTvNAM9m>!ztdEIeOO)ox!kv@F@4aBt9Ip#BUzW1x6Q|XhiB@UJ z^DFE#fcPwMMkIeCA}!Ml+7%pXl9rfI8(5CgTn-pkV|k`Z<#=J0&#e|XsS6Ve*cMIr zrfe~+JaoSH88)WV&N4F*kNEvRKX+57YZ!LMUz2MjEWyjbzx!3L^c{H<yx=oIKB%Ed z1=6>zA-u*WjhBcS4Xtc<1=U0|I}@YUkm|6jTZymMT*>w37D#+K+;cW~6FnY*GsZ~5 zG>0QG#bxztnDMJdq^VM6LtHfSXkGl~=QG0KyZ)b4%2Z7AN9lCZ(ih_B($iFzMnht+ O_`8xR!i0j!=MS()!{J{5 literal 0 HcmV?d00001 diff --git a/man/multi_tox.Rd b/man/multi_tox.Rd index ebd17c4..7abbb9f 100644 --- a/man/multi_tox.Rd +++ b/man/multi_tox.Rd @@ -40,37 +40,55 @@ details. } \details{ The predictions are symmetric, i.e. it does not matter which of the toxicant -models is 1 or 2 as long as the concentration arguments are supplied in the -right order. See the example below. +models is a or b as long as the concentration arguments are supplied in the +same order. -This method is only suitable for experiments without environmental stress. -Any environmental stress in \code{model_a} or \code{model_b} is ignored. +This method is only suitable for experiments without or with low +environmental stress. Any environmental stress supplied as arguments to +\code{\link{ecxsys}} in \code{model_a} or \code{model_b} is ignored. } \examples{ -toxicant_a <- ecxsys( - concentration = c(0, 0.05, 0.5, 5, 30), - hormesis_concentration = 0.5, - survival_tox_observed = c(90, 81, 92, 28, 0), +# Using a data set which is included in this package. See ?multiple_stress +ms <- multiple_stress +esfen <- ms[ms$food == 1 & ms$prochloraz == 0, ] +proch <- ms[ms$food == 1 & ms$esfenvalerate == 0, ] + +model_esfen <- ecxsys( + concentration = esfen$esfenvalerate, + survival_tox_observed = esfen$survival, + hormesis_concentration = 0.1 ) -toxicant_b <- ecxsys( - concentration = c(0, 0.1, 1, 10, 100, 1000), - hormesis_concentration = 10, - survival_tox_observed = c(26, 25, 24, 27, 5, 0), - survival_max = 30 +model_proch <- ecxsys( + concentration = proch$prochloraz, + survival_tox_observed = proch$survival, + hormesis_concentration = 100 ) -multi_tox( - toxicant_a , - toxicant_b , - c(0, 0.02, 0.2, 2, 20), - c(0, 0.03, 0.4, 5, 10) + +# Predict the survival at 8 different esfenvalerate concentrations +# but keep the prochloraz concentration constant at 32: +mt <- multi_tox( + model_esfen, + model_proch, + c(0, 0.0001, 0.001, 0.01, 0.1, 0.316, 1, 3.16), + 32, + sa_contribution = 0.8 ) +mt[1:3] # The concentrations and survival of the 8 mixtures. + +# Predict the survival at 4 different combinations +# of esfenvalerate and prochloraz: +mt <- multi_tox( + model_esfen, + model_proch, + c(0.1, 0.2, 0.3, 0.4), + c(0, 1, 32, 100), + sa_contribution = 0.8 +) +mt[1:3] # The concentrations and survival of the 4 mixtures. -# Example of symmetric prediction: -conc_a <- c(0, 0.03, 0.3, 3) -conc_b <- 5.5 -sa_contrib <- 0.75 -mix_a <- multi_tox(toxicant_a , toxicant_b , conc_a, conc_b, sa_contrib) -mix_b <- multi_tox(toxicant_b , toxicant_a , conc_b, conc_a, sa_contrib) -identical(mix_a$survival, mix_b$survival) } +\references{ +Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic + effects of toxicant mixtures. Manuscript submitted for publication. +} diff --git a/man/multiple_stress.Rd b/man/multiple_stress.Rd new file mode 100644 index 0000000..8667622 --- /dev/null +++ b/man/multiple_stress.Rd @@ -0,0 +1,34 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/data.R +\docType{data} +\name{multiple_stress} +\alias{multiple_stress} +\title{Survival of Daphnia magna exposed to multiple stressors} +\format{ +A data frame with 58 rows and 4 variables: +\describe{ + \item{food}{The amount of food in percent.} + \item{esfenvalerate}{The concentration of esfenvalerate in µg/L.} + \item{prochloraz}{The concentration of prochloraz in µg/L.} + \item{survival}{The mean survival in percent.} +} +} +\source{ +\href{https://doi.org/10.1021/acs.est.9b04293}{Shahid, N., Liess, M., + Knillmann, S., 2019. Environmental Stress Increases Synergistic Effects + of Pesticide Mixtures on Daphnia magna. Environ. Sci. Technol. 53, + 12586–12593.} + + Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects + of toxicant mixtures. Manuscript submitted for publication. +} +\usage{ +multiple_stress +} +\description{ +A data set of mixture toxicity experiments. Individuals of Daphnia magna were +exposed to combinations of food stress and multiple concentrations of the +pesticides esfenvalerate and prochloraz. The survival was recorded at 21 days +after contamination. +} +\keyword{datasets} diff --git a/man/plot_ecxsys.Rd b/man/plot_ecxsys.Rd index 0cb3187..d5c09b0 100644 --- a/man/plot_ecxsys.Rd +++ b/man/plot_ecxsys.Rd @@ -26,7 +26,7 @@ plot_survival( ) } \arguments{ -\item{model}{The list returned from \code{\link{ecxsys}}.} +\item{model}{The object returned from \code{\link{ecxsys}}.} \item{which}{A vector of names to plot. Allowed are the column names of the \code{model$curves} data frame. There is also diff --git a/man/stressaddition-package.Rd b/man/stressaddition-package.Rd index 9908cea..795b1ee 100644 --- a/man/stressaddition-package.Rd +++ b/man/stressaddition-package.Rd @@ -4,30 +4,30 @@ \name{stressaddition-package} \alias{stressaddition} \alias{stressaddition-package} -\title{stressaddition: Modeling Tri-Phasic Concentration-Response Relationships} +\title{stressaddition: Modelling Tri-Phasic Concentration-Response Relationships} \description{ -The implementation of the tri-phasic concentration-response model - introduced in "Predicting low-concentration effects of pesticides" by - Liess, Henz and Knillmann (2019) <doi: 10.1038/s41598-019-51645-4>. It - allows modeling of ecotoxicological experiments where the response shows - signs of hormesis. +The stress addition approach is an alternative to the traditional + concentration addition or effect addition models. It allows the modelling + of tri-phasic concentration-response relationships either as single toxicant + experiments, in combination with an environmental stressor or as mixtures of + two toxicants. } \details{ -See the publication linked below for more information including - equations. +This R package contains the definitions of the + \link[=ecxsys]{ECx-SyS} model and the \link[=multi_tox]{Multi-TOX} model. + See the publications linked below for more information including equations. - In the paper the model is introduced in the context of survival - experiments. However, we conjecture that it can also be used to model - other concentration dependent responses. - - This package is not on CRAN. It is hosted on the UFZ GitLab server. Visit - the repository linked below for the newest version. See the readme file - there for instructions on how to install and update the package. + Author contributions: M. Liess conceived the ECx-SyS and Multi-TOX models. + S. Henz developed this R package. N. Shahid contributed to the optimization + of the Multi-TOX model and provided the \link{multiple_stress} data set. } \references{ \href{https://doi.org/10.1038/s41598-019-51645-4}{Liess, M., Henz, S. & Knillmann, S. Predicting low-concentration effects of pesticides. Sci Rep 9, 15248 (2019).} + + Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects of + toxicant mixtures. Manuscript submitted for publication. } \seealso{ Useful links: @@ -46,7 +46,8 @@ Authors: Other contributors: \itemize{ - \item Helmholtz-Zentrum fuer Umweltforschung GmbH - UFZ [copyright holder, funder] + \item Naeem Shahid (\href{https://orcid.org/0000-0001-6581-1654}{ORCID}) [contributor] + \item Helmholtz-Zentrum fuer Umweltforschung GmbH - UFZ [copyright holder] } } -- GitLab