Skip to content
Snippets Groups Projects
Commit e8a0f4dc authored by Stefano Serafin's avatar Stefano Serafin
Browse files

fixed parameter histogram plots

parent 51152354
No related branches found
No related tags found
No related merge requests found
......@@ -167,12 +167,14 @@ def plot_CBL_PE(exp,figtitle,parameter_id=0,plot_spread=False,ax=None):
ax1.set_ylabel('probability density')
#
ax2 = fig.add_subplot(1,2,2)
ax2.hist(par_tran[0,:], bins=np.linspace(np.log(0.01),np.log(5),51), density=True, histtype='step', color='blue', label='initial')
#ax2.hist(par_tran[0,:], bins=np.linspace(np.log(0.01),np.log(5),51), density=True, histtype='step', color='blue', label='initial')
ax1.hist(par_tran[0,:], bins=np.linspace(-20,20,51), density=True, histtype='step', color='blue', label='initial')
ax2.axvspan(p10beg_tran,p90beg_tran,color='blue',lw=0,alpha=0.1)
ax2.hist(par_tran[-1,:], bins=np.linspace(np.log(0.01),np.log(5),51), density=True, histtype='step', color='red', label='final')
#ax2.hist(par_tran[-1,:], bins=np.linspace(np.log(0.01),np.log(5),51), density=True, histtype='step', color='red', label='final')
ax1.hist(par_tran[-1,:], bins=np.linspace(-20,20,51), density=True, histtype='step', color='red', label='final')
ax2.axvspan(p10end_tran,p90end_tran,color='red',lw=0,alpha=0.1)
ax2.axvline(x=pt[parameter_id](exp.nr.parameter_true,kind='dir'),color='black',dashes=[3,1])
ax2.set_xlabel(r'log($p$)')
ax2.set_xlabel(r"$p^\prime")
ax2.set_ylabel('probability density')
#
fig.tight_layout()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment