被紅色框線(xiàn)圈住的散點(diǎn)圖

 

基礎(chǔ)繪圖函數(shù)箱線(xiàn)圖

有人留言說(shuō) 和原圖不是很像，因?yàn)榕渖珱](méi)有按照論文中提供的代碼來(lái)。 下面是完全重復(fù)論文中的代碼

cols <- c("#E69F00", "#56B4E9", "#009E73")
boxplot(log10(rel_crAss)~country,data=HMP,col=cols,
       axes=F,xlab=NULL,ylab=NULL,
       horizontal = T)
axis(2,at=c(1,2,3),labels=c("China", "Europe", "US"),las=1)
title("a",adj=0,line=0)

 

• 讀入數(shù)據(jù)

HMP<-read.table("data/HMP.txt")

• 最基本的散點(diǎn)圖

plot(rel_res~rel_crAss,data=HMP)

 

畫(huà)圖用plot()函數(shù)，需要指定畫(huà)圖用到的變量y和x，還有畫(huà)圖用到的數(shù)據(jù)data

原始代碼分別對(duì) rel_res 和 rel_crAss取了log10

plot(log10(rel_res)~log10(rel_crAss),data=HMP)

 

取log10以后可以看到散點(diǎn)分布的更加均勻了。

• 按照國(guó)家分組填充顏色

cols <- c("#E69F00", "#56B4E9", "#009E73", "#F0E442", "#0072B2", "#D55E00", "#CC79A7")
plot(log10(rel_res)~log10(rel_crAss), data=HMP,
    bg=cols[as.factor(HMP$country)],pch=21)

 

• 更改點(diǎn)的大小

plot(log10(rel_res)~log10(rel_crAss), data=HMP,
    bg=cols[as.factor(HMP$country)],pch=21,cex=2)

 

• 更改x軸和y軸的標(biāo)簽

plot(log10(rel_res)~log10(rel_crAss), data=HMP, bg=cols[as.factor(HMP$country)], pch=21,
    ylab = "Normalized ARG abundance (log10)",
    xlab="Normalized crAssphage abundance (log10)", cex=2)

 

• 更改坐標(biāo)軸的范圍

plot(log10(rel_res)~log10(rel_crAss), data=HMP,
    bg=cols[as.factor(HMP$country)], pch=21,
    ylab = "Normalized ARG abundance (log10)",
    xlab="Normalized crAssphage abundance (log10)",
    cex=2,
    ylim=c(2.5, 4.5))

 

接下來(lái)將箱線(xiàn)圖和散點(diǎn)圖按照上下拼接到一起，用到的是par(fig=c(a,b,c,d)),這里需要滿(mǎn)足 a<b,c<d

具體可以參考鏈接 https://blog.csdn.net/qingchongxinshuru/article/details/52004182

par(fig=c(0,1,0,0.75))
plot(log10(rel_res)~log10(rel_crAss), data=HMP,
    bg=cols[as.factor(HMP$country)], pch=21,
    ylab = "Normalized ARG abundance (log10)",
    xlab="Normalized crAssphage abundance (log10)",
    cex=2,
    ylim=c(2.5, 4.5))
par(fig=c(0,1,0.5,1),new=T)
boxplot(log10(rel_crAss)~country,data=HMP,col=cols,
       axes=F,xlab=NULL,ylab=NULL,
       horizontal = T)
axis(2,at=c(1,2,3),labels=c("China", "Europe", "US"),las=1)
title("a",adj=0,line=0)