24小时热门版块排行榜    

查看: 380  |  回复: 3

xiezhancnu

木虫 (小有名气)

[求助] perl程序求助~ 已有2人参与

chr        start        end
chr10        2        4
chr1        2        3
chr3        2        4
chr1        1        2
chr2        1        3
chr10        3        5

排序成:
chr        start        end
chr1        1        2
chr1        2        3
chr2        1        3
chr3        2        4
chr10        2        4
chr10        3        5
首先按照第一列染色体顺序排chr1 chr2 ……chr9,chr10,chr11……chr21,chr22,chrX,chrY
相同染色体则按照位置的前后顺序进行排列。

先谢过!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

退化or进化

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
为什么要自己写perl脚本呢?Linux下的sort不就可以完成?或者更专业的用bedtools啊...
2楼2015-05-08 08:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上官神速

新虫 (初入文坛)

【答案】应助回帖

@ARGV==1 || die "No way";
my ($text) = @ARGV;
my @tmp;
my $hd;
open($hd,"< $text" || die "can not open it";
while(<$hd>
{
      chomp;
          s/\s+//g;
          s/chr//;
          push(@tmp,$_);
}
my @sorted = sort {$a <=> $b } @tmp;
foreach(@sorted){s/(.*)([0-9])([0-9])$/$1 $2 $3/;print "chr$_\n";}
结果:
chr1 1 2
chr1 2 3
chr2 1 3
chr3 2 4
chr10 2 4
chr10 3 5
3楼2015-06-12 16:22:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上官神速

新虫 (初入文坛)

【答案】应助回帖

@ARGV==1 || die "No way";
my ($text) = @ARGV;
my @tmp;
my $hd;
open($hd,"< $text" || die "can not open it";
while(<$hd>
{
      chomp;
          s/^\s+//g;
          s/chr//;
          s/\s+/-/;
          push(@tmp,$_);
}
my @sorted = sort {$a <=> $b } @tmp;
foreach(@sorted){s/\-/ /g;print "chr$_\n";}


结果:chr1 2        3
chr1 1        2
chr2 1        3
chr3 2        4
chr10 2        3
chr10 3        5
4楼2015-06-12 16:39:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiezhancnu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 材料调剂 5+3 想要一壶桃花水 2026-03-25 6/300 2026-03-25 18:20 by xcjcqu
[考研] 材料与化工考研调剂 +9 孅華 2026-03-22 9/450 2026-03-25 13:09 by cmz0325
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 07化学280分求调剂 +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[考研] B区考研调剂 +4 yqdszhdap- 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 环境学硕288求调剂 +8 皮皮皮123456 2026-03-22 8/400 2026-03-23 23:47 by 热情沙漠
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
信息提示
请填处理意见