Changeset 665 for plugins/base1/se.lu.thep.wenni/trunk/base/base1/base_plugin_script/wenni.pl

Timestamp:

Apr 17, 2008, 12:46:13 AM (15 years ago)

Author:

Jari Häkkinen

Message:

Added support for using configurable columns. Useful if BASE extra values are used for intensity and background SD.

File:

• plugins/base1/se.lu.thep.wenni/trunk/base/base1/base_plugin_script/wenni.pl (modified) (9 diffs)

plugins/base1/se.lu.thep.wenni/trunk/base/base1/base_plugin_script/wenni.pl

@@ -29 +29 @@
 sub basefile_options($$);
 sub BaseFileConverter($);
-sub cleanup($);
+sub cleanup($$);
 sub cmdline_options($);
 sub nni($);
-sub NNIFileConverter($$);
+sub NNIFileConverter($);
 sub result($$);
 
@@ -49 +49 @@
 # determine beta and number of neighbours. Presedence order is (low to
 # high priority): default value, basefile value, command line value.
-my %option=( beta => 0.6, datatype => "derived",
-             neighbours => 4, nodelete => 0 );
+my %option=(
+    beta => 0.6,
+    datatype => "derived",
+    intensity1 => "intensity1",
+    intensity2 => "intensity2",
+    neighbours => 4,
+    nodelete => 0,
+    BCh1SD => "BCh1SD",
+    BCh2SD => "BCh2SD"
+    );
 basefile_options(\%option,$basefile);
 foreach my $key (keys %cmdline_option) {
@@ -59 +67 @@
 }
 
-BaseFileConverter($option{'datatype'});
-NNIFileConverter($option{'beta'},$option{'datatype'});
+BaseFileConverter(\%option);
+NNIFileConverter(\%option);
 nni($option{'neighbours'});
-result($basefile,$option{'datatype'});
+result($basefile,\%option);
 if (!$option{'nodelete'}) {
-    cleanup($basefile);
+    cleanup($basefile,\%option);
 }
 
@@ -92 +100 @@
 
 sub BaseFileConverter($) {
-    my $datatype=$_[0];
+    my $option=$_[0];
     my $command="$FindBin::Bin/BaseFileConverter";
     $command.=" BASEfile.data";
     $command.=" wenni_";
-    if ($datatype eq "derived") {
-      $command.=" -assayFields intensity1";
-      $command.=" -assayFields intensity2";
+    if ($option{"datatype"} eq "derived") {
+      $command.=" -assayFields $$option{intensity1}";
+      $command.=" -assayFields $$option{intensity2}";
     }
     else {
@@ -106 +114 @@
       $command.=" -assayFields BCh2Mean";
     }
-    $command.=" -assayFields BCh1SD";
-    $command.=" -assayFields BCh2SD";
+    $command.=" -assayFields $$option{BCh1SD}";
+    $command.=" -assayFields $$option{BCh2SD}";
     return system("$command");
 }
 
-sub cleanup($) {
+sub cleanup($$) {
     my $basefile=$_[0];
+    my $option=$_[1];
     unlink($basefile);
     unlink("wenni_BCh1Mean.data");
-    unlink("wenni_BCh1SD.data");
+    unlink("wenni_$$option{BCh1SD}.data");
     unlink("wenni_BCh2Mean.data");
-    unlink("wenni_BCh2SD.data");
+    unlink("wenni_$$option{BCh2SD}.data");
     unlink("wenni_FCh1Mean.data");
     unlink("wenni_FCh2Mean.data");
     unlink("wenni_imputed.data");
-    unlink("wenni_intensity1.data");
-    unlink("wenni_intensity2.data");
+    unlink("wenni_$$option{intensity1}.data");
+    unlink("wenni_$$option{intensity2}.data");
     unlink("wenni_logratio.data");
     unlink("wenni_weight.data");
@@ -187 +196 @@
 }
 
-sub NNIFileConverter($$) {
-    my $beta=$_[0];
-    my $datatype=$_[1];
+sub NNIFileConverter($) {
+    my $option=$_[0];
     my $command="$FindBin::Bin/NNIFileConverter";
-    $command.=" -beta $beta";
-    if ($datatype eq "derived") {
-      $command.=" -fg1 wenni_intensity1.data";
-      $command.=" -fg2 wenni_intensity2.data";
+    $command.=" -beta $$option{'beta'}";
+    if ($$option{'datatype'} eq "derived") {
+      $command.=" -fg1 wenni_$$option{intensity1}.data";
+      $command.=" -fg2 wenni_$$option{intensity2}.data";
       $command.=" -datatype derived";
     }
@@ -204 +212 @@
         $command.=" -datatype raw";
     }
-    $command.=" -bgstd1 wenni_BCh1SD.data";
-    $command.=" -bgstd2 wenni_BCh2SD.data";
+    $command.=" -bgstd1 wenni_$$option{BCh1SD}.data";
+    $command.=" -bgstd2 wenni_$$option{BCh2SD}.data";
     $command.=" -logratio wenni_logratio.data";
     $command.=" -weight wenni_weight.data";
@@ -213 +221 @@
 sub result($$) {
     my $basefile=$_[0];
-    my $datatype=$_[1];
+    my $option=$_[1];
+    my $datatype=$$option{'datatype'};
     print "BASEfile\n";
     print "section\tspots\n";
@@ -255 +264 @@
     }
     else {
-        open(F1FILE,"<wenni_intensity1.data");
-        open(F2FILE,"<wenni_intensity2.data");
+        open(F1FILE,"<wenni_$$option{intensity1}.data");
+        open(F2FILE,"<wenni_$$option{intensity2}.data");
     }
     while (<IMPUTEFILE>) {