如何定量、提取表达某基因的细胞（及比例）

一，提取表达某基因的所有细胞

expr <- sce_endo@assays$RNA
gene_expression <- expr %>% .['PLAT',] %>% as.data.frame() %>% t()
gene_expression <- as.data.frame(gene_expression)
colnames(gene_expression) <- 'PLAT'
gene_expression$cell <- rownames(gene_expression)
gene_expression_sel <- gene_expression[which(gene_expression$PLAT > 0),]
sce_select <- sce[,rownames(gene_expression_sel)]
dim(gene_expression_sel) / dim(gene_expression)

二，计算每个cluster中表达某基因的细胞比例

for (i in c(0:15)){rate = length(WhichCells(object = sce_endo,expression = VCAM1 > 0,idents = i)) /table(sce_endo$seurat_clusters)[[i+1]]print(rate)
}

如何定量、提取表达某基因的细胞（及比例）相关推荐

seurat提取表达矩阵_Hemberg-lab单细胞转录组数据分析
单细胞RNA-seq简介混合RNA-seq2000年末的重大技术突破,取代微阵列表达芯片被广泛使用通过混合大量细胞获取足够RNA用于建库测序,来定量每个基因的平均表达水平用于比较转录组,例如比较 ...
seurat提取表达矩阵_单细胞数据分析神器——Seurat
近年来,单细胞技术日益火热,并且有着愈演愈烈的趋势.在2015年至2017年,甚至对某细胞群体或组织进行单细胞测序,解析其细胞成分就能发一篇CNS级别的文章.近两三年,单细胞技术从最开始的基因组,转录 ...
seurat提取表达矩阵_Seurat
library(dplyr) library(Seurat) # 10X的数据可以使用Read10X这个函数,会返回一个UMI count矩阵,其中的每个值表示每个基因(行)在每个细胞(列)的分子数量 ...
seurat提取表达矩阵_10X scRNA免疫治疗学习笔记-3-走Seurat标准流程
刘小泽写于19.10.15 笔记目的:根据生信技能树的单细胞转录组课程探索10X Genomics技术相关的分析课程链接在:http://jm.grazy.cn/index/mulitcourse/ ...
寻找U2OS中表达的基因及其promoter并用于后续annotation
方法1.RNA-seq得到不同表达程度基因方法2. 直接download U2OS_gene.csv https://cancer.sanger.ac.uk/cell_lines/download ...
seurat提取表达矩阵_如何改造Seurat包的DoHeatmap函数？
刘小泽写于19.12.4 分析过单细胞数据的小伙伴应该都使用过Seurat包,其中有个函数叫DoHeatmap,具体操作可以看: 单细胞转录组学习笔记-17-用Seurat包分析文章数据前言走完S ...
MPB：林科院袁志林组-提取杨树人工林土壤微生物菌体细胞的4种方法
为进一步提高<微生物组实验手册>稿件质量,本项目新增大众评审环节.文章在通过同行评审后,采用公众号推送方式分享全文,任何人均可在线提交修改意见.公众号格式显示略有问题,建议电脑端点击文末阅 ...
人体肠道细菌与自身细胞的比例究竟是多少？
关于人体细菌与人体细胞的比例,有众多文章和学术报告反复引用1972年关于细菌数量是人体细胞10倍的结论来突出细菌的重要性.然而这一46年前的粗糙结论显然不符合最新研究的结果. 以色列魏茨曼科学研究所的 ...
seurat提取表达矩阵_本周最新文献速递20200719
本周最新文献速递一 (将不完全相似的表型进行荟萃分析,找出新的遗传位点,这是我第二次见到文献进行这种操作,一个新思路,值得借鉴) 文献题目: Multivariate genomic scan im ...

如何定量、提取表达某基因的细胞（及比例）

如何定量、提取表达某基因的细胞（及比例）相关推荐

最新文章

热门文章