aboutsummaryrefslogblamecommitdiff
path: root/chm/calpredict.html
blob: dff8ec412f7e66fc9f25d6c6b751e02a3ea6fb00 (plain) (tree)









































































































                                                                                                                                                                                                    
<html><head><title>Estimate measurement results including confidence intervals</title>
<link rel="stylesheet" type="text/css" href="Rchm.css">
</head>
<body>

<table width="100%"><tr><td>calpredict(chemCal)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
<param name="keyword" value="R:   calpredict">
<param name="keyword" value=" Estimate measurement results including confidence intervals">
</object>


<h2>Estimate measurement results including confidence intervals</h2>


<h3>Description</h3>

<p>
This function generates estimates for x values from y values, including
a confidence interval for the x values.  The formulas in this function used
for prediction of concentrations from (replicate) measurements are taken from 
the "Handbook of Chemometrics and Qualimetrics Part A" by D. L. Massart, 
Vandeginste, B. G. M., Buydens, L. M. C., De Jong, S., Lewi, P. J. and
Smeyers-Verbeke, J, Elsevier, Amsterdam, 1997 and from the EURACHEM/CITAC
report on "Quantifying uncertainty in analytical measurement", 2000, 
pp. 111f.
</p>


<h3>Usage</h3>

<pre>
  calpredict(yobs,xi,yi,xunit="",level=0.95,intercept=FALSE,syobs=FALSE)
</pre>


<h3>Arguments</h3>

<table summary="R argblock">
<tr valign="top"><td><code>yobs</code></td>
<td>
A numeric vector containing the observed data.
</td></tr>
<tr valign="top"><td><code>xi</code></td>
<td>
A vector of x values of the calibration. 
</td></tr>
<tr valign="top"><td><code>yi</code></td>
<td>
A vector of y values of the calibration. 
</td></tr>
<tr valign="top"><td><code>xunit</code></td>
<td>
The unit of the given values on the x axis as a character string.
</td></tr>
<tr valign="top"><td><code>level</code></td>
<td>
The desired confidence level for the confidence interval of the
estimates. Defaults to 0.95.
</td></tr>
<tr valign="top"><td><code>intercept</code></td>
<td>
Logical value determining if an intercept is to be fitted or not. 
Default is FALSE.
</td></tr>
<tr valign="top"><td><code>syobs</code></td>
<td>
If TRUE, a standard deviation for the given y values is 
calculated, and the resulting confidence interval will
include this variability (not validated yet). If FALSE (default), this
standard deviation is not included in the 
confidence interval. If a numeric value is given, 
it is used for the standard deviation of "real samples",
in addition to the standard deviation of the y values
in the calibration (also not validated yet).
</td></tr>
</table>

<h3>Value</h3>

<p>
A list containing the estimate, its standard deviation and its 
confidence interval.</p>

<h3>Author(s)</h3>

<p>
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>
</p>


<h3>Examples</h3>

<pre>
data(pahCalibration)
attach(pahCalibration)
y &lt;- c(51.2,51.4,51.1,51.8)
estimate &lt;- calpredict(y,conc,acenaphthene,xunit="mg/L")
</pre>



<hr><div align="center"><a href="00Index.html">[Package Contents]</a></div>

</body></html>

Contact - Imprint