41 | | <a name="l00050"></a><a class="code" href="classemix.html#469e910479b3402589104ef3bb1b5741">00050</a> <span class="keywordtype">double</span> <a class="code" href="classemix.html#469e910479b3402589104ef3bb1b5741" title="Compute log-probability of argument val.">evalpdflog</a> ( <span class="keyword">const</span> vec &val )<span class="keyword"> const </span>{<span class="keywordtype">int</span> i; <span class="keywordtype">double</span> sum=0.0; <span class="keywordflow">for</span> ( i=0;i<<a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a>.length();i++ ) {sum+=<a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a> ( i ) *<a class="code" href="classemix.html#7155d827417b53fb6aa5cff1459e6eb4" title="Component (epdfs).">Coms</a> ( i )->evalpdflog ( val );} <span class="keywordflow">return</span> log ( sum );}; |
42 | | <a name="l00051"></a>00051 |
43 | | <a name="l00052"></a>00052 <span class="comment">//Access methods</span> |
44 | | <a name="l00054"></a><a class="code" href="classemix.html#35baf0f62307455cebfbf462144b4dd8">00054</a> <span class="comment"></span> vec& <a class="code" href="classemix.html#35baf0f62307455cebfbf462144b4dd8" title="returns a pointer to the internal mean value. Use with Care!">_w</a>() {<span class="keywordflow">return</span> <a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a>;} |
45 | | <a name="l00055"></a>00055 }; |
| 41 | <a name="l00050"></a><a class="code" href="classemix.html#469e910479b3402589104ef3bb1b5741">00050</a> <span class="keywordtype">double</span> <a class="code" href="classemix.html#469e910479b3402589104ef3bb1b5741" title="Compute log-probability of argument val.">evalpdflog</a> ( <span class="keyword">const</span> vec &val )<span class="keyword"> const </span>{ |
| 42 | <a name="l00051"></a>00051 <span class="keywordtype">int</span> i; |
| 43 | <a name="l00052"></a>00052 <span class="keywordtype">double</span> sum=0.0; |
| 44 | <a name="l00053"></a>00053 <span class="keywordflow">for</span> ( i=0;i<<a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a>.length();i++ ) {sum+=<a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a> ( i ) *<a class="code" href="classemix.html#7155d827417b53fb6aa5cff1459e6eb4" title="Component (epdfs).">Coms</a> ( i )->evalpdflog ( val );} |
| 45 | <a name="l00054"></a>00054 <span class="keywordflow">return</span> log ( sum ); |
| 46 | <a name="l00055"></a>00055 }; |
47 | | <a name="l00057"></a>00057 |
48 | | <a name="l00058"></a>00058 <span class="preprocessor">#endif //MX_H</span> |
| 48 | <a name="l00057"></a>00057 <span class="comment">//Access methods</span> |
| 49 | <a name="l00059"></a><a class="code" href="classemix.html#35baf0f62307455cebfbf462144b4dd8">00059</a> <span class="comment"></span> vec& <a class="code" href="classemix.html#35baf0f62307455cebfbf462144b4dd8" title="returns a pointer to the internal mean value. Use with Care!">_w</a>() {<span class="keywordflow">return</span> <a class="code" href="classemix.html#cdc8f4e2a5e43595ac5450d8a7f9d1ab" title="weights of the components">w</a>;} |
| 50 | <a name="l00060"></a>00060 }; |
| 51 | <a name="l00061"></a>00061 |
| 52 | <a name="l00066"></a><a class="code" href="classeprod.html">00066</a> <span class="keyword">class </span><a class="code" href="classeprod.html" title="Chain rule decomposition of epdf.">eprod</a>: <span class="keyword">public</span> <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> { |
| 53 | <a name="l00067"></a>00067 <span class="keyword">protected</span>: |
| 54 | <a name="l00068"></a>00068 Array<epdf*> epdfs; |
| 55 | <a name="l00069"></a>00069 Array<mpdf*> mpdfs; |
| 56 | <a name="l00070"></a>00070 <span class="keyword">public</span>: |
| 57 | <a name="l00071"></a>00071 |
| 58 | <a name="l00072"></a>00072 |
| 59 | <a name="l00073"></a>00073 }; |
| 60 | <a name="l00074"></a>00074 |
| 61 | <a name="l00078"></a><a class="code" href="classmmix.html">00078</a> <span class="keyword">class </span><a class="code" href="classmmix.html" title="Mixture of mpdfs with constant weights.">mmix</a> : <span class="keyword">public</span> <a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> { |
| 62 | <a name="l00079"></a>00079 <span class="keyword">protected</span>: |
| 63 | <a name="l00081"></a><a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f">00081</a> Array<mpdf*> <a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a>; |
| 64 | <a name="l00083"></a><a class="code" href="classmmix.html#507a5e6df53e6e9a6a330f1271cf0914">00083</a> <a class="code" href="classemix.html" title="Mixture of epdfs.">emix</a> <a class="code" href="classmmix.html#507a5e6df53e6e9a6a330f1271cf0914" title="Internal epdf.">Epdf</a>; |
| 65 | <a name="l00084"></a>00084 <span class="keyword">public</span>: |
| 66 | <a name="l00086"></a><a class="code" href="classmmix.html#6769a962f1f4d0ecc5b903b5f7d26ad3">00086</a> <a class="code" href="classmmix.html#6769a962f1f4d0ecc5b903b5f7d26ad3" title="Default constructor.">mmix</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>, <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> ) : <a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> ( rv, rvc ), <a class="code" href="classmmix.html#507a5e6df53e6e9a6a330f1271cf0914" title="Internal epdf.">Epdf</a> ( rv ) {<a class="code" href="classmpdf.html#7aa894208a32f3487827df6d5054424c" title="pointer to internal epdf">ep</a>=&<a class="code" href="classmmix.html#507a5e6df53e6e9a6a330f1271cf0914" title="Internal epdf.">Epdf</a>;}; |
| 67 | <a name="l00088"></a><a class="code" href="classmmix.html#e7681490a5c5b2d0a43f4e5842e10a9d">00088</a> <span class="keywordtype">void</span> <a class="code" href="classmmix.html#e7681490a5c5b2d0a43f4e5842e10a9d" title="Set weights w and components R.">set_parameters</a> ( <span class="keyword">const</span> vec &w, <span class="keyword">const</span> Array<mpdf*> &<a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a> ) { |
| 68 | <a name="l00089"></a>00089 Array<epdf*> Eps ( Coms.length()); |
| 69 | <a name="l00090"></a>00090 |
| 70 | <a name="l00091"></a>00091 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i<Coms.length();i++ ) {<a class="code" href="classmpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>* Ci=<a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a>(i); |
| 71 | <a name="l00092"></a>00092 Eps ( i ) =& ( <a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a> ( i )->_epdf() );} |
| 72 | <a name="l00093"></a>00093 <a class="code" href="classmmix.html#507a5e6df53e6e9a6a330f1271cf0914" title="Internal epdf.">Epdf</a>.<a class="code" href="classemix.html#6ce82b24d7fe540ddbf44b472a005555" title="Set weights w and components R.">set_parameters</a> ( w,Eps ); |
| 73 | <a name="l00094"></a>00094 }; |
| 74 | <a name="l00095"></a>00095 |
| 75 | <a name="l00096"></a><a class="code" href="classmmix.html#1a66629ec5c2b6f4740606cf737f0193">00096</a> <span class="keywordtype">void</span> <a class="code" href="classmmix.html#1a66629ec5c2b6f4740606cf737f0193" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( <span class="keyword">const</span> vec &cond ) { |
| 76 | <a name="l00097"></a>00097 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i<<a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a>.length();i++ ) {<a class="code" href="classmmix.html#460632c7a0afd41b7840902181a7d85f" title="Component (epdfs).">Coms</a> ( i )->condition ( cond );} |
| 77 | <a name="l00098"></a>00098 }; |
| 78 | <a name="l00099"></a>00099 }; |
| 79 | <a name="l00100"></a>00100 <span class="preprocessor">#endif //MX_H</span> |