This is a very old post, but I just figured out what’s going on. Carlos got wrong answers in his MCMC because with the diffuse prior there is a massive phase transition and you need Nested Sampling. There is nothing “wrong” about the diffuse prior in this problem.
Cheers,
Brendon
google reader seems to have a cache of its posts but that ended at the end of 2009. So sometime around the 1st of the year 2010, WP started producing garbage RSS.
]]>Also thanks for the other comments. I will probably include those comments by moving this page off to gWave shortly — it is much easier to update and make changes there than in the blog format.
]]>Have you looked at the Sherpa scripts page. For instance, some of the plot options are discussed in the plotting thread, although I couldn’t see the error bars explicitly addressed. Try
p = get_data_plot_prefs()
print(p)
p["yerrorbars"] = False
To change the units of the y error bars, use the set_analysis() routine, for instance
set_analysis("wave", type="counts")
Also, for error estimation, we suggest you use
conf()
rather than projection nowadays, and note that the error analysis will be done in parallel (so expect all your cores to get used up).
And finally,
,ahelp set_analysis
can also be used to avoid quoting the first argument by using a comma before the command. This is an ipython-ism, and can be used for other commands, e.g.
,load_pha foo.pha
Doug
]]>