Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

is de_genes keyed on condition only or cell_drug_dose_comb #144

Open
bhomass opened this issue Sep 12, 2023 · 1 comment
Open

is de_genes keyed on condition only or cell_drug_dose_comb #144

bhomass opened this issue Sep 12, 2023 · 1 comment

Comments

@bhomass
Copy link

bhomass commented Sep 12, 2023

in train.py, it says de_genes is keyed on cell_drug_dose_comb

        # genes for every cell_drug_dose combination.
        bool_de = dataset.var_names.isin(
            np.array(dataset.de_genes[cell_drug_dose_comb])
        )

But, in lincs.py, de_genes is keyed on eval_category
adata.uns['rank_genes_groups_cov'] = {cat: de_genes_quick[extract_drug(cat)] for cat in adata.obs.eval_category.unique() if extract_drug(cat) != 'DMSO'}

where eval_category is
adata.obs["eval_category"] = adata.obs["cov_drug_name"]
so, no dosage included

checking lincs_full_smiles_sciplex_genes.h5ad though, indeed it is keyed on cell-line_drug.

which convention did you ultimately choose? without code modification in train.py, the code crashes.

@bhomass
Copy link
Author

bhomass commented Sep 15, 2023

when it comes to sciplex_complete_middle_subset_lincs_genes.h5ad, we are back to cell-line_drug_dose again.

There has to be a set convention, or else the code in train.py can not work for all datasets.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant