Skip to content

Commit

Permalink
enable filter peak by qvalue in integrate module
Browse files Browse the repository at this point in the history
  • Loading branch information
wbaopaul committed May 1, 2022
1 parent e3d05c9 commit 4564c50
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions scripts/integrate_peak.sh
Original file line number Diff line number Diff line change
Expand Up @@ -35,11 +35,13 @@ ABS_PATH=`cd "$OUTPUT_DIR"; pwd`
## reconstruct the peak-cell matrix
## not re-call cells
mtx_files='TMP'
re='^[0-9]+([.][0-9]+)?$'
for pk0 in "${peaks[@]}"
do
if [ "$pk0" -eq "$pk0" ] 2>/dev/null;then
continue ## the gap parameter
if [[ $pk0 =~ $re ]] 2>/dev/null;then
continue ## the gap and qvalue parameter
fi
echo "Reconstruct matrix for sample related to $pk0: "
sample0=$(basename $pk0)
sample0=`echo $sample0 | awk -F. '{print $1}'`
sample0=${sample0/_features_BlacklistRemoved/}
Expand Down

0 comments on commit 4564c50

Please sign in to comment.