File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change 11PBMC10k\_ multiome data prepare for cellxgene VIP
22================
33
4+ To generate a multiome instance in cellxgene VIP, three items are required. You can generate them by following this notebook.
5+
6+ + One h5ad file. Please follows ** part 1** .
7+ + Three .rds file (annotation.rds, links.rds, peaks.rds). Please follows ** part 1** .
8+ + bigwig files (.bw files) and one index file (named bw.cluster). Please follows ** part 2** .
9+
10+
11+ ## Part 1
412Following R code is based on
513< https://satijalab.org/signac/articles/pbmc_multiomic.html >
614
@@ -252,7 +260,7 @@ saveRDS(Links(pbmc[["peaks"]], file='links.rds'))
252260SaveH5Seurat(pbmc , filename = ' pbmc10k_multiome.h5Seurat' )
253261Convert(' pbmc10k_multiome.h5Seurat' , dest = ' h5ad' )
254262```
255-
263+ ## Part 2
256264Generate celltype.txt for bam splitting (In Python)
257265
258266``` python
@@ -332,3 +340,11 @@ ls *bw | sed -e 's/\.bam$//' > bw.cluster_c2
332340sed -i ' s/./ /g' bw.cluster_c2
333341paste bw.cluster_c1 bw.cluster_c2 > bw.cluster
334342```
343+ Example output of bw.cluster:
344+ CellType | |
345+ ------------- | -------------
346+ ASDC.bw | ASDC
347+ B.intermediate.bw | B intermediate
348+ B.memory.bw | B memory
349+ B.naive.bw | B naive
350+
You can’t perform that action at this time.
0 commit comments