diff --git a/genect_rnaseq/fpkm_subset_panel.pl b/genect_rnaseq/fpkm_subset_panel.pl index 9b611ec953bb256f0c831d2761c0bd31803c3ccc..4c37da57d943a8f4d9740700b24e011fc1af6b9b 100755 --- a/genect_rnaseq/fpkm_subset_panel.pl +++ b/genect_rnaseq/fpkm_subset_panel.pl @@ -9,14 +9,14 @@ my ($genelist,$fpkm_file) = ("") x 2; GetOptions( 'g|genes=s' => \$genelist, 'f|fpkm=s' => \$fpkm_file, - 'i|refinfo=s'=> \$refinfo); + 'i|geneinfo=s'=> \$refinfo); my $fpkm_out = $fpkm_file; $fpkm_out =~ s/\.fpkm.+txt/\.fpkm.capture.txt/; open OUT, ">$fpkm_out" or die $!; #gene_info.human.txt file is required for gene alias information -open GENEINFO, "<$refinfo/gene_info.human.txt" or die $!; +open GENEINFO, "<$geneinfo" or die $!; my %geneinfo; while(my $gline = <GENEINFO>){ chomp $gline; diff --git a/genect_rnaseq/geneabundance.sh b/genect_rnaseq/geneabundance.sh index 7eb6d3ff2b8a89cf0c77df26d9c37720bfb56eb1..8f6f6d7282efde84c5d139f9e0e6bad7f7851ebc 100644 --- a/genect_rnaseq/geneabundance.sh +++ b/genect_rnaseq/geneabundance.sh @@ -19,6 +19,7 @@ do p) pair_id=$OPTARG;; s) stranded=$OPTARG;; f) filter=$OPTARG;; + i) ginfo=$OPTARG;; h) usage;; esac done @@ -54,6 +55,6 @@ if [[ -f $filter ]] then cd .. mv ${pair_id}.fpkm.txt ${prefix}.fpkm.ori.txt - perl ${baseDir}/fpkm_subset_panel.pl -f ${prefix}.fpkm.ori.txt -g $filter + perl ${baseDir}/fpkm_subset_panel.pl -f ${prefix}.fpkm.ori.txt -g $filter -i $ginfo mv ${prefix}.fpkm.capture.txt ${pair_id}.fpkm.txt fi