Learn R Programming

drc (version 1.2-0)

CRS.4a: The Cedergreen-Ritz-Streibig model

Description

'CRS.4a', 'CRS.4b' and 'CRS.4c' provide the Cedergreen-Ritz-Streibig modified log-logistic model for describing hormesis with the lower limit equal to 0. 'UCRS.4a', 'UCRS.4b' and 'UCRS.4c' provide the Cedergreen-Ritz-Streibig modified log-logistic model for describing u-shaped hormesis with the lower limit equal to 0.

Usage

CRS.4a(names = c("b", "d", "e", "f"))

  UCRS.4a(names = c("b", "d", "e", "f"))

Arguments

names
a vector of character strings giving the names of the parameters. The default is reasonable (see above).

Value

Details

The model is given by the expression $$f(x) = 0 + \frac{d-0+f \exp(-1/x)}{1+\exp(b(\log(x)-\log(e)))}$$ which is a five-parameter model. It is a modification of the four-parameter logistic curve to take hormesis into account. The u-shaped model is given by the expression $$f(x) = 0 + d - \frac{d-0+f \exp(-1/x^{\alpha})}{1+\exp(b(\log(x)-\log(e)))}$$ The a,b,c models are obtained by setting alpha equal to 1, 0.5 and 0.25, respectively.

References

See the reference under cedergreen.

See Also

Similar functions are CRS.5a and UCRS.5a, but with an extra parameter for the lower limit.

Examples

Run this code
## Modified logistic model
model1 <- multdrc(hormesis[,c(2,1)], fct=CRS.4a())
summary(model1)
ED(model1, c(50))

model2 <- multdrc(hormesis[,c(2,1)], fct=CRS.4b())
summary(model2)
ED(model2, c(50))

model3 <- multdrc(hormesis[,c(2,1)], fct=CRS.4c())
summary(model3)
ED(model3, c(50))

rm(model1, model2, model3)

Run the code above in your browser using DataLab