Skip to content

Error on scglue.genomics.rna_anchored_guidance_graph #128

Description

@peiruzhang

Hello, I have the error when running guidance = scglue.genomics.rna_anchored_guidance_graph(rna, atac) when using the data provided in the tutorial. The error message is as follows:

ValueError Traceback (most recent call last)
Cell In[58], line 2
1 #print(rna.var)
----> 2 guidance = scglue.genomics.rna_anchored_guidance_graph(rna, atac)
3 #print(len(atac))
4 #guidance

File ~/anaconda3/envs/glue/lib/python3.10/site-packages/scglue/genomics.py:528, in rna_anchored_guidance_graph(rna, gene_region, promoter_len, extend_range, extend_fn, signs, propagate_highly_variable, corrupt_rate, random_state, *others)
526 raise RuntimeError("Length of others and signs must match!")
527 if set(signs).difference({-1, 1}):
--> 528 raise RuntimeError("signs can only contain {-1, 1}!")
530 rna_bed = Bed(rna.var.assign(name=rna.var_names))
531 other_beds = [Bed(other.var.assign(name=other.var_names)) for other in others]

File ~/anaconda3/envs/glue/lib/python3.10/site-packages/scglue/utils.py:496, in ConstrainedDataFrame.init(self, *args, **kwargs)
494 df = pd.DataFrame(*args, **kwargs)
495 df = self.rectify(df)
--> 496 self.verify(df)
497 super().init(df)

File ~/anaconda3/envs/glue/lib/python3.10/site-packages/scglue/genomics.py:65, in Bed.verify(cls, df)
63 super(Bed, cls).verify(df)
64 print("当前列名:", df.columns)
---> 65 print("数据内容示例:", df.head().to_csv(sep='\t', na_rep='nan'))
66 if len(df.columns) != len(cls.COLUMNS) or np.any(df.columns != cls.COLUMNS):
67 raise ValueError("Invalid BED format!")

ValueError: Invalid BED format!

When I deleted the duplicate columns in rna.var, it showed an error indicating that highly_variable was missing. Then I tried to remove the column name consistency check, and the following error message appeared.

Image

I will appreciate if you could help me to figure it out. Thanks!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions