Skip to content
Snippets Groups Projects
Commit 86e04962 authored by Tim O'Donnell's avatar Tim O'Donnell
Browse files

update

parent 1c47306e
No related branches found
No related tags found
No related merge requests found
...@@ -89,11 +89,16 @@ fi ...@@ -89,11 +89,16 @@ fi
# ******************************************************** # ********************************************************
echo "Beginning production run" echo "Beginning production run"
time python make_multiallelic_training_data.py \ if [ -f "$SCRIPT_DIR/train.multiallelic.csv" ]; then
--hits "$(mhcflurry-downloads path data_mass_spec_annotated)/annotated_ms.csv.bz2" \ echo "Using existing multiallelic train data."
--expression "$(mhcflurry-downloads path data_curated)/rna_expression.csv.bz2" \ cp "$SCRIPT_DIR/train.multiallelic.csv" .
--decoys-per-hit 1 \ else
--out train.multiallelic.csv time python make_multiallelic_training_data.py \
--hits "$(mhcflurry-downloads path data_mass_spec_annotated)/annotated_ms.csv.bz2" \
--expression "$(mhcflurry-downloads path data_curated)/rna_expression.csv.bz2" \
--decoys-per-hit 1 \
--out train.multiallelic.csv
fi
ALLELE_LIST=$(bzcat "$MONOALLELIC_TRAIN" | cut -f 1 -d , | grep -v allele | uniq | sort | uniq) ALLELE_LIST=$(bzcat "$MONOALLELIC_TRAIN" | cut -f 1 -d , | grep -v allele | uniq | sort | uniq)
ALLELE_LIST+=$(cat train.multiallelic.csv | cut -f 7 -d , | gerp -v hla | uniq | tr ' ' '\n' | sort | uniq) ALLELE_LIST+=$(cat train.multiallelic.csv | cut -f 7 -d , | gerp -v hla | uniq | tr ' ' '\n' | sort | uniq)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment