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

fix

parent 1541ad78
No related branches found
No related tags found
No related merge requests found
...@@ -92,9 +92,7 @@ mhcflurry-class1-train-presentation-models \ ...@@ -92,9 +92,7 @@ mhcflurry-class1-train-presentation-models \
--affinity-predictor "$(mhcflurry-downloads path models_class1_pan)/models.combined" \ --affinity-predictor "$(mhcflurry-downloads path models_class1_pan)/models.combined" \
--cleavage-predictor-with-flanks "$(mhcflurry-downloads path models_class1_cleavage)/models.selected" \ --cleavage-predictor-with-flanks "$(mhcflurry-downloads path models_class1_cleavage)/models.selected" \
--cleavage-predictor-without-flanks "$(mhcflurry-downloads path models_class1_cleavage_variants)/models.selected.no_flank" \ --cleavage-predictor-without-flanks "$(mhcflurry-downloads path models_class1_cleavage_variants)/models.selected.no_flank" \
--out-models-dir "$(pwd)/models" \ --out-models-dir "$(pwd)/models"
--worker-log-dir "$SCRATCH_DIR/$DOWNLOAD_NAME" \
$PARALLELISM_ARGS
cp $SCRIPT_ABSOLUTE_PATH . cp $SCRIPT_ABSOLUTE_PATH .
bzip2 -f "$LOG" bzip2 -f "$LOG"
......
...@@ -53,6 +53,7 @@ parser.add_argument( ...@@ -53,6 +53,7 @@ parser.add_argument(
def run(): def run():
args = parser.parse_args(sys.argv[1:]) args = parser.parse_args(sys.argv[1:])
hit_df = pandas.read_csv(args.hits) hit_df = pandas.read_csv(args.hits)
original_sample_ids = hit_df.sample_id.unique()
numpy.testing.assert_equal(hit_df.hit_id.nunique(), len(hit_df)) numpy.testing.assert_equal(hit_df.hit_id.nunique(), len(hit_df))
hit_df = hit_df.loc[ hit_df = hit_df.loc[
(hit_df.mhc_class == "I") & (hit_df.mhc_class == "I") &
...@@ -74,8 +75,7 @@ def run(): ...@@ -74,8 +75,7 @@ def run():
assert not (args.only_pmid and args.exclude_pmid) assert not (args.only_pmid and args.exclude_pmid)
pmids = list(args.only_pmid) + list(args.exclude_pmid) pmids = list(args.only_pmid) + list(args.exclude_pmid)
hit_pmids = hit_df.pmid.unique() missing = [pmid for pmid in pmids if pmid not in original_sample_ids]
missing = [pmid for pmid in pmids if pmid not in hit_pmids]
assert not missing, missing assert not missing, missing
mask = hit_df.pmid.isin(pmids) mask = hit_df.pmid.isin(pmids)
......
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