Skip to content

'inf' in marginalized log likelihood

posterior plots fail due to this error

Traceback (most recent call last):
  File "/home/vinaya.valsan/virtualenvs/igwn-py39-20221118-rapidpe-rift-dev/bin/compute_posterior.py", line 7, in <module>
    exec(compile(f.read(), __file__, 'exec'))
  File "/home/vinaya.valsan/rapidPE/RapidPE_RIFT_Developments/RapidPE-RIFT_fork/rapidpe-rift-pipe/bin/postscripts/compute_posterior.py", line 227, in <module>
    sample_dict = postutils.get_posterior_samples(
  File "/home/vinaya.valsan/rapidPE/RapidPE_RIFT_Developments/RapidPE-RIFT_fork/rapidpe-rift-pipe/src/rapidpe_rift_pipe/postscript_utils.py", line 225, in get_posterior_samples
    N = N_mn.rvs(1)[0]
  File "/cvmfs/software.igwn.org/conda/envs/igwn-py39-20221118/lib/python3.9/site-packages/scipy/stats/_multivariate.py", line 3257, in rvs
    return self._dist.rvs(self.n, self.p, size, random_state)
  File "/cvmfs/software.igwn.org/conda/envs/igwn-py39-20221118/lib/python3.9/site-packages/scipy/stats/_multivariate.py", line 3206, in rvs
    return random_state.multinomial(n, p, size)
  File "mtrand.pyx", line 4249, in numpy.random.mtrand.RandomState.multinomial
  File "_common.pyx", line 376, in numpy.random._common.check_array_constraint
  File "_common.pyx", line 362, in numpy.random._common._check_array_cons_bounded_0_1
ValueError: pvals < 0, pvals > 1 or pvals contains NaNs

from postutils.get_posterior_samples

Margll_sel = [       inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
        inf  4.6343778        inf  4.5907099  4.678034
 20.046766         inf  4.5894783        inf        inf
        inf        inf        inf        inf        inf
        inf        inf        inf 13.745892         inf
        inf        inf        inf        inf        inf
        inf  4.5876511        inf        inf        inf
        inf        inf        inf        inf        inf
        inf        inf  4.5898155        inf        inf
        inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
        inf  4.651869         inf  4.5930583        inf
        inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
        inf        inf        inf        inf        inf
 35.795673         inf        inf        inf        inf
        inf        inf        inf        inf]
margL_normed  = [nan nan nan nan nan nan nan nan nan nan nan  0. nan  0.  0.  0.
 nan  0. nan nan nan nan nan nan nan nan nan nan  0. nan nan nan
 nan nan nan nan  0. nan nan nan nan nan nan nan nan nan nan  0.
 nan nan nan nan nan nan nan nan nan nan nan nan nan  0. nan  0.
 nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
 nan nan nan nan nan nan nan nan nan nan  0. nan nan nan nan nan
 nan nan nan]

for this run /home/vinaya.valsan/rapidPE/RapidPE_RIFT_Developments/RapidPE-RIFT_fork/pycbc_pastro/output/pycbc_pastro/S230715bw