Skip to content

Commit

Permalink
Fix omitted variants, correct test cases. (#84)
Browse files Browse the repository at this point in the history
* addert assert to make sure no variants are omitted when merging annotations

* changed expected files, sorts input files in aggregate_abscores

* fixup! Format Python code with psf/black pull_request

---------

Co-authored-by: “Marcel-Mueck” <“mueckm1@gmail.com”>
Co-authored-by: PMBio <PMBio@users.noreply.github.com>
  • Loading branch information
3 people authored May 14, 2024
1 parent 6f04801 commit 5715e99
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion deeprvat/annotations/annotations.py
Original file line number Diff line number Diff line change
Expand Up @@ -1082,7 +1082,7 @@ def aggregate_abscores(
delayed(process_chunk)(
i, abs_splice_res_dir, tissues_to_exclude, tissue_agg_function, ca_shortened
)
for i in tqdm(os.listdir(abs_splice_res_dir))
for i in tqdm(sorted(os.listdir(abs_splice_res_dir)))
)
all_absplice_scores = list(output_generator)

Expand Down Expand Up @@ -1751,6 +1751,10 @@ def process_vep(
assert vep_file.pos.isna().sum() == 0
assert vep_file.ref.isna().sum() == 0
assert vep_file.alt.isna().sum() == 0
assert (
vep_file[["chrom", "pos", "ref", "alt"]].drop_duplicates().shape
== vcf_df[["chrom", "pos", "ref", "alt"]].drop_duplicates().shape
)

if "pos" in vep_file.columns:
vep_file["pos"] = vep_file["pos"].astype(int)
Expand Down
Binary file not shown.
Binary file not shown.

0 comments on commit 5715e99

Please sign in to comment.