view Perl/ConvertOrthofinderMatrix.pl @ 13:152d7c43478b draft default tip

Uploaded

#!/usr/bin/perl

use strict;

my $indir = $ARGV[0];
my $matrix = $ARGV[1];
my $out = $ARGV[2];
my $strain_names = $ARGV[3];

my %strains_of_gb;
open(F,$strain_names);
while(<F>){
	my $line = $_;
	$line =~s/\n//g;$line =~s/\r//g;
	my ($gb,$strain) = split(/\t/,$line);
	$strains_of_gb{$gb} = $strain;
}
close(F);



my $cl_num = 0;
my $nb_strains = 1;
open(O,">$out");
open(U,">$out.upsetr.txt");
open(M,">$out.accessory_01matrix.txt");
open(F,$matrix);
my $firstline = <F>;
$firstline =~s/\n//g;$firstline =~s/\r//g;
my @infos = split(/\t/,$firstline);
print O "ClutserID";
print U "ClutserID";
print M "Gene";
for (my $j=1; $j <= $#infos; $j++){
        my $gbfile = $infos[$j];
        $gbfile =~s/\"//g;
        $gbfile =~s/\.gb\.filt//g;
	$gbfile =~s/\.gb\.rmdup//g;
        my $strain = $strains_of_gb{$gbfile};
        print O "\t".$strain;
        print U "\t".$strain;
        print M "\t".$strain;
        $nb_strains++;
}
print O "\n";
print U "\n";
print M "\n";
while(<F>){
        $cl_num++;
        my @infos = split(/\t/,$_);
        print O $cl_num;
        print U $cl_num;
        my $concat_accessory = "";
        for (my $i = 1; $i <= $#infos; $i++){
                my $val = $infos[$i];
                $val =~s/\"//g;
		$val =~s/\n//g;$val =~s/\r//g;
                if ($val =~/\w+/){
                        print U "\t1";
                        $concat_accessory .= "\t1";
                }
                else{
                        print U "\t0";
                        $concat_accessory .= "\t0";
                }
                my @genes = split(/, /,$val);
                my $concat = "";
                foreach my $gene(@genes){
                        my $prot_id = $gene;
                        $concat .= "$prot_id,"
                }
                chop($concat);
                if (scalar @genes == 0){
                        $concat = "-";
                }
                print O "\t".$concat;
        }
        if ($concat_accessory =~/0/){
                print M $cl_num.$concat_accessory."\n";
        }
        print O "\n";
        print U "\n";
}
close(F);
close(O);
close(U);
close(M);