In [24]:

    

import matplotlib
import Bio.SeqIO
import Bio.Restriction
import scipy
import matplotlib.pyplot as plt
%matplotlib inline
import math
import random
import numpy as np


    

In [33]:

    

TCF = [x for x in np.arange(0.65, 1.36, 0.01) ]
FP = [np.ceil(y*63) for y in TCF]

plt.plot(TCF,FP)
plt.xlabel("TCF")
plt.ylabel("FP")
plt.title("For given example, UFC = 63")


    

    
Out[33]:





<matplotlib.text.Text at 0x7f42c682c990>






    

In [47]:

    

TCF = [x for x in np.arange(0.65, 1.36, 0.01) ]
FP = [np.ceil(y*63) for y in TCF]
plt.plot(TCF,FP, label="UFC (ex.) = 63")

TCF = [x for x in np.arange(0.65, 1.36, 0.01) ]
FP = [np.ceil(y*43) for y in TCF]
plt.plot(TCF,FP, label="UFC (min) = 43")

TCF = [x for x in np.arange(0.65, 1.36, 0.01) ]
FP = [np.ceil(y*85) for y in TCF]
plt.plot(TCF,FP, label='UFC (max) = 85')

plt.xlabel("TCF")
plt.ylabel("FP")
plt.legend(loc = 'upper left')
plt.show()


    

In [27]:

    

len(TCF)


    

    
Out[27]:





71

In [29]:

    

len(FP)


    

    
Out[29]:





71

In [ ]: