The third top-level HDF group, named [`likelihood_stats'], contains quantities to obtain the prior \(p(\vec{\vartheta }|H)\)  and likelihood \(p(\vec{d}(t)|\vec{\vartheta },H)\)  from Eq. REF  for each independent sample.
In order to obtain the prior for each independent sample, the [`likelihood_stats'] HDF group contains a dataset of the natural logarithm of the prior probabilities called [`likelihood_stats/prior']. The datasets in the [`likelihood_stats'] HDF group are one-dimensional arrays indexed by the independent sample (eg. the \(i\) -th element corresponds to the prior probability of the \(i\) -th independent sample) as well. In order to obtain the likelihood for each independent sample, there is a dataset containing the natural logarithm of the likelihood ratio \(\Lambda \)  called [`likelihood_stats/loglr']. The likelihood ratio \(\Lambda \)  is defined as [1]}
\(\log \Lambda = \log \frac{p(\vec{d}(t)|\vec{\vartheta },H)}{p(\vec{d}(t)|\vec{n})}\) 
