Matthew Pitkin edited To_draw_points_from_the__.tex  almost 8 years ago

Commit id: e2914013e276ac9a50195893f616db2b6d75a4ac

deletions | additions      

       

To draw points from the Fermi-Dirac-like distribution in Eq.~\ref{eq:fermidirac} we would draw values uniformly between 0 and 1 and then for each value invert the cdf. The cdf function can be inverted to give  \begin{equation}  H_0 = -\sigma \log{}\left(-e^{-\mu/\sigma} + \left(1+e^{\mu/\sigma} \right)^{-C(H_0)} + e^{1-\mu/\sigma}\left(1+e^{\mu/sigma} e^{1-\mu/\sigma}\left(1+e^{\mu/\sigma}  \right)^{-C(H_0)} \right). \end{equation}  It is definitely worth testing out this as a prior for the searches. I should also look into the Jacobian between this and a prior distribution as I might need that when thinking about using posterior samples of $h_0$ to work out the underlying distributions of $h_0$ (or $\varepsilon$/$Q_{22}$).