aboutsummaryrefslogblamecommitdiff
path: root/docs/reference/drcfit.html
blob: 75a168ac6d22ee07c02a8192a6ef663875a62c29 (plain) (tree)


















































































































































































































                                                                                                                                                                                                                                                                                                                                                                                                                                                     









                                                                                                                                                                                                                                                                                                                                                      






































                                                                                                     
<!-- Generated by pkgdown: do not edit by hand -->
<!DOCTYPE html>
<html>
  <head>
  <meta charset="utf-8">
<meta http-equiv="X-UA-Compatible" content="IE=edge">
<meta name="viewport" content="width=device-width, initial-scale=1.0">

<title>Fit dose-response models using the drc package — drcfit • drfit</title>

<!-- jquery -->
<script src="https://code.jquery.com/jquery-3.1.0.min.js" integrity="sha384-nrOSfDHtoPMzJHjVTdCopGqIqeYETSXhZDFyniQ8ZHcVy08QesyHcnOUpMpqnmWq" crossorigin="anonymous"></script>
<!-- Bootstrap -->

<link href="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.7/css/bootstrap.min.css" rel="stylesheet" integrity="sha384-BVYiiSIFeK1dGmJRAkycuHAHRg32OmUcww7on3RYdg4Va+PmSTsz/K68vbdEjh4u" crossorigin="anonymous">
<script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.7/js/bootstrap.min.js" integrity="sha384-Tc5IQib027qvyjSMfHjOMaLkfuWVxZxUPnCJA7l2mCWNIpG9mGCD8wGNIcPD7Txa" crossorigin="anonymous"></script>

<!-- Font Awesome icons -->
<link href="https://maxcdn.bootstrapcdn.com/font-awesome/4.6.3/css/font-awesome.min.css" rel="stylesheet" integrity="sha384-T8Gy5hrqNKT+hzMclPo118YTQO6cYprQmhrYwIiQ/3axmI1hQomh7Ud2hPOy8SP1" crossorigin="anonymous">


<!-- pkgdown -->
<link href="../pkgdown.css" rel="stylesheet">
<script src="../jquery.sticky-kit.min.js"></script>
<script src="../pkgdown.js"></script>

<!-- mathjax -->
<script src='https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML'></script>

<!--[if lt IE 9]>
<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
<![endif]-->


  </head>

  <body>
    <div class="container template-reference-topic">
      <header>
      <div class="navbar navbar-default navbar-fixed-top" role="navigation">
  <div class="container">
    <div class="navbar-header">
      <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar">
        <span class="icon-bar"></span>
        <span class="icon-bar"></span>
        <span class="icon-bar"></span>
      </button>
      <a class="navbar-brand" href="../index.html">drfit</a>
    </div>
    <div id="navbar" class="navbar-collapse collapse">
      <ul class="nav navbar-nav">
        <li>
  <a href="../reference/index.html">Reference</a>
</li>
      </ul>
      
      <ul class="nav navbar-nav navbar-right">
        
      </ul>
    </div><!--/.nav-collapse -->
  </div><!--/.container -->
</div><!--/.navbar -->

      
      </header>

      <div class="row">
  <div class="col-md-9 contents">
    <div class="page-header">
    <h1>Fit dose-response models using the drc package</h1>
    </div>

    
    <p>Fit dose-response relationships to dose-response data and calculate
  biometric results for (eco)toxicity evaluation using the drc package</p>
    

    <pre><span class='fu'>drcfit</span>(<span class='no'>data</span>, <span class='kw'>chooseone</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='kw'>probit</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='kw'>logit</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>,
  <span class='kw'>weibull</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='kw'>linlogit</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='kw'>level</span> <span class='kw'>=</span> <span class='fl'>0.95</span>,
  <span class='kw'>showED50</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='kw'>EDx</span> <span class='kw'>=</span> <span class='kw'>NULL</span>)</pre>
    
    <h2 class="hasAnchor" id="arguments"><a class="anchor" href="#arguments"></a> Arguments</h2>
    <table class="ref-arguments">
    <colgroup><col class="name" /><col class="desc" /></colgroup>
    <tr>
      <th>data</th>
      <td><p>A data frame containing dose-response data. The data frame has to contain
    at least a factor called &#8220;substance&#8221;, a numeric vector &#8220;dose&#8221;
    with the dose values, a vector called &#8220;unit&#8221; containing the unit
    used for the dose and a numeric vector &#8220;response&#8221; with the response
    values of the test system normalized between 0 and 1. Such a data frame can
    be easily obtained if a compliant RODBC data source is available for use in
    conjunction with the function <code><a href='drdata.html'>drdata</a></code>.</p>
<p>If there is a column called &#8220;ok&#8221; and it is set to &#8220;no fit&#8221; in
    a specific line, then the corresponding data point will be excluded from
    the fitting procedure, although it will be plotted.</p></td>
    </tr>
    <tr>
      <th>probit</th>
      <td><p>A boolean defining if cumulative density curves of normal distributions
    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><a href='drfit.html'>drfit</a></code>. Parameter e from <code>LN.2</code> is listed
    as a here, and parameter b from there is listed as b.</p></td>
    </tr>
    <tr>
      <th>logit</th>
      <td><p>A boolean defining if cumulative density curves of logistic distributions
    <code>plogis</code> are fitted to the decadic logarithm of the dose.
    Default is FALSE.
    Again the parameter definitions used in the model are different to the
    ones used in <code><a href='drfit.html'>drfit</a></code>. Parameter e from <code>LL.2</code> is listed
    as a here, and parameter b from LL.2 is listed as b.</p></td>
    </tr>
    <tr>
      <th>weibull</th>
      <td><p>A boolean defining if Weibull dose-response models
    (<code>W1.2</code> are fitted to the untransformed dose. Default is FALSE.
    Note that the results differ from the ones obtained with
    <code><a href='drfit.html'>drfit</a></code>, due to a different model specification.</p></td>
    </tr>
    <tr>
      <th>linlogit</th>
      <td><p>A boolean defining if the linear-logistic function
    <code><a href='linlogitf.html'>linlogitf</a></code> as defined by van Ewijk and Hoekstra 1993 is
    fitted to the data. Default is FALSE. Obtaining the ED50 (and EDx values
    in general) uses <code>ED</code> internally and does not always give a
    result.</p></td>
    </tr>
    <tr>
      <th>level</th>
      <td><p>The level for the confidence interval listed for the log ED50.</p></td>
    </tr>
    <tr>
      <th>chooseone</th>
      <td><p>If TRUE (default), the models are tried in the order linlogit, probit,
    logit, weibull, and the first model that produces a valid fit is used.
    If FALSE, all models that are set to TRUE and that can be fitted will be
    reported.</p></td>
    </tr>
    <tr>
      <th>EDx</th>
      <td><p>A vector of inhibition values x in percent for which the corresponding doses
    EDx should be reported.</p></td>
    </tr>
    <tr>
      <th>showED50</th>
      <td><p>If set to TRUE, the ED50 and its confidence interval on the original dose
    scale (not log scale) is included in the output.</p></td>
    </tr>
    </table>
    
    <h2 class="hasAnchor" id="value"><a class="anchor" href="#value"></a>Value</h2>

    <p>A dataframe with the attribute <code>models</code> holding a list of the fitted
  dose-response models of class <code>nls</code>. The dataframe has at least
  one line for each substance.</p>
<p>The following variables are in the dataframe:</p>
<dt>Substance</dt><dd><p>The name of the substance</p></dd>
  <dt>ndl</dt><dd><p>The number of dose levels in the raw data</p></dd>
  <dt>n</dt><dd><p>The total number of data points in the raw data used for the fit</p></dd>
  <dt>lld</dt><dd><p>The decadic logarithm of the lowest dose</p></dd>
  <dt>lhd</dt><dd><p>The total number of data points in the raw data used for the fit</p></dd>
  <dt>mtype</dt><dd><p>If the data did not show a mean response &lt; 0.5 at the highest dose level,
    the modeltype is set to &#8220;inactive&#8221;. If the mean response at the
    lowest dose is smaller than 0.5, the modeltype is set to &#8220;active&#8221;.
    In both cases, no fitting procedure is carried out. If the fitted ED50 
    is higher than the highest dose, &#8220;no fit&#8221; is given here.</p></dd>
  <dt>logED50</dt><dd><p>The decadic logarithm of the ED50</p></dd>
  <dt>low %</dt><dd><p>The lower bound of the confidence interval of log ED50. 
    The name of the column depends on the requested confidence <code>level</code>.</p></dd>
  <dt>high %</dt><dd><p>The higher bound of the confidence interval of log ED50.
    The name of the column depends on the requested confidence <code>level</code>.</p></dd>
  <dt>unit</dt><dd><p>The unit used for the dose levels in the dose-response data</p></dd>
  <dt>sigma</dt><dd><p>The square root of the estimated variance of the random error as returned
    by <code>summary.drc</code>.</p></dd>
  <dt>a</dt><dd><p>For the linlogit model, this is the parameter e from <code>BC.4</code>.
    For the probit and the logit model, this is the ED50. For the weibull 
    model, this is parameter e from <code>W1.2</code>. Note that the Weibull 
    model is fitted to the untransformed data.</p></dd>
  <dt>b</dt><dd><p>For the linlogit, probit, logit and weibull models, these are the
    parameters b from <code>BC.4</code>, <code>LN.2</code>, 
    <code>LL.2</code> and <code>W1.2</code>, respectively.
    Note that the parameter definitions (and in the case of Weibull, the model
    used) are different to the ones used in <code>drfit</code>.</p></dd>
  <dt>c</dt><dd><p>Only the &#8220;linlogit&#8221; fit produces a third parameter <code>c</code>, which is
    the parameter f from the <code>BC.4</code> function.</p></dd>
  If the parameter 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 EDx,
  the corresponding dose levels are given in addition, together with their
  confidence intervals as calculated by ED from the drc package.

    
    <h2 class="hasAnchor" id="note"><a class="anchor" href="#note"></a>Note</h2>

    <p>There is a demo for each dataset that can be accessed by
  <code>demo(dataset)</code></p>
    
    <h2 class="hasAnchor" id="see-also"><a class="anchor" href="#see-also"></a>See also</h2>

    <p>Further examples are given in help pages to the datasets
  <code><a href='antifoul.html'>antifoul</a></code>, <code><a href='IM1xIPC81.html'>IM1xIPC81</a></code> and
  <code><a href='IM1xVibrio.html'>IM1xVibrio</a></code>.</p>
    

    <h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
    <pre class="examples"><div class='input'><span class='fu'>data</span>(<span class='no'>antifoul</span>)
<span class='no'>r</span> <span class='kw'>&lt;-</span> <span class='fu'>drcfit</span>(<span class='no'>antifoul</span>, <span class='kw'>showED50</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='kw'>EDx</span> <span class='kw'>=</span> <span class='fu'>c</span>(<span class='fl'>5</span>, <span class='fl'>10</span>, <span class='fl'>20</span>))</div><div class='output co'>#&gt; <span class='message'></span>
#&gt; <span class='message'>TBT: Fitting data...</span></div><div class='output co'>#&gt; <span class='message'></span>
#&gt; <span class='message'>Zn Pyrithion: Fitting data...</span></div><div class='input'><span class='fu'>format</span>(<span class='no'>r</span>, <span class='kw'>digits</span> <span class='kw'>=</span> <span class='fl'>2</span>)</div><div class='output co'>#&gt;      Substance ndl   n  lld lhd  mtype logED50  2.5%  97.5%   unit sigma    a
#&gt; 1          TBT  38 135 -2.7 2.4 probit   -0.16 -0.28 -0.072 microM  0.19 0.68
#&gt; 2 Zn Pyrithion  27  81 -2.1 2.0 probit   -0.40 -0.52 -0.303 microM  0.23 0.40
#&gt;       b ED50 ED50 2.5% ED50 97.5%  EDx5 EDx5 2.5% EDx5 97.5% EDx10 EDx10 2.5%
#&gt; 1 -0.64 0.68      0.52       0.85 0.053     0.015      0.091 0.093      0.040
#&gt; 2 -1.04 0.40      0.30       0.50 0.082     0.023      0.142 0.117      0.048
#&gt;   EDx10 97.5% EDx20 EDx20 2.5% EDx20 97.5%
#&gt; 1        0.15  0.18       0.11        0.26
#&gt; 2        0.19  0.18       0.10        0.26</div></pre>
  </div>
  <div class="col-md-3 hidden-xs hidden-sm" id="sidebar">
    <h2>Contents</h2>
    <ul class="nav nav-pills nav-stacked">
      <li><a href="#arguments">Arguments</a></li>
      
      <li><a href="#value">Value</a></li>

      <li><a href="#note">Note</a></li>

      <li><a href="#see-also">See also</a></li>
      
      <li><a href="#examples">Examples</a></li>
    </ul>

    <h2>Author</h2>
    
  Johannes Ranke <a href='mailto:jranke@uni-bremen.de'>jranke@uni-bremen.de</a> 
  <a href = 'http://www.uft.uni-bremen.de/chemie/ranke'>http://www.uft.uni-bremen.de/chemie/ranke</a>
  The functionality of the drc package used under the hood in this function
  was written by Christian Ritz.

  </div>
</div>

      <footer>
      <div class="copyright">
  <p>Developed by Johannes Ranke.</p>
</div>

<div class="pkgdown">
  <p>Site built with <a href="http://hadley.github.io/pkgdown/">pkgdown</a>.</p>
</div>

      </footer>
   </div>

  </body>
</html>

Contact - Imprint