aboutsummaryrefslogtreecommitdiff
path: root/man/drcfit.Rd
diff options
context:
space:
mode:
authorranke <ranke@d1b72e20-2ee0-0310-a1c4-ad5adbbefcdc>2014-02-25 21:26:05 +0000
committerranke <ranke@d1b72e20-2ee0-0310-a1c4-ad5adbbefcdc>2014-02-25 21:26:05 +0000
commit189aad5e39df513075fc7c7a64e646b1a354bc4c (patch)
tree39d44b4495eab9893e4c90eaefb78d4ff96ea6b5 /man/drcfit.Rd
parent35362cfc7c02f12fd9689f68a772a804d895535a (diff)
More detailed documentation of the result dataframe returned
by drfit and drcfit. git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/drfit@99 d1b72e20-2ee0-0310-a1c4-ad5adbbefcdc
Diffstat (limited to 'man/drcfit.Rd')
-rw-r--r--man/drcfit.Rd95
1 files changed, 71 insertions, 24 deletions
diff --git a/man/drcfit.Rd b/man/drcfit.Rd
index 2da0418..b10418a 100644
--- a/man/drcfit.Rd
+++ b/man/drcfit.Rd
@@ -25,11 +25,10 @@
the fitting procedure, although it will be plotted.}
\item{probit}{
A boolean defining if cumulative density curves of normal distributions
- \code{\link{pnorm}} are fitted against the decadic logarithm of the dose.
- Default ist TRUE.
+ are fitted against the decadic logarithm of the dose. Default ist TRUE.
Note that the parameter definitions used in the model are different to the
ones used in \code{\link{drfit}}. Parameter e from \code{\link{LN.2}} is listed
- as a here, and parameter b from LN.2 is listed as b.}
+ as a here, and parameter b from there is listed as b.}
\item{logit}{
A boolean defining if cumulative density curves of logistic distributions
\code{\link{plogis}} are fitted to the decadic logarithm of the dose.
@@ -66,27 +65,74 @@
}
}
\value{
- \item{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 \dQuote{inactive}. If the mean response at the lowest
- dose is smaller than 0.5, the modeltype is set to \dQuote{active}. In
- both cases, no fitting procedure is carried out. Every successful fit is
- reported in one line. Parameters of the fitted curves are only reported
- if the fitted ED50 is not higher than the highest dose.
+ A dataframe with the attribute \code{models} holding a list of the fitted
+ dose-response models of class \code{\link{nls}}. The dataframe has at least
+ one line for each substance.
- \code{ndl} is the number of dose levels in the raw data, \code{n} is the
- total number of data points in the raw data used for the fit
- \code{lld} is the decadic logarithm of the lowest dose and
- \code{lhd} is the decadic logarithm of the highest dose.
-
- If the parameter \code{showED50} was set to TRUE, the ED50 values and their
- confidence intervals are also included on the original dose scale.
-
- If one or more response leves were specified in the argument \code{EDx},
- the corresponding dose levels are given, together with their confidence
- intervals.
+ The following variables are in the dataframe:
+ \item{Substance}{
+ The name of the substance
+ }
+ \item{ndl}{
+ The number of dose levels in the raw data
+ }
+ \item{n}{
+ The total number of data points in the raw data used for the fit
+ }
+ \item{lld}{
+ The decadic logarithm of the lowest dose
+ }
+ \item{lhd}{
+ The total number of data points in the raw data used for the fit
+ }
+ \item{mtype}{
+ If the data did not show a mean response < 0.5 at the highest dose level,
+ the modeltype is set to \dQuote{inactive}. If the mean response at the
+ lowest dose is smaller than 0.5, the modeltype is set to \dQuote{active}.
+ In both cases, no fitting procedure is carried out. If the fitted ED50
+ is higher than the highest dose, \dQuote{no fit} is given here.
+ }
+ \item{logED50}{
+ The decadic logarithm of the ED50
+ }
+ \item{low \%}{
+ The lower bound of the confidence interval of log ED50.
+ The name of the column depends on the requested confidence \code{level}.
}
+ \item{high \%}{
+ The higher bound of the confidence interval of log ED50.
+ The name of the column depends on the requested confidence \code{level}.
+ }
+ \item{unit}{
+ The unit used for the dose levels in the dose-response data
+ }
+ \item{sigma}{
+ The square root of the estimated variance of the random error as returned
+ by \code{\link{summary.drc}}.
+ }
+ \item{a}{
+ For the linlogit model, this is the parameter e from \code{\link{BC.4}}.
+ For the probit and the logit model, this is the ED50. For the weibull
+ model, this is parameter e from \code{\link{W1.2}}. Note that the Weibull
+ model is fitted to the untransformed data.
+ }
+ \item{b}{
+ For the linlogit, probit, logit and weibull models, these are the
+ parameters b from \code{\link{BC.4}}, \code{\link{LN.2}},
+ \code{\link{LL.2}} and \code{\link{W1.2}}, respectively.
+ Note that the parameter definitions (and in the case of Weibull, the model
+ used) are different to the ones used in \code{\link{drfit}}.
+ }
+ \item{c}{
+ Only the \dQuote{linlogit} fit produces a third parameter \code{c}, which is
+ the parameter f from the \code{\link{BC.4}} function.
+ }
+ If the parameter \code{showED50} was set to TRUE, the ED50 values and their
+ confidence intervals are also included on the original dose scale.
+
+ If one or more response leves were specified in the argument \code{EDx},
+ the corresponding dose levels are given in addition, together with their
+ confidence intervals as calculated by \code{\link{ED}} from the drc package.
}
\examples{
data(antifoul)
@@ -101,9 +147,10 @@ format(r, digits = 2)
\code{\link{IM1xVibrio}}.
}
\author{
- Johannes Ranke
- \email{jranke@uni-bremen.de}
+ Johannes Ranke \email{jranke@uni-bremen.de}
\url{http://www.uft.uni-bremen.de/chemie/ranke}
+ The functionality of the drc package used under the hood in this function
+ was written by Christian Ritz.
}
\keyword{models}
\keyword{regression}

Contact - Imprint