专栏首页生物信息学、python、R、linuxseurat包分析多组对比单细胞数据

seurat包分析多组对比单细胞数据

library(Seurat)
#import data
#C_data T_data 为要分析的data.frame

Control<-CreateSeuratObject(counts =C_data,min.cells = 5, min.features = 10,project = "control")
Treat<-CreateSeuratObject(counts =T_data,min.cells = 5, min.features = 10,project = "treat")
#将多个数据合成一个list
T_C<-list(Treat,Control)
names(T_C)<-c("T","C")
#分别对每个找Variable Features
for (i in 1:length(T_C)){
  T_C[[i]]<-FindVariableFeatures(T_C[[i]], selection.method = "vst", 
                       nfeatures = 2000, verbose = FALSE)
}
#找到交集的feature
T_C<- FindIntegrationAnchors(object.list = T_C, dims = 1:20)
#整合数据
T_C <- IntegrateData(anchorset = T_C, dims = 1:20)
DefaultAssay(T_C) <- "integrated"
# Run the standard workflow for visualization and clustering
T_C <- ScaleData(T_C, verbose = FALSE)
T_C <- RunPCA(T_C, npcs = 30, verbose = FALSE)
# t-SNE and Clustering
T_C <- RunUMAP(T_C, reduction = "pca", dims = 1:20)
T_C <- FindNeighbors(T_C, reduction = "pca", dims = 1:20)
T_C <- FindClusters(T_C, resolution = 0.5)
# Visualization
p1 <- DimPlot(T_C, reduction = "umap", group.by = "stim")
p2 <- DimPlot(T_C, reduction = "umap", label = TRUE)
plot_grid(p1, p2)

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • pandas的一些小知识

    生信编程日常
  • R plot / ggplot2加透明度

    2.plot修改 plot画图的时候,直接修改col就可以,一个比较方便的方法是用rgb()函数。这个函数提供三个值,默认情况下是0到1的值,可以修改maxC...

    生信编程日常
  • Seaborn入门(一): 实现Boxplot

    与之相关的函数seaborn.boxplot的具体参数为: seaborn.boxplot(x=None, y=None, hue=None, data=No...

    生信编程日常
  • 【融职培训】Web前端学习 第3章 JavaScript基础教程4 条件语句

    例如:我们登录一个软件,输入用户名和密码的时候,输入正确则登录成功,输入错误,会弹出相关提示,这就是一个典型的分支语句(在实际开发中,这个功能或设计多个分支语句...

    学习猿地
  • Web前端学习 第3章 JavaScript基础教程4 条件语句

    例如:我们登录一个软件,输入用户名和密码的时候,输入正确则登录成功,输入错误,会弹出相关提示,这就是一个典型的分支语句(在实际开发中,这个功能或设计多个分支语句...

    学习猿地
  • 【译】Understanding SOLID Principles - Single Responsibility

    之前的第一篇文章阐述了依赖倒置原则(DIP)能够使我们编写的代码变得低耦合,同时具有很好的可测试性,接下来我们来简单了解下单一职责原则的基本概念:

    littlelyon
  • 七夕情人节,看 ---大数据时代里的爱情!

    从前,在西雅图的一家Pony Expresso咖啡店里,一个男人与一个女人开始了对这个绵长而又神秘的事物的体验,这个事物已得到了愈来愈多科学研究,而我们称其为爱...

    小莹莹
  • 小程序之保存图片到相册的几种方法

    honey缘木鱼
  • centos6下ActiveMQ+Zookeeper消息中间件集群部署记录

    由于最近一个项目并发请求压力比较大,所以考虑改进架构,引入消息中间件集群作为一个缓冲消息队列,具体需求: 1)将大量的WebService请求报文发送到mq集群...

    洗尽了浮华
  • task8 GAN text-to-image

    https://www.tensorflow.org/api_docs/python/tf/layers/batch_normalization https:...

    平凡的学生族

扫码关注云+社区

领取腾讯云代金券