43 | | <div class="dynheader"> |
44 | | Collaboration diagram for bdm::ARX:</div> |
45 | | <div class="dynsection"> |
46 | | <p><center><img src="classbdm_1_1ARX__coll__graph.png" border="0" usemap="#bdm_1_1ARX__coll__map" alt="Collaboration graph"></center> |
47 | | <map name="bdm_1_1ARX__coll__map"> |
48 | | <area shape="rect" href="classbdm_1_1BMEF.html" title="Estimator for Exponential family." alt="" coords="15,332,108,359"><area shape="rect" href="classbdm_1_1BM.html" title="Bayesian Model of a system, i.e. all uncertainty is modeled by probabilities." alt="" coords="5,199,83,225"><area shape="rect" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects." alt="" coords="55,7,167,33"><area shape="rect" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables." alt="" coords="73,87,148,113"><area shape="rect" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density." alt="" coords="107,199,192,225"><area shape="rect" href="classldmat.html" title="Matrix stored in LD form, (commonly known as UD)." alt="" coords="213,279,272,305"><area shape="rect" href="classbdm_1_1egiw.html" title="Gauss-inverse-Wishart density stored in LD form." alt="" coords="108,385,193,412"><area shape="rect" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form..." alt="" coords="216,199,277,225"><area shape="rect" href="classbdm_1_1eEF.html" title="General conjugate exponential family posterior density." alt="" coords="109,279,189,305"></map> |
49 | | <center><font size="2">[<a href="graph_legend.html">legend</a>]</font></center></div> |
50 | | |
51 | | <p> |
52 | | <a href="classbdm_1_1ARX-members.html">List of all members.</a><table border="0" cellpadding="0" cellspacing="0"> |
| 40 | |
| 41 | <p> |
| 42 | <a href="classbdm_1_1ARX-members.html">List of all members.</a><hr><a name="_details"></a><h2>Detailed Description</h2> |
| 43 | Linear Autoregressive model with Gaussian noise. |
| 44 | <p> |
| 45 | Regression of the following kind: <p class="formulaDsp"> |
| 46 | <img class="formulaDsp" alt="\[ y_t = \theta_1 \psi_1 + \theta_2 + \psi_2 +\ldots + \theta_n \psi_n + r e_t \]" src="form_61.png"> |
| 47 | <p> |
| 48 | where unknown parameters <code>rv</code> are <img class="formulaInl" alt="$[\theta r]$" src="form_51.png">, regression vector <img class="formulaInl" alt="$\psi=\psi(y_{1:t},u_{1:t})$" src="form_52.png"> is a known function of past outputs and exogeneous variables <img class="formulaInl" alt="$u_t$" src="form_53.png">. Distrubances <img class="formulaInl" alt="$e_t$" src="form_54.png"> are supposed to be normally distributed: <p class="formulaDsp"> |
| 49 | <img class="formulaDsp" alt="\[ e_t \sim \mathcal{N}(0,1). \]" src="form_62.png"> |
| 50 | <p> |
| 51 | <p> |
| 52 | Extension for time-variant parameters <img class="formulaInl" alt="$\theta_t,r_t$" src="form_55.png"> may be achived using exponential forgetting (Kulhavy and Zarrop, 1993). In such a case, the forgetting factor <code>frg</code> <img class="formulaInl" alt="$\in <0,1>$" src="form_56.png"> should be given in the constructor. Time-invariant parameters are estimated for <code>frg</code> = 1. <table border="0" cellpadding="0" cellspacing="0"> |
55 | | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="44914d0b259204b3446db82b989bd626"></a><!-- doxytag: member="bdm::ARX::ARX" ref="44914d0b259204b3446db82b989bd626" args="(const RV &rv, const mat &V0, const double &nu0, const double frg0=1.0)" --> |
56 | | </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#44914d0b259204b3446db82b989bd626">ARX</a> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &<a class="el" href="classbdm_1_1BM.html#18d6db4af8ee42077741d9e3618153ca">rv</a>, const mat &V0, const double &nu0, const double frg0=1.0)</td></tr> |
| 55 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="43ed6114f04a3a8756fe2b42eaa35f98"></a><!-- doxytag: member="bdm::ARX::ARX" ref="43ed6114f04a3a8756fe2b42eaa35f98" args="(const double frg0=1.0)" --> |
| 56 | </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#43ed6114f04a3a8756fe2b42eaa35f98">ARX</a> (const double frg0=1.0)</td></tr> |
96 | | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7c8d1fe774fe1da50293d50ad8aef43d"></a><!-- doxytag: member="bdm::ARX::predictor" ref="7c8d1fe774fe1da50293d50ad8aef43d" args="(const RV &rv0, const vec &rgr) const " --> |
97 | | <a class="el" href="classbdm_1_1enorm.html">enorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#7c8d1fe774fe1da50293d50ad8aef43d">predictor</a> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &rv0, const vec &rgr) const </td></tr> |
| 95 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="749827323c034f11bec61b6e2fc3d42a"></a><!-- doxytag: member="bdm::ARX::epredictor" ref="749827323c034f11bec61b6e2fc3d42a" args="(const vec &rgr) const " --> |
| 96 | <a class="el" href="classbdm_1_1enorm.html">enorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#749827323c034f11bec61b6e2fc3d42a">epredictor</a> (const vec &rgr) const </td></tr> |
100 | | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5b73b70457f49ce4ad8660d729172dfd"></a><!-- doxytag: member="bdm::ARX::predictor" ref="5b73b70457f49ce4ad8660d729172dfd" args="(const RV &rv0) const " --> |
101 | | <a class="el" href="classbdm_1_1enorm.html">enorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#5b73b70457f49ce4ad8660d729172dfd">predictor</a> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &rv0) const </td></tr> |
102 | | |
103 | | <tr><td class="mdescLeft"> </td><td class="mdescRight">Constructs a predictive density (marginal density on data). <br></td></tr> |
104 | | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="02d9e91f21a700947a7b7eec1beed956"></a><!-- doxytag: member="bdm::ARX::predictor" ref="02d9e91f21a700947a7b7eec1beed956" args="(const RV &rv0, const RV &rvc0) const " --> |
105 | | <a class="el" href="classbdm_1_1mlnorm.html">mlnorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#02d9e91f21a700947a7b7eec1beed956">predictor</a> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &rv0, const <a class="el" href="classbdm_1_1RV.html">RV</a> &rvc0) const </td></tr> |
106 | | |
107 | | <tr><td class="mdescLeft"> </td><td class="mdescRight">conditional version of the predictor <br></td></tr> |
108 | | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="2ce8c6599497ffb94dfcb66d1fe7aca6"></a><!-- doxytag: member="bdm::ARX::predictor_student" ref="2ce8c6599497ffb94dfcb66d1fe7aca6" args="(const RV &rv0, const RV &rvc0) const " --> |
109 | | <a class="el" href="classbdm_1_1mlstudent.html">mlstudent</a> * </td><td class="memItemRight" valign="bottom"><b>predictor_student</b> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &rv0, const <a class="el" href="classbdm_1_1RV.html">RV</a> &rvc0) const </td></tr> |
| 99 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="4cdf5e2a7d3480ec31f6247ed4289b15"></a><!-- doxytag: member="bdm::ARX::epredictor" ref="4cdf5e2a7d3480ec31f6247ed4289b15" args="() const " --> |
| 100 | <a class="el" href="classbdm_1_1enorm.html">enorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#4cdf5e2a7d3480ec31f6247ed4289b15">epredictor</a> () const </td></tr> |
| 101 | |
| 102 | <tr><td class="mdescLeft"> </td><td class="mdescRight">Constructs a predictive density <img class="formulaInl" alt="$ f(d_{t+1} |d_{t}, \ldots d_{0}) $" src="form_112.png">. <br></td></tr> |
| 103 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="el" href="classbdm_1_1mlnorm.html">mlnorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#74fe8ae2d88bee8639510fd0eaf73513">predictor</a> () const </td></tr> |
| 104 | |
| 105 | <tr><td class="mdescLeft"> </td><td class="mdescRight">conditional version of the predictor <a href="#74fe8ae2d88bee8639510fd0eaf73513"></a><br></td></tr> |
| 106 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c6a2428a46407fe45b4c7a99069c0801"></a><!-- doxytag: member="bdm::ARX::predictor_student" ref="c6a2428a46407fe45b4c7a99069c0801" args="() const " --> |
| 107 | <a class="el" href="classbdm_1_1mlstudent.html">mlstudent</a> * </td><td class="memItemRight" valign="bottom"><b>predictor_student</b> () const </td></tr> |
115 | | const <a class="el" href="classbdm_1_1egiw.html">egiw</a> * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1ARX.html#ab2c55205a324e9d698fbd8ac229ad4f">_e</a> () const </td></tr> |
116 | | |
117 | | <tr><td class="mdescLeft"> </td><td class="mdescRight">Returns a pointer to the <a class="el" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> representing posterior density on parameters. Use with care! <br></td></tr> |
| 113 | const <a class="el" href="classbdm_1_1egiw.html">egiw</a> * </td><td class="memItemRight" valign="bottom"><b>_e</b> () const </td></tr> |
| 114 | |
| 115 | <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5912dbcf28ae711e30b08c2fa766a3e6"></a><!-- doxytag: member="bdm::ARX::_copy_" ref="5912dbcf28ae711e30b08c2fa766a3e6" args="(bool changerv=false)" --> |
| 116 | <a class="el" href="classbdm_1_1BMEF.html">BMEF</a> * </td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BMEF.html#5912dbcf28ae711e30b08c2fa766a3e6">_copy_</a> (bool changerv=false)</td></tr> |
| 117 | |
| 118 | <tr><td class="mdescLeft"> </td><td class="mdescRight">Flatten the posterior as if to keep nu0 data. <br></td></tr> |
| 119 | <tr><td colspan="2"><div class="groupHeader">Mathematical operations</div></td></tr> |
184 | | <hr><a name="_details"></a><h2>Detailed Description</h2> |
185 | | Linear Autoregressive model with Gaussian noise. |
186 | | <p> |
187 | | Regression of the following kind: <p class="formulaDsp"> |
188 | | <img class="formulaDsp" alt="\[ y_t = \theta_1 \psi_1 + \theta_2 + \psi_2 +\ldots + \theta_n \psi_n + r e_t \]" src="form_61.png"> |
189 | | <p> |
190 | | where unknown parameters <code>rv</code> are <img class="formulaInl" alt="$[\theta r]$" src="form_51.png">, regression vector <img class="formulaInl" alt="$\psi=\psi(y_{1:t},u_{1:t})$" src="form_52.png"> is a known function of past outputs and exogeneous variables <img class="formulaInl" alt="$u_t$" src="form_53.png">. Distrubances <img class="formulaInl" alt="$e_t$" src="form_54.png"> are supposed to be normally distributed: <p class="formulaDsp"> |
191 | | <img class="formulaDsp" alt="\[ e_t \sim \mathcal{N}(0,1). \]" src="form_62.png"> |
192 | | <p> |
193 | | <p> |
194 | | Extension for time-variant parameters <img class="formulaInl" alt="$\theta_t,r_t$" src="form_55.png"> may be achived using exponential forgetting (Kulhavy and Zarrop, 1993). In such a case, the forgetting factor <code>frg</code> <img class="formulaInl" alt="$\in <0,1>$" src="form_56.png"> should be given in the constructor. Time-invariant parameters are estimated for <code>frg</code> = 1. <hr><h2>Member Function Documentation</h2> |
| 175 | <hr><h2>Member Function Documentation</h2> |
246 | | <p>References <a class="el" href="libEF_8h-source.html#l00206">bdm::egiw::_nu()</a>, <a class="el" href="libEF_8h-source.html#l00202">bdm::egiw::_V()</a>, <a class="el" href="arx_8h-source.html#l00038">est</a>, <a class="el" href="libBM_8h-source.html#l00510">bdm::BM::evalll</a>, <a class="el" href="libEF_8h-source.html#l00077">bdm::BMEF::frg</a>, <a class="el" href="libEF_8h-source.html#l00079">bdm::BMEF::last_lognc</a>, <a class="el" href="libEF_8cpp-source.html#l00044">bdm::egiw::lognc()</a>, <a class="el" href="arx_8h-source.html#l00042">nu</a>, <a class="el" href="libDC_8cpp-source.html#l00058">ldmat::opupdt()</a>, <a class="el" href="libEF_8h-source.html#l00208">bdm::egiw::pow()</a>, and <a class="el" href="arx_8h-source.html#l00040">V</a>.</p> |
| 227 | <p>References <a class="el" href="libEF_8h-source.html#l00204">bdm::egiw::_nu()</a>, <a class="el" href="libEF_8h-source.html#l00202">bdm::egiw::_V()</a>, <a class="el" href="arx_8h-source.html#l00038">est</a>, <a class="el" href="libBM_8h-source.html#l00583">bdm::BM::evalll</a>, <a class="el" href="libEF_8h-source.html#l00077">bdm::BMEF::frg</a>, <a class="el" href="libEF_8h-source.html#l00079">bdm::BMEF::last_lognc</a>, <a class="el" href="libEF_8cpp-source.html#l00044">bdm::egiw::lognc()</a>, <a class="el" href="arx_8h-source.html#l00042">nu</a>, <a class="el" href="libDC_8cpp-source.html#l00058">ldmat::opupdt()</a>, <a class="el" href="libEF_8h-source.html#l00197">bdm::egiw::pow()</a>, and <a class="el" href="arx_8h-source.html#l00040">V</a>.</p> |
| 228 | |
| 229 | </div> |
| 230 | </div><p> |
| 231 | <a class="anchor" name="74fe8ae2d88bee8639510fd0eaf73513"></a><!-- doxytag: member="bdm::ARX::predictor" ref="74fe8ae2d88bee8639510fd0eaf73513" args="() const " --> |
| 232 | <div class="memitem"> |
| 233 | <div class="memproto"> |
| 234 | <table class="memname"> |
| 235 | <tr> |
| 236 | <td class="memname"><a class="el" href="classbdm_1_1mlnorm.html">mlnorm</a>< <a class="el" href="classldmat.html">ldmat</a> > * bdm::ARX::predictor </td> |
| 237 | <td>(</td> |
| 238 | <td class="paramname"> </td> |
| 239 | <td> ) </td> |
| 240 | <td> const<code> [virtual]</code></td> |
| 241 | </tr> |
| 242 | </table> |
| 243 | </div> |
| 244 | <div class="memdoc"> |
| 245 | |
| 246 | <p> |
| 247 | conditional version of the predictor |
| 248 | <p> |
| 249 | |
| 250 | <p> |
| 251 | <----------- TODO |
| 252 | <p>Reimplemented from <a class="el" href="classbdm_1_1BM.html#598b25e3f3d96a5bc00a5faeb5b3c912">bdm::BM</a>.</p> |
| 253 | |
| 254 | <p>References <a class="el" href="libBM_8h-source.html#l00272">bdm::epdf::dimension()</a>, <a class="el" href="arx_8h-source.html#l00038">est</a>, <a class="el" href="libEF_8cpp-source.html#l00107">bdm::egiw::mean_mat()</a>, <a class="el" href="libDC_8h-source.html#l00289">ldmat::rows()</a>, <a class="el" href="libEF_8h-source.html#l00762">bdm::mlnorm< sq_T >::set_parameters()</a>, and <a class="el" href="arx_8h-source.html#l00040">V</a>.</p> |