eggplantbren/DNest4

Python postprocessing functions yielding different results

tripathi opened this issue · 7 comments

tripathi commented

Hi!

When running the Gaussian python example from commit 190180.. on a Mac (OSX10.10) with Python 2.7, it runs into NaNs for log_x_diff, and thus log_z when calculated by analysis.py's postprocess->compute_stats function. See
analysispy.txt

However, if I manually postprocess the output using the postprocess function in deprecated.py, I get a reasonable value for log_z. See
deprecatedpy.txt

Should I continue to use the deprecated function?

Note I've used MACOSX_DEPLOYMENT_TARGET=10.10 python setup.py install for builds (although I get the same result with TARGET=10.9).

dfm commented

I can't reproduce this even with the same commit. What version of numpy? It's doing something weird with the way it handles underflow...

tripathi commented

Thanks for the quick response! I have numpy 1.9.2 (installed using canopy).

dfm commented

Can you also attach the text files that the Gaussian example saves?

tripathi commented

posterior_sample.txt
weights.txt
levels.txt
sample_info.txt
sample.txt
sampler_state.txt
sample_log_X.txt (You can see all the infs)
stats.txt

Note these are from a different run as that mentioned above, so I'm also attaching the new log file
newlog.txt

Thanks!

omaclaren commented

(I just came across the same issue btw - I'll try to post a more informative comment later)

eggplantbren commented

The workaround is to always use

import dnest4.classic
dnest4.classic.postprocess()
omaclaren commented

Awesome, thanks, this worked for me.

And, I just realised that this is also discussed under #10

(sorry!). So to confirm - classic works for me but the default doesn't (i.e. it gives nans). I'm running Python 2.7.11.