【原】如何利用perl快速筛选差异基因？

我们只要将我们的input.txt文件和perl脚本放到同一个文件下，然后双击相应的脚本就可以得到我们的结果文件。

下面我们就将代码放在这里：

筛选上调差异基因的代码：

#!/usr/bin/perl -w

open IN,"<input.txt";

open OUT,">up.txt";

while (<IN>) {

    chomp; 

    my @gene=split /\s+/,$_; 

    if ($gene[1] <=0.05 && $gene[3] >=2) {

        print OUT "$_\n";

    } 

    else {

        next;

    }

}

close IN;

close OUT;

筛选下调差异基因的代码：

#!/usr/bin/perl -w

open IN,"<input.txt";

open OUT,">up.txt";

while (<IN>) {

    chomp; 

    my @gene=split /\s+/,$_; 

    if ($gene[1] <=0.05 && $gene[3] <=-2) {

        print OUT "$_\n";

    } 

    else {

        next;

    }

}

close IN;

close OUT;