strawlab / best Goto Github PK

Hello,

when I run example/smart_drug.py, it gives

[-----------------100%-----------------] 110000 of 110000 complete in 11.1 sec
Traceback (most recent call last):
File "smart_drug.py", line 27, in
fig = best.plot.make_figure(M)
File "/usr/local/lib/python2.7/site-packages/best/plot.py", line 130, in make_figure
posterior_std1 = M.trace('group1_std')[:]
File "/usr/local/lib/python2.7/site-packages/pymc/Model.py", line 834, in trace
return self.db.trace(name, chain)
File "/usr/local/lib/python2.7/site-packages/pymc/database/base.py", line 373, in trace
trace = copy.copy(self._traces[name])
KeyError: 'group1_std'

and if I do the M.summary(), it shows the correct means (~ 101) but huge SD (~ hundreds). I am using python 2.7 on OSX. I am wondering what should I do to fix/local the problem?

Sincerely,
David Yu

Hi,
I have installed the software on python 2.7 and all the modules load up OK. The module runs properly but the plot routine does not work. See attached the screen shot of the output. It is in the Matplot.pyplot code when it tries to build the sub-plots is has a problem with the attribute "axisbg". It is not able to determine what it is. Any idea? Charles [email protected].

work.

Thank you for the excellent package. Can you mention what version of pymc were you using? I couldn't find the NoncentralT on any of the pymc installation or codebases I searched. Thank you

Hello,

in the function hdi_of_mcmc, could it be that the standard configuration actually calculates a slightly off HDI?

let's assume we want to calculate the 0.95 interval, then the two lines would subtract the remaining 5% two times and thus calculating the 0.9 interval?

ci_idx_inc = int(np.floor( cred_mass*len(sorted_pts) ))
n_cis = len(sorted_pts) - ci_idx_inc

or am I mistaken somewhere?

Best regards,
Clemens Sauerzopf