root/doc/html/classbdm_1_1migamma.html @ 270

Revision 270, 18.3 kB (checked in by smidl, 16 years ago)

Changes in the very root classes!
* rv and rvc are no longer compulsory,
* samplecond does not return ll
* BM has drv

Line 
1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
3<title>mixpp: bdm::migamma Class Reference</title>
4<link href="doxygen.css" rel="stylesheet" type="text/css">
5<link href="tabs.css" rel="stylesheet" type="text/css">
6</head><body>
7<!-- Generated by Doxygen 1.5.6 -->
8<div class="navigation" id="top">
9  <div class="tabs">
10    <ul>
11      <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
12      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
13      <li><a href="modules.html"><span>Modules</span></a></li>
14      <li><a href="namespaces.html"><span>Namespaces</span></a></li>
15      <li class="current"><a href="classes.html"><span>Classes</span></a></li>
16      <li><a href="files.html"><span>Files</span></a></li>
17    </ul>
18  </div>
19  <div class="tabs">
20    <ul>
21      <li><a href="classes.html"><span>Alphabetical&nbsp;List</span></a></li>
22      <li><a href="annotated.html"><span>Class&nbsp;List</span></a></li>
23      <li><a href="hierarchy.html"><span>Class&nbsp;Hierarchy</span></a></li>
24      <li><a href="functions.html"><span>Class&nbsp;Members</span></a></li>
25    </ul>
26  </div>
27  <div class="navpath"><a class="el" href="namespacebdm.html">bdm</a>::<a class="el" href="classbdm_1_1migamma.html">migamma</a>
28  </div>
29</div>
30<div class="contents">
31<h1>bdm::migamma Class Reference</h1><!-- doxytag: class="bdm::migamma" --><!-- doxytag: inherits="bdm::mEF" --><code>#include &lt;<a class="el" href="libEF_8h-source.html">libEF.h</a>&gt;</code>
32<p>
33<div class="dynheader">
34Inheritance diagram for bdm::migamma:</div>
35<div class="dynsection">
36<p><center><img src="classbdm_1_1migamma__inherit__graph.png" border="0" usemap="#bdm_1_1migamma__inherit__map" alt="Inheritance graph"></center>
37<map name="bdm_1_1migamma__inherit__map">
38<area shape="rect" href="classbdm_1_1migamma__fix.html" title="Inverse&#45;Gamma random walk around a fixed point." alt="" coords="5,316,149,343"><area shape="rect" href="classbdm_1_1mEF.html" title="Exponential family model." alt="" coords="35,161,120,188"><area shape="rect" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies." alt="" coords="32,84,123,111"><area shape="rect" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects." alt="" coords="21,7,133,33"></map>
39<center><font size="2">[<a href="graph_legend.html">legend</a>]</font></center></div>
40
41<p>
42<a href="classbdm_1_1migamma-members.html">List of all members.</a><hr><a name="_details"></a><h2>Detailed Description</h2>
43Inverse-Gamma random walk.
44<p>
45Mean value, <img class="formulaInl" alt="$ \mu $" src="form_113.png">, of this density is given by <code>rvc</code> . Standard deviation of the random walk is proportional to one <img class="formulaInl" alt="$ k $" src="form_114.png">-th the mean. This is achieved by setting <img class="formulaInl" alt="$ \alpha=\mu/k^2+2 $" src="form_115.png"> and <img class="formulaInl" alt="$ \beta=\mu(\alpha-1)$" src="form_116.png">.<p>
46The standard deviation of the walk is then: <img class="formulaInl" alt="$ \mu/\sqrt(k)$" src="form_117.png">. <table border="0" cellpadding="0" cellspacing="0">
47<tr><td></td></tr>
48<tr><td colspan="2"><br><h2>Public Member Functions</h2></td></tr>
49<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="8b10ab922e2a7bae2fb6bb3efc7b6151"></a><!-- doxytag: member="bdm::migamma::set_parameters" ref="8b10ab922e2a7bae2fb6bb3efc7b6151" args="(int len, double k0)" -->
50void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#8b10ab922e2a7bae2fb6bb3efc7b6151">set_parameters</a> (int len, double k0)</td></tr>
51
52<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Set value of <code>k</code>. <br></td></tr>
53<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7a34b1e2e3aa2250d7c0ed7df1665b8c"></a><!-- doxytag: member="bdm::migamma::condition" ref="7a34b1e2e3aa2250d7c0ed7df1665b8c" args="(const vec &amp;val)" -->
54void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#7a34b1e2e3aa2250d7c0ed7df1665b8c">condition</a> (const vec &amp;val)</td></tr>
55
56<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Update <code>ep</code> so that it represents this <a class="el" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> conditioned on <code>rvc</code> = cond. <br></td></tr>
57<tr><td colspan="2"><div class="groupHeader">Constructors</div></td></tr>
58<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="b109fa502a9ab521dfb48412fd45fca7"></a><!-- doxytag: member="bdm::migamma::migamma" ref="b109fa502a9ab521dfb48412fd45fca7" args="()" -->
59&nbsp;</td><td class="memItemRight" valign="bottom"><b>migamma</b> ()</td></tr>
60
61<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="a0126f741c6d2b6016df95a2410071e3"></a><!-- doxytag: member="bdm::migamma::migamma" ref="a0126f741c6d2b6016df95a2410071e3" args="(const migamma &amp;m)" -->
62&nbsp;</td><td class="memItemRight" valign="bottom"><b>migamma</b> (const <a class="el" href="classbdm_1_1migamma.html">migamma</a> &amp;m)</td></tr>
63
64<tr><td colspan="2"><div class="groupHeader">Matematical operations</div></td></tr>
65<tr><td class="memItemLeft" nowrap align="right" valign="top">virtual vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#f0c1db6fcbb3aae2dd6123884457a367">samplecond</a> (const vec &amp;cond)</td></tr>
66
67<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Returns a sample from the density conditioned on <code>cond</code>, <img class="formulaInl" alt="$x \sim epdf(rv|cond)$" src="form_10.png"><a href="#f0c1db6fcbb3aae2dd6123884457a367"></a><br></td></tr>
68<tr><td class="memItemLeft" nowrap align="right" valign="top">virtual mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#afe4185b26baeb03688202e254d3b005">samplecond_m</a> (const vec &amp;cond, int N)</td></tr>
69
70<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Returns.  <a href="#afe4185b26baeb03688202e254d3b005"></a><br></td></tr>
71<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="6336a8a72462e2a56a3989a220f18b1b"></a><!-- doxytag: member="bdm::migamma::evallogcond" ref="6336a8a72462e2a56a3989a220f18b1b" args="(const vec &amp;dt, const vec &amp;cond)" -->
72virtual double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#6336a8a72462e2a56a3989a220f18b1b">evallogcond</a> (const vec &amp;dt, const vec &amp;cond)</td></tr>
73
74<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Shortcut for conditioning and evaluation of the internal <a class="el" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>. In some cases, this operation can be implemented efficiently. <br></td></tr>
75<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0b0ed1ed663071bb7cf4a1349eb94fcb"></a><!-- doxytag: member="bdm::migamma::evallogcond_m" ref="0b0ed1ed663071bb7cf4a1349eb94fcb" args="(const mat &amp;Dt, const vec &amp;cond)" -->
76virtual vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#0b0ed1ed663071bb7cf4a1349eb94fcb">evallogcond_m</a> (const mat &amp;Dt, const vec &amp;cond)</td></tr>
77
78<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix version of evallogcond. <br></td></tr>
79<tr><td colspan="2"><div class="groupHeader">Access to attributes</div></td></tr>
80<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5571482d150fbcb72cc36f6694ce1a10"></a><!-- doxytag: member="bdm::migamma::_rv" ref="5571482d150fbcb72cc36f6694ce1a10" args="()" -->
81<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><b>_rv</b> ()</td></tr>
82
83<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="26001264236846897bd11e4baad47245"></a><!-- doxytag: member="bdm::migamma::_rvc" ref="26001264236846897bd11e4baad47245" args="()" -->
84<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><b>_rvc</b> ()</td></tr>
85
86<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="1c2bae3e1e90874e72941863974ec0ed"></a><!-- doxytag: member="bdm::migamma::dimension" ref="1c2bae3e1e90874e72941863974ec0ed" args="()" -->
87int&nbsp;</td><td class="memItemRight" valign="bottom"><b>dimension</b> ()</td></tr>
88
89<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="35e135910aed187b7290742f50e61bc8"></a><!-- doxytag: member="bdm::migamma::dimensionc" ref="35e135910aed187b7290742f50e61bc8" args="()" -->
90int&nbsp;</td><td class="memItemRight" valign="bottom"><b>dimensionc</b> ()</td></tr>
91
92<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="1892fe3933488942253679f068e9e7f6"></a><!-- doxytag: member="bdm::migamma::_epdf" ref="1892fe3933488942253679f068e9e7f6" args="()" -->
93<a class="el" href="classbdm_1_1epdf.html">epdf</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><b>_epdf</b> ()</td></tr>
94
95<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="05e843fd11c410a99dad2b88c55aca80"></a><!-- doxytag: member="bdm::migamma::_e" ref="05e843fd11c410a99dad2b88c55aca80" args="()" -->
96<a class="el" href="classbdm_1_1epdf.html">epdf</a> *&nbsp;</td><td class="memItemRight" valign="bottom"><b>_e</b> ()</td></tr>
97
98<tr><td colspan="2"><div class="groupHeader">Connection to other objects</div></td></tr>
99<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7631a5570e4ade1420065e8df78f4401"></a><!-- doxytag: member="bdm::migamma::set_rvc" ref="7631a5570e4ade1420065e8df78f4401" args="(const RV &amp;rvc0)" -->
100void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_rvc</b> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;rvc0)</td></tr>
101
102<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="18ac26bc2f96ae01ef4eb06178abbd75"></a><!-- doxytag: member="bdm::migamma::set_rv" ref="18ac26bc2f96ae01ef4eb06178abbd75" args="(const RV &amp;rv0)" -->
103void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_rv</b> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;rv0)</td></tr>
104
105<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="f8e3798150b42fd1f3e16ddbbe0e7045"></a><!-- doxytag: member="bdm::migamma::isnamed" ref="f8e3798150b42fd1f3e16ddbbe0e7045" args="()" -->
106bool&nbsp;</td><td class="memItemRight" valign="bottom"><b>isnamed</b> ()</td></tr>
107
108<tr><td colspan="2"><br><h2>Protected Attributes</h2></td></tr>
109<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="a31b39d4179551b593c9e0d7d756783a"></a><!-- doxytag: member="bdm::migamma::epdf" ref="a31b39d4179551b593c9e0d7d756783a" args="" -->
110<a class="el" href="classbdm_1_1eigamma.html">eigamma</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#a31b39d4179551b593c9e0d7d756783a">epdf</a></td></tr>
111
112<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Internal <a class="el" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> that arise by conditioning on <code>rvc</code>. <br></td></tr>
113<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="dc56bc9da542e0103ec16b9be8e5e38c"></a><!-- doxytag: member="bdm::migamma::k" ref="dc56bc9da542e0103ec16b9be8e5e38c" args="" -->
114double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#dc56bc9da542e0103ec16b9be8e5e38c">k</a></td></tr>
115
116<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Constant <img class="formulaInl" alt="$k$" src="form_41.png">. <br></td></tr>
117<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c9847093da59a9ba0ebb68d2c592f5dc"></a><!-- doxytag: member="bdm::migamma::_alpha" ref="c9847093da59a9ba0ebb68d2c592f5dc" args="" -->
118vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#c9847093da59a9ba0ebb68d2c592f5dc">_alpha</a></td></tr>
119
120<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of epdf.alpha <br></td></tr>
121<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0d854c047001b5465cf1ba21f52904b5"></a><!-- doxytag: member="bdm::migamma::_beta" ref="0d854c047001b5465cf1ba21f52904b5" args="" -->
122vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1migamma.html#0d854c047001b5465cf1ba21f52904b5">_beta</a></td></tr>
123
124<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of epdf.beta <br></td></tr>
125<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7c1900976ff13dbc09c9729b3bbff9e6"></a><!-- doxytag: member="bdm::migamma::dimc" ref="7c1900976ff13dbc09c9729b3bbff9e6" args="" -->
126int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6">dimc</a></td></tr>
127
128<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">dimension of the condition <br></td></tr>
129<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5a5f08950daa08b85b01ddf4e1c36288"></a><!-- doxytag: member="bdm::migamma::rvc" ref="5a5f08950daa08b85b01ddf4e1c36288" args="" -->
130<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288">rvc</a></td></tr>
131
132<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">random variable in condition <br></td></tr>
133<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5eea43c56d38e4441bfb30270db949c0"></a><!-- doxytag: member="bdm::migamma::ep" ref="5eea43c56d38e4441bfb30270db949c0" args="" -->
134<a class="el" href="classbdm_1_1epdf.html">epdf</a> *&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0">ep</a></td></tr>
135
136<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">pointer to internal <a class="el" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> <br></td></tr>
137</table>
138<hr><h2>Member Function Documentation</h2>
139<a class="anchor" name="f0c1db6fcbb3aae2dd6123884457a367"></a><!-- doxytag: member="bdm::migamma::samplecond" ref="f0c1db6fcbb3aae2dd6123884457a367" args="(const vec &amp;cond)" -->
140<div class="memitem">
141<div class="memproto">
142      <table class="memname">
143        <tr>
144          <td class="memname">virtual vec bdm::mpdf::samplecond           </td>
145          <td>(</td>
146          <td class="paramtype">const vec &amp;&nbsp;</td>
147          <td class="paramname"> <em>cond</em>          </td>
148          <td>&nbsp;)&nbsp;</td>
149          <td><code> [inline, virtual, inherited]</code></td>
150        </tr>
151      </table>
152</div>
153<div class="memdoc">
154
155<p>
156Returns a sample from the density conditioned on <code>cond</code>, <img class="formulaInl" alt="$x \sim epdf(rv|cond)$" src="form_10.png">.
157<p>
158<dl compact><dt><b>Parameters:</b></dt><dd>
159  <table border="0" cellspacing="2" cellpadding="0">
160    <tr><td valign="top"></td><td valign="top"><em>cond</em>&nbsp;</td><td>is numeric value of <code>rv</code> </td></tr>
161  </table>
162</dl>
163
164<p>References <a class="el" href="libBM_8h-source.html#l00315">bdm::mpdf::condition()</a>, <a class="el" href="libBM_8h-source.html#l00288">bdm::mpdf::ep</a>, and <a class="el" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10">bdm::epdf::sample()</a>.</p>
165
166<p>Referenced by <a class="el" href="libPF_8h-source.html#l00143">bdm::MPF&lt; BM_T &gt;::bayes()</a>, <a class="el" href="libPF_8cpp-source.html#l00007">bdm::PF::bayes()</a>, and <a class="el" href="libDS_8cpp-source.html#l00044">bdm::ArxDS::step()</a>.</p>
167
168</div>
169</div><p>
170<a class="anchor" name="afe4185b26baeb03688202e254d3b005"></a><!-- doxytag: member="bdm::migamma::samplecond_m" ref="afe4185b26baeb03688202e254d3b005" args="(const vec &amp;cond, int N)" -->
171<div class="memitem">
172<div class="memproto">
173      <table class="memname">
174        <tr>
175          <td class="memname">virtual mat bdm::mpdf::samplecond_m           </td>
176          <td>(</td>
177          <td class="paramtype">const vec &amp;&nbsp;</td>
178          <td class="paramname"> <em>cond</em>, </td>
179        </tr>
180        <tr>
181          <td class="paramkey"></td>
182          <td></td>
183          <td class="paramtype">int&nbsp;</td>
184          <td class="paramname"> <em>N</em></td><td>&nbsp;</td>
185        </tr>
186        <tr>
187          <td></td>
188          <td>)</td>
189          <td></td><td></td><td><code> [inline, virtual, inherited]</code></td>
190        </tr>
191      </table>
192</div>
193<div class="memdoc">
194
195<p>
196Returns.
197<p>
198<dl compact><dt><b>Parameters:</b></dt><dd>
199  <table border="0" cellspacing="2" cellpadding="0">
200    <tr><td valign="top"></td><td valign="top"><em>N</em>&nbsp;</td><td>samples from the density conditioned on <code>cond</code>, <img class="formulaInl" alt="$x \sim epdf(rv|cond)$" src="form_10.png">. </td></tr>
201    <tr><td valign="top"></td><td valign="top"><em>cond</em>&nbsp;</td><td>is numeric value of <code>rv</code> </td></tr>
202    <tr><td valign="top"></td><td valign="top"><em>ll</em>&nbsp;</td><td>is a return value of log-likelihood of the sample. </td></tr>
203  </table>
204</dl>
205
206<p>References <a class="el" href="libBM_8h-source.html#l00315">bdm::mpdf::condition()</a>, <a class="el" href="libBM_8h-source.html#l00272">bdm::epdf::dimension()</a>, <a class="el" href="libBM_8h-source.html#l00288">bdm::mpdf::ep</a>, and <a class="el" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10">bdm::epdf::sample()</a>.</p>
207
208</div>
209</div><p>
210<hr>The documentation for this class was generated from the following file:<ul>
211<li><a class="el" href="libEF_8h-source.html">libEF.h</a></ul>
212</div>
213<hr size="1"><address style="text-align: right;"><small>Generated on Wed Feb 11 23:34:07 2009 for mixpp by&nbsp;
214<a href="http://www.doxygen.org/index.html">
215<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address>
216</body>
217</html>
Note: See TracBrowser for help on using the browser.