Skip to content

Commit 971ac82

Browse files
committed
Fix report polishing bug
1 parent 7646229 commit 971ac82

File tree

2 files changed

+5
-6
lines changed

2 files changed

+5
-6
lines changed

PostProcess/add_risk_score.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -82,7 +82,7 @@ def compute_risk_score(report_fname: str, report_outfname: str, alternative: boo
8282
fout.write("\t".join(header) + "\n") # write header to out put report
8383
for line in fin:
8484
score_diff, score_diff_abs = _compute_risk_score(line)
85-
fout.write(f"{line}\t{score_diff}\t{score_diff_abs}\n")
85+
fout.write(f"{line.strip()}\t{score_diff}\t{score_diff_abs}\n")
8686
except (IOError, Exception) as e:
8787
raise OSError(
8888
f"An error occurred while computing risk scores for {report_fname}"

PostProcess/remove_n_and_dots.py

Lines changed: 4 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -52,10 +52,9 @@ def polish_report(report_chunks: TextFileReader, report_fname: str) -> None:
5252
header = True # only for first iteration
5353
for chunk in report_chunks:
5454
assert "rsID" in chunk.columns.tolist()
55-
chunk: pd.DataFrame = chunk.replace("n", "NA") # replace ns with NAs
56-
chunk["rsID"] = chunk["rsID"].str.replace(
57-
".", "NA"
58-
) # replace . in rsids with NAs
55+
chunk: pd.DataFrame = chunk.replace({"n": "NA"}) # replace ns with NAs
56+
# replace . in rsids with NAs
57+
chunk["rsID"] = chunk["rsID"].str.replace(".", "NA", regex=False)
5958
chunk.to_csv(
6059
f"{report_fname}.tmp", header=header, mode="a", sep="\t", index=False
6160
)
@@ -79,7 +78,7 @@ def remove_n_dots() -> None:
7978
polish_report(read_report_chunks(report_fname), report_fname)
8079
if not os.path.isfile(f"{report_fname}.tmp"):
8180
raise FileNotFoundError(f"Polished report {report_fname}.tmp not created?")
82-
code = subprocess.call(f"mv {report_fname}.tmp {report_fname}")
81+
code = subprocess.call(f"mv {report_fname}.tmp {report_fname}", shell=True)
8382
if code != 0:
8483
raise subprocess.SubprocessError(f"Failed renaming {report_fname}.tmp")
8584

0 commit comments

Comments
 (0)