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