Skip to content

Commit 4de5421

Browse files
committed
fix: convert identifiers to lowercase
Signed-off-by: Josh Loecker <joshloecker@icloud.com>
1 parent 666fb84 commit 4de5421

1 file changed

Lines changed: 2 additions & 3 deletions

File tree

main/como/rnaseq_preprocess.py

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -359,12 +359,12 @@ async def _write_counts_matrix(
359359
"""
360360
study_metrics = _organize_gene_counts_files(data_dir=como_context_dir)
361361
counts: list[pd.DataFrame] = await asyncio.gather(*[_create_sample_counts_matrix(metric) for metric in study_metrics])
362-
rna_specific_sample_names = set(config_df.loc[config_df["library_prep"] == rna.value, "sample_name"].tolist())
362+
rna_specific_sample_names = set(config_df.loc[config_df["library_prep"].str.lower() == rna.value.lower(), "sample_name"].tolist())
363363

364364
final_matrix: pd.DataFrame = functools.reduce(lambda left, right: pd.merge(left, right, on="ensembl_gene_id", how="outer"), counts)
365365
final_matrix.fillna(value=0, inplace=True)
366366
final_matrix.iloc[:, 1:] = final_matrix.iloc[:, 1:].astype(np.uint64)
367-
final_matrix = final_matrix[["ensembl_gene_id", *rna_specific_sample_names]]
367+
final_matrix = cast(pd.DataFrame, final_matrix[["ensembl_gene_id", *rna_specific_sample_names]])
368368

369369
output_counts_matrix_filepath.parent.mkdir(parents=True, exist_ok=True)
370370
final_matrix.to_csv(output_counts_matrix_filepath, index=False)
@@ -459,7 +459,6 @@ async def _create_config_df( # noqa: C901
459459

460460
fragment_label = f"{context_name}_{label}_fragment_size.txt"
461461
frag_paths = [p for p in aux_lookup["fragment"].values() if p.name == fragment_label]
462-
print(f"HERE: {prep}")
463462
if not frag_paths and prep.lower() != RNAType.TRNA.value.lower():
464463
logger.warning(f"No fragment file for '{label}'; defaulting to 100 bp (needed for zFPKM).")
465464
mean_frag = 100.0

0 commit comments

Comments
 (0)