def plotLabelKinsPage(Proteins, protTODict, medDict10, medDict1000, fig, name): time = numpy.linspace(0, 140, 200) k10 = qMS.growthRate(92) k1000 = qMS.growthRate(47) Large1Ax = [] a = 0 for prot in Proteins: d_10 = protTODict[prot]['10'] #p_10String = qMS.calcPercent(p_10, sigfig=2).split('.')[0]+'%' d_1000 = protTODict[prot]['1000'] #p_1000String = qMS.calcPercent(p_1000, sigfig=2).split('.')[0]+'% ax = fig.add_subplot(5,2,a+1) ax.scatter([i[0] for i in medDict10[prot]], [i[1] for i in medDict10[prot]], c = 'r') ax.scatter([i[0] for i in medDict1000[prot]], [i[1] for i in medDict1000[prot]], c = 'b') ax.plot(time, qMS.maxLabFunc(k10, time), c='r') ax.plot(time, qMS.maxLabFunc(k1000, time), c='b') ax.plot(time, qMS.overLabelingFunc(k10, time, d_10), c='r', ls='--') ax.plot(time, qMS.overLabelingFunc(k1000, time, d_1000), c='b', ls='--') ax.set_xlim([0,140]) ax.set_ylim([0,1.00]) ax.set_yticks([0,.25,.50,.75,1.00]) #if a < 7: # plt.setp(ax.get_xticklabels(), visible=False) #if a%2==1: # plt.setp(ax.get_yticklabels(), visible=False) ax.set_xticks([0,20,40,60,80,100,120,140]) ax.yaxis.tick_left() ax.xaxis.tick_bottom() ax.text(0.5, 0.95, prot[4:], transform=ax.transAxes, fontsize=12, verticalalignment='top', horizontalalignment='center') Large1Ax.append(ax) a = a+1 fig.text(0.005, 0.1165, "fraction labeled", fontsize=12, rotation='vertical', verticalalignment='center') fig.text(0.005, 0.317, "fraction labeled", fontsize=12, rotation='vertical', verticalalignment='center') fig.text(0.005, 0.505, "fraction labeled", fontsize=12, rotation='vertical', verticalalignment='center') fig.text(0.005, 0.7, "fraction labeled", fontsize=12, rotation='vertical', verticalalignment='center') fig.text(0.005, 0.9, "fraction labeled", fontsize=12, rotation='vertical', verticalalignment='center') fig.text(0.25, 0.0025, "time (mins)", fontsize=12, horizontalalignment='center') fig.text(0.75, 0.0025, "time (mins)", fontsize=12, horizontalalignment='center') fig.tight_layout() #pylab.savefig(name) return Large1Ax
def turnoverResid10(d, y, t): k = qMS.growthRate(92) #return y - qMS.overLabelingFunc(k, t, d) return qMS.overLabelingFunc(k, t, d) - y