def _plotMTCombined(conf, data, subject='__'):
mt, time, idS, id, p = data['all']
time2 = numpy.linspace(time.min(),time.max(),30000)
figname = os.path.join(grPath,subject+'_'+conf[0]+'-'+conf[1]+'_'+conf[2]+'_'+'_MTGraph_All'+grExt)
figtitle = subject+' '+conf[0]+'-'+conf[1]+' '+conf[2]+'Combined replications'
fig = pylab.figure()
host = HostAxes(fig, [0.15, 0.1, 0.65, 0.8])
par1 = ParasiteAxes(host, sharex=host)
host.parasites.append(par1)
host.set_ylabel("MT (ms)")
host.set_xlabel("Time (ms)")
host.axis["right"].set_visible(False)
par1.axis["right"].set_visible(True)
par1.set_ylabel("ID")
par1.axis["right"].major_ticklabels.set_visible(True)
par1.axis["right"].label.set_visible(True)
fig.add_axes(host)
sigmoidArr = sigmoid(p,time2)
t2 = time2[numpy.where(abs(sigmoidArr - (sigmoidArr.min() + sigmoidArr.max())/2) < eps2)[0]]
#t50 = abs(t2-time2.max())
t50 = t2
if isinstance(t50, numpy.ndarray):
t50 = t50.mean()
#host.set_ylim(mt.min(),mt.max())
if conf[0] == '4':
host.set_ylim(200,900)
host.set_xlim(0, 70000)
else:
host.set_ylim(200,1600)
host.set_xlim(0, 90000)
host.scatter(time,mt,s=0.1,marker='x',linewidth=0.1,label="MT")
host.plot(time2,sigmoid(p,time2),'-r',linewidth=1.5,label="Fitted MT")
par1.scatter(time, id, s=1, c='r',label="ID", linewidth=0.5)
msg = "MT50 = %f" % t50
host.text( 0.3 , 0.95, msg, fontsize = 8,\
horizontalalignment='left', verticalalignment='center',\
transform = host.transAxes)
host.vlines(t50,mt.min(),mt.max())
#host.set_ylim(mt.min(),mt.max())
if conf[0] == '4':
host.set_ylim(200,900)
host.set_xlim(0, 70000)
else:
host.set_ylim(200,1600)
host.set_xlim(0, 90000)
par1.set_ylim(min(int(conf[0]),int(conf[1])),max(int(conf[0]),int(conf[1])))
#host.legend()
host.set_title(figtitle)
fig.savefig(figname)
pylab.close()
from mpl_toolkits.axes_grid.parasite_axes import HostAxes, ParasiteAxes
import matplotlib.pyplot as plt
if __name__ == "__main__":
fig = plt.figure(1)
host = HostAxes(fig, [0.15, 0.1, 0.65, 0.8])
par1 = ParasiteAxes(host, sharex=host)
par2 = ParasiteAxes(host, sharex=host)
host.parasites.append(par1)
host.parasites.append(par2)
host.set_ylabel("Density")
host.set_xlabel("Distance")
host.axis["right"].set_visible(False)
par1.axis["right"].set_visible(True)
par1.set_ylabel("Temperature")
par1.axis["right"].major_ticklabels.set_visible(True)
par1.axis["right"].label.set_visible(True)
par2.set_ylabel("Velocity")
offset = (60, 0)
new_axisline = par2._grid_helper.new_fixed_axis
par2.axis["right2"] = new_axisline(loc="right", axes=par2, offset=offset)
fig.add_axes(host)
host.set_xlim(0, 2)
host.set_ylim(0, 2)
from mpl_toolkits.axes_grid.parasite_axes import HostAxes, ParasiteAxes
import matplotlib.pyplot as plt
if __name__ == "__main__":
fig = plt.figure(1)
host = HostAxes(fig, [0.15, 0.1, 0.65, 0.8])
par1 = ParasiteAxes(host, sharex=host)
par2 = ParasiteAxes(host, sharex=host)
host.parasites.append(par1)
host.parasites.append(par2)
host.set_ylabel("Density")
host.set_xlabel("Distance")
host.axis["right"].set_visible(False)
par1.axis["right"].set_visible(True)
par1.set_ylabel("Temperature")
par1.axis["right"].major_ticklabels.set_visible(True)
par1.axis["right"].label.set_visible(True)
par2.set_ylabel("Velocity")
offset = (60, 0)
new_axisline = par2._grid_helper.new_fixed_axis
par2.axis["right2"] = new_axisline(loc="right",
axes=par2,
offset=offset)
fig.add_axes(host)
