用python获取复杂度的P值

#coding:utf-8
import os
import matplotlib.pyplot as plt
plt.switch_backend('agg')

path = os.getcwd()
f_fasta = open(path+'/fac_cpx.txt','r')
list_data = []
for line in f_fasta.readlines():
	line = line.strip('\n')
	list_data.append(line)
f_fasta.close()

#print(len(list_data))
List_Num = []
for i in range(len(list_data)):
	a_seq = list_data[i]
	A_seq = a_seq.split(',')
	A_seq = list(map(int,A_seq))
	List_Num.append(A_seq)

list_box_1_5 = []
for j in range(5):
	list_box = []
	for i in range(len(List_Num)):
		list_box.append(List_Num[i][j])
	list_box_1_5.append(list_box)

#list_box_1_5 includes 5 types ,and each has 4226

#LCSD = [16,66,91,93,92]
LCSD = [16,49,70,78,81]
#a = [[1,1,3,4,5],[6,1,1,1,5],[3,2,7,8,9]]
def P_value(a):
    P = []
    for j in range(len(LCSD)):
        a_1 =[]
        for i in range(len(a[j])):
            if LCSD[j] >= a[j][i]:
                a_1.append(a[j][i])
        p = round(len(a_1)/len(a[j]),3)
        P.append(p)
    return P


P = P_value(list_box_1_5)
print(P)

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值