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: work/git/mixpp/bdm/stat/libBM.h Source File</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 Page</span></a></li> |
---|
12 | <li><a href="pages.html"><span>Related Pages</span></a></li> |
---|
13 | <li><a href="modules.html"><span>Modules</span></a></li> |
---|
14 | <li><a href="annotated.html"><span>Classes</span></a></li> |
---|
15 | <li class="current"><a href="files.html"><span>Files</span></a></li> |
---|
16 | </ul> |
---|
17 | </div> |
---|
18 | <h1>work/git/mixpp/bdm/stat/libBM.h</h1><a href="libBM_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 |
---|
19 | <a name="l00017"></a>00017 <span class="preprocessor">#ifndef BM_H</span> |
---|
20 | <a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define BM_H</span> |
---|
21 | <a name="l00019"></a>00019 <span class="preprocessor"></span> |
---|
22 | <a name="l00020"></a>00020 <span class="preprocessor">#include <itpp/itbase.h></span> |
---|
23 | <a name="l00021"></a>00021 <span class="preprocessor">#include "../itpp_ext.h"</span> |
---|
24 | <a name="l00022"></a>00022 <span class="comment">//#include <std></span> |
---|
25 | <a name="l00023"></a>00023 |
---|
26 | <a name="l00024"></a>00024 <span class="keyword">using namespace </span>itpp; |
---|
27 | <a name="l00025"></a>00025 |
---|
28 | <a name="l00027"></a><a class="code" href="classstr.html">00027</a> <span class="keyword">class </span><a class="code" href="classstr.html" title="Structure of RV (used internally), i.e. expanded RVs.">str</a> { |
---|
29 | <a name="l00028"></a>00028 <span class="keyword">public</span>: |
---|
30 | <a name="l00030"></a><a class="code" href="classstr.html#8c4e7241505999c272eb15e586ba7923">00030</a> ivec <a class="code" href="classstr.html#8c4e7241505999c272eb15e586ba7923" title="vector id ids (non-unique!)">ids</a>; |
---|
31 | <a name="l00032"></a><a class="code" href="classstr.html#a55be419d33638dd0b2de3ecd64a4f27">00032</a> ivec <a class="code" href="classstr.html#a55be419d33638dd0b2de3ecd64a4f27" title="vector of times">times</a>; |
---|
32 | <a name="l00034"></a><a class="code" href="classstr.html#11ba7967bca045fe27fab8e70b6c0598">00034</a> <a class="code" href="classstr.html#11ba7967bca045fe27fab8e70b6c0598" title="Default constructor.">str</a> ( ivec ids0, ivec times0 ) :<a class="code" href="classstr.html#8c4e7241505999c272eb15e586ba7923" title="vector id ids (non-unique!)">ids</a> ( ids0 ),<a class="code" href="classstr.html#a55be419d33638dd0b2de3ecd64a4f27" title="vector of times">times</a> ( times0 ) { |
---|
33 | <a name="l00035"></a>00035 it_assert_debug ( times0.length() ==ids0.length(),<span class="stringliteral">"Incompatible input"</span> ); |
---|
34 | <a name="l00036"></a>00036 }; |
---|
35 | <a name="l00037"></a>00037 }; |
---|
36 | <a name="l00038"></a>00038 |
---|
37 | <a name="l00045"></a><a class="code" href="classRV.html">00045</a> <span class="keyword">class </span><a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> { |
---|
38 | <a name="l00046"></a>00046 <span class="keyword">protected</span>: |
---|
39 | <a name="l00048"></a><a class="code" href="classRV.html#f2396efe2e6d656ebe9e6f06018cbe25">00048</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#f2396efe2e6d656ebe9e6f06018cbe25" title="size = sum of sizes">tsize</a>; |
---|
40 | <a name="l00050"></a><a class="code" href="classRV.html#0d7b36e2bbccf880c8fcf1e8cc43c1a9">00050</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#0d7b36e2bbccf880c8fcf1e8cc43c1a9" title="len = number of individual rvs">len</a>; |
---|
41 | <a name="l00052"></a><a class="code" href="classRV.html#1bd7165140f4b880a7f344bbb1c433f9">00052</a> ivec <a class="code" href="classRV.html#1bd7165140f4b880a7f344bbb1c433f9" title="Vector of unique IDs.">ids</a>; |
---|
42 | <a name="l00054"></a><a class="code" href="classRV.html#c878aab13f34b420e1eb5b485563682b">00054</a> ivec <a class="code" href="classRV.html#c878aab13f34b420e1eb5b485563682b" title="Vector of sizes.">sizes</a>; |
---|
43 | <a name="l00056"></a><a class="code" href="classRV.html#646e530c47a2dd38254b778d9f11ac89">00056</a> ivec <a class="code" href="classRV.html#646e530c47a2dd38254b778d9f11ac89" title="Vector of shifts from current time.">times</a>; |
---|
44 | <a name="l00058"></a><a class="code" href="classRV.html#df5d0030b277a1db2f1fd5fb79152acb">00058</a> Array<std::string> <a class="code" href="classRV.html#df5d0030b277a1db2f1fd5fb79152acb" title="Array of names.">names</a>; |
---|
45 | <a name="l00059"></a>00059 |
---|
46 | <a name="l00060"></a>00060 <span class="keyword">private</span>: |
---|
47 | <a name="l00062"></a>00062 <span class="keywordtype">void</span> init ( ivec in_ids, Array<std::string> in_names, ivec in_sizes, ivec in_times ); |
---|
48 | <a name="l00063"></a>00063 <span class="keyword">public</span>: |
---|
49 | <a name="l00065"></a>00065 <a class="code" href="classRV.html#6a22f5c06d6403870d4e414b1993bf9e" title="Constructor of empty RV.">RV</a> ( Array<std::string> in_names, ivec in_sizes, ivec in_times ); |
---|
50 | <a name="l00067"></a>00067 <a class="code" href="classRV.html#6a22f5c06d6403870d4e414b1993bf9e" title="Constructor of empty RV.">RV</a> ( Array<std::string> in_names, ivec in_sizes ); |
---|
51 | <a name="l00069"></a>00069 <a class="code" href="classRV.html#6a22f5c06d6403870d4e414b1993bf9e" title="Constructor of empty RV.">RV</a> ( Array<std::string> in_names ); |
---|
52 | <a name="l00071"></a>00071 <a class="code" href="classRV.html#6a22f5c06d6403870d4e414b1993bf9e" title="Constructor of empty RV.">RV</a> (); |
---|
53 | <a name="l00072"></a>00072 |
---|
54 | <a name="l00074"></a>00074 <span class="keyword">friend</span> std::ostream &<a class="code" href="classRV.html#cad6f15535d10437f6690cb4492176a8" title="Printing output e.g. for debugging.">operator<< </a>( std::ostream &os, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv ); |
---|
55 | <a name="l00075"></a>00075 |
---|
56 | <a name="l00077"></a><a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0">00077</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classRV.html#f2396efe2e6d656ebe9e6f06018cbe25" title="size = sum of sizes">tsize</a>;} ; |
---|
57 | <a name="l00079"></a><a class="code" href="classRV.html#c114a6f3ff06796cc2f4dacba74291eb">00079</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#c114a6f3ff06796cc2f4dacba74291eb" title="Return length (number of entries) of the RV.">length</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classRV.html#0d7b36e2bbccf880c8fcf1e8cc43c1a9" title="len = number of individual rvs">len</a>;} ; |
---|
58 | <a name="l00080"></a>00080 |
---|
59 | <a name="l00081"></a>00081 <span class="comment">//TODO why not inline and later??</span> |
---|
60 | <a name="l00082"></a>00082 |
---|
61 | <a name="l00084"></a>00084 ivec <a class="code" href="classRV.html#0c493a7a897f4d75682f4264e40779a1" title="Find indexes of self in another rv,.">findself</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2 ) <span class="keyword">const</span>; |
---|
62 | <a name="l00086"></a>00086 <span class="keywordtype">bool</span> <a class="code" href="classRV.html#82cc7fa22ce91cf461a522ad62014c65" title="Compare if rv2 is identical to this RV.">equal</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2 ) <span class="keyword">const</span>; |
---|
63 | <a name="l00088"></a>00088 <span class="keywordtype">bool</span> <span class="keyword">add</span> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2 ); |
---|
64 | <a name="l00090"></a>00090 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classRV.html#b9d175c327c21488b1e2fb756a84e149" title="Subtract another variable from the current one.">subt</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2 ) <span class="keyword">const</span>; |
---|
65 | <a name="l00092"></a>00092 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classRV.html#b26d25a122be3303806af7726a625a0f" title="Select only variables at indeces ind.">subselect</a> ( <span class="keyword">const</span> ivec &ind ) <span class="keyword">const</span>; |
---|
66 | <a name="l00094"></a>00094 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classRV.html#248dfb2387976db2658489347669278c" title="Select only variables at indeces ind.">operator() </a>( <span class="keyword">const</span> ivec &ind ) <span class="keyword">const</span>; |
---|
67 | <a name="l00096"></a>00096 <span class="keywordtype">void</span> <a class="code" href="classRV.html#e0a158d53ab9430627241f844f4c9a60" title="Shift time shifted by delta.">t</a> ( <span class="keywordtype">int</span> delta ); |
---|
68 | <a name="l00098"></a>00098 <a class="code" href="classstr.html" title="Structure of RV (used internally), i.e. expanded RVs.">str</a> <a class="code" href="classRV.html#52f3ade796814da627ff5a44d00986be" title="generate str from rv, by expanding sizes">tostr</a>() <span class="keyword">const</span>; |
---|
69 | <a name="l00101"></a>00101 ivec <a class="code" href="classRV.html#bb724fa4e2d9ed7bfd0993b5975018a4">dataind</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &crv ) <span class="keyword">const</span>; |
---|
70 | <a name="l00104"></a>00104 <span class="keywordtype">void</span> <a class="code" href="classRV.html#bb724fa4e2d9ed7bfd0993b5975018a4">dataind</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2, ivec &selfi, ivec &rv2i ) <span class="keyword">const</span>; |
---|
71 | <a name="l00105"></a>00105 |
---|
72 | <a name="l00107"></a><a class="code" href="classRV.html#b88ef5124713c9e30fd9446846b3f044">00107</a> Array<std::string>& <a class="code" href="classRV.html#b88ef5124713c9e30fd9446846b3f044" title="access function">_names</a>() {<span class="keywordflow">return</span> <a class="code" href="classRV.html#df5d0030b277a1db2f1fd5fb79152acb" title="Array of names.">names</a>;}; |
---|
73 | <a name="l00108"></a>00108 |
---|
74 | <a name="l00110"></a><a class="code" href="classRV.html#4a0069bc69f5a98d7ce9ed59077f2e6c">00110</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#4a0069bc69f5a98d7ce9ed59077f2e6c" title="access function">id</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classRV.html#1bd7165140f4b880a7f344bbb1c433f9" title="Vector of unique IDs.">ids</a> ( at );}; |
---|
75 | <a name="l00112"></a><a class="code" href="classRV.html#20db23a5585d64a6e24404cc8379bc28">00112</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#20db23a5585d64a6e24404cc8379bc28" title="access function">size</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classRV.html#c878aab13f34b420e1eb5b485563682b" title="Vector of sizes.">sizes</a> ( at );}; |
---|
76 | <a name="l00114"></a><a class="code" href="classRV.html#af068cf478864fe1a4fbfc6b4693b088">00114</a> <span class="keywordtype">int</span> <a class="code" href="classRV.html#af068cf478864fe1a4fbfc6b4693b088" title="access function">time</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classRV.html#646e530c47a2dd38254b778d9f11ac89" title="Vector of shifts from current time.">times</a> ( at );}; |
---|
77 | <a name="l00116"></a><a class="code" href="classRV.html#3ae4dd6e368f814a02710636ac8e1721">00116</a> std::string <a class="code" href="classRV.html#3ae4dd6e368f814a02710636ac8e1721" title="access function">name</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classRV.html#df5d0030b277a1db2f1fd5fb79152acb" title="Array of names.">names</a> ( at );}; |
---|
78 | <a name="l00117"></a>00117 |
---|
79 | <a name="l00119"></a><a class="code" href="classRV.html#f58db96c9f5b8a538e26f154b1047c62">00119</a> <span class="keywordtype">void</span> <a class="code" href="classRV.html#f58db96c9f5b8a538e26f154b1047c62" title="access function">set_id</a> ( <span class="keywordtype">int</span> at, <span class="keywordtype">int</span> id0 ) {<a class="code" href="classRV.html#1bd7165140f4b880a7f344bbb1c433f9" title="Vector of unique IDs.">ids</a> ( at ) =id0;}; |
---|
80 | <a name="l00121"></a><a class="code" href="classRV.html#70b24c39c5130b1e4753fa2eef495433">00121</a> <span class="keywordtype">void</span> <a class="code" href="classRV.html#70b24c39c5130b1e4753fa2eef495433" title="access function">set_size</a> ( <span class="keywordtype">int</span> at, <span class="keywordtype">int</span> size0 ) {<a class="code" href="classRV.html#c878aab13f34b420e1eb5b485563682b" title="Vector of sizes.">sizes</a> ( at ) =size0; <a class="code" href="classRV.html#f2396efe2e6d656ebe9e6f06018cbe25" title="size = sum of sizes">tsize</a>=sum ( <a class="code" href="classRV.html#c878aab13f34b420e1eb5b485563682b" title="Vector of sizes.">sizes</a> );}; |
---|
81 | <a name="l00123"></a><a class="code" href="classRV.html#af668289c016d12ac2a9ddb85416006b">00123</a> <span class="keywordtype">void</span> <a class="code" href="classRV.html#af668289c016d12ac2a9ddb85416006b" title="access function">set_time</a> ( <span class="keywordtype">int</span> at, <span class="keywordtype">int</span> time0 ) {<a class="code" href="classRV.html#646e530c47a2dd38254b778d9f11ac89" title="Vector of shifts from current time.">times</a> ( at ) =time0;}; |
---|
82 | <a name="l00124"></a>00124 |
---|
83 | <a name="l00126"></a>00126 <span class="keywordtype">void</span> <a class="code" href="classRV.html#fc4a3175ff678cba7ae4109654f66ce5" title="Assign unused ids to this rv.">newids</a>(); |
---|
84 | <a name="l00127"></a>00127 }; |
---|
85 | <a name="l00128"></a>00128 |
---|
86 | <a name="l00130"></a>00130 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="group__core.html#g33c114e83980d883c5b211c47d5322a4" title="Concat two random variables.">concat</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv1, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv2 ); |
---|
87 | <a name="l00131"></a>00131 |
---|
88 | <a name="l00133"></a>00133 <span class="keyword">extern</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="group__core.html#g9ea0562597470f6058ec209ee72db5fa" title="Default empty RV that can be used as default argument.">RV0</a>; |
---|
89 | <a name="l00134"></a>00134 |
---|
90 | <a name="l00136"></a>00136 |
---|
91 | <a name="l00137"></a><a class="code" href="classfnc.html">00137</a> <span class="keyword">class </span><a class="code" href="classfnc.html" title="Class representing function of variable represented by rv.">fnc</a> { |
---|
92 | <a name="l00138"></a>00138 <span class="keyword">protected</span>: |
---|
93 | <a name="l00140"></a><a class="code" href="classfnc.html#22d51d10a7901331167f64f80d1af8e9">00140</a> <span class="keywordtype">int</span> <a class="code" href="classfnc.html#22d51d10a7901331167f64f80d1af8e9" title="Length of the output vector.">dimy</a>; |
---|
94 | <a name="l00141"></a>00141 <span class="keyword">public</span>: |
---|
95 | <a name="l00143"></a><a class="code" href="classfnc.html#8a1da32233c3006a4239cf601c413538">00143</a> <a class="code" href="classfnc.html#8a1da32233c3006a4239cf601c413538" title="default constructor">fnc</a> ( <span class="keywordtype">int</span> dy ) :<a class="code" href="classfnc.html#22d51d10a7901331167f64f80d1af8e9" title="Length of the output vector.">dimy</a> ( dy ) {}; |
---|
96 | <a name="l00145"></a><a class="code" href="classfnc.html#7fad2e6cb35909a165ec5a53def4a18c">00145</a> <span class="keyword">virtual</span> vec <a class="code" href="classfnc.html#7fad2e6cb35909a165ec5a53def4a18c" title="function evaluates numerical value of at cond ">eval</a> ( <span class="keyword">const</span> vec &cond ) { |
---|
97 | <a name="l00146"></a>00146 <span class="keywordflow">return</span> vec ( 0 ); |
---|
98 | <a name="l00147"></a>00147 }; |
---|
99 | <a name="l00148"></a>00148 |
---|
100 | <a name="l00150"></a><a class="code" href="classfnc.html#d28babf615a7240c34ee57779250f602">00150</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classfnc.html#d28babf615a7240c34ee57779250f602" title="function substitutes given value into an appropriate position">condition</a>(<span class="keyword">const</span> vec &val){}; |
---|
101 | <a name="l00151"></a>00151 |
---|
102 | <a name="l00153"></a><a class="code" href="classfnc.html#a8891973d0ca48ce38e1886df45ca298">00153</a> <span class="keywordtype">int</span> <a class="code" href="classfnc.html#a8891973d0ca48ce38e1886df45ca298" title="access function">_dimy</a>()<span class="keyword"> const</span>{<span class="keywordflow">return</span> <a class="code" href="classfnc.html#22d51d10a7901331167f64f80d1af8e9" title="Length of the output vector.">dimy</a>;} |
---|
103 | <a name="l00154"></a>00154 |
---|
104 | <a name="l00156"></a><a class="code" href="classfnc.html#17164c202f6feee3d708b8caab6306ab">00156</a> <span class="keyword">virtual</span> <a class="code" href="classfnc.html#17164c202f6feee3d708b8caab6306ab" title="Destructor for future use;.">~fnc</a>() {}; |
---|
105 | <a name="l00157"></a>00157 }; |
---|
106 | <a name="l00158"></a>00158 |
---|
107 | <a name="l00159"></a>00159 <span class="keyword">class </span><a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>; |
---|
108 | <a name="l00160"></a>00160 |
---|
109 | <a name="l00162"></a>00162 |
---|
110 | <a name="l00163"></a><a class="code" href="classepdf.html">00163</a> <span class="keyword">class </span><a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> { |
---|
111 | <a name="l00164"></a>00164 <span class="keyword">protected</span>: |
---|
112 | <a name="l00166"></a><a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9">00166</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>; |
---|
113 | <a name="l00167"></a>00167 <span class="keyword">public</span>: |
---|
114 | <a name="l00169"></a><a class="code" href="classepdf.html#d8eb760037b3bad5a0d64081606697cd">00169</a> <a class="code" href="classepdf.html#d8eb760037b3bad5a0d64081606697cd" title="default constructor">epdf</a>() :<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> ( ) {}; |
---|
115 | <a name="l00170"></a>00170 |
---|
116 | <a name="l00172"></a><a class="code" href="classepdf.html#c95b1a27a8dd9507bb9a5a3cb2809c7a">00172</a> <a class="code" href="classepdf.html#d8eb760037b3bad5a0d64081606697cd" title="default constructor">epdf</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv0 ) :<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> ( rv0 ) {}; |
---|
117 | <a name="l00173"></a>00173 |
---|
118 | <a name="l00174"></a>00174 <span class="comment">// //! Returns the required moment of the epdf</span> |
---|
119 | <a name="l00175"></a>00175 <span class="comment">// virtual vec moment ( const int order = 1 );</span> |
---|
120 | <a name="l00176"></a>00176 |
---|
121 | <a name="l00178"></a>00178 <span class="keyword">virtual</span> vec <a class="code" href="classepdf.html#8019654e494bf5e458f6fb947e11b262" title="Returns a sample, from density .">sample</a> () <span class="keyword">const</span> =0; |
---|
122 | <a name="l00180"></a>00180 <span class="keyword">virtual</span> mat <a class="code" href="classepdf.html#76608914c3b19e150292d5c56e93e508" title="Returns N samples from density .">sample_m</a> ( <span class="keywordtype">int</span> N ) <span class="keyword">const</span>; |
---|
123 | <a name="l00181"></a>00181 |
---|
124 | <a name="l00183"></a>00183 <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classepdf.html#e1996af1da1fa1214270066a96ca113e" title="Compute log-probability of argument val.">evallog</a> ( <span class="keyword">const</span> vec &val ) <span class="keyword">const</span> =0; |
---|
125 | <a name="l00184"></a>00184 |
---|
126 | <a name="l00186"></a><a class="code" href="classepdf.html#2495a04bbacb9b55fe5a3a59b78affca">00186</a> <span class="keyword">virtual</span> vec <a class="code" href="classepdf.html#2495a04bbacb9b55fe5a3a59b78affca" title="Compute log-probability of multiple values argument val.">evallog_m</a> ( <span class="keyword">const</span> mat &Val )<span class="keyword"> const </span>{ |
---|
127 | <a name="l00187"></a>00187 vec x ( Val.cols() ); |
---|
128 | <a name="l00188"></a>00188 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i<Val.cols();i++ ) {x ( i ) =<a class="code" href="classepdf.html#e1996af1da1fa1214270066a96ca113e" title="Compute log-probability of argument val.">evallog</a> ( Val.get_col ( i ) ) ;} |
---|
129 | <a name="l00189"></a>00189 <span class="keywordflow">return</span> x; |
---|
130 | <a name="l00190"></a>00190 } |
---|
131 | <a name="l00192"></a><a class="code" href="classepdf.html#e87dc8260a5c37bc1b03eb66174741a0">00192</a> <span class="keyword">virtual</span> <a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>* <a class="code" href="classepdf.html#e87dc8260a5c37bc1b03eb66174741a0" title="Return conditional density on the given RV, the remaining rvs will be in conditioning...">condition</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> )<span class="keyword"> const </span>{it_warning ( <span class="stringliteral">"Not implemented"</span> ); <span class="keywordflow">return</span> NULL;} |
---|
132 | <a name="l00194"></a><a class="code" href="classepdf.html#38de9f59b65ee06028554f3f74b66025">00194</a> <span class="keyword">virtual</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classepdf.html#38de9f59b65ee06028554f3f74b66025" title="Return marginal density on the given RV, the remainig rvs are intergrated out.">marginal</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> )<span class="keyword"> const </span>{it_warning ( <span class="stringliteral">"Not implemented"</span> ); <span class="keywordflow">return</span> NULL;} |
---|
133 | <a name="l00195"></a>00195 |
---|
134 | <a name="l00197"></a>00197 <span class="keyword">virtual</span> vec <a class="code" href="classepdf.html#bf0a070a8f3e67a93604ca724638b870" title="return expected value">mean</a>() <span class="keyword">const</span> =0; |
---|
135 | <a name="l00198"></a>00198 |
---|
136 | <a name="l00200"></a>00200 <span class="keyword">virtual</span> vec <a class="code" href="classepdf.html#29f17cea338692f4cb6d4ac9f7d89e1f" title="return expected variance (not covariance!)">variance</a>() <span class="keyword">const</span> = 0; |
---|
137 | <a name="l00201"></a>00201 |
---|
138 | <a name="l00203"></a><a class="code" href="classepdf.html#0a322dd106f04c0a2915e3d4f4227396">00203</a> <span class="keyword">virtual</span> <a class="code" href="classepdf.html#0a322dd106f04c0a2915e3d4f4227396" title="Destructor for future use;.">~epdf</a>() {}; |
---|
139 | <a name="l00205"></a><a class="code" href="classepdf.html#ca0d32aabb4cbba347e0c37fe8607562">00205</a> <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classepdf.html#ca0d32aabb4cbba347e0c37fe8607562" title="access function, possibly dangerous!">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>;} |
---|
140 | <a name="l00207"></a><a class="code" href="classepdf.html#7fb94ce90d1ac7077d29f7d6a6c3e0a5">00207</a> <span class="keywordtype">void</span> <a class="code" href="classepdf.html#7fb94ce90d1ac7077d29f7d6a6c3e0a5" title="modifier function - useful when copying epdfs">_renewrv</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &in_rv ) {<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>=in_rv;} |
---|
141 | <a name="l00209"></a>00209 }; |
---|
142 | <a name="l00210"></a>00210 |
---|
143 | <a name="l00211"></a>00211 |
---|
144 | <a name="l00213"></a>00213 <span class="comment">//TODO Samplecond can be generalized</span> |
---|
145 | <a name="l00214"></a>00214 |
---|
146 | <a name="l00215"></a><a class="code" href="classmpdf.html">00215</a> <span class="keyword">class </span><a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> { |
---|
147 | <a name="l00216"></a>00216 <span class="keyword">protected</span>: |
---|
148 | <a name="l00218"></a><a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb">00218</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb" title="modeled random variable">rv</a>; |
---|
149 | <a name="l00220"></a><a class="code" href="classmpdf.html#acb7dda792b3cd5576f39fa3129abbab">00220</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classmpdf.html#acb7dda792b3cd5576f39fa3129abbab" title="random variable in condition">rvc</a>; |
---|
150 | <a name="l00222"></a><a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c">00222</a> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>; |
---|
151 | <a name="l00223"></a>00223 <span class="keyword">public</span>: |
---|
152 | <a name="l00224"></a>00224 |
---|
153 | <a name="l00226"></a><a class="code" href="classmpdf.html#3f172b79ec4a5ebc87898a5381141f1b">00226</a> <span class="keyword">virtual</span> vec <a class="code" href="classmpdf.html#3f172b79ec4a5ebc87898a5381141f1b" title="Returns a sample from the density conditioned on cond, .">samplecond</a> ( <span class="keyword">const</span> vec &cond, <span class="keywordtype">double</span> &ll ) { |
---|
154 | <a name="l00227"></a>00227 this-><a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond ); |
---|
155 | <a name="l00228"></a>00228 vec temp= <a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#8019654e494bf5e458f6fb947e11b262" title="Returns a sample, from density .">sample</a>(); |
---|
156 | <a name="l00229"></a>00229 ll=<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#e1996af1da1fa1214270066a96ca113e" title="Compute log-probability of argument val.">evallog</a> ( temp );<span class="keywordflow">return</span> temp; |
---|
157 | <a name="l00230"></a>00230 }; |
---|
158 | <a name="l00232"></a><a class="code" href="classmpdf.html#b1dae6171ee39a6a05976c7b1007a3c5">00232</a> <span class="keyword">virtual</span> mat <a class="code" href="classmpdf.html#b1dae6171ee39a6a05976c7b1007a3c5" title="Returns.">samplecond_m</a> ( <span class="keyword">const</span> vec &cond, vec &ll, <span class="keywordtype">int</span> N ) { |
---|
159 | <a name="l00233"></a>00233 this-><a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond ); |
---|
160 | <a name="l00234"></a>00234 mat temp ( <a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb" title="modeled random variable">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>(),N ); vec smp ( <a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb" title="modeled random variable">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>() ); |
---|
161 | <a name="l00235"></a>00235 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i<N;i++ ) {smp=<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#8019654e494bf5e458f6fb947e11b262" title="Returns a sample, from density .">sample</a>() ;temp.set_col ( i, smp );ll ( i ) =<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#e1996af1da1fa1214270066a96ca113e" title="Compute log-probability of argument val.">evallog</a> ( smp );} |
---|
162 | <a name="l00236"></a>00236 <span class="keywordflow">return</span> temp; |
---|
163 | <a name="l00237"></a>00237 }; |
---|
164 | <a name="l00239"></a><a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83">00239</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( <span class="keyword">const</span> vec &cond ) {it_error ( <span class="stringliteral">"Not implemented"</span> );}; |
---|
165 | <a name="l00240"></a>00240 |
---|
166 | <a name="l00242"></a><a class="code" href="classmpdf.html#2ef8a6374029d990a678782f6decebbe">00242</a> <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classmpdf.html#2ef8a6374029d990a678782f6decebbe" title="Shortcut for conditioning and evaluation of the internal epdf. In some cases, this...">evallogcond</a> ( <span class="keyword">const</span> vec &dt, <span class="keyword">const</span> vec &cond ) {<span class="keywordtype">double</span> tmp; this-><a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );tmp = <a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#e1996af1da1fa1214270066a96ca113e" title="Compute log-probability of argument val.">evallog</a> ( dt ); it_assert_debug(std::isfinite(tmp),<span class="stringliteral">"Infinite value"</span>); <span class="keywordflow">return</span> tmp; |
---|
167 | <a name="l00243"></a>00243 }; |
---|
168 | <a name="l00244"></a>00244 |
---|
169 | <a name="l00246"></a><a class="code" href="classmpdf.html#95fcff214848f66f1b489459370573fa">00246</a> <span class="keyword">virtual</span> vec <a class="code" href="classmpdf.html#95fcff214848f66f1b489459370573fa" title="Matrix version of evallogcond.">evallogcond_m</a> ( <span class="keyword">const</span> mat &Dt, <span class="keyword">const</span> vec &cond ) {this-><a class="code" href="classmpdf.html#0f95a0cc6ab40611f46804682446ed83" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );<span class="keywordflow">return</span> <a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>-><a class="code" href="classepdf.html#2495a04bbacb9b55fe5a3a59b78affca" title="Compute log-probability of multiple values argument val.">evallog_m</a> ( Dt );}; |
---|
170 | <a name="l00247"></a>00247 |
---|
171 | <a name="l00249"></a><a class="code" href="classmpdf.html#6788be9f3a888796499c5293a318fcfb">00249</a> <span class="keyword">virtual</span> <a class="code" href="classmpdf.html#6788be9f3a888796499c5293a318fcfb" title="Destructor for future use;.">~mpdf</a>() {}; |
---|
172 | <a name="l00250"></a>00250 |
---|
173 | <a name="l00252"></a><a class="code" href="classmpdf.html#581ecf362185d37c08bb31cb9d046d6f">00252</a> <a class="code" href="classmpdf.html#581ecf362185d37c08bb31cb9d046d6f" title="Default constructor.">mpdf</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv0, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rvc0 ) :<a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb" title="modeled random variable">rv</a> ( rv0 ),<a class="code" href="classmpdf.html#acb7dda792b3cd5576f39fa3129abbab" title="random variable in condition">rvc</a> ( rvc0 ) {}; |
---|
174 | <a name="l00254"></a><a class="code" href="classmpdf.html#15ef062183b1ccdf794732d5fa0b77cd">00254</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classmpdf.html#15ef062183b1ccdf794732d5fa0b77cd" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classmpdf.html#acb7dda792b3cd5576f39fa3129abbab" title="random variable in condition">rvc</a>;} |
---|
175 | <a name="l00256"></a><a class="code" href="classmpdf.html#71256ffb5fbd08f41d650e606a5bd585">00256</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classmpdf.html#71256ffb5fbd08f41d650e606a5bd585" title="access function">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classmpdf.html#f6687c07ff07d47812dd565368ca59eb" title="modeled random variable">rv</a>;} |
---|
176 | <a name="l00258"></a><a class="code" href="classmpdf.html#e17780ee5b2cfe05922a6c56af1462f8">00258</a> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>& <a class="code" href="classmpdf.html#e17780ee5b2cfe05922a6c56af1462f8" title="access function">_epdf</a>() {<span class="keywordflow">return</span> *<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>;} |
---|
177 | <a name="l00260"></a><a class="code" href="classmpdf.html#75ded3b0f657cd7da6590691a810963c">00260</a> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classmpdf.html#75ded3b0f657cd7da6590691a810963c" title="access function">_e</a>() {<span class="keywordflow">return</span> <a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>;} |
---|
178 | <a name="l00261"></a>00261 }; |
---|
179 | <a name="l00262"></a>00262 |
---|
180 | <a name="l00265"></a><a class="code" href="classdatalink__e2e.html">00265</a> <span class="keyword">class </span><a class="code" href="classdatalink__e2e.html">datalink_e2e</a> { |
---|
181 | <a name="l00266"></a>00266 <span class="keyword">protected</span>: |
---|
182 | <a name="l00268"></a><a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216">00268</a> <span class="keywordtype">int</span> <a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216" title="Remember how long val should be.">valsize</a>; |
---|
183 | <a name="l00270"></a><a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8">00270</a> <span class="keywordtype">int</span> <a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8" title="Remember how long val of &quot;Up&quot; should be.">valupsize</a>; |
---|
184 | <a name="l00272"></a><a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae">00272</a> ivec <a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a>; |
---|
185 | <a name="l00273"></a>00273 <span class="keyword">public</span>: |
---|
186 | <a name="l00275"></a><a class="code" href="classdatalink__e2e.html#366d0aca52e56d4108f85390f87c1b9f">00275</a> <a class="code" href="classdatalink__e2e.html#366d0aca52e56d4108f85390f87c1b9f" title="Constructor.">datalink_e2e</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up ) : |
---|
187 | <a name="l00276"></a>00276 <a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216" title="Remember how long val should be.">valsize</a> ( rv.count() ), <a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8" title="Remember how long val of &quot;Up&quot; should be.">valupsize</a> ( rv_up.count() ), <a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a> ( rv.dataind ( rv_up ) ) { |
---|
188 | <a name="l00277"></a>00277 it_assert_debug ( <a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a>.length() ==<a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216" title="Remember how long val should be.">valsize</a>,<span class="stringliteral">"rv is not fully in rv_up"</span> ); |
---|
189 | <a name="l00278"></a>00278 } |
---|
190 | <a name="l00280"></a><a class="code" href="classdatalink__e2e.html#1b114a76be7aa45bf33697e521fd041d">00280</a> vec <a class="code" href="classdatalink__e2e.html#1b114a76be7aa45bf33697e521fd041d" title="Get val for myself from val of &quot;Up&quot;.">get_val</a> ( <span class="keyword">const</span> vec &val_up ) {it_assert_debug ( <a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8" title="Remember how long val of &quot;Up&quot; should be.">valupsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> ); <span class="keywordflow">return</span> get_vec ( val_up,<a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a> );} |
---|
191 | <a name="l00282"></a><a class="code" href="classdatalink__e2e.html#1ee38edb1778e8900a8d7535c3b3f053">00282</a> <span class="keywordtype">void</span> <a class="code" href="classdatalink__e2e.html#1ee38edb1778e8900a8d7535c3b3f053" title="Fill val of &quot;Up&quot; by my pieces.">fill_val</a> ( vec &val_up, <span class="keyword">const</span> vec &val ) { |
---|
192 | <a name="l00283"></a>00283 it_assert_debug ( <a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216" title="Remember how long val should be.">valsize</a>==val.length(),<span class="stringliteral">"Wrong val"</span> ); |
---|
193 | <a name="l00284"></a>00284 it_assert_debug ( <a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8" title="Remember how long val of &quot;Up&quot; should be.">valupsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> ); |
---|
194 | <a name="l00285"></a>00285 set_subvector ( val_up, <a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a>, val ); |
---|
195 | <a name="l00286"></a>00286 } |
---|
196 | <a name="l00287"></a>00287 }; |
---|
197 | <a name="l00288"></a>00288 |
---|
198 | <a name="l00290"></a><a class="code" href="classdatalink__m2e.html">00290</a> <span class="keyword">class </span><a class="code" href="classdatalink__m2e.html" title="data link between">datalink_m2e</a>: <span class="keyword">public</span> <a class="code" href="classdatalink__e2e.html">datalink_e2e</a> { |
---|
199 | <a name="l00291"></a>00291 <span class="keyword">protected</span>: |
---|
200 | <a name="l00293"></a><a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3">00293</a> <span class="keywordtype">int</span> <a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3" title="Remember how long cond should be.">condsize</a>; |
---|
201 | <a name="l00295"></a><a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a">00295</a> ivec <a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a>; |
---|
202 | <a name="l00297"></a><a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28">00297</a> ivec <a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>; |
---|
203 | <a name="l00298"></a>00298 |
---|
204 | <a name="l00299"></a>00299 <span class="keyword">public</span>: |
---|
205 | <a name="l00301"></a><a class="code" href="classdatalink__m2e.html#ff8f6e8a974e22ced461b232fadbdce0">00301</a> <a class="code" href="classdatalink__m2e.html#ff8f6e8a974e22ced461b232fadbdce0" title="Constructor.">datalink_m2e</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rvc, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up ) : |
---|
206 | <a name="l00302"></a>00302 <a class="code" href="classdatalink__e2e.html">datalink_e2e</a> ( rv,rv_up ), <a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3" title="Remember how long cond should be.">condsize</a> ( rvc.count() ) { |
---|
207 | <a name="l00303"></a>00303 <span class="comment">//establish v2c connection</span> |
---|
208 | <a name="l00304"></a>00304 rvc.<a class="code" href="classRV.html#bb724fa4e2d9ed7bfd0993b5975018a4">dataind</a> ( rv_up, <a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>, <a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a> ); |
---|
209 | <a name="l00305"></a>00305 } |
---|
210 | <a name="l00307"></a><a class="code" href="classdatalink__m2e.html#6db351e5e2158791e97144d212844084">00307</a> vec <a class="code" href="classdatalink__m2e.html#6db351e5e2158791e97144d212844084" title="Construct condition.">get_cond</a> ( <span class="keyword">const</span> vec &val_up ) { |
---|
211 | <a name="l00308"></a>00308 vec tmp ( <a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3" title="Remember how long cond should be.">condsize</a> ); |
---|
212 | <a name="l00309"></a>00309 set_subvector ( tmp,<a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>,val_up ( <a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a> ) ); |
---|
213 | <a name="l00310"></a>00310 <span class="keywordflow">return</span> tmp; |
---|
214 | <a name="l00311"></a>00311 } |
---|
215 | <a name="l00312"></a>00312 <span class="keywordtype">void</span> fill_val_cond ( vec &val_up, <span class="keyword">const</span> vec &val, <span class="keyword">const</span> vec &cond ) { |
---|
216 | <a name="l00313"></a>00313 it_assert_debug ( <a class="code" href="classdatalink__e2e.html#3ce1d13efa0b133a4d202c11bd57d216" title="Remember how long val should be.">valsize</a>==val.length(),<span class="stringliteral">"Wrong val"</span> ); |
---|
217 | <a name="l00314"></a>00314 it_assert_debug ( <a class="code" href="classdatalink__e2e.html#db962bc6201d1ec12dce00d7c1a8b4f8" title="Remember how long val of &quot;Up&quot; should be.">valupsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> ); |
---|
218 | <a name="l00315"></a>00315 set_subvector ( val_up, <a class="code" href="classdatalink__e2e.html#d4c3a76395b46d9d901eb232aacdafae" title="val-to-val link, indeces of the upper val">v2v_up</a>, val ); |
---|
219 | <a name="l00316"></a>00316 set_subvector ( val_up, <a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a>, cond ); |
---|
220 | <a name="l00317"></a>00317 } |
---|
221 | <a name="l00318"></a>00318 }; |
---|
222 | <a name="l00321"></a><a class="code" href="classdatalink__m2m.html">00321</a> <span class="keyword">class </span><a class="code" href="classdatalink__m2m.html">datalink_m2m</a>: <span class="keyword">public</span> <a class="code" href="classdatalink__m2e.html" title="data link between">datalink_m2e</a> { |
---|
223 | <a name="l00322"></a>00322 <span class="keyword">protected</span>: |
---|
224 | <a name="l00324"></a><a class="code" href="classdatalink__m2m.html#d1044bf4846f878290100429600273b2">00324</a> ivec <a class="code" href="classdatalink__m2m.html#d1044bf4846f878290100429600273b2" title="cond-to-cond link, indeces of the upper cond">c2c_up</a>; |
---|
225 | <a name="l00326"></a><a class="code" href="classdatalink__m2m.html#056193d459cd5e5fd323b5cbc6a55bad">00326</a> ivec <a class="code" href="classdatalink__m2m.html#056193d459cd5e5fd323b5cbc6a55bad" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>; |
---|
226 | <a name="l00327"></a>00327 <span class="keyword">public</span>: |
---|
227 | <a name="l00329"></a><a class="code" href="classdatalink__m2m.html#bd64962e779b6dd89529046245a8a98d">00329</a> <a class="code" href="classdatalink__m2m.html#bd64962e779b6dd89529046245a8a98d" title="Constructor.">datalink_m2m</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rvc, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up, <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rvc_up ) : |
---|
228 | <a name="l00330"></a>00330 <a class="code" href="classdatalink__m2e.html" title="data link between">datalink_m2e</a> ( rv, rvc, rv_up) { |
---|
229 | <a name="l00331"></a>00331 <span class="comment">//establish c2c connection</span> |
---|
230 | <a name="l00332"></a>00332 rvc.<a class="code" href="classRV.html#bb724fa4e2d9ed7bfd0993b5975018a4">dataind</a> ( rvc_up, <a class="code" href="classdatalink__m2m.html#056193d459cd5e5fd323b5cbc6a55bad" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>, <a class="code" href="classdatalink__m2m.html#d1044bf4846f878290100429600273b2" title="cond-to-cond link, indeces of the upper cond">c2c_up</a> ); |
---|
231 | <a name="l00333"></a>00333 it_assert_debug(<a class="code" href="classdatalink__m2m.html#056193d459cd5e5fd323b5cbc6a55bad" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>.length()+<a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>.length()==<a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3" title="Remember how long cond should be.">condsize</a>, <span class="stringliteral">"cond is not fully given"</span>); |
---|
232 | <a name="l00334"></a>00334 } |
---|
233 | <a name="l00336"></a><a class="code" href="classdatalink__m2m.html#51dd05760f99282a1da5f51eb88e7e42">00336</a> vec <a class="code" href="classdatalink__m2m.html#51dd05760f99282a1da5f51eb88e7e42" title="Get cond for myself from val and cond of &quot;Up&quot;.">get_cond</a> ( <span class="keyword">const</span> vec &val_up, <span class="keyword">const</span> vec &cond_up ) { |
---|
234 | <a name="l00337"></a>00337 vec tmp ( <a class="code" href="classdatalink__m2e.html#b3927998b71fce3ab030eba1c3c3cec3" title="Remember how long cond should be.">condsize</a> ); |
---|
235 | <a name="l00338"></a>00338 set_subvector ( tmp,<a class="code" href="classdatalink__m2e.html#29604fb622b6a319d2d5d2709fe01f28" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>,val_up ( <a class="code" href="classdatalink__m2e.html#6b46f983c46621b6bfe2b5da3d467d2a" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a> ) ); |
---|
236 | <a name="l00339"></a>00339 set_subvector ( tmp,<a class="code" href="classdatalink__m2m.html#056193d459cd5e5fd323b5cbc6a55bad" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>,cond_up ( <a class="code" href="classdatalink__m2m.html#d1044bf4846f878290100429600273b2" title="cond-to-cond link, indeces of the upper cond">c2c_up</a> ) ); |
---|
237 | <a name="l00340"></a>00340 <span class="keywordflow">return</span> tmp; |
---|
238 | <a name="l00341"></a>00341 } |
---|
239 | <a name="l00343"></a>00343 |
---|
240 | <a name="l00344"></a>00344 }; |
---|
241 | <a name="l00345"></a>00345 |
---|
242 | <a name="l00349"></a><a class="code" href="classmepdf.html">00349</a> <span class="keyword">class </span><a class="code" href="classmepdf.html" title="Unconditional mpdf, allows using epdf in the role of mpdf.">mepdf</a> : <span class="keyword">public</span> <a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> { |
---|
243 | <a name="l00350"></a>00350 <span class="keyword">public</span>: |
---|
244 | <a name="l00352"></a><a class="code" href="classmepdf.html#0f39973ffd94c1e20e2d8aacc4adc6fa">00352</a> <a class="code" href="classmepdf.html#0f39973ffd94c1e20e2d8aacc4adc6fa" title="Default constructor.">mepdf</a> (<span class="keyword">const</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* em ) :<a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> ( em-><a class="code" href="classmpdf.html#71256ffb5fbd08f41d650e606a5bd585" title="access function">_rv</a>(),<a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a>() ) {<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>=<span class="keyword">const_cast<</span><a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>*<span class="keyword">></span>(em);}; |
---|
245 | <a name="l00353"></a><a class="code" href="classmepdf.html#bc87d89dbc3cf5eb6a0ad7d6d572488f">00353</a> <span class="keywordtype">void</span> <a class="code" href="classmepdf.html#bc87d89dbc3cf5eb6a0ad7d6d572488f" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( <span class="keyword">const</span> vec &cond ) {} |
---|
246 | <a name="l00354"></a>00354 }; |
---|
247 | <a name="l00355"></a>00355 |
---|
248 | <a name="l00358"></a><a class="code" href="classcompositepdf.html">00358</a> <span class="keyword">class </span><a class="code" href="classcompositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</a> { |
---|
249 | <a name="l00359"></a>00359 <span class="keyword">protected</span>: |
---|
250 | <a name="l00361"></a><a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f">00361</a> <span class="keywordtype">int</span> <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a>; |
---|
251 | <a name="l00363"></a><a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d">00363</a> Array<mpdf*> <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a>; |
---|
252 | <a name="l00364"></a>00364 <span class="keyword">public</span>: |
---|
253 | <a name="l00365"></a>00365 <a class="code" href="classcompositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</a> ( Array<mpdf*> A0 ) : <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a> ( A0.length() ), <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a> ( A0 ) {}; |
---|
254 | <a name="l00367"></a>00367 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classcompositepdf.html#635d219fb3e32852400d6f98aa4bdc93" title="find common rv, flag">getrv</a> ( <span class="keywordtype">bool</span> checkoverlap=<span class="keyword">false</span> ); |
---|
255 | <a name="l00369"></a>00369 <span class="keywordtype">void</span> <a class="code" href="classcompositepdf.html#5a4479785c1fdb2cc351fbea07b899ac" title="common rvc of all mpdfs is written to rvc">setrvc</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rvc ); |
---|
256 | <a name="l00370"></a>00370 }; |
---|
257 | <a name="l00371"></a>00371 |
---|
258 | <a name="l00379"></a><a class="code" href="classDS.html">00379</a> <span class="keyword">class </span><a class="code" href="classDS.html" title="Abstract class for discrete-time sources of data.">DS</a> { |
---|
259 | <a name="l00380"></a>00380 <span class="keyword">protected</span>: |
---|
260 | <a name="l00382"></a><a class="code" href="classDS.html#1012025a1a3a1fdd24aa48698a3d3454">00382</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classDS.html#1012025a1a3a1fdd24aa48698a3d3454" title="Observed variables, returned by getdata().">Drv</a>; |
---|
261 | <a name="l00384"></a><a class="code" href="classDS.html#0a88ffaeb3c7c2cf30c69bc41a000629">00384</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classDS.html#0a88ffaeb3c7c2cf30c69bc41a000629" title="Action variables, accepted by write().">Urv</a>; <span class="comment">//</span> |
---|
262 | <a name="l00385"></a>00385 <span class="keyword">public</span>: |
---|
263 | <a name="l00387"></a>00387 <span class="keywordtype">void</span> <a class="code" href="classDS.html#db2dacc9e71a36eeb9c5c2ee402eeeb1" title="Returns full vector of observed data.">getdata</a> ( vec &dt ); |
---|
264 | <a name="l00389"></a>00389 <span class="keywordtype">void</span> <a class="code" href="classDS.html#db2dacc9e71a36eeb9c5c2ee402eeeb1" title="Returns full vector of observed data.">getdata</a> ( vec &dt, ivec &indeces ); |
---|
265 | <a name="l00391"></a>00391 <span class="keywordtype">void</span> <a class="code" href="classDS.html#0dfb6788b70c0d48a1f5276ed6c93039" title="Accepts action variable and schedule it for application.">write</a> ( vec &ut ); |
---|
266 | <a name="l00393"></a>00393 <span class="keywordtype">void</span> <a class="code" href="classDS.html#0dfb6788b70c0d48a1f5276ed6c93039" title="Accepts action variable and schedule it for application.">write</a> ( vec &ut, ivec &indeces ); |
---|
267 | <a name="l00399"></a>00399 <span class="keywordtype">void</span> <a class="code" href="classDS.html#9a35ca9c9321dce2bf63bf668f785b75" title="Method that assigns random variables to the datasource. Typically, the datasource...">linkrvs</a> ( <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &drv, <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &urv ); |
---|
268 | <a name="l00400"></a>00400 |
---|
269 | <a name="l00402"></a>00402 <span class="keywordtype">void</span> <a class="code" href="classDS.html#470ba89e8e7c3c4d830803221b4e9be6" title="Moves from to , i.e. perfroms the actions and reads response of the system.">step</a>(); |
---|
270 | <a name="l00403"></a>00403 |
---|
271 | <a name="l00404"></a>00404 }; |
---|
272 | <a name="l00405"></a>00405 |
---|
273 | <a name="l00410"></a><a class="code" href="classBM.html">00410</a> <span class="keyword">class </span><a class="code" href="classBM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a> { |
---|
274 | <a name="l00411"></a>00411 <span class="keyword">protected</span>: |
---|
275 | <a name="l00413"></a><a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88">00413</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a>; |
---|
276 | <a name="l00415"></a><a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979">00415</a> <span class="keywordtype">double</span> <a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979" title="Logarithm of marginalized data likelihood.">ll</a>; |
---|
277 | <a name="l00417"></a><a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129">00417</a> <span class="keywordtype">bool</span> <a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>; |
---|
278 | <a name="l00418"></a>00418 <span class="keyword">public</span>: |
---|
279 | <a name="l00419"></a>00419 |
---|
280 | <a name="l00421"></a><a class="code" href="classBM.html#43d2ad6a6c9f1bcd22fe8906c477e50e">00421</a> <a class="code" href="classBM.html#43d2ad6a6c9f1bcd22fe8906c477e50e" title="Default constructor.">BM</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv0, <span class="keywordtype">double</span> ll0=0,<span class="keywordtype">bool</span> evalll0=<span class="keyword">true</span> ) :<a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a> ( rv0 ), <a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979" title="Logarithm of marginalized data likelihood.">ll</a> ( ll0 ),<a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a> ( evalll0 ) {<span class="comment">//Fixme: test rv</span> |
---|
281 | <a name="l00422"></a>00422 }; |
---|
282 | <a name="l00424"></a><a class="code" href="classBM.html#e8fc5fae042f2b46d4d3827f9587d47d">00424</a> <a class="code" href="classBM.html#43d2ad6a6c9f1bcd22fe8906c477e50e" title="Default constructor.">BM</a> ( <span class="keyword">const</span> <a class="code" href="classBM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a> &B ) : <a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a> ( B.<a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a> ), <a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979" title="Logarithm of marginalized data likelihood.">ll</a> ( B.<a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979" title="Logarithm of marginalized data likelihood.">ll</a> ), <a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a> ( B.<a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a> ) {} |
---|
283 | <a name="l00425"></a>00425 |
---|
284 | <a name="l00429"></a>00429 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classBM.html#a892eff438aab2dd1a9e2efcb7fb5bdf" title="Incremental Bayes rule.">bayes</a> ( <span class="keyword">const</span> vec &dt ) = 0; |
---|
285 | <a name="l00431"></a>00431 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classBM.html#0186270f75189677f390fe088a9947e9" title="Batch Bayes rule (columns of Dt are observations).">bayesB</a> ( <span class="keyword">const</span> mat &Dt ); |
---|
286 | <a name="l00433"></a>00433 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>& <a class="code" href="classBM.html#8554809fec160bee279bcd5fa14ead5c" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>() <span class="keyword">const</span> =0; |
---|
287 | <a name="l00434"></a>00434 |
---|
288 | <a name="l00436"></a>00436 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classBM.html#db500c4836f883c175a3cf3ab5a3ce24" title="Returns a pointer to the epdf representing posterior density on parameters. Use with...">_e</a>() <span class="keyword">const</span> =0; |
---|
289 | <a name="l00437"></a>00437 |
---|
290 | <a name="l00440"></a><a class="code" href="classBM.html#8a8ce6df431689964c41cc6c849cfd06">00440</a> <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classBM.html#8a8ce6df431689964c41cc6c849cfd06">logpred</a> ( <span class="keyword">const</span> vec &dt )<span class="keyword"> const</span>{it_error ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> 0.0;} |
---|
291 | <a name="l00442"></a><a class="code" href="classBM.html#cd0660f2a1a344b56ac39802708ff165">00442</a> vec <a class="code" href="classBM.html#cd0660f2a1a344b56ac39802708ff165" title="Matrix version of logpred.">logpred_m</a> ( <span class="keyword">const</span> mat &dt )<span class="keyword"> const</span>{vec tmp ( dt.cols() );<span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i<dt.cols();i++ ) {tmp ( i ) =<a class="code" href="classBM.html#8a8ce6df431689964c41cc6c849cfd06">logpred</a> ( dt.get_col ( i ) );}<span class="keywordflow">return</span> tmp;} |
---|
292 | <a name="l00443"></a>00443 |
---|
293 | <a name="l00445"></a><a class="code" href="classBM.html#5594d68ee9aa6fc8c1e79019da5c9d56">00445</a> <span class="keyword">virtual</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classBM.html#5594d68ee9aa6fc8c1e79019da5c9d56" title="Constructs a predictive density (marginal density on data).">predictor</a> ( <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &<a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a> )<span class="keyword"> const </span>{it_error ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> NULL;}; |
---|
294 | <a name="l00446"></a>00446 |
---|
295 | <a name="l00448"></a><a class="code" href="classBM.html#ca0f02b3b4144e0895cc14f7e0374bdd">00448</a> <span class="keyword">virtual</span> <a class="code" href="classBM.html#ca0f02b3b4144e0895cc14f7e0374bdd" title="Destructor for future use;.">~BM</a>() {}; |
---|
296 | <a name="l00450"></a><a class="code" href="classBM.html#126bd2595c48e311fc2a7ab72876092a">00450</a> <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classBM.html#126bd2595c48e311fc2a7ab72876092a" title="access function">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classBM.html#af00f0612fabe66241dd507188cdbf88" title="Random variable of the posterior.">rv</a>;} |
---|
297 | <a name="l00452"></a><a class="code" href="classBM.html#87f4a547d2c29180be88175e5eab9c88">00452</a> <span class="keywordtype">double</span> <a class="code" href="classBM.html#87f4a547d2c29180be88175e5eab9c88" title="access function">_ll</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classBM.html#5623fef6572a08c2b53b8c87b82dc979" title="Logarithm of marginalized data likelihood.">ll</a>;} |
---|
298 | <a name="l00454"></a><a class="code" href="classBM.html#1ffa9f23669aabecc3760c06c6987522">00454</a> <span class="keywordtype">void</span> <a class="code" href="classBM.html#1ffa9f23669aabecc3760c06c6987522" title="access function">set_evalll</a> ( <span class="keywordtype">bool</span> evl0 ) {<a class="code" href="classBM.html#bf6fb59b30141074f8ee1e2f43d03129" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>=evl0;} |
---|
299 | <a name="l00455"></a>00455 |
---|
300 | <a name="l00458"></a><a class="code" href="classBM.html#eb58c81d6a7b75b05fc6f276eed78887">00458</a> <span class="keyword">virtual</span> <a class="code" href="classBM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a>* <a class="code" href="classBM.html#eb58c81d6a7b75b05fc6f276eed78887">_copy_</a> ( <span class="keywordtype">bool</span> changerv=<span class="keyword">false</span> ) {it_error ( <span class="stringliteral">"function _copy_ not implemented for this BM"</span> ); <span class="keywordflow">return</span> NULL;}; |
---|
301 | <a name="l00459"></a>00459 }; |
---|
302 | <a name="l00460"></a>00460 |
---|
303 | <a name="l00470"></a><a class="code" href="classBMcond.html">00470</a> <span class="keyword">class </span><a class="code" href="classBMcond.html" title="Conditional Bayesian Filter.">BMcond</a> { |
---|
304 | <a name="l00471"></a>00471 <span class="keyword">protected</span>: |
---|
305 | <a name="l00473"></a><a class="code" href="classBMcond.html#9ba793c8ec453f04d372d17195ed8dec">00473</a> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classBMcond.html#9ba793c8ec453f04d372d17195ed8dec" title="Identificator of the conditioning variable.">rvc</a>; |
---|
306 | <a name="l00474"></a>00474 <span class="keyword">public</span>: |
---|
307 | <a name="l00476"></a>00476 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classBMcond.html#3543cfe274c35f9033c24d8036a96c4d" title="Substitute val for rvc.">condition</a> ( <span class="keyword">const</span> vec &val ) =0; |
---|
308 | <a name="l00478"></a><a class="code" href="classBMcond.html#9a38dcdabef902fd91da3c9fb6ac27f7">00478</a> <a class="code" href="classBMcond.html#9a38dcdabef902fd91da3c9fb6ac27f7" title="Default constructor.">BMcond</a> ( <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> &rv0 ) :<a class="code" href="classBMcond.html#9ba793c8ec453f04d372d17195ed8dec" title="Identificator of the conditioning variable.">rvc</a> ( rv0 ) {}; |
---|
309 | <a name="l00480"></a><a class="code" href="classBMcond.html#d40c72250cf93640f982d06de20815b2">00480</a> <span class="keyword">virtual</span> <a class="code" href="classBMcond.html#d40c72250cf93640f982d06de20815b2" title="Destructor for future use.">~BMcond</a>() {}; |
---|
310 | <a name="l00482"></a><a class="code" href="classBMcond.html#3fa60348b2da6b4208bb95b8d146900a">00482</a> <span class="keyword">const</span> <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classBMcond.html#3fa60348b2da6b4208bb95b8d146900a" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classBMcond.html#9ba793c8ec453f04d372d17195ed8dec" title="Identificator of the conditioning variable.">rvc</a>;} |
---|
311 | <a name="l00483"></a>00483 }; |
---|
312 | <a name="l00484"></a>00484 |
---|
313 | <a name="l00486"></a>00486 <span class="preprocessor">#endif // BM_H</span> |
---|
314 | </pre></div></div> |
---|
315 | <hr size="1"><address style="text-align: right;"><small>Generated on Thu Jan 15 10:50:25 2009 for mixpp by |
---|
316 | <a href="http://www.doxygen.org/index.html"> |
---|
317 | <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address> |
---|
318 | </body> |
---|
319 | </html> |
---|