image.png 今天试着重复的图片对应着的是论文中的Figure3a,,曼哈顿图展示的是XP-CLR score,这个分支具体是什么含义还得仔细学习,群体基因组类的论文通常都有这个 数据可以直接在论文的补充材料下载 首先是读入数据 df<-readxl::read_excel("NG/41588_2021_831_MOESM6_ESM.xlsx", sheet = "Fig3a") head(df)...
这里用到了一个R包ggh4x,可以很方便的操作坐标轴,具体可以参考推文通过ggh4x调整ggplot2图形刻度 完整的作图代码 library(ggplot2) #install.packages("ggh4x") library(ggh4x) p1<-ggplot(data=df,aes(x=x,y=`XP-CLR score`))+ geom_jitter(aes(color=Chromosome),width = 0.5)+ scale_y_continuous(#e...