diff --git a/PostProcess/add_risk_score.py b/PostProcess/add_risk_score.py index 125d7e2..5f88505 100755 --- a/PostProcess/add_risk_score.py +++ b/PostProcess/add_risk_score.py @@ -82,7 +82,7 @@ def compute_risk_score(report_fname: str, report_outfname: str, alternative: boo fout.write("\t".join(header) + "\n") # write header to out put report for line in fin: score_diff, score_diff_abs = _compute_risk_score(line) - fout.write(f"{line}\t{score_diff}\t{score_diff_abs}\n") + fout.write(f"{line.strip()}\t{score_diff}\t{score_diff_abs}\n") except (IOError, Exception) as e: raise OSError( f"An error occurred while computing risk scores for {report_fname}" diff --git a/PostProcess/remove_n_and_dots.py b/PostProcess/remove_n_and_dots.py index 1e26b89..04a917f 100755 --- a/PostProcess/remove_n_and_dots.py +++ b/PostProcess/remove_n_and_dots.py @@ -52,10 +52,9 @@ def polish_report(report_chunks: TextFileReader, report_fname: str) -> None: header = True # only for first iteration for chunk in report_chunks: assert "rsID" in chunk.columns.tolist() - chunk: pd.DataFrame = chunk.replace("n", "NA") # replace ns with NAs - chunk["rsID"] = chunk["rsID"].str.replace( - ".", "NA" - ) # replace . in rsids with NAs + chunk: pd.DataFrame = chunk.replace({"n": "NA"}) # replace ns with NAs + # replace . in rsids with NAs + chunk["rsID"] = chunk["rsID"].str.replace(".", "NA", regex=False) chunk.to_csv( f"{report_fname}.tmp", header=header, mode="a", sep="\t", index=False ) @@ -79,7 +78,7 @@ def remove_n_dots() -> None: polish_report(read_report_chunks(report_fname), report_fname) if not os.path.isfile(f"{report_fname}.tmp"): raise FileNotFoundError(f"Polished report {report_fname}.tmp not created?") - code = subprocess.call(f"mv {report_fname}.tmp {report_fname}") + code = subprocess.call(f"mv {report_fname}.tmp {report_fname}", shell=True) if code != 0: raise subprocess.SubprocessError(f"Failed renaming {report_fname}.tmp")