前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >单细胞数据拟时序分析-destiny

单细胞数据拟时序分析-destiny

作者头像
生信编程日常
发布2020-04-01 15:37:32
3.5K0
发布2020-04-01 15:37:32
举报
文章被收录于专栏:生物信息学、python、R、linux

单细胞数据分析常用到建立trajectory和pseudoTime,拟时序分析可以用 Diffusion( Destiny R package)

代码语言:javascript
复制
#Diffusion PseudoTime  Analysis
library(destiny) # 加载 destiny...
data(guo_norm) # 测试用的data
class(guo_norm)

image.png

代码语言:javascript
复制
dm <- DiffusionMap(ct,k = 3)
plot(dm)

image.png

给每个细胞添加注释信息,如这个细胞的类型或者属于的类群

代码语言:javascript
复制
palette(cube_helix(6)) #用cube_helix创建连续的颜色
#palette(hue_pal()(6))#也可以用ggplot2里面的默认颜色
plot(dm, pch = 20, # pch for prettier points
col_by = 'num_cells', # or “col” with a vector or one color
legend_main = 'Cell stage')

image.png

代码语言:javascript
复制
#2D plot
plot(dm, 1:2, pch = 20, col_by = 'num_cells',
legend_main = 'Cell stage')

image.png

代码语言:javascript
复制
#3D plot
library(rgl)
plot3d(eigenvectors(dm)[, 1:3],
col = log2(guo_norm$num_cells),
type = 's', radius = .01)
view3d(theta = 10, phi = 30, zoom = .8)
# now use your mouse to rotate the plot in the window
rgl.close()

image.png

同样可以用ggplot画出来

代码语言:javascript
复制
qplot(DC1, DC2, data = dm, colour = factor(num_cells)) +
scale_color_cube_helix()

image.png

代码语言:javascript
复制
# or alternatively:
dif<-fortify(dm)#转化为data.frame
ggplot(dif, aes(DC1, DC2, color = factor(num_cells)))+geom_point()

image.png

代码语言:javascript
复制
#plot 特征值
plot(eigenvalues(dm), ylim = 0:1, pch = 20,
xlab = 'Diffusion component (DC)', ylab = 'Eigenvalue')

image.png

detiny的数据输入格式为Biobase包建立的ExpressionSet格式的文件,如果我们的数据是表达矩阵,则数据需要转化成这个格式,如seurat包里面的数据Seurat.object可以这样转化:

代码语言:javascript
复制
library(Biobase)
ct <-GetAssayData(object = Seurat.object)
ct<-ct[VariableFeatures(Seurat.object),]
ct <- as.ExpressionSet(as.data.frame(t(ct)))
#添加注释信息
#. Annotations can be accessed directly via ct$column and ct[['column']]. 
ct$celltype <- DPT@meta.data[,c("integrated_merge_cluster")]
dm <- DiffusionMap(ct,k = 10)
palette(cube_helix(4)) # configure color palette
plot(dm, pch = 20, # pch for prettier points
col_by = "celltype")

image.png

参考: [http://10.30.30.253:8787/help/library/destiny/doc/Diffusion-Maps.pdf] [https://bioconductor.org/packages/release/bioc/vignettes/destiny/inst/doc/DPT.pdf] https://broadinstitute.github.io/2019_scWorkshop/functional-pseudotime-analysis.html

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档