seurat 单细胞数据分析中 VizDimLoadings 函数
2022/8/27 23:25:55
本文主要是介绍seurat 单细胞数据分析中 VizDimLoadings 函数,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
前期处理:https://www.jianshu.com/p/fef17a1babc2
#可视化对每个主成分影响比较大的基因集
001、
dat <- pbmc[["pca"]]@feature.loadings ## 数据来源 dat[1:3, 1:3] dat <- dat[order(-dat[,1]),][1:29,1] dat <- as.data.frame(dat) dat <- cbind(gene = rownames(dat), dat) dat <- rbind(dat, MALAT1 = c("MALAT1", -0.10)) ## ???? dat$num <- nrow(dat):1 par(mai = c(1, 1, 1, 1),mgp = c(2.5,0.7,0)) plot(dat$dat,pch = 19, dat$num, yaxt ="n", xlab = "PC_1", ylab = "") axis(2, 1:30, rev(dat$gene), col.axis = "black", las = 2) ## 绘图
标准结果:
VizDimLoadings(pbmc, dims = 1, reduction = "pca")
这篇关于seurat 单细胞数据分析中 VizDimLoadings 函数的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-15鸿蒙生态设备数量超8亿台
- 2024-05-13TiDB + ES:转转业财系统亿级数据存储优化实践
- 2024-05-09“2024鸿蒙零基础快速实战-仿抖音App开发(ArkTS版)”实战课程已上线
- 2024-05-09聊聊如何通过arthas-tunnel-server来远程管理所有需要arthas监控的应用
- 2024-05-09log4j2这么配就对了
- 2024-05-09nginx修改Content-Type
- 2024-05-09Redis多数据源,看这篇就够了
- 2024-05-09Google Chrome驱动程序 124.0.6367.62(正式版本)去哪下载?
- 2024-05-09有没有大佬知道这种数据应该怎么抓取呀?
- 2024-05-09这种运行结果里的10.100000001,怎么能最快改成10.1?