【perl脚本】把01数据转化为基因型数据

最新推荐文章于 2022-04-13 17:06:58 发布

weixin_30436101

最新推荐文章于 2022-04-13 17:06:58 发布

阅读量481

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/liulele/p/8964632.html

本文介绍了一个使用Perl脚本处理生物信息数据的例子，展示了如何读取和解析文本文件，并根据特定规则转换数据。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

use strict;
use warnings;

open DATA, "<data.txt" or die ("cannot open file:$!");
open OUT, ">out.txt";

my @num = (5,6,6,5,5,5);
my @type = ("A","B","C","D","E","F","G");

foreach my $tmp1 (<DATA>){
    chomp $tmp1;
    my @arr;
    @arr = split(/    /,$tmp1);
    print OUT $arr[0]."\t";
    my $order;
    my $o = 0;
    foreach my $tmp2 (@num){
        my $gen = "";
        foreach my $tmp3 (1...$tmp2){
            if ($arr[$o+$tmp3]){
                $gen = $gen.$type[$tmp3-1]
            }
        }
        if (length($gen) == 1){
            print OUT $gen."    ".$gen;
        }elsif(length($gen) == 2){
            print OUT substr($gen,0,1)."    ".substr($gen,1,1);
        }else{
        print OUT "error";
        }
        print OUT "\t";
        $o += $tmp2;
    }
    print OUT "\n";

}

close(DATA);
close(OUT);