大橙子网站建设,新征程启航
为企业提供网站建设、域名注册、服务器等服务
怎么用R语言ggplot2散点图并添加拟合曲线和置信区间,相信很多没有经验的人对此束手无策,为此本文总结了问题出现的原因和解决方法,通过这篇文章希望你能解决这个问题。
目前创新互联已为上千余家的企业提供了网站建设、域名、虚拟主机、成都网站托管、企业网站设计、临泉网站维护等服务,公司将坚持客户导向、应用为本的策略,正道将秉承"和谐、参与、激情"的文化,与客户和合作伙伴齐心协力一起成长,共同发展。
今天的内容是补充材料Fig4中的小a中的第三个小图
数据准备好是如下格式
需要示例数据可以直接在文末留言
mydf<-read.csv("Single_Cell/covid_metadata_c.csv",header=T,row.names = 1)
dim(mydf)
mydf
这里作图用到的是ggscatter()
函数,这个函数来自ggpubr
这个包。
library(ggpubr)
ggscatter(mydf,x="DPS",y="ISG")
这里需要注意的是想x,y后面跟着的值需要加双引号
p1<-ggscatter(mydf,x="DPS",y="ISG",
add = "reg.line")
p2<-ggscatter(mydf,x="DPS",y="ISG",
add = "reg.line",
conf.int = T)
p3<-ggscatter(mydf,x="DPS",y="ISG",
add = "reg.line",
conf.int = T)+
labs(y = "Mean ISG module score",
x = "Days post-symptom onset")
p4<-ggscatter(mydf,x="DPS",y="ISG",
add = "reg.line",
conf.int = T)+
labs(y = "Mean ISG module score",
x = "Days post-symptom onset")+
stat_cor(label.x = 8, label.y = 0.2)
ggarrange(p1,p2,p3,p4,
ncol = 2,
nrow = 2,
labels = c("p1","p2","p3","p4"))
最终的结果如下