# 导出SingleR需要的数据

    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12

    tp_samples <- list.files('~/GEO/GSE193337')
    tp_dir <- file.path('~/GEO/GSE193337', tp_samples)
    names(tp_dir) <- tp_samples
    counts <- Seurat::Read10X(data.dir = tp_dir)
    sce <- Seurat::CreateSeuratObject(counts, project = 'prostate',
                                min.cells = 3, min.features = 200)
    rm(counts)
    gc()
    table(sce$orig.ident)
    sce <- Seurat::NormalizeData(sce)
    logumi <- Seurat::GetAssayData(sce, slot="data")
    saveRDS(logumi, 'logumi.rds')

注释
导出SingleR需要的数据
https://occdn.limour.top/2394.html
Author
Limour
Posted on
October 4, 2022
Licensed under
【转载】货拉拉案当事人对《底线》“帮你拉案”的声明 Previous
构造Pseudo-bulk的表达量矩阵加速SingleR运行 Next

#### Search

Keyword