@@ -22,26 +22,25 @@ in regular intervals to check for updates from those sources.
Please cite this package if you use it in your analysis. See `citation("stressaddition")` for details.
## Example
In the paper we use the model in the context of survival experiments. However, it can also be applicable in modeling other concentration or dose dependent responses. For this reason the more general term "effect" instead of "survival" is used throughout the package.
```r
library(stressaddition)
model<-ecxsys(
concentration=c(0,0.05,0.5,5,30),
hormesis_concentration=0.5,
effect_tox_observed=c(90,81,92,28,0),
effect_tox_env_observed=c(29,27,33,5,0)
survival_tox_observed=c(90,81,92,28,0),
survival_tox_env_observed=c(29,27,33,5,0)
)
# Plot the effect and the system stress:
par(mfrow=c(2,1))
plot_effect(model)
plot_survival(model)
plot_stress(model)
# The LC50 of the effect under the influence of toxicant and system tress:
ec(model,"effect_tox_sys",50)
# The LC50 under the influence of toxicant and system tress:
lc(model,"survival_tox_sys",50)
# The LC10 of the effect under the influence of toxicant, system and environmental tress:
ec(model,"effect_tox_env_sys",10)
# The LC10 under the influence of toxicant, system and environmental tress: