Show
Ignore:
Timestamp:
01/28/09 10:01:09 (15 years ago)
Author:
smidl
Message:

doc

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • doc/html/merger_8h-source.html

    r234 r255  
    2323<a name="l00017"></a>00017 <span class="preprocessor">#include "<a class="code" href="mixef_8h.html" title="Bayesian Filtering for mixtures of exponential family (EF) members.">mixef.h</a>"</span> 
    2424<a name="l00018"></a>00018  
    25 <a name="l00019"></a>00019 <span class="keyword">using namespace </span>itpp; 
     25<a name="l00019"></a>00019 <span class="keyword">namespace </span>bdm{ 
    2626<a name="l00020"></a>00020 <span class="keyword">using</span> std::string; 
    2727<a name="l00021"></a>00021  
    28 <a name="l00028"></a><a class="code" href="classmerger.html">00028</a> <span class="keyword">class </span><a class="code" href="classmerger.html" title="Function for general combination of pdfs.">merger</a> : <span class="keyword">public</span> <a class="code" href="classcompositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</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> { 
     28<a name="l00028"></a><a class="code" href="classbdm_1_1merger.html">00028</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1merger.html" title="Function for general combination of pdfs.">merger</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1compositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</a>, <span class="keyword">public</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> { 
    2929<a name="l00029"></a>00029 <span class="keyword">protected</span>: 
    30 <a name="l00031"></a><a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f">00031</a>         <a class="code" href="classMixEF.html" title="Mixture of Exponential Family Densities.">MixEF</a> <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>; 
    31 <a name="l00033"></a><a class="code" href="classmerger.html#01e222f5df660e1b25c712e82da8df8f">00033</a>         Array&lt;datalink_m2e*&gt; <a class="code" href="classmerger.html#01e222f5df660e1b25c712e82da8df8f" title="Data link for each mpdf in mpdfs.">dls</a>; 
    32 <a name="l00035"></a><a class="code" href="classmerger.html#e741e6096a7acec3dfa477efa1d8160e">00035</a>         Array&lt;RV&gt; <a class="code" href="classmerger.html#e741e6096a7acec3dfa477efa1d8160e" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a>; 
    33 <a name="l00037"></a><a class="code" href="classmerger.html#d0ff2953f904f12a355f72f51813ca7a">00037</a>         Array&lt;datalink_m2e*&gt; <a class="code" href="classmerger.html#d0ff2953f904f12a355f72f51813ca7a" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a>; 
     30<a name="l00031"></a><a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81">00031</a>         <a class="code" href="classbdm_1_1MixEF.html" title="Mixture of Exponential Family Densities.">MixEF</a> <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>; 
     31<a name="l00033"></a><a class="code" href="classbdm_1_1merger.html#ea3c9f4739f857877fd5285b3973289c">00033</a>         Array&lt;datalink_m2e*&gt; <a class="code" href="classbdm_1_1merger.html#ea3c9f4739f857877fd5285b3973289c" title="Data link for each mpdf in mpdfs.">dls</a>; 
     32<a name="l00035"></a><a class="code" href="classbdm_1_1merger.html#135b938f841cb88ad00a8bbded818703">00035</a>         Array&lt;RV&gt; <a class="code" href="classbdm_1_1merger.html#135b938f841cb88ad00a8bbded818703" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a>; 
     33<a name="l00037"></a><a class="code" href="classbdm_1_1merger.html#6c48bf36ec1f3f7b52e3159cf428a193">00037</a>         Array&lt;datalink_m2e*&gt; <a class="code" href="classbdm_1_1merger.html#6c48bf36ec1f3f7b52e3159cf428a193" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a>; 
    3434<a name="l00038"></a>00038  
    35 <a name="l00040"></a><a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69">00040</a>         <span class="keywordtype">int</span> <a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69" title="Number of samples used in approximation.">Ns</a>; 
    36 <a name="l00042"></a><a class="code" href="classmerger.html#6acf33e1c64aaaee20689725e8110f65">00042</a>         <span class="keywordtype">int</span> <a class="code" href="classmerger.html#6acf33e1c64aaaee20689725e8110f65" title="Number of components in a mixture.">Nc</a>; 
    37 <a name="l00044"></a><a class="code" href="classmerger.html#e68a72e2740007ae68545e378c7701ad">00044</a>         <span class="keywordtype">double</span> <a class="code" href="classmerger.html#e68a72e2740007ae68545e378c7701ad" title="Prior on the log-normal merging model.">beta</a>; 
    38 <a name="l00046"></a><a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801">00046</a>         <a class="code" href="classeEmp.html" title="Weighted empirical density.">eEmp</a> <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>; 
     35<a name="l00040"></a><a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109">00040</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109" title="Number of samples used in approximation.">Ns</a>; 
     36<a name="l00042"></a><a class="code" href="classbdm_1_1merger.html#214ca7c566e8f75d45bb243d7c74580d">00042</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1merger.html#214ca7c566e8f75d45bb243d7c74580d" title="Number of components in a mixture.">Nc</a>; 
     37<a name="l00044"></a><a class="code" href="classbdm_1_1merger.html#1a85d299f202a70afe0cdc1045c26a5d">00044</a>         <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1merger.html#1a85d299f202a70afe0cdc1045c26a5d" title="Prior on the log-normal merging model.">beta</a>; 
     38<a name="l00046"></a><a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94">00046</a>         <a class="code" href="classbdm_1_1eEmp.html" title="Weighted empirical density.">eEmp</a> <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>; 
    3939<a name="l00047"></a>00047  
    4040<a name="l00048"></a>00048 <span class="keyword">public</span>: 
    41 <a name="l00050"></a><a class="code" href="classmerger.html#3f9a955329eab65026e8322e9027d469">00050</a>         <a class="code" href="classmerger.html#3f9a955329eab65026e8322e9027d469" title="Default constructor.">merger</a> ( <span class="keyword">const</span> Array&lt;mpdf*&gt; &amp;S ) : 
    42 <a name="l00051"></a>00051                         <a class="code" href="classcompositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</a> ( S ), <a class="code" href="classepdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> ( <a class="code" href="classcompositepdf.html#635d219fb3e32852400d6f98aa4bdc93" title="find common rv, flag">getrv</a> ( false ) ), 
    43 <a name="l00052"></a>00052                         <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a> ( Array&lt;<a class="code" href="classBMEF.html" title="Estimator for Exponential family.">BMEF</a>*&gt; ( 0 ),vec ( 0 ) ), <a class="code" href="classmerger.html#01e222f5df660e1b25c712e82da8df8f" title="Data link for each mpdf in mpdfs.">dls</a> ( <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classmerger.html#e741e6096a7acec3dfa477efa1d8160e" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classmerger.html#d0ff2953f904f12a355f72f51813ca7a" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>(<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>,0) { 
    44 <a name="l00053"></a>00053                 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> ztmp; 
     41<a name="l00050"></a><a class="code" href="classbdm_1_1merger.html#30ae0fa8a88737185f598bf63fec0569">00050</a>         <a class="code" href="classbdm_1_1merger.html#30ae0fa8a88737185f598bf63fec0569" title="Default constructor.">merger</a> ( <span class="keyword">const</span> Array&lt;mpdf*&gt; &amp;S ) : 
     42<a name="l00051"></a>00051                         <a class="code" href="classbdm_1_1compositepdf.html" title="Abstract composition of pdfs, a base for specific classes this abstract class is...">compositepdf</a> ( S ), <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> ( <a class="code" href="classbdm_1_1compositepdf.html#c73d39acc4378eee6a63155c3517d3c9" title="find common rv, flag">getrv</a> ( false ) ), 
     43<a name="l00052"></a>00052                         <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a> ( Array&lt;<a class="code" href="classbdm_1_1BMEF.html" title="Estimator for Exponential family.">BMEF</a>*&gt; ( 0 ),vec ( 0 ) ), <a class="code" href="classbdm_1_1merger.html#ea3c9f4739f857877fd5285b3973289c" title="Data link for each mpdf in mpdfs.">dls</a> ( <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classbdm_1_1merger.html#135b938f841cb88ad00a8bbded818703" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classbdm_1_1merger.html#6c48bf36ec1f3f7b52e3159cf428a193" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ), <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>(<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>,0) { 
     44<a name="l00053"></a>00053                 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> ztmp; 
    4545<a name="l00054"></a>00054                 <span class="comment">// Extend rv by rvc!</span> 
    46 <a name="l00055"></a>00055                 <a class="code" href="classRV.html" title="Class representing variables, most often random variables.">RV</a> rvc; <a class="code" href="classcompositepdf.html#5a4479785c1fdb2cc351fbea07b899ac" title="common rvc of all mpdfs is written to rvc">setrvc</a> ( <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>,rvc ); 
    47 <a name="l00056"></a>00056                 <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#935790526b630dec4ffefa9ec0c2b6fb" title="Add (concat) another variable to the current one,.">add</a> ( rvc ); 
    48 <a name="l00057"></a>00057                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;<a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a>;i++ ) { 
     46<a name="l00055"></a>00055                 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> rvc; <a class="code" href="classbdm_1_1compositepdf.html#c71fc2c51f49d797e61f479f543c75ce" title="common rvc of all mpdfs is written to rvc">setrvc</a> ( <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>,rvc ); 
     47<a name="l00056"></a>00056                 <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#87841b5ee43997b79789c0c22047e224" title="Add (concat) another variable to the current one,.">add</a> ( rvc ); 
     48<a name="l00057"></a>00057                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;<a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a>;i++ ) { 
    4949<a name="l00058"></a>00058                         <span class="comment">//Establich connection between mpdfs and merger</span> 
    50 <a name="l00059"></a>00059                         <a class="code" href="classmerger.html#01e222f5df660e1b25c712e82da8df8f" title="Data link for each mpdf in mpdfs.">dls</a> ( i ) = <span class="keyword">new</span> <a class="code" href="classdatalink__m2e.html" title="data link between">datalink_m2e</a> ( <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a> ( i )-&gt;<a class="code" href="classepdf.html#ca0d32aabb4cbba347e0c37fe8607562" title="access function, possibly dangerous!">_rv</a>(), <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rvc(), <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> ); 
     50<a name="l00059"></a>00059                         <a class="code" href="classbdm_1_1merger.html#ea3c9f4739f857877fd5285b3973289c" title="Data link for each mpdf in mpdfs.">dls</a> ( i ) = <span class="keyword">new</span> <a class="code" href="classbdm_1_1datalink__m2e.html" title="data link between">datalink_m2e</a> ( <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( i )-&gt;<a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21" title="access function, possibly dangerous!">_rv</a>(), <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rvc(), <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a> ); 
    5151<a name="l00060"></a>00060                         <span class="comment">// find out what is missing in each mpdf</span> 
    52 <a name="l00061"></a>00061                         ztmp= <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rv(); 
    53 <a name="l00062"></a>00062                         ztmp.<a class="code" href="classRV.html#935790526b630dec4ffefa9ec0c2b6fb" title="Add (concat) another variable to the current one,.">add</a> ( <a class="code" href="classcompositepdf.html#637877061b917a16a0421bb5c4cd480d" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rvc() ); 
    54 <a name="l00063"></a>00063                         <a class="code" href="classmerger.html#e741e6096a7acec3dfa477efa1d8160e" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( i ) =<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#b9d175c327c21488b1e2fb756a84e149" title="Subtract another variable from the current one.">subt</a> ( ztmp ); 
    55 <a name="l00064"></a>00064                         <a class="code" href="classmerger.html#d0ff2953f904f12a355f72f51813ca7a" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( i ) = <span class="keyword">new</span> <a class="code" href="classdatalink__m2e.html" title="data link between">datalink_m2e</a> ( <a class="code" href="classmerger.html#e741e6096a7acec3dfa477efa1d8160e" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( i ), ztmp, <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a> ) ; 
     52<a name="l00061"></a>00061                         ztmp= <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rv(); 
     53<a name="l00062"></a>00062                         ztmp.<a class="code" href="classbdm_1_1RV.html#87841b5ee43997b79789c0c22047e224" title="Add (concat) another variable to the current one,.">add</a> ( <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( i )-&gt;_rvc() ); 
     54<a name="l00063"></a>00063                         <a class="code" href="classbdm_1_1merger.html#135b938f841cb88ad00a8bbded818703" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( i ) =<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#aec44dabdf0a6d90fbae95e1356eda39" title="Subtract another variable from the current one.">subt</a> ( ztmp ); 
     55<a name="l00064"></a>00064                         <a class="code" href="classbdm_1_1merger.html#6c48bf36ec1f3f7b52e3159cf428a193" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( i ) = <span class="keyword">new</span> <a class="code" href="classbdm_1_1datalink__m2e.html" title="data link between">datalink_m2e</a> ( <a class="code" href="classbdm_1_1merger.html#135b938f841cb88ad00a8bbded818703" title="Array of rvs that are not modelled by mpdfs at all (aux).">rvzs</a> ( i ), ztmp, <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a> ) ; 
    5656<a name="l00065"></a>00065                 }; 
    5757<a name="l00066"></a>00066                 <span class="comment">//Set Default values of parameters</span> 
    58 <a name="l00067"></a>00067                 <a class="code" href="classmerger.html#e68a72e2740007ae68545e378c7701ad" title="Prior on the log-normal merging model.">beta</a>=2.0; 
    59 <a name="l00068"></a>00068                 <a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69" title="Number of samples used in approximation.">Ns</a>=100; 
    60 <a name="l00069"></a>00069                 <a class="code" href="classmerger.html#6acf33e1c64aaaee20689725e8110f65" title="Number of components in a mixture.">Nc</a>=10; 
    61 <a name="l00070"></a>00070                 <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classMixEF.html#6576024e16523da5cbaaf233512c53dc" title="Set which method is to be used.">set_method</a> ( EM ); 
     58<a name="l00067"></a>00067                 <a class="code" href="classbdm_1_1merger.html#1a85d299f202a70afe0cdc1045c26a5d" title="Prior on the log-normal merging model.">beta</a>=2.0; 
     59<a name="l00068"></a>00068                 <a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109" title="Number of samples used in approximation.">Ns</a>=100; 
     60<a name="l00069"></a>00069                 <a class="code" href="classbdm_1_1merger.html#214ca7c566e8f75d45bb243d7c74580d" title="Number of components in a mixture.">Nc</a>=10; 
     61<a name="l00070"></a>00070                 <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classbdm_1_1MixEF.html#664529d52cc667383b39eeb440ccd577" title="Set which method is to be used.">set_method</a> ( EM ); 
    6262<a name="l00071"></a>00071         } 
    63 <a name="l00073"></a><a class="code" href="classmerger.html#79420e251308df4502c834171d518ad9">00073</a>         <span class="keywordtype">void</span> <a class="code" href="classmerger.html#79420e251308df4502c834171d518ad9" title="Set internal parameters used in approximation.">set_parameters</a> ( <span class="keywordtype">double</span> beta0, <span class="keywordtype">int</span> Ns0, <span class="keywordtype">int</span> Nc0 ) {<a class="code" href="classmerger.html#e68a72e2740007ae68545e378c7701ad" title="Prior on the log-normal merging model.">beta</a>=beta0;<a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69" title="Number of samples used in approximation.">Ns</a>=Ns0;<a class="code" href="classmerger.html#6acf33e1c64aaaee20689725e8110f65" title="Number of components in a mixture.">Nc</a>=Nc0;<a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#a4215f6a5a04d07b43f7ebaa942e15f1" title="Set sample.">set_n</a>(Ns0,<span class="keyword">false</span>);} 
    64 <a name="l00075"></a><a class="code" href="classmerger.html#02d2a445e44a87f37ee3a2b9b40cffd1">00075</a>         <span class="keywordtype">void</span> <a class="code" href="classmerger.html#02d2a445e44a87f37ee3a2b9b40cffd1" title="Initialize the proposal density. This function must be called before merge()!">init</a>() {  
    65 <a name="l00076"></a>00076                 Array&lt;vec&gt; Smps ( <a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a> ); 
     63<a name="l00073"></a><a class="code" href="classbdm_1_1merger.html#bc867c8e2b4709041bf8802ae9c78e37">00073</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1merger.html#bc867c8e2b4709041bf8802ae9c78e37" title="Set internal parameters used in approximation.">set_parameters</a> ( <span class="keywordtype">double</span> beta0, <span class="keywordtype">int</span> Ns0, <span class="keywordtype">int</span> Nc0 ) {<a class="code" href="classbdm_1_1merger.html#1a85d299f202a70afe0cdc1045c26a5d" title="Prior on the log-normal merging model.">beta</a>=beta0;<a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109" title="Number of samples used in approximation.">Ns</a>=Ns0;<a class="code" href="classbdm_1_1merger.html#214ca7c566e8f75d45bb243d7c74580d" title="Number of components in a mixture.">Nc</a>=Nc0;<a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#dccd02eaa4c45e858a6351723686ac85" title="Set sample.">set_n</a>(Ns0,<span class="keyword">false</span>);} 
     64<a name="l00075"></a><a class="code" href="classbdm_1_1merger.html#8ae79a651e93d6e2d611556001af8aa8">00075</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1merger.html#8ae79a651e93d6e2d611556001af8aa8" title="Initialize the proposal density. This function must be called before merge()!">init</a>() {  
     65<a name="l00076"></a>00076                 Array&lt;vec&gt; Smps ( <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ); 
    6666<a name="l00077"></a>00077                 <span class="comment">//Gibbs sampling</span> 
    67 <a name="l00078"></a>00078                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;<a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a>;i++ ) {Smps ( i ) =zeros ( 0 );} 
     67<a name="l00078"></a>00078                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;<a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a>;i++ ) {Smps ( i ) =zeros ( 0 );} 
    6868<a name="l00079"></a>00079         } 
    69 <a name="l00081"></a>00081         <span class="keywordtype">void</span> <a class="code" href="classmerger.html#bc23accc77ec766bce7202934e2c180c" title="Create a mixture density, make sure to call init() before the first call.">merge</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>* g0 ); 
    70 <a name="l00083"></a><a class="code" href="classmerger.html#bc23accc77ec766bce7202934e2c180c">00083</a>         <span class="keywordtype">void</span> <a class="code" href="classmerger.html#bc23accc77ec766bce7202934e2c180c" title="Create a mixture density, make sure to call init() before the first call.">merge</a> () {<a class="code" href="classmerger.html#bc23accc77ec766bce7202934e2c180c" title="Create a mixture density, make sure to call init() before the first call.">merge</a> ( &amp; ( <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classMixEF.html#efb3e20c2151d91c4fc080b7722a2069" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>() ) );}; 
     69<a name="l00081"></a>00081         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1merger.html#b72f4662c9ac760a8b17ebe5d697f3c2" title="Create a mixture density, make sure to call init() before the first call.">merge</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* g0 ); 
     70<a name="l00083"></a><a class="code" href="classbdm_1_1merger.html#b72f4662c9ac760a8b17ebe5d697f3c2">00083</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1merger.html#b72f4662c9ac760a8b17ebe5d697f3c2" title="Create a mixture density, make sure to call init() before the first call.">merge</a> () {<a class="code" href="classbdm_1_1merger.html#b72f4662c9ac760a8b17ebe5d697f3c2" title="Create a mixture density, make sure to call init() before the first call.">merge</a> ( &amp; ( <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classbdm_1_1MixEF.html#33d0b3da1d10bf149d41ee74f6284a19" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>() ) );}; 
    7171<a name="l00084"></a>00084  
    72 <a name="l00086"></a>00086         vec <a class="code" href="classmerger.html#217c65f335c603033fa09a2ecc64bc05" title="Merge log-likelihood values.">lognorm_merge</a> ( mat &amp;lW ); 
    73 <a name="l00089"></a><a class="code" href="classmerger.html#379198c3d2063bfa63f5d1245a2511ba">00089</a>         vec <a class="code" href="classmerger.html#379198c3d2063bfa63f5d1245a2511ba">sample</a> ( )<span class="keyword"> const </span>{ <span class="keywordflow">return</span> <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classMixEF.html#efb3e20c2151d91c4fc080b7722a2069" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>().<a class="code" href="classepdf.html#8019654e494bf5e458f6fb947e11b262" title="Returns a sample,  from density .">sample</a>();} 
    74 <a name="l00090"></a><a class="code" href="classmerger.html#632cd7e0bcd149c0dc85042063364f6b">00090</a>         <span class="keywordtype">double</span> <a class="code" href="classmerger.html#632cd7e0bcd149c0dc85042063364f6b" title="Compute log-probability of argument val.">evallog</a> ( <span class="keyword">const</span> vec &amp;dt )<span class="keyword"> const </span>{ 
     72<a name="l00086"></a>00086         vec <a class="code" href="classbdm_1_1merger.html#b90f1b8a21a693176dd953a3b983aa0b" title="Merge log-likelihood values.">lognorm_merge</a> ( mat &amp;lW ); 
     73<a name="l00089"></a><a class="code" href="classbdm_1_1merger.html#0ef73e76eeee44c8e3896e3a271ad974">00089</a>         vec <a class="code" href="classbdm_1_1merger.html#0ef73e76eeee44c8e3896e3a271ad974">sample</a> ( )<span class="keyword"> const </span>{ <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classbdm_1_1MixEF.html#33d0b3da1d10bf149d41ee74f6284a19" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>().<a class="code" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10" title="Returns a sample,  from density .">sample</a>();} 
     74<a name="l00090"></a><a class="code" href="classbdm_1_1merger.html#c47543661897638cab9f506c76a1f92f">00090</a>         <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1merger.html#c47543661897638cab9f506c76a1f92f" title="Compute log-probability of argument val.">evallog</a> ( <span class="keyword">const</span> vec &amp;dt )<span class="keyword"> const </span>{ 
    7575<a name="l00091"></a>00091                 vec dtf=ones ( dt.length() +1 ); 
    7676<a name="l00092"></a>00092                 dtf.set_subvector ( 0,dt ); 
    77 <a name="l00093"></a>00093                 <span class="keywordflow">return</span> <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classMixEF.html#424ca64f36d4e41de7a7e7ae921d35ea">logpred</a> ( dtf ); 
     77<a name="l00093"></a>00093                 <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>.<a class="code" href="classbdm_1_1MixEF.html#da724da464a75e07521941e430929efa">logpred</a> ( dtf ); 
    7878<a name="l00094"></a>00094         } 
    79 <a name="l00095"></a><a class="code" href="classmerger.html#3a46999cd77047d5d4ac6c1cb735a83e">00095</a>         vec <a class="code" href="classmerger.html#3a46999cd77047d5d4ac6c1cb735a83e" title="return expected value">mean</a>()<span class="keyword"> const </span>{ 
    80 <a name="l00096"></a>00096                 <span class="keyword">const</span> Vec&lt;double&gt; &amp;w = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b2bfb73b72486a5c89f2ab850c7a9b" title="Potentially dangerous, use with care.">_w</a>(); 
    81 <a name="l00097"></a>00097                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b747eca73b16f30370827ba4cc3575" title="access function">_samples</a>(); 
    82 <a name="l00098"></a>00098                 vec tmp=zeros ( <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>() );  
    83 <a name="l00099"></a>00099                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69" title="Number of samples used in approximation.">Ns</a>; i++ ) { 
     79<a name="l00095"></a><a class="code" href="classbdm_1_1merger.html#82c5ad7677315fde8e9fc9381c17ece5">00095</a>         vec <a class="code" href="classbdm_1_1merger.html#82c5ad7677315fde8e9fc9381c17ece5" title="return expected value">mean</a>()<span class="keyword"> const </span>{ 
     80<a name="l00096"></a>00096                 <span class="keyword">const</span> Vec&lt;double&gt; &amp;w = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#d7f83cc0415cd44ae7cc8b4bdad93aef" title="Potentially dangerous, use with care.">_w</a>(); 
     81<a name="l00097"></a>00097                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#c24966b0aaeb767bc8a6b4fd60931be2" title="access function">_samples</a>(); 
     82<a name="l00098"></a>00098                 vec tmp=zeros ( <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>() );  
     83<a name="l00099"></a>00099                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109" title="Number of samples used in approximation.">Ns</a>; i++ ) { 
    8484<a name="l00100"></a>00100                         tmp+=w ( i ) *S ( i ); 
    8585<a name="l00101"></a>00101                 } 
    … …  
    8787<a name="l00103"></a>00103         } 
    8888<a name="l00104"></a>00104         mat covariance()<span class="keyword"> const </span>{ 
    89 <a name="l00105"></a>00105                 <span class="keyword">const</span> vec &amp;w = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b2bfb73b72486a5c89f2ab850c7a9b" title="Potentially dangerous, use with care.">_w</a>(); 
    90 <a name="l00106"></a>00106                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b747eca73b16f30370827ba4cc3575" title="access function">_samples</a>(); 
     89<a name="l00105"></a>00105                 <span class="keyword">const</span> vec &amp;w = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#d7f83cc0415cd44ae7cc8b4bdad93aef" title="Potentially dangerous, use with care.">_w</a>(); 
     90<a name="l00106"></a>00106                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#c24966b0aaeb767bc8a6b4fd60931be2" title="access function">_samples</a>(); 
    9191<a name="l00107"></a>00107                  
    92 <a name="l00108"></a>00108                 vec mea = <a class="code" href="classmerger.html#3a46999cd77047d5d4ac6c1cb735a83e" title="return expected value">mean</a>();  
     92<a name="l00108"></a>00108                 vec mea = <a class="code" href="classbdm_1_1merger.html#82c5ad7677315fde8e9fc9381c17ece5" title="return expected value">mean</a>();  
    9393<a name="l00109"></a>00109                  
    9494<a name="l00110"></a>00110                 cout &lt;&lt; sum(w) &lt;&lt; <span class="stringliteral">","</span> &lt;&lt; w*w &lt;&lt;endl; 
    9595<a name="l00111"></a>00111                  
    96 <a name="l00112"></a>00112                 mat Tmp=zeros(<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>(), <a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>()); 
    97 <a name="l00113"></a>00113                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classmerger.html#ce84de1824875b4262eefb640e8c6d69" title="Number of samples used in approximation.">Ns</a>; i++ ) { 
     96<a name="l00112"></a>00112                 mat Tmp=zeros(<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>(), <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>()); 
     97<a name="l00113"></a>00113                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classbdm_1_1merger.html#e91b83e6d9076848c71ae27664072109" title="Number of samples used in approximation.">Ns</a>; i++ ) { 
    9898<a name="l00114"></a>00114                         Tmp+=w ( i ) *outer_product(S ( i ), S(i)); 
    9999<a name="l00115"></a>00115                 } 
    100100<a name="l00116"></a>00116                 <span class="keywordflow">return</span> Tmp-outer_product(mea,mea); 
    101101<a name="l00117"></a>00117         } 
    102 <a name="l00118"></a><a class="code" href="classmerger.html#d6481bfb85df3c3773cee6c83687790a">00118</a>         vec <a class="code" href="classmerger.html#d6481bfb85df3c3773cee6c83687790a" title="return expected variance (not covariance!)">variance</a>()<span class="keyword"> const </span>{ 
    103 <a name="l00119"></a>00119                 <span class="keyword">const</span> vec &amp;w = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b2bfb73b72486a5c89f2ab850c7a9b" title="Potentially dangerous, use with care.">_w</a>(); 
    104 <a name="l00120"></a>00120                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classeEmp.html#31b747eca73b16f30370827ba4cc3575" title="access function">_samples</a>(); 
     102<a name="l00118"></a><a class="code" href="classbdm_1_1merger.html#bcef70ece256b1ab1fa7a3a354eb1cc4">00118</a>         vec <a class="code" href="classbdm_1_1merger.html#bcef70ece256b1ab1fa7a3a354eb1cc4" title="return expected variance (not covariance!)">variance</a>()<span class="keyword"> const </span>{ 
     103<a name="l00119"></a>00119                 <span class="keyword">const</span> vec &amp;w = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#d7f83cc0415cd44ae7cc8b4bdad93aef" title="Potentially dangerous, use with care.">_w</a>(); 
     104<a name="l00120"></a>00120                 <span class="keyword">const</span> Array&lt;vec&gt; &amp;S = <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>.<a class="code" href="classbdm_1_1eEmp.html#c24966b0aaeb767bc8a6b4fd60931be2" title="access function">_samples</a>(); 
    105105<a name="l00121"></a>00121                  
    106 <a name="l00122"></a>00122                 vec tmp=zeros(<a class="code" href="classepdf.html#74da992e3f5d598da8850b646b79b9d9" title="Identified of the random variable.">rv</a>.<a class="code" href="classRV.html#f5c7b8bd589eef09ccdf3329a0addea0" title="Return number of scalars in the RV.">count</a>()); 
     106<a name="l00122"></a>00122                 vec tmp=zeros(<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>()); 
    107107<a name="l00123"></a>00123                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;Ns; i++ ) { 
    108108<a name="l00124"></a>00124                         tmp+=w ( i ) *pow(S ( i ),2); 
    109109<a name="l00125"></a>00125                 } 
    110 <a name="l00126"></a>00126                 <span class="keywordflow">return</span> tmp-pow(<a class="code" href="classmerger.html#3a46999cd77047d5d4ac6c1cb735a83e" title="return expected value">mean</a>(),2); 
     110<a name="l00126"></a>00126                 <span class="keywordflow">return</span> tmp-pow(<a class="code" href="classbdm_1_1merger.html#82c5ad7677315fde8e9fc9381c17ece5" title="return expected value">mean</a>(),2); 
    111111<a name="l00127"></a>00127         } 
    112 <a name="l00129"></a><a class="code" href="classmerger.html#51e425ee1dd46d73353526e6a1c24ca3">00129</a>         <span class="keyword">virtual</span> <a class="code" href="classmerger.html#51e425ee1dd46d73353526e6a1c24ca3" title="for future use">~merger</a>() { 
    113 <a name="l00130"></a>00130                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classcompositepdf.html#c1646fe41ff5934d010e54ccaa46133f" title="Number of mpdfs in the composite.">n</a>; i++ ) { 
    114 <a name="l00131"></a>00131                         <span class="keyword">delete</span> <a class="code" href="classmerger.html#01e222f5df660e1b25c712e82da8df8f" title="Data link for each mpdf in mpdfs.">dls</a> ( i ); 
    115 <a name="l00132"></a>00132                         <span class="keyword">delete</span> <a class="code" href="classmerger.html#d0ff2953f904f12a355f72f51813ca7a" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( i ); 
     112<a name="l00129"></a><a class="code" href="classbdm_1_1merger.html#febef179bb5e074ae6e34d15b4596273">00129</a>         <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1merger.html#febef179bb5e074ae6e34d15b4596273" title="for future use">~merger</a>() { 
     113<a name="l00130"></a>00130                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a>; i++ ) { 
     114<a name="l00131"></a>00131                         <span class="keyword">delete</span> <a class="code" href="classbdm_1_1merger.html#ea3c9f4739f857877fd5285b3973289c" title="Data link for each mpdf in mpdfs.">dls</a> ( i ); 
     115<a name="l00132"></a>00132                         <span class="keyword">delete</span> <a class="code" href="classbdm_1_1merger.html#6c48bf36ec1f3f7b52e3159cf428a193" title="Data Links of rv0 mpdfs - these will be conditioned the (rv,rvc) of mpdfs.">zdls</a> ( i ); 
    116116<a name="l00133"></a>00133                 } 
    117117<a name="l00134"></a>00134         }; 
    118118<a name="l00135"></a>00135  
    119 <a name="l00137"></a><a class="code" href="classmerger.html#6e430acd4eb7eea6078df0eb0a429821">00137</a>         <a class="code" href="classMixEF.html" title="Mixture of Exponential Family Densities.">MixEF</a>&amp; <a class="code" href="classmerger.html#6e430acd4eb7eea6078df0eb0a429821" title="Access function.">_Mix</a>() {<span class="keywordflow">return</span> <a class="code" href="classmerger.html#021d1522cea70db8545d708afd32ec2f" title="Internal mixture of EF models.">Mix</a>;} 
    120 <a name="l00139"></a><a class="code" href="classmerger.html#f79677063ec72cdc175280a3fae8ebc2">00139</a>         <a class="code" href="classeEmp.html" title="Weighted empirical density.">eEmp</a>&amp; <a class="code" href="classmerger.html#f79677063ec72cdc175280a3fae8ebc2" title="Access function.">_Smp</a>() {<span class="keywordflow">return</span> <a class="code" href="classmerger.html#b13c9e8921e6bc6b8dd7bbafb54fa801" title="Projection to empirical density.">eSmp</a>;} 
     119<a name="l00137"></a><a class="code" href="classbdm_1_1merger.html#13f8534f626a37bed91a6d9db54d6d2f">00137</a>         <a class="code" href="classbdm_1_1MixEF.html" title="Mixture of Exponential Family Densities.">MixEF</a>&amp; <a class="code" href="classbdm_1_1merger.html#13f8534f626a37bed91a6d9db54d6d2f" title="Access function.">_Mix</a>() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1merger.html#67ae71cbf0686d7c7bf116fa5db86d81" title="Internal mixture of EF models.">Mix</a>;} 
     120<a name="l00139"></a><a class="code" href="classbdm_1_1merger.html#23f5f8b80df579b05681006fbcfe3699">00139</a>         <a class="code" href="classbdm_1_1eEmp.html" title="Weighted empirical density.">eEmp</a>&amp; <a class="code" href="classbdm_1_1merger.html#23f5f8b80df579b05681006fbcfe3699" title="Access function.">_Smp</a>() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1merger.html#522d43f7aa13fe6701e93fba174dec94" title="Projection to empirical density.">eSmp</a>;} 
    121121<a name="l00140"></a>00140 }; 
    122122<a name="l00141"></a>00141  
    123 <a name="l00142"></a>00142  
     123<a name="l00142"></a>00142 } 
    124124<a name="l00143"></a>00143  
    125125<a name="l00144"></a>00144 <span class="preprocessor">#endif // MER_H</span> 
    126126</pre></div></div> 
    127 <hr size="1"><address style="text-align: right;"><small>Generated on Thu Jan 15 10:50:25 2009 for mixpp by&nbsp; 
     127<hr size="1"><address style="text-align: right;"><small>Generated on Tue Jan 27 16:29:52 2009 for mixpp by&nbsp; 
    128128<a href="http://www.doxygen.org/index.html"> 
    129129<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address>