gene annotation for unpaired data #208
Unanswered
ruohuchengxhe
asked this question in
Q&A
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
Uh oh!
There was an error while loading. Please reload this page.
-
Hi,
I'm using unpaired data and currently processing public scATACseq data using pycisTopic. I use peak.bed as input and have the following files available:
barcodes.tsv
peak_annotation.tsv
peaks.bed
matrix.mtx
fragments.tsv.gz
I have also obtained tss.bed file as the tutorial suggested.

then created cistopic_obj using the following code:
what I have right now in the cistopic_obj are:
list(vars(cistopic_obj).keys())['fragment_matrix',
'binary_matrix',
'cell_names',
'region_names',
'cell_data',
'region_data',
'project',
'path_to_fragments',
'selected_model',
'projections',
'gene']
I manually refered peaks to get gene names.
I would like to know should I annotate the cells at this moment by skiping DARs and do gene activities?
Any help is appreciated!
Beta Was this translation helpful? Give feedback.
All reactions