MJさむらいの実験室

この間も書いた海外生活で重宝するインスタント味噌汁なのであるが、

またやってしまった。

味噌と具が別々の袋に入っているちょっと高級なインスタント味噌汁を買ってきたのだが、

具、味噌と入れないといけないところを

具、具と入れてしまった。

味噌を付け足さないと飲めないので

具、具、味噌汁というコンビネーションが出来上がった。


長ネギとアゲの超豪華味噌汁だ。


今回はこれでよいのだが、いつのひか味噌だけを飲まないといけない日が来る。

ブラストクラストというプログラムを使って、ＥＳＴデーターなどのたくさんの配列の冗長を取り除く事ができます。

目で確認できるほどなら問題ないのですが、データーがたくさんあると数えるのが大変なのでパールでプログラムを作ってみました。

このプログラムを使うと、ブラストクラストの吐き出すデーターのクラスターされた数をカウントできます。

僕の場合は１２０９３７の配列をクラスタリングされてできた結果が６７６９４にクラスタリングされました。

下に張ったプログラムを実行すると

ぶわーっと各クラスターに何個のシークエンスがクラスタリングされたかを吐き出した後に

２０個以上クラスタリングされたクラスターの数、
１９個から６個クラスタリングされたクラスターの数、
５個クラスタリングされたクラスターの数、
４個クラスタリングされたクラスターの数、
３個クラスタリングされたクラスターの数、
２個クラスタリングされたクラスターの数、
１個（クラスタリングされなかったもの）クラスタリングされたクラスターの数、
を吐き出します。

クラスターの数は最後の方にプリントされます。

20= 532
19-5>1069
5=253
4=349
3=571
2=1035
1=63885
67694
total 67694

こんな感じです。
これは例えば３つのシークエンスがクラスタリングされたものが５７１あるということを表しています。
同様に１９から５よりも上の数のシークエンスがクラスタリングされたものは１０６９あるという意味です。

こういう処理をやろうとしている人にしか
このプログラムの意味はわからないと思いますが丁寧に説明するのは面倒なのであんまり説明しません。

きっとこんな事をしている人たちはバイオインフォマジシャン達なので
こんなしょぼいプログラムは必要ないでしょう。

一応、独学で奮闘している人にとって参考になるかなあと思ってのっけておきます。
（僕も独学ですので適当なプログラミングです。）


プログラムはここから。

#!/usr/bin/perl -w
use strict;
use warnings;
use POSIX;

my $counter_20 = 0;
my $counter_5_20= 0;
my $counter_5= 0;
my $counter_4= 0;
my $counter_3= 0;
my $counter_2= 0;
my $counter_1= 0;
my $counter= 0;

my $name = <STDIN>;

open (FASTANAME, $name);

while (<FASTANAME>)
{

my @filename = split( " ", $_);
 
my $clusternumber = @filename;
print $clusternumber, "\n";
if ($clusternumber == 1)
  {$counter_1++
}elsif ($clusternumber == 2)
  {$counter_2++
}elsif ($clusternumber == 3)
  {$counter_3++
}elsif ($clusternumber == 4)
  {$counter_4++
}elsif ($clusternumber == 5)
  {$counter_5++
}elsif ($clusternumber >= 20)
  {$counter_20++
}else
  {$counter_5_20++
}
$counter++;
}
close (FASTANAME);

print "20\= $counter_20\n";
print '19>5',"$counter_5_20\n";
print "5\=$counter_5\n";
print "4\=$counter_4\n";
print "3\=$counter_3\n";
print "2\=$counter_2\n";
print "1\=$counter_1\n";
print "$counter\n";
my $total = $counter_1 + $counter_2 + $counter_3 + $counter_4 + $counter_5 + $counter_5_20 + $counter_20;
print "total $total\n";

こちらでの仕事ではないのですが、永らくくすぶっていた仕事をようやく投稿にまでこぎつけました。

まあ、だらだらしていた自分が悪いのですが、ここまで来たらさっさとけりをつけたいところです。

すんなりと終わる事を期待しています。

ＦＩＧのフォーマットやらを直すのにえらく時間がかかってしまって疲れました。