...都是准确的
这个代码不报错
for i in {2..6};do Rscript $scriptdir/meta_stat.r -i ../05.export_data/taxa_summary/absolute_abundance/feature_table_tax_L$i.txt \ -m $metadata -t Group \ -a WSD_Lig_PPI -b WSD_PPI \ -o meta_stat --name WSD_Ligatur_vsKD_Ligatur_Lppi$idone#...
use strict;
use warnings;
print STDERR "gene symbol column number: ";
my $geneSymbolCol=<STDIN>;
chomp($geneSymbolCol);
$geneSymbolCol--;
my $expFile="probeMatrix.txt";
my $gplFile="ann.txt";
my $expFileWF="geneMatrix.txt";
my %hash=();
my @sampleName=();
open(EXP,"$expFile") or d...