1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
<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 <- c(51.2,51.4,51.1,51.8)
estimate <- calpredict(y,conc,acenaphthene,xunit="mg/L")
</pre>
<hr><div align="center"><a href="00Index.html">[Package Contents]</a></div>
</body></html>
|
