How label Group list after integratin scRNA-seq and scATAC-seq data
0
0
Entering edit mode
4 days ago

Hello everyOne i have a question related to created Group list after integrating scRNA-seq data with scATAC-seq data we have identify totaly 6 clusters from both scRNA and scATAC-seq whats the logic behind this to generate group listenter image description here

groupList <- SimpleList(
    G1 = SimpleList(
        ATAC = Merged.proj2.SPG$cellNames[
            Merged.proj2.SPG$Clusters %in% c("C4")  
        ],
        RNA = colnames(SPG.seRNA)[
            SPG.seRNA$idents.rename2 %in% c("G1") 
        ]
    ),
    preL = SimpleList(
        ATAC = Merged.proj2.SPc$cellNames[
            Merged.proj2.SPc$Clusters %in% c("C3")  
        ],
        RNA = colnames(SPc.seRNA)[
            SPc.seRNA$idents.rename2 %in% c("preL") 
        ]
    ),
    ePL = SimpleList(
        ATAC = Merged.proj2.SPc$cellNames[
            Merged.proj2.SPc$Clusters %in% c("C5")  
        ],
        RNA = colnames(SPc.seRNA)[
            SPc.seRNA$idents.rename2 %in% c("ePL")
        ]
    ),
    IPL = SimpleList(
        ATAC = Merged.proj2.SPc$cellNames[
            Merged.proj2.SPc$Clusters %in% c("C2")  
        ],
        RNA = colnames(SPc.seRNA)[
            SPc.seRNA$idents.rename2 %in% c("IPL")
        ]
    )
    mPL = SimpleList(
        ATAC = Merged.proj2.SPc$cellNames[
            Merged.proj2.SPc$Clusters %in% c("C6")  
        ],
        RNA = colnames(SPc.seRNA)[
            SPc.seRNA$idents.rename2 %in% c("mPL")
        ]
    ),
    L = SimpleList(
        ATAC = Merged.proj2.SPc$cellNames[
            Merged.proj2.SPc$Clusters %in% c("C1")  
        ],
        RNA = colnames(SPc.seRNA)[
            SPc.seRNA$idents.rename2 %in% c("L")
        ]
    ),)
ArchR • 150 views
ADD COMMENT

Login before adding your answer.

Traffic: 1643 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6