美国加州大学圣地亚哥分校生信课程练习题week1.1

问题需求:在text中统计pattern重复次数,需要注意的是TGGACTCTG和TGGACTCTGGACTCTG都符合pattern。

text:GGGTTATGGACTCTGGACTCTCAATGGACTCTGGACTCACTGGACTCTGGACTCGGCAATGGACTCATGGACTCTTGGACTCTTGGACTCTACTCGTGGACTCTGGACTCTAATGGACTCAATGGACTCTGGACTCATTTGGACTCGTTGGACTCCTGGACTCGTTGGACTCTGGACTCTTGGACTCCACTCTGGACTCACTGGACTCTGGACTCTTTGGACTCACCGCTGGACTCCTATGGACTCAGTGGACTCGAGTGGACTCCAGTAGTAGGATGTGGACTCTTGGACTCCGGGGACTTAAGCCTGGACTCGGTGGACTCACGATGGACTCTGGACTCTAACTGGACTCGGTGGACTCTGGACTCACCGACAGCGTCTGGACTCCACTGCAGAATGGACTCTGTGGACTCGAACTGGACTCGGCGGGAGATTTGGACTCATGGACTCTGGACTCACTGGACTCCATGGACTCAATTTTGGACTCTGGACTCTCGTGGACTCTGGACTCCCCTGGACTCCTATGGACTCCGTGGACTCTGGACTCTGGACTCTGGACTCTGGACTCGGACTGGACTCATGGACTCTGGACTCTGGACTCCGTTGGACTCTTGGACTCCTCTGGACTCTGTGGACTCAGTGGACTCTTGGACTCGTGGACTCATGGACTCGGTTTTAGTGGACTCGGTATGGACTCTCTGGACTCTGGACTCAGCATGGACTCATGGACTCCCGAGTCATGGACTCTGGACTCGGACGCTGGACTCTGGACTCTATGGACTCATTGGACTCTGGACTCTTTGGACTCAAATGGACTCACTGGACTCTGGACTCCATGGACTCGTGGACTCACTGGACTCGATCCGAATGGACTCTGGACTCACTCCTGGACTCCATGGACTCTTGGACTCGGGTGGACTCGTGGACTCCTCTGGACTCAGATGGACTCAGACCCTGGACTCGTTGGACTCTCTTTGGACTC

pattren:TGGACTCTG

代码实现:

python代码实现:

txt, i= {},1path ="dataset_2_6.txt"
f =open(path,encoding='utf-8')
for line in f:    
  txt[i] = line
  i +=1    
  txt.update(txt)
f.close() 
text = str(txt[1].strip())
pattern = str(txt[2].strip())
def PatternCount(text, pattern):
count=0
for i in range(len(text)-len(pattern)+1):
  if text[i:i+len(pattern)] == pattern:
    count=count+1
return count 
print(PatternCount(text, pattern))

R代码实现:

library(stringr)
data = read.table('dataset_2_6.txt', stringsAsFactors = F)
str = str_split(data[1,1],'') %>%
        as.data.frame()

ncount = 0
for (i in 1:(nrow(str) - nchar(data[2,1]) + 1)) {
        j = i+nchar(data[2,1]) - 1
        rep.1 = paste(str[i:j,1], sep = '')
        
        rep.2 = ''
        
        for (m in 1:length(rep.1)) {
                rep.2 = paste(rep.2,rep.1[m], sep = '')
        }
        
        if (rep.2 == data[2,1]) {
                ncount = ncount + 1
        }
        
        if (i == (nrow(str) - nchar(data[2,1]) + 1)) {
               print(ncount) 
        }
}

shell代码实现:

text=$(sed -n '1p' dataset_2_6.txt)
pattern=$(sed -n '2p' dataset_2_6.txt)
len_text=echo ${#text}
len_pattern=echo ${#pattern}
len_range=$(($len_text-$len_pattern))
count=0
for((i=0;i<=$len_range;i++))
  do
    temp=$(echo ${text:$i:$len_pattern})
    if [ $temp == $pattern ]; then
      count=$(($count+1))
    fi
done
echo $count
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 美国加州大学圣地亚哥分校Python实战项目总结 课程内容和老师教学情况 课程内容当然就是讲Python啦!Pyt...
    哲哲99阅读 1,846评论 0 0
  • 表情是什么,我认为表情就是表现出来的情绪。表情可以传达很多信息。高兴了当然就笑了,难过就哭了。两者是相互影响密不可...
    Persistenc_6aea阅读 127,941评论 2 7
  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 11,321评论 0 4

友情链接更多精彩内容