免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
12下一页
最近访问板块 发新帖
查看: 10310 | 回复: 10
打印 上一主题 下一主题

请问如何用perl统计一个字符串里面特定字符的个数? [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2011-09-09 16:34 |只看该作者 |倒序浏览
本帖最后由 gigigou 于 2011-09-09 21:56 编辑

比如说有一个文件,每一行都是类似于以下的字符串(其实是DNA序列)
TCTCATGTGAAAAACTATATCAATAATATAAAAACA
TTGCTGGCTCCGGTGCGCTTGTGACGGCCCGTGGAA
TGTTCAATCGTCTTGAATGCAAAAATTGTTTTATTA
ACATGNNANNNNNNNNNNNNNNNNNNNNNNNNNNNN

用perl处理的时候先打开文件,然后while(<IN>
但是这样每次$_得到的都是整个一行的字符,也就是第一次  $_='TCTCATGTGAAAAACTATATCAATAATATAAAAACA'
我想要统计其中某个特定字符的个数(比如G的个数),想用shift,pop,但是变量类型不对
我刚开始学perl,请教各位应该怎么做?
谢谢!

谢谢各位的帮助!
祝大家中秋快乐!

论坛徽章:
0
2 [报告]
发表于 2011-09-09 16:51 |只看该作者
$_='TCTCATGTGAAAAACTATATCAATAATATAAAAACA';
my $count = () = /G/g;
print $count, "\n";

论坛徽章:
0
3 [报告]
发表于 2011-09-09 16:52 |只看该作者
C:\>perl -e "@a=shift=~/a/g;print scalar(@a)" asdaad
3

论坛徽章:
0
4 [报告]
发表于 2011-09-09 16:58 |只看该作者
剛好之前有寫過~XD
  1. While(<IN>)
  2. {
  3.   $all_length=length($_);
  4.   $_=~s/G//g;
  5.   $no_G_length=length($_);
  6.   $G_count=$all_length-$no_G_length;
  7.   $sum_G_count+=$G_count;
  8. }
  9. print $sum_G_count;
复制代码

论坛徽章:
0
5 [报告]
发表于 2011-09-09 17:33 |只看该作者
本帖最后由 thomaslyang 于 2011-09-09 18:19 编辑
$_='TCTCATGTGAAAAACTATATCAATAATATAAAAACA';
my $count = () = /G/g;
print $count, "\n";
3P主义 发表于 2011-09-09 16:51
  1. #!/usr/bin/python
  2.   print( "TCTCATGTGAAAAACTATATCAATAATATAAAAACA".count("G"), end="\n");
复制代码
  1. #!/usr/bin/perl -w

  2. my $string = "TCTCATGTGAAAAACTATATCAATAATATAAAAACA";

  3. my $count = 0;
  4. ++$count while($string =~ m/A/g);
复制代码

论坛徽章:
0
6 [报告]
发表于 2011-09-09 17:37 |只看该作者
  1. perl -0777 -wnl -e 'print s/a//g, "\n"' <filename>
复制代码

论坛徽章:
0
7 [报告]
发表于 2011-09-09 17:51 |只看该作者
再來一個~XD
  1. $/="G";
  2. @G_splice=<DATA>;
  3. print $G_splice[$#G_splice]=~/G/ ? scalar(@G_splice) : $#G_splice;
复制代码
TIMTOWTDI

论坛徽章:
0
8 [报告]
发表于 2011-09-09 18:16 |只看该作者
本帖最后由 thomaslyang 于 2011-09-09 18:28 编辑
zzmmgjj 发表于 2011-09-09 17:37
  1. $_ = "TCTCATGTGAAAAACTATATCAATAATATAAAAACA";

  2. print s/(A)/$1/g, "\n";
复制代码

论坛徽章:
0
9 [报告]
发表于 2011-09-09 21:55 |只看该作者
回复 5# thomaslyang


    谢谢!
实在太感谢了!

论坛徽章:
0
10 [报告]
发表于 2011-09-09 21:55 |只看该作者
谢谢以上各位!
非常感谢你们的热情帮助!
祝大家中秋快乐!
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP