如果第一次閱讀，請(qǐng)查看寫(xiě)在前面

import math#讀取文件 with open('../examples/ros_bio28_PROB.txt') as f:file = f.readlines()dna = file[0].rstrip() probability = file[1] probability = probability.split(' ')#根據(jù)提供的各個(gè)概率計(jì)算序列相加的對(duì)數(shù) all_lo = [] for i in probability:i = float(i)gpro = i/2apro = (1-i)/2loga = 0for j in dna:if j == 'A' or j == 'T':loga += math.log(apro, 10)else:loga += math.log(gpro, 10)all_lo.append(loga)#格式化為三位小數(shù) result = '' for line in all_lo:line = format(line, '.3f')result += line + ' ' print(result)

總結(jié)

以上是生活随笔為你收集整理的Rosalind第28题——ros_bio28_PROB的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。