drfit (version 0.03-13)

drfit: Fit dose-response models

Description

Fit dose-response relationships to dose-response data and calculate biometric results for (eco)toxicity evaluation

Usage

drfit(data, startlogEC50 = NA, chooseone = TRUE, lognorm = TRUE, logis = FALSE, 
    linearlogis = FALSE, b0 = 2, f0 = 0)

Arguments

data

A data frame as returned from drdata. The data frame has to contain at least a factor called "substance", a vector called "unit" containing the unit used for the dose, a column "response" with the

startlogEC50

Especially for the linearlogis model, a suitable log10 of the EC50 has to be given by the user, since it is not correctly estimated for data showing hormesis with the default estimation method.

lognorm

A boolean defining if cumulative density curves of normal distributions are fitted to the data. Default ist TRUE.

logis

A boolean defining if cumulative densitiy curves of logistic distributions are fitted to the data. Default is FALSE.

linearlogis

A boolean defining if the linear-logistic function as defined by van Ewijk and Hoekstra 1993 is fitted to the data. Default is FALSE.

chooseone

If TRUE (default), the models are tried in the order linearlogis, logis and lognorm, and the first model that produces a valid fit is used. Usually this will be the one with the lowest residual standard deviation. If FALSE, all models that are se

b0,f0

If the linearlogistic model is fitted, b0 and f0 give the possibility to adapt the starting values for the parameters b and f.

Value

• results
  A data frame containing at least one line for each substance. If the data did not show a mean response < 0.5 at the highest dose level, the modeltype is set to "none". Every successful fit is reported in one line. Parameters of the fitted curves are only reported if the fitted EC50 is not higher than the highest dose.

Examples

data(antifoul)
r <- drfit(antifoul)
format(r,digits=2)