Skip to content

Commit

Permalink
coemstic
Browse files Browse the repository at this point in the history
  • Loading branch information
vbrover committed Mar 15, 2021
1 parent 708be6c commit 8b4329e
Show file tree
Hide file tree
Showing 2 changed files with 18 additions and 18 deletions.
12 changes: 6 additions & 6 deletions phylogeny/distTree_inc.sh
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ QC=1


if [ $QC == 1 ]; then
section "QC ..."
section "QC"
$INC/qc.sh go
fi

Expand Down Expand Up @@ -91,17 +91,17 @@ mv $INC/tree.new $INC/tree
#tail -n +5 leaf_errors.dm | sort -k 2 -g -r > leaf_errors.txt

if [ -e $INC/outlier-genogroup ]; then
section "Database: genogroup outliers ..."
section "Database: genogroup outliers"
$INC/objects_in_tree.sh $INC/outlier-genogroup null
mv $INC/outlier-genogroup $INC/hist/outlier-genogroup.$VER
fi


if [ $QC == 1 ]; then
section "QC ..."
section "QC"
$INC/qc.sh go
fi
section "Tree QC ..."
section "Tree QC"
$THIS/makeDistTree -threads 15 -data $INC/ -variance $VARIANCE -qc -noqual > $INC/hist/makeDistTree-qc.$VER
else
VER=`cat $INC/version`
Expand All @@ -125,7 +125,7 @@ if [ -e $INC/phen ]; then
OLD_ROOT=`grep '^Old root: ' $INC/hist/tree_quality_phen.$VER | sed 's/^Old root: //1'`
NEW_ROOT=`grep '^New root: ' $INC/hist/tree_quality_phen.$VER | sed 's/^New root: //1'`

section "Setting root and sorting ..."
section "Setting root and sorting"
if [ ! $NEW_ROOT ]; then
NEW_ROOT=$OLD_ROOT
fi
Expand All @@ -138,7 +138,7 @@ if [ -e $INC/phen ]; then

if [ -n "$RELDIR" ]; then
echo ""
section "Release ..."
section "Release"
if [ ! -e $RELDIR ]; then
error "$RELDIR does not exist"
fi
Expand Down
24 changes: 12 additions & 12 deletions phylogeny/distTree_inc_complete.sh
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,7 @@ OBJS=$2


#if false; then
section "QC $INC/ ..."
section "QC $INC/"
if [ -s $INC/tree ]; then
echo "$INC/tree must be empty"
exit 1
Expand All @@ -28,21 +28,21 @@ fi

$THIS/../sort.sh $OBJS

section "Computing dissimilarities ..."
section "Computing dissimilarities"
$THIS/../list2pairs $OBJS > $INC/dissim_request
$THIS/distTree_inc_request2dissim.sh $INC $INC/dissim_request $INC/dissim.raw
rm $INC/dissim_request
cat $INC/dissim.raw | grep -vwi "nan" | grep -vwi "inf" > $INC/dissim
rm $INC/dissim.raw

section "data.dm ..."
section "data.dm"
$THIS/../dm/pairs2dm $INC/dissim 1 "dissim" 6 -distance > $INC/../data.dm
echo "nan:"
set +o errexit
grep -wic 'nan' $INC/../data.dm
set -o errexit

section "Outliers ..."
section "Outliers"
$THIS/../dm/dm2objs $INC/../data -noprogress | sort > $INC/tree.list

$THIS/../setMinus $OBJS $INC/tree.list > $INC/outlier-alien
Expand All @@ -53,7 +53,7 @@ rm $INC/outlier-alien

HYBRIDNESS_MIN=`cat $INC/hybridness_min`
if [ $HYBRIDNESS_MIN != 0 ]; then
section "distTriangle ..."
section "distTriangle"
#cat data.dm | sed 's/nan/inf/g' > $INC/data1.dm
mkdir $INC/clust
$THIS/../dm/distTriangle $INC/../data "dissim" -clustering_dir $INC/clust -hybridness_min $HYBRIDNESS_MIN -hybrid $INC/hybrid.new
Expand All @@ -65,14 +65,14 @@ if [ $HYBRIDNESS_MIN != 0 ]; then
mv $INC/clust/1/data.dm $INC/../data.dm
rm -r $INC/clust/
$THIS/../dm/dm2objs $INC/../data | sort > $INC/tree.list
section "Hybrid ..."
section "Hybrid"
$THIS/distTree_inc_hybrid.sh $INC
#echo "Unhybrid ..."
#echo "Unhybrid"
#$THIS/distTree_inc_unhybrid.sh $INC
fi


section "Tree ..."
section "Tree"
HYBRID=""
if [ $HYBRIDNESS_MIN != 0 ]; then
DISSIM_BOUNDARY=`cat $INC/dissim_boundary`
Expand All @@ -81,12 +81,12 @@ fi
VARIANCE=`cat $INC/variance`
$THIS/makeDistTree -threads 5 -data $INC/../data -dissim_attr "dissim" -variance $VARIANCE -optimize -subgraph_iter_max 10 $HYBRID -output_tree $INC/tree > $INC/hist/makeDistTree-complete.1

section "Database ..."
section "Database"
$INC/objects_in_tree.sh $INC/tree.list 1
if [ $HYBRIDNESS_MIN != 0 ]; then
section "Hybrid ..."
section "Hybrid"
$THIS/distTree_inc_hybrid.sh $INC
#echo "Unhybrid ..."
#echo "Unhybrid"
#$THIS/distTree_inc_unhybrid.sh $INC
fi

Expand All @@ -96,7 +96,7 @@ cp $INC/tree $INC/hist/tree.1


if [ -e $INC/phen ]; then
section "Quality ..."
section "Quality"
LARGE=0
if [ -e $INC/large ]; then
LARGE=1
Expand Down

0 comments on commit 8b4329e

Please sign in to comment.