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

[results] different input files or loom versions having different auccell loom results #544

Open
lexuziandyu opened this issue May 11, 2024 · 1 comment
Labels
question Further information is requested

Comments

@lexuziandyu
Copy link

lexuziandyu commented May 11, 2024

My version is:

pyscenic==0.12.1 loompy==2.0.17 anndata==0.9.2 scanpy==1.9.8

Compared with my previous pyscenic output file, the loom lacks the "MetaData" info.

Previous:

Class: H5Group
Filename: count.loom
Group: /attrs
Listing:
name obj_type dataset.dims dataset.type_class
CreationDate H5I_DATASET 0 H5T_STRING
LOOM_SPEC_VERSION H5I_DATASET 0 H5T_STRING
MetaData H5I_DATASET 0 H5T_STRING

Now:

Class: H5Group
Filename: celltype1.loom
Group: /attrs
Listing:
name obj_type dataset.dims dataset.type_class
LOOM_SPEC_VERSION H5I_DATASET 1 H5T_STRING

These two process are also different in: the previous one using a large SeuratObj, converted to a loom file, containing all celltypes, while the new one is subseted and converted to loom file only 1 celltype within.
The new output is a little bit different. When I perform

regulonThresholds = get_regulon_thresholds(loom)

It shows that:

Error in [[.H5Group(loom[["attrs"]], "MetaData") :
An object with name MetaData does not exist in this group

However, I can perform get_regulons_AUC and other functions on the Large aucellRseults properly. Is this possible that the input only-1-celltype loom file is too small for aucell step to compute the thresholds, or I need change my version to get the thresholds?

@lexuziandyu lexuziandyu added the question Further information is requested label May 11, 2024
@lexuziandyu
Copy link
Author

Update: after re-performing on the all-celltype-data, there is still a lack of "MetaData" slot. So it is believed that the problem is about the version.

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

No branches or pull requests

1 participant