ggplot画 ump 和tsne 从seurat中使用addmodule得到的umap 使用ggplot画图
ggplot画 ump 和tsne 从seurat中使用addmodule得到的umap 使用ggplot画图 subset_data
inflammatory_gene=read.xlsx(“G:/silicosis/sicosis/silicosis_ST/yll/0214/harmony_cluster/HALLMARK_INFLAMMATORY_RESPONSE.xlsx”)
head(inflammatory_gene)
inflammatory_gene=capitalize(tolower(inflammatory_gene$gene_symbol))
inflammatory_gene=list(inflammatory_gene)
load(file = “G:/silicosis/sicosis/yll/macrophage/no cluster2/macrophage_clean.rds”)
#对给定的基因集合进行打分
subset_data=AddModuleScore(subset_data,
features = inflammatory_gene,
name = “inflammatory_gene”)
Warning: The following features are not present in the object: Ccl20, Cxcl11, Cxcl6, Icoslg, Il4r, Cxcl8, Ros1, not searching for symbol synonyms
#结果保存在这里
colnames(subset_data@meta.data)
mydata<- FetchData(subset_data,vars = c(“UMAP_1”,“UMAP_2”,“inflammatory_gene1”,“stim”))
head(mydata)
a <- ggplot(mydata,aes(x = UMAP_1,y =UMAP_2,colour = inflammatory_gene1))+
geom_point(size = 1)+
scale_color_gradientn(values = seq(0,1,0.2),colours = c(‘blue’,‘cyan’,‘green’,‘yellow’,‘orange’,‘red’))
a+ theme_bw() + theme(panel.border = element_blank(), panel.grid.major = element_blank(),
panel.grid.minor = element_blank(), axis.line = element_line(colour = “black”))
mydata<- FetchData(subset_data,vars = c(“UMAP_1”,“UMAP_2”,“inflammatory_gene1”,“stim”,“cell.type”))
head(mydata)
class_avg <- mydata %>%
group_by(cell.type) %>%
summarise(
UMAP_1 = median(UMAP_1),
UMAP_2 = median(UMAP_2),
inflammatory_gene1=mean(inflammatory_gene1)
)
head(class_avg)
A tibble: 2 x 4
cell.type UMAP_1 UMAP_2 inflammatory_gene1
1 Macrophage-1 -3.66 -0.153 0.0131
2 Macrophage-2 11.2 1.77 0.0837
ggplot(mydata,aes(x = UMAP_1,y =UMAP_2,
colour = inflammatory_gene1))+geom_point(size = 1)+
scale_color_gradientn(values = seq(0,1,0.2),colours = c(‘blue’,‘cyan’,‘green’,‘yellow’,‘orange’,‘red’))+ #face = c(“plain”, “bold”, “italic”)
geom_text(aes(label = cell.type,fontface=“plain”,family=“sans”), data = class_avg,colour = “red”)+ # ggrepel::geom_text_repel(aes(label = cell.type), data = class_avg) +#添加标签
theme(text=element_text(family=“sans”,size=18)) +
facet_wrap(~stim)
getwd()
[1] “G:/silicosis/sicosis/silicosis_ST/yll/0214/harmony_cluster”
在这里插入图片描述
ggplot画 ump 和tsne 从seurat中使用addmodule得到的umap 使用ggplot画图相关推荐
- R语言ggplot画世界地图并根据条件给国家上色
期末结束,来分享一下 Exploration Data Analysis这门课我做的一个图. 语言依旧是 R, 用的Rstudio,这玩意儿真的拉闸,作为一个数据分析主要语言的软件,我数据量一大它就各 ...
- 函数图像轻松画:教你用永中图象
函数图像轻松画:教你用永中图象 函数图像轻松画:教你用永中图象 转载于:https://blog.51cto.com/premium/933220
- 命名空间system中不存在data_patternplot包:用ggplot解决你对线性填充,不!所有填充的全部幻想。...
patternplot包:用ggplot解决你对线性填充,不!所有填充的全部幻想.mp.weixin.qq.com 写在前面 patternplot包,提供了丰度的图形可视化填充选项,但是目前我尽然 ...
- VC++ 利用MFC的CWindowDC类实现画线功能 在桌面窗口中画线 绘制彩色线条 CPen nPenStyle nWidth crColor
目录 利用MFC的CWindowDC类实现画线功能 在桌面窗口中画线 绘制彩色线条 CPen nPenStyle nWidth crColor 接上:VC++ 绘制线条 OnLButtonDown函数 ...
- R语言图层—利用ggplot画蝙蝠侠logo及分层
R语言图层-利用ggplot画蝙蝠侠logo及分层 batman logo是由六个函数构成的. 利用ggplot2将多个线画到同一幅图案里面. 以下是r代码: require(ggplot2)f1 & ...
- python ggplot画等值线图_传说中的画图神器Plotnine,Python中的ggplot2
R的忠实用户知道,ggplot2可以使您在处理探索性数据分析和数据可视化时更加简单.它使得创建优雅而强大的情节变得如此容易,从而有助于解读数据中的潜在关系. 那么Python用户是否也有类似的画图工具 ...
- r语言ggplot画两条曲线_R语言作图——Line plot with error
原创:黄小仙 为了画今天的这个图,小仙决定凭空想象一台可以实时监控基因表达水平的设备,成功得到了这么一组数据. 想要画的图是这样子滴. Step1.绘图数据的准备 首先要把你想要绘图的数据调整成R语言 ...
- python画三维立体图-如何在论文中画出漂亮的插图?
----2020.08.07增---- 看到评论区有人说"没代码没教程所以没帮助"-- 好吧,我寻思链接.参考资料都在回答中给出来了呀(可能不够明显?) 于是,重新整理.注释了一下 ...
- 为什么python画不了图-解决python中使用plot画图,图不显示的问题
python matplotlib画的图怎么显示不出来 初学python,matplotlib库画图不显示分享助别熬夜了,你等不来的人,等到多晚都不会来的. 见图.python matplotlib ...
- 将不同数据来源的ggplot图绘制到同一张图中,并添加统一的图例
在R语言中使用ggplot绘制多幅图时,如果每幅图的数据来源或者风格不同时很难用分面绘制.这时可以使用ggpubr包的ggarrange函数将多幅图汇总到一张画布上. 还可以通过common.lege ...
最新文章
- 为反对种族歧视,GitHub正在把master/slave等术语替换掉
- 在Linux环境下使用OpenSSL对消息和文件进行加密(转载)
- 谷歌发安全警告:社交网络威胁用户隐私
- LeetCode算法入门- Longest Common Prefix -day13
- Scintilla 3 24在MFC中的使用 动态 静态
- opencv+python机读卡识别
- The Python Challenge 题解
- Linux之crontab(计划任务)
- Java项目集成永中文档在线预览
- 利用Python+xarray实现遥感数据——海表温度的经验正交函数(EOF)分解——xarray学习文档02
- 如何安装数据库和数据库安装不了如何解决
- ffmpeg 推流 m3u8 播放 EasyDarwin 存储录像和播放
- 到全国各省市关于公务员考试的信息
- 最简单的rar压缩包文档解密方法,rar压缩包不能复制打印怎么办?
- 关于召回率和hit rate
- 最新parity 以太坊客户端2.7.2稳定版配置参数说明
- 编写函数比较字符串的大小
- FTP、FTPS和SFTP
- MAC远程桌面控制Ubuntu
- iOS程序员眼中的首次使用产品体验
热门文章
- 应用系统安全检测技术
- Unity打包exe生成的文件夹合并成单独的exe
- spring boot + Spring Security 入门
- C++:Timer类实现
- cf 581A— Vasya the Hipster
- 如何合理选择AI加速器?
- Cisco ASA 9.16(Adaptive Security Appliance (ASA) Software)
- 简单实现antd的表单设计
- python 移动平均线_Python中的SMA(短期移动平均线)
- 访问页面出错:The requested URL / was not found on this server. 网站打开出现404