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
|
<html lang="en">
<head>
<title>polyconf - optim_doc</title>
<meta http-equiv="Content-Type" content="text/html">
<meta name="description" content="optim_doc">
<meta name="generator" content="makeinfo 4.13">
<link title="Top" rel="start" href="index.html#Top">
<link rel="up" href="Residual-optimization.html#Residual-optimization" title="Residual optimization">
<link rel="prev" href="wpolyfit.html#wpolyfit" title="wpolyfit">
<link rel="next" href="LinearRegression.html#LinearRegression" title="LinearRegression">
<link href="http://www.gnu.org/software/texinfo/" rel="generator-home" title="Texinfo Homepage">
<!--
Additional documentation for the optim package for Octave.
Copyright (C) <Olaf Till <i7tiol@t-online.de>>
You can redistribute this documentation and/or modify it under the terms
of the GNU General Public License as published by the Free Software
Foundation; either version 3 of the License, or (at your option) any
later version.
This documentation is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
Public License for more details.
You should have received a copy of the GNU General Public License along
with this documentation; if not, see <http://www.gnu.org/licenses/>.-->
<meta http-equiv="Content-Style-Type" content="text/css">
<style type="text/css"><!--
pre.display { font-family:inherit }
pre.format { font-family:inherit }
pre.smalldisplay { font-family:inherit; font-size:smaller }
pre.smallformat { font-family:inherit; font-size:smaller }
pre.smallexample { font-size:smaller }
pre.smalllisp { font-size:smaller }
span.sc { font-variant:small-caps }
span.roman { font-family:serif; font-weight:normal; }
span.sansserif { font-family:sans-serif; font-weight:normal; }
--></style>
</head>
<body>
<div class="node">
<a name="polyconf"></a>
<p>
Next: <a rel="next" accesskey="n" href="LinearRegression.html#LinearRegression">LinearRegression</a>,
Previous: <a rel="previous" accesskey="p" href="wpolyfit.html#wpolyfit">wpolyfit</a>,
Up: <a rel="up" accesskey="u" href="Residual-optimization.html#Residual-optimization">Residual optimization</a>
<hr>
</div>
<h3 class="section">2.11 Confidence and prediction intervals for polynomial fitting</h3>
<p><a name="index-polyconf-59"></a>
<h4 class="subheading">Helptext:</h4>
<p><a name="XREFpolyconf"></a>
<pre class="verbatim">[y,dy] = polyconf(p,x,s)
Produce prediction intervals for the fitted y. The vector p
and structure s are returned from polyfit or wpolyfit. The
x values are where you want to compute the prediction interval.
polyconf(...,['ci'|'pi'])
Produce a confidence interval (range of likely values for the
mean at x) or a prediction interval (range of likely values
seen when measuring at x). The prediction interval tells
you the width of the distribution at x. This should be the same
regardless of the number of measurements you have for the value
at x. The confidence interval tells you how well you know the
mean at x. It should get smaller as you increase the number of
measurements. Error bars in the physical sciences usually show
a 1-alpha confidence value of erfc(1/sqrt(2)), representing
one standandard deviation of uncertainty in the mean.
polyconf(...,1-alpha)
Control the width of the interval. If asking for the prediction
interval 'pi', the default is .05 for the 95% prediction interval.
If asking for the confidence interval 'ci', the default is
erfc(1/sqrt(2)) for a one standard deviation confidence interval.
Example:
[p,s] = polyfit(x,y,1);
xf = linspace(x(1),x(end),150);
[yf,dyf] = polyconf(p,xf,s,'ci');
plot(xf,yf,'g-;fit;',xf,yf+dyf,'g.;;',xf,yf-dyf,'g.;;',x,y,'xr;data;');
plot(x,y-polyval(p,x),';residuals;',xf,dyf,'g-;;',xf,-dyf,'g-;;');
</pre>
<!-- -->
</body></html>
|
