Description Usage Format Details Source References Examples

Data consist of 5 mixture, 6 dilutions, three replicates, and 12 common controls; in total 102 onservations.

1 |

A data frame with 102 observations on the following 3 variables.

`dose`

a numeric vector of dose values

`pct`

a numeric vector denoting the grouping according to the mixtures percentages

`rgr`

a numeric vector of response values (relative growth rates)

The dataset is analysed in Soerensen et al (2007). The asymmetric Voelund model is appropriate, whereas the symmetric Hewlett model is not.

The dataset is kindly provided by Nina Cedergreen, Department of Agricultural Sciences, Royal Veterinary and Agricultural University, Denmark.

Soerensen, H. and Cedergreen, N. and Skovgaard, I. M. and Streibig, J. C. (2007)
An isobole-based statistical model and test for synergism/antagonism in binary mixture toxicity experiments,
*Environmental and Ecological Statistics*, **14**, 383–397.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | ```
## Fitting the model with freely varying ED50 values
mecter.free <- drm(rgr ~ dose, pct, data = mecter,
fct = LL.4(), pmodels = list(~1, ~1, ~1, ~factor(pct) - 1))
## Lack-of-fit test
modelFit(mecter.free) # not really acceptable
summary(mecter.free)
## Plotting isobole structure
isobole(mecter.free, exchange = 0.02)
## Fitting the concentration addition model
mecter.ca <- mixture(mecter.free, model = "CA")
## Comparing to model with freely varying e parameter
anova(mecter.ca, mecter.free) # rejected
## Plotting isobole based on concentration addition
isobole(mecter.free, mecter.ca, exchange = 0.02) # poor fit
## Fitting the Hewlett model
mecter.hew <- mixture(mecter.free, model = "Hewlett")
## Comparing to model with freely varying e parameter
anova(mecter.hew, mecter.free) # rejected
## Plotting isobole based on the Hewlett model
isobole(mecter.free, mecter.hew, exchange = 0.02) # poor fit
## Fitting the Voelund model
mecter.voe<-mixture(mecter.free, model = "Voelund")
## Comparing to model with freely varying e parameter
anova(mecter.voe, mecter.free) # accepted
## Plotting isobole based on the Voelund model
isobole(mecter.free, mecter.voe, exchange = 0.02) # good fit
``` |

drc documentation built on May 29, 2017, 11:31 a.m.

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.