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画图相关推荐

  1. R语言ggplot画世界地图并根据条件给国家上色

    期末结束,来分享一下 Exploration Data Analysis这门课我做的一个图. 语言依旧是 R, 用的Rstudio,这玩意儿真的拉闸,作为一个数据分析主要语言的软件,我数据量一大它就各 ...

  2. 函数图像轻松画:教你用永中图象

    函数图像轻松画:教你用永中图象 函数图像轻松画:教你用永中图象 转载于:https://blog.51cto.com/premium/933220

  3. 命名空间system中不存在data_patternplot包:用ggplot解决你对线性填充,不!所有填充的全部幻想。...

    patternplot包:用ggplot解决你对线性填充,不!所有填充的全部幻想.​mp.weixin.qq.com 写在前面 patternplot包,提供了丰度的图形可视化填充选项,但是目前我尽然 ...

  4. VC++ 利用MFC的CWindowDC类实现画线功能 在桌面窗口中画线 绘制彩色线条 CPen nPenStyle nWidth crColor

    目录 利用MFC的CWindowDC类实现画线功能 在桌面窗口中画线 绘制彩色线条 CPen nPenStyle nWidth crColor 接上:VC++ 绘制线条 OnLButtonDown函数 ...

  5. R语言图层—利用ggplot画蝙蝠侠logo及分层

    R语言图层-利用ggplot画蝙蝠侠logo及分层 batman logo是由六个函数构成的. 利用ggplot2将多个线画到同一幅图案里面. 以下是r代码: require(ggplot2)f1 & ...

  6. python ggplot画等值线图_传说中的画图神器Plotnine,Python中的ggplot2

    R的忠实用户知道,ggplot2可以使您在处理探索性数据分析和数据可视化时更加简单.它使得创建优雅而强大的情节变得如此容易,从而有助于解读数据中的潜在关系. 那么Python用户是否也有类似的画图工具 ...

  7. r语言ggplot画两条曲线_R语言作图——Line plot with error

    原创:黄小仙 为了画今天的这个图,小仙决定凭空想象一台可以实时监控基因表达水平的设备,成功得到了这么一组数据. 想要画的图是这样子滴. Step1.绘图数据的准备 首先要把你想要绘图的数据调整成R语言 ...

  8. python画三维立体图-如何在论文中画出漂亮的插图?

    ----2020.08.07增---- 看到评论区有人说"没代码没教程所以没帮助"-- 好吧,我寻思链接.参考资料都在回答中给出来了呀(可能不够明显?) 于是,重新整理.注释了一下 ...

  9. 为什么python画不了图-解决python中使用plot画图,图不显示的问题

    python matplotlib画的图怎么显示不出来 初学python,matplotlib库画图不显示分享助别熬夜了,你等不来的人,等到多晚都不会来的. 见图.python matplotlib ...

  10. 将不同数据来源的ggplot图绘制到同一张图中,并添加统一的图例

    在R语言中使用ggplot绘制多幅图时,如果每幅图的数据来源或者风格不同时很难用分面绘制.这时可以使用ggpubr包的ggarrange函数将多幅图汇总到一张画布上. 还可以通过common.lege ...

最新文章

  1. 为反对种族歧视,GitHub正在把master/slave等术语替换掉
  2. 在Linux环境下使用OpenSSL对消息和文件进行加密(转载)
  3. 谷歌发安全警告:社交网络威胁用户隐私
  4. LeetCode算法入门- Longest Common Prefix -day13
  5. Scintilla 3 24在MFC中的使用 动态 静态
  6. opencv+python机读卡识别
  7. The Python Challenge 题解
  8. Linux之crontab(计划任务)
  9. Java项目集成永中文档在线预览
  10. 利用Python+xarray实现遥感数据——海表温度的经验正交函数(EOF)分解——xarray学习文档02
  11. 如何安装数据库和数据库安装不了如何解决
  12. ffmpeg 推流 m3u8 播放 EasyDarwin 存储录像和播放
  13. 到全国各省市关于公务员考试的信息
  14. 最简单的rar压缩包文档解密方法,rar压缩包不能复制打印怎么办?
  15. 关于召回率和hit rate
  16. 最新parity 以太坊客户端2.7.2稳定版配置参数说明
  17. 编写函数比较字符串的大小
  18. FTP、FTPS和SFTP
  19. MAC远程桌面控制Ubuntu
  20. iOS程序员眼中的首次使用产品体验

热门文章

  1. 应用系统安全检测技术
  2. Unity打包exe生成的文件夹合并成单独的exe
  3. spring boot + Spring Security 入门
  4. C++:Timer类实现
  5. cf 581A— Vasya the Hipster
  6. 如何合理选择AI加速器?
  7. Cisco ASA 9.16(Adaptive Security Appliance (ASA) Software)
  8. 简单实现antd的表单设计
  9. python 移动平均线_Python中的SMA(短期移动平均线)
  10. 访问页面出错:The requested URL / was not found on this server. 网站打开出现404