Rosalind第28题——ros_bio28_PROB
生活随笔
收集整理的這篇文章主要介紹了
Rosalind第28题——ros_bio28_PROB
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
如果第一次閱讀,請查看寫在前面
import math#讀取文件 with open('../examples/ros_bio28_PROB.txt') as f:file = f.readlines()dna = file[0].rstrip() probability = file[1] probability = probability.split(' ')#根據提供的各個概率計算序列相加的對數 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)#格式化為三位小數 result = '' for line in all_lo:line = format(line, '.3f')result += line + ' ' print(result)總結
以上是生活随笔為你收集整理的Rosalind第28题——ros_bio28_PROB的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: uniapp 微信登入+tp6 获取微信
- 下一篇: 有哪些你看过五遍以上的电影?