Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Fixed a wrong reference to a bed file
  • Loading branch information
jenzopr committed Oct 30, 2015
1 parent c3e0bd6 commit 7dc9a05
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/admire
Expand Up @@ -436,7 +436,7 @@ for a; do
bedtools sort -i $RESULTS/$N/$g.pvals.bed > $RESULTS/$N/$g.pvals.sorted.bed
fi
bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$g.pvals.passed_fdr_$Q.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.passed_fdr_$Q.abs.bed
bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$p.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.all.abs.bed
bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$g.pvals.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.all.abs.bed
awk -v g="$g" 'BEGIN{FS="\t";OFS=",";split(g,c,"-");print "Feature","Gene","GeneID","Type","Chr","Start","Stop","Strand","P-val "c[1],"Q-val "c[1],"P-val "c[2],"Q-val "c[2],"Higher in","with Q-val","Coverage","Min.abs.diff","Median.abs.diff","Max.abs.diff";}
{n = split($4,a,";");for(i=1; i<=n; i++){
split(a[i],b,"=");
Expand Down

0 comments on commit 7dc9a05

Please sign in to comment.