应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 生物科学 » 求助 一个perl代码修改~~看看吧...
51/1返回列表
查看: 2479  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

landspace

新虫 (初入文坛)

[求助] 求助 一个perl代码修改~~看看吧...

请教,


这是我的一个代码 我想的是利用blast出来的结果（s. start, s. end,） 提取下面这个基因组中有关的序列。

我不知道哪里错了还请大家帮忙指正！！

O(∩_∩)O谢谢！！麻烦了！
CODE:
01.#!/user/bin/perl -w

02.

03.$ARGV[0]= “H:\Pro\blast\bin\test1.2”; # blast文件路径赋值给ARGV[0]

04.open(FIN,"$ARGV[0]")||die; #打开blast文件

05.open(FIN1,"H:\Pro\blast\bin\data\ testdatabase.fa ")||die; #打开基因组文件

06.open(FOUT,"> “H:\Pro\perl\bin\out]")||die; # 输出文件

07.

08.

09.while(){ #逐行读取文件

10.chomp; #并加换行符

11.

12.$l=$_; #将每一行的文件赋值给$|

13.if($l=~/^>(\S+)/){ #匹配以>的开头

14.$id=$1; #提取序列名称

15.$sq{$id}=""; #构建哈希表%sq，以序列名称$id为key

16.}

17.else{

18.$sq{$id}.=$l; #逐行将序列写入字符串

19.}

20.}

21.

22.while(){ #逐行读取文件

23.chomp; #去掉换行符

24.next if/^#/; #去掉注释行 跳过 以#开头的行

25.@a=split("\t"); #以空格为边界打散成数据组，分割文件

26.next if(@a < 8);#跳过元素小于8的行 此步可以省略吗？

27.$st=$a[8]; $len=$a[9]-$a[8];#计算序列长度 s.end –s.start？

28.if($a[8] eq "-"){ #表示$a[8]”-”符号

29.#$st=$a[7]; $len=-$len; #基因不存在？

30.}

31.$seq=substr($sq{$a[0]},$st-3000,$len+3000);#将基因组中$st即s.start为起始的$len个碱基写入子串seq 上游加3000bp

32.if($a[3] eq "-"){

33.$seq=reverse($seq); #倒置序列，得到互补序列

34.$seq=lc($seq);

35.$seq=~s/a/T/g;

36.$seq=~s/t/A/g;

37.$seq=~s/c/G/g;

38.$seq=~s/g/C/g;

39.}

40.#$a[-1]=~s/intron Parent=//; #序列名称输出 对应关系不确定？

41.if($a[8]=~/^ID=([^;]+);/){

42.$name=$1;

43.}

44.else{

45.print "$a[8]\n";

46.}

47.print FOUT ">$name\n$seq\n";#输出

test1.2文件是：
# Fields: Query id, Subject id, % identity, alignment length, mismatches, gap openings, q. start, q. end, s. start, s. end, e-value, bit score

CeraIN-CO3_12S_12546bp      gi|154076492|ref    99.43        4377         4       1       1       4356         2321         6697         0.0   8510

CeraIN-CO3_12S_12546bp      gi|154076401|ref 100.00        21     0       0       10946       10966       4091         4071         0.64 42.1

……


testdatabase.fa 文件是：

>gi|154076413|ref| (based on CSol2.0 SCAFFOLD6002)

TACACCATTATCTTCCATAGTTTTAAATATAGACATGTAATACTCAGCAAAATGGTTCCGCAGTTGCCTC

CAATATTTCTCTATTCTCTCGTTAGAAGTACTCTTACCCTTTATGAAACTTTTTTCTCCAGCTTTCTCAT

CTGTATGATGGAAACGTAAAGCCATCTGAAGTAATTCTACAACTGTATTTTCAGTCCCTTTGTCTGATCG


>gi|154074987|ref| (based on CSol2.0 SCAFFOLD5715)

TAATTGAGAACGGGCGCAGTAGTTTTCCCGTAACCTTTATTGCAGTGGCTCGGCTACTTTCGTAGGTTTG

GCGTATTTAGTTTCCGTCCGATTTTCGTACGAGCTCTCCAAGCGCCGCGTGTTAAATTCGATCGCCTTTT

TTTGCGCCTATTTTCTTTCAAATAATTGCGGAAGCTCTGATGGCTAGTTAATTCTAACTCTGACGCGTTT

TTTGCACTGCGCGTTGTTTCCAAATAACCACGGAGCTTAATTTCCAAATTGACGACGATTTTCTCGCTGC

CGGGTCCATTGTCCGATATGGTGTCGCATCCGTAATTTGCCTCAGAGTATGTATCGACCTCTGACTCAGA

GCTCAGTGCAGAATCCCCATCGGCGCTCATGGTTAAGTCGTGCGACGTTTGCCACTCAAAAATATATGCA

……

错误代码：Unrecognized escape \P passed through at test.pl line 3.
Unrecognized escape \P passed through at test.pl line 5.
Unrecognized escape \d passed through at test.pl line 5.
Unrecognized escape \P passed through at test.pl line 6.
Unrecognized escape \p passed through at test.pl line 6.
Missing braces on \o{} at test.pl line 6, within string
Missing right curly or square bracket at test.pl line 48, at end of line
syntax error at test.pl line 48, at EOF
Execution of test.pl aborted due to compilation errors.

谢谢，麻烦了！！
金币只有一个了- -
回复此楼

» 收录本帖的淘帖专辑推荐

perl语言专栏

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-01-11 21:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

landspace

新虫 (初入文坛)
我更改了一下 后来就出现了这个问题
CODE:
#!/user/bin/perl -w

$ARGV[0]= "H:\\Pro\\blast\\bin\\test1.2"; # blast文件路径赋值给ARGV[0]
open(FIN,"$ARGV[0]")||die; #打开blast文件
open(FIN1,"H:\\Pro\\blast\\bin\\data\\testdatabase.fa")||die; #打开基因组文件
open(FOUT, "> H:\\Pro\\perl\\bin\\out")||die; # 输出文件


while(){ #逐行读取文件
        chomp; #并加换行符

        $l=$_; #将每一行的文件赋值给$|
        if($l=~/^>(\S+)/){ #匹配以>的开头
                $id=$1; #提取序列名称
                $sq{$id}=""; #构建哈希表%sq，以序列名称$id为key
        }
        else{
                $sq{$id}.=$l; #逐行将序列写入字符串
        }
}
while(){ #逐行读取文件
        chomp; #去掉换行符
next if/^#/;  #去掉注释行 跳过 以#开头的行
        @a=split("\t"); #以空格为边界打散成数据组，分割文件
        next if(@a < 8);#跳过元素小于8的行 此步可以省略吗？
        $st=$a[8]; $len=$a[9]-$a[8];#计算序列长度 s.end –s.start？
        if($a[8] eq "-"){ #表示$a[8]”-”符号
                #$st=$a[7]; $len=-$len; #基因不存在？
        }
        $seq=substr($sq{$a[0]},$st-3000,$len+3000);#将基因组中$st即s.start为起始的$len个碱基写入子串seq 上游加3000bp
        if($a[3] eq "-"){
                $seq=reverse($seq); #倒置序列，得到互补序列
                $seq=lc($seq);
                $seq=~s/a/T/g;
                $seq=~s/t/A/g;
                $seq=~s/c/G/g;
                $seq=~s/g/C/g;
        }
        #$a[-1]=~s/intron Parent=//; #序列名称输出 对应关系不确定？
        if($a[8]=~/^ID=([^;]+);/){
                $name=$1;
        }
        else{
                print "$a[8]\n";
        }
        print FOUT ">$name\n$seq\n";#输出

}

Use of uninitialized value in substr at test2.pl line 30, line 5.
897
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 5.
Use of uninitialized value in substr at test2.pl line 30, line 6.
896
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 6.
Use of uninitialized value in substr at test2.pl line 30, line 7.
947
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 7.
Use of uninitialized value in substr at test2.pl line 30, line 8.
896
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 8.
Use of uninitialized value in substr at test2.pl line 30, line 9.
947
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 9.
Use of uninitialized value in substr at test2.pl line 30, line 10.
900
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 10.
Use of uninitialized value in substr at test2.pl line 30, line 11.
951
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 11.
Use of uninitialized value in substr at test2.pl line 30, line 12.
896
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 12.
Use of uninitialized value in substr at test2.pl line 30, line 13.
896
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 13.
Use of uninitialized value in substr at test2.pl line 30, line 14.
946
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 14.
Use of uninitialized value in substr at test2.pl line 30, line 15.
185
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 15.
Use of uninitialized value in substr at test2.pl line 30, line 16.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 16.
Use of uninitialized value in substr at test2.pl line 30, line 17.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 17.
Use of uninitialized value in substr at test2.pl line 30, line 18.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 18.
Use of uninitialized value in substr at test2.pl line 30, line 19.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 19.
Use of uninitialized value in substr at test2.pl line 30, line 20.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 20.
Use of uninitialized value in substr at test2.pl line 30, line 21.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 21.
Use of uninitialized value in substr at test2.pl line 30, line 22.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 22.
Use of uninitialized value in substr at test2.pl line 30, line 23.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 23.
Use of uninitialized value in substr at test2.pl line 30, line 24.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 24.
Use of uninitialized value in substr at test2.pl line 30, line 25.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 25.
Use of uninitialized value in substr at test2.pl line 30, line 26.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 26.
Use of uninitialized value in substr at test2.pl line 30, line 27.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 27.
Use of uninitialized value in substr at test2.pl line 30, line 28.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 28.
Use of uninitialized value in substr at test2.pl line 30, line 29.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 29.
Use of uninitialized value in substr at test2.pl line 30, line 30.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 30.
Use of uninitialized value in substr at test2.pl line 30, line 31.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 31.
Use of uninitialized value in substr at test2.pl line 30, line 32.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 32.
Use of uninitialized value in substr at test2.pl line 30, line 33.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 33.
Use of uninitialized value in substr at test2.pl line 30, line 34.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 34.
Use of uninitialized value in substr at test2.pl line 30, line 35.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 35.
Use of uninitialized value in substr at test2.pl line 30, line 36.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 36.
Use of uninitialized value in substr at test2.pl line 30, line 37.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 37.
Use of uninitialized value in substr at test2.pl line 30, line 38.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 38.
Use of uninitialized value in substr at test2.pl line 30, line 39.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 39.
Use of uninitialized value in substr at test2.pl line 30, line 40.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 40.
Use of uninitialized value in substr at test2.pl line 30, line 41.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 41.
Use of uninitialized value in substr at test2.pl line 30, line 42.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 42.
Use of uninitialized value in substr at test2.pl line 30, line 43.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 43.
Use of uninitialized value in substr at test2.pl line 30, line 44.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 44.
Use of uninitialized value in substr at test2.pl line 30, line 45.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 45.
Use of uninitialized value in substr at test2.pl line 30, line 46.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 46.
Use of uninitialized value in substr at test2.pl line 30, line 47.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 47.
Use of uninitialized value in substr at test2.pl line 30, line 48.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 48.
Use of uninitialized value in substr at test2.pl line 30, line 49.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 49.
Use of uninitialized value in substr at test2.pl line 30, line 50.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 50.
Use of uninitialized value in substr at test2.pl line 30, line 51.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 51.
Use of uninitialized value in substr at test2.pl line 30, line 52.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 52.
Use of uninitialized value in substr at test2.pl line 30, line 53.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 53.
Use of uninitialized value in substr at test2.pl line 30, line 54.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 54.
Use of uninitialized value in substr at test2.pl line 30, line 55.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 55.
Use of uninitialized value in substr at test2.pl line 30, line 56.
832
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 56.
Use of uninitialized value in substr at test2.pl line 30, line 57.
362
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 57.
Use of uninitialized value in substr at test2.pl line 30, line 58.
510
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 58.
Use of uninitialized value in substr at test2.pl line 30, line 59.
752
Use of uninitialized value $name in concatenation (.) or string at test2.pl line
46, line 59.

请问这个是什么原因？
赞一下 回复此楼
3楼2013-01-13 12:13:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

wizardfan

至尊木虫 (著名写手)

优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
提示的很清楚 Unrecognized escape \P passed through at test.pl line 3.
$ARGV[0]= “H:\Pro\blast\bin\test1.2”;
在perl里\是特殊字符，表示转换成特定符号，比如最常用的\n，就是回车。如果你要用目录的话，要用两个\\而不是单个\
赞一下 回复此楼
2楼2013-01-12 07:43:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wizardfan

至尊木虫 (著名写手)

优秀版主

【答案】应助回帖

你自己的分析是什么？知道use strict;有什么用？
赞一下 回复此楼
4楼2013-01-14 07:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hydraphenix

银虫 (小有名气)
太菜了，版主都不屑于回答了。
赞一下 回复此楼
5楼2013-01-14 17:16:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +15 koalala 2026-08-24 17/850 2026-08-24 14:16 by hhs666
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭