美文网首页
2020-04-17

2020-04-17

作者: 吹不散的烽烟 | 来源:发表于2020-04-17 23:04 被阅读0次

自学单细胞笔记

认识sce object

[参考文献]https://www.researchgate.net/publication/332045971_Orchestrating_Single-Cell_Analysis_with_Bioconductor
image.png

(图片引自Robert A. Amezquita et al, 2019 )

学习sce Object

环境变量设置

rm(list = ls())
library(BiocManager)
# BiocManager::install("SingleCellExperiment")
# BiocManager::install("scran")
# BiocManager::install("scater")
library(SingleCellExperiment)
library(scran)
library(scater)

构建一个表达矩阵

counts_matrix <- data.frame(cell_1 = rpois(10, 10),
                            cell_2 = rpois(10, 10),
                            cell_3 = rpois(10, 30))
rownames(counts_matrix) <- paste0("gene_", 1:10)
counts_matrix <- as.matrix(counts_matrix) # must be a matrix object!

构建sce对象

sce <- SingleCellExperiment(assays = list(counts = counts_matrix))
class(sce)
sce

取出表达矩阵counts

assay(sce,"counts") #方法1
counts(sce) #方法2
#metadata(1): log.exprs.offset
#assays(2): counts logcounts

添加矩阵

sce <- scran::computeSumFactors(sce)
sce
sce <- scater::normalize(sce)
sce

提取添加的矩阵

assay(sce,"logcounts") #方法1
logcounts(sce)  #方法2
assays(sce) ### 查看assays列表

添加矩阵counts_100

counts_100 <- assay(sce, "counts") + 100
counts_100
assay(sce, "counts_100") <- counts_100 # assign a new entry to assays slot
sce
assays(sce)

构建sample metadata:batch

cell_metadata <- data.frame(batch = c(1, 1, 2))
rownames(cell_metadata) <- paste0("cell_", 1:3) ##paste和paste0区别sep="是否有空格"
## From scratch:
sce <- SingleCellExperiment(assays = list(counts = counts_matrix),
                            colData = cell_metadata) ##添加sample metadata: batch effect通过colData
## Appending to existing object (requires DataFrame() coercion)
## colData(sce) <- DataFrame(cell_metadata)
sce
colData(sce) #获取batch effect

利用scater添加calculateQCMetrics 查看colData

sce <- scater::calculateQCMetrics(sce)
colData(sce)[, 1:5]
sce$batch
## colData(sce)$batch # same as above

构建子集

sce[, sce$batch == 1] ##构建子集batch=1

添加feature metadata

rowRanges(sce)  # empty
rowData(sce)[, 1:3]
sce[c("gene_1", "gene_4"), ]
## sce[c(1, 4), ] # same as above in this case

添加sizeFactors

sce <- scran::computeSumFactors(sce)
sce <- scater::normalize(sce)
sizeFactors(sce)

添加reducedDim:PCA,TSNE,UMAP

sce <- scater::runPCA(sce)#降维分析
reducedDim(sce, "PCA")
sce <- scater::runTSNE(sce, perplexity = 0.1)
reducedDim(sce, "TSNE")
reducedDims(sce)
u <- uwot::umap(t(logcounts(sce)), n_neighbors = 2)
reducedDim(sce, "UMAP_uwot") <- u
reducedDim(sce, "UMAP_uwot")
reducedDims(sce)

添加others metadata

my_genes <- c("gene_1", "gene_5")
metadata(sce) <- list(favorite_genes = my_genes)
metadata(sce)
your_genes <- c("gene_4", "gene_8")
metadata(sce)$your_genes <- your_genes
metadata(sce)

相关文章

  • 2020-04-19

    2020-04-17 菜菜_d868 字数 348 · 阅读 0 2020-04-17 22:06 姓名:邢彩颜 ...

  • 2020-04-18

    2020-04-17 2020-04-17 成长日志第303天 家名:温暖有爱之家 家规:真诚待人不自欺,学会看见...

  • 【日更4】读书笔记,一起学习,一起进步

    《这样读书就够了》读书笔记(四) 2020-04-17 今天将...

  • 2020-07-13 再不学会,我就老了

    http://tech.sina.com.cn/csj/2020-04-17/doc-iirczymi684830...

  • 2021-10-11

    iOS布局-autoresizingMask Code&Ocean 2020-04-17 14:56:58 22...

  • 学习小组54期 1st day--忍冬

    2020-04-17 今天的内容: https://m.umu.cn/session/article/29dNxa...

  • 晚间日记4

    文|董大壳 今天2020-04-17 起床:06:00 就寝:24:00 天气:大风 心情:一般般 任务清单 今天...

  • 复制粘贴

    2020-04-17 抢购的活动 !!! 继续尝试海报 我的早餐推荐吧!!!(长图制作) 选品: 玩具的留到6/1...

  • 2020-04-19

    再跑 默蝶一 字数 127 · 阅读 1 2020-04-17 22:44 跑步连续打卡第四天了吧?自己无法有规定...

  • 人 育

    时间:2011年9月21日(星期三) 中午1:27纯文本|2020-04-17 收件人: 黑龙江-13校园文化惠自...

网友评论

      本文标题:2020-04-17

      本文链接:https://www.haomeiwen.com/subject/bukpvhtx.html