extract specific gene expression data from the ABA RNA-Seq BrainSpan Matrix

extract.data.cmd

$ perl extract.data.pl GENESYMBL 2 |perl -ne '{chomp; $row++; @cols = split (/\t/, $_); $colNo = 0; foreach $data (@cols) {$colNo++;  $h{$colNo}{$row} = $data; } if (eof) { foreach $colNo (sort {$a<=>$b} keys %h) { $newRow =""; foreach $row (sort {$a<=>$b} keys %{$h{$colNo}} ) {$data = $h{$colNo}{$row}; $newRow .="$data\t"; } $newRow =~ s/\t$//; print "$newRow\n"; } } }' |sed 's/"//g' >gene.data/GENESYMBL.txt; paste sample.info.txt gene.data/GENESYMBL.txt >gene.data/GENESYMBL.data; rm gene.data/GENESYMBL.txt; cut -f4 gene.data/GENESYMBL.data >gene.data/GENESYMBL.val;