【画图】与SARS-CoV-2病毒结合的ACE2受体在TCGA肺癌数据的表达情况（请不要过度解读这个图的结果！）
【画图】与SARS-CoV-2病毒结合的ACE2受体在TCGA肺癌数据的表达情况（请不要过度解读这个图的结果！）
Chris生命科学小站
发布于 2023-02-28 19:10:17
发布于 2023-02-28 19:10:17
1、获得分组信息
load("~/TCGAOS/TCGAOS/TCGA-LUAD_data.rda")
TCGA_LUADrace<-data.frame(row.names = rownames(coldata),race=coldata@listData[["race"]])2、获得ACE2表达信息
load("~/TCGAOS/TCGAOS/TCGA-LUAD_expr_data.rda")
colnames(vsd.TCGA.expr)<-str_replace_all(colnames(vsd.TCGA.expr), "[.]", "-")
ACE2_TCGA_LUAD_expr<-t(vsd.TCGA.expr["ENSG00000130234",])
colnames(ACE2_TCGA_LUAD_expr)<-"ACE2_expression"3、把上面两个信息揉在一起
ACE2_TCGA_LUAD_expr_m<-merge(ACE2_TCGA_LUAD_expr,TCGA_LUADrace,by=0)4、ggplot2做图
p2<-ggplot(ACE2_TCGA_LUAD_expr_m,aes(x=race,y=ACE2_expression))+
geom_boxplot(fill="#FF3D2E",colour="#1F3552",alpha=0.7)+
geom_jitter()+
ggtitle("TCGA_LUAD ACE2 expression")+
theme_bw()+
theme(axis.text.x=element_text(colour="black", size = 11))5、拼图输出pdf
p<-cowplot::plot_grid(p1,p2,align = "h")
ggsave(p,filename = "p.pdf",width = 23.5,height = 10)简单的Boxplot图是最常见的，供大家参考！
上面这个代码复制粘贴，是不能出图的，因为分组和表达的数据是直接加载之前保存的信息。
大家可以再文末赞赏，可以获得出图代码和清洗好的数据！
注意，不要过度解读这个这个图结果的意义！这仅仅是一个画图代码分享！
本文参与 腾讯云自媒体分享计划，分享自微信公众号。
原始发表：2020-02-20，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读