跟著Nature Medicine學作圖--箱線圖+散點圖

封面

從這個系列開始,師兄就帶著大家從各大頂級期刊中的Figuer入手症概,從仿照別人的作圖風格到最后實現(xiàn)自己游刃有余的套用在自己的分析數(shù)據(jù)上蕾额!這一系列絕對是高質(zhì)量!還不趕緊點贊+在看彼城,學起來诅蝶!

參考文獻

本期分享的是昨天更新的Nature Medicine上面那篇文章中的一個箱線圖和散點圖。

之所以講這張圖是因為這張圖既通過散點圖描繪了兩組樣本點的分布差異募壕,同時在X和Y兩個尺度上分別繪制箱線圖调炬,既反映了兩組在X和Y兩個尺度上的顯著差異,也更好的闡述了X和Y的相關(guān)關(guān)系舱馅。此外箱線圖還加上了連線缰泡,顯示出了從術(shù)后1天,到術(shù)后3個月的變化趨勢代嗤,非常好的呈現(xiàn)方式棘钞,值得大家學習缠借!

示例數(shù)據(jù)和代碼獲取

https://mp.weixin.qq.com/s/sPi2GxZeDpDjLpcXlPzsvw

讀圖

原圖

這張圖通過廣義估計方程(GEE)分析,發(fā)現(xiàn)在袖狀胃切除(SG)治療后B. thetaiotaommicron豐度的增加與BMI的下降相關(guān)宜猜。

效果預覽

復現(xiàn)效果

示例數(shù)據(jù)和R包載入

# 加載包:
library(ggplot2)

# 讀取數(shù)據(jù):
new_data <- read.csv("data01.csv")

> new_data
   Bacteroides_thetaiotaomicron      BMI    group   group2
1                     -13.36254 3.941414 group_0M  sample1
2                     -16.49327 3.927730 group_0M  sample2
3                     -15.85773 3.692952 group_0M  sample3
4                     -16.86678 3.948474 group_0M  sample4
5                     -15.57506 3.579352 group_0M  sample5
6                     -12.02062 3.759882 group_0M  sample6
7                     -16.75177 4.155833 group_0M  sample7
8                     -15.62603 4.045753 group_0M  sample8
9                     -13.00599 3.570529 group_0M  sample9
10                    -15.37026 3.522664 group_0M sample10
11                    -14.68648 3.931393 group_0M sample11
12                    -13.89668 3.971936 group_0M sample12
13                    -16.94041 3.505565 group_0M sample13
14                    -15.68604 4.082791 group_0M sample14
15                    -13.40720 4.124388 group_0M sample15

繪制

  • 首先繪制最中間的散點圖:
# 繪制散點圖:
ggplot(new_data)+
  # 繪制基本散點圖
  geom_point(aes(Bacteroides_thetaiotaomicron, BMI,color=group))+
  # 設(shè)置顏色
  scale_color_manual(values = c(group_0M="#ff00ff", group_3M="#8ac53e"))+
  # 設(shè)置主題
  theme_classic()+
  # 設(shè)置坐標軸范圍
  scale_x_continuous(breaks = c(-17:11))+
  scale_y_continuous(breaks = seq(3.2,4.0,0.2))+
  # 設(shè)置x軸和y軸標簽
  xlab("Bacteroides thetaiotaomicron (11021)")+
  ylab("BMI(lg)")+
  # 去掉圖例:
  theme(legend.position = 'none')

ggsave(filename = "scatter_plot.pdf",height=5,width = 5)
散點圖
  • 再繪制左邊的箱線圖
# 繪制左邊的箱線圖:
ggplot(new_data,aes(group, BMI))+
  # 加上誤差棒泼返;由于自帶的箱形圖沒有胡須末端沒有短橫線,使用誤差條的方式補上
  stat_boxplot(geom = "errorbar",width=0.15,aes(color=group))+ 
  # 繪制基本箱線圖
  geom_boxplot(aes(color=group),fill="white")+
  # 加上散點之間的連線
  geom_line(aes(group=group2), color="black",linetype="dashed",size=0.2, alpha=0.8)+
  # 箱線圖加散點
  geom_jitter(aes(color=group,fill=group),width =0.05,shape = 21)+ #設(shè)置為向水平方向抖動的散點圖姨拥,width指定了向水平方向抖動绅喉,不改變縱軸的值
  # 設(shè)置顏色
  scale_color_manual(values = c(group_0M="#ff00ff", group_3M="#8ac53e"))+
  scale_fill_manual(values = c(group_0M="#ff00ff", group_3M="#8ac53e"))+
  # 設(shè)置主題
  theme_classic()+
  # 設(shè)置坐標軸范圍
  scale_x_discrete(labels=c("0M","3M"))+
  scale_y_continuous(breaks = seq(3.2,4.0,0.2))+
  # 設(shè)置x軸和y軸標簽
  xlab("")+
  ylab("")+
  # 去掉圖例:
  theme(legend.position = 'none')

ggsave(filename = "left_boxplot.pdf",height=5,width = 2)
左側(cè)箱線圖
  • 最后繪制下方的箱線圖
# 繪制下面的箱線圖:
ggplot(new_data,aes(group, Bacteroides_thetaiotaomicron))+
  # 加上誤差棒;由于自帶的箱形圖沒有胡須末端沒有短橫線叫乌,使用誤差條的方式補上
  stat_boxplot(geom = "errorbar",width=0.15,aes(color=group))+ 
  # 繪制基本箱線圖
  geom_boxplot(aes(color=group),fill="white")+
  # 加上散點之間的連線
  geom_line(aes(group=group2), color="black", linetype="dashed", size=0.2, alpha=0.8)+
  # 箱線圖加散點
  geom_jitter(aes(color=group,fill=group),width =0.05,shape = 21)+ #設(shè)置為向水平方向抖動的散點圖柴罐,width指定了向水平方向抖動,不改變縱軸的值
  # 設(shè)置顏色
  scale_color_manual(values = c(group_0M="#ff00ff", group_3M="#8ac53e"))+
  scale_fill_manual(values = c(group_0M="#ff00ff", group_3M="#8ac53e"))+
  # 設(shè)置主題
  theme_classic()+
  # 設(shè)置坐標軸范圍
  scale_x_discrete(labels=c("0M","3M"))+
  scale_y_continuous(breaks = c((-17):(-11)))+
  # 設(shè)置x軸和y軸標簽
  xlab("")+
  ylab("")+
  # 去掉圖例
  theme(legend.position = 'none')+
  # 旋轉(zhuǎn)坐標軸
  coord_flip()

ggsave(filename = "down_boxplot.pdf",height=2,width = 5)
下方箱線圖
  • 最后综芥,用AI將三個圖拼接在一起就完事啦丽蝎!

結(jié)果展示

復現(xiàn)效果

示例數(shù)據(jù)和代碼獲取

https://mp.weixin.qq.com/s/sPi2GxZeDpDjLpcXlPzsvw

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市膀藐,隨后出現(xiàn)的幾起案子屠阻,更是在濱河造成了極大的恐慌,老刑警劉巖额各,帶你破解...
    沈念sama閱讀 217,277評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件国觉,死亡現(xiàn)場離奇詭異,居然都是意外死亡虾啦,警方通過查閱死者的電腦和手機麻诀,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,689評論 3 393
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來傲醉,“玉大人蝇闭,你說我怎么就攤上這事∮脖希” “怎么了呻引?”我有些...
    開封第一講書人閱讀 163,624評論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長吐咳。 經(jīng)常有香客問我逻悠,道長,這世上最難降的妖魔是什么韭脊? 我笑而不...
    開封第一講書人閱讀 58,356評論 1 293
  • 正文 為了忘掉前任童谒,我火速辦了婚禮,結(jié)果婚禮上沪羔,老公的妹妹穿的比我還像新娘饥伊。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 67,402評論 6 392
  • 文/花漫 我一把揭開白布撵渡。 她就那樣靜靜地躺著融柬,像睡著了一般。 火紅的嫁衣襯著肌膚如雪趋距。 梳的紋絲不亂的頭發(fā)上粒氧,一...
    開封第一講書人閱讀 51,292評論 1 301
  • 那天,我揣著相機與錄音节腐,去河邊找鬼外盯。 笑死,一個胖子當著我的面吹牛翼雀,可吹牛的內(nèi)容都是我干的饱苟。 我是一名探鬼主播,決...
    沈念sama閱讀 40,135評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼狼渊,長吁一口氣:“原來是場噩夢啊……” “哼箱熬!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起狈邑,我...
    開封第一講書人閱讀 38,992評論 0 275
  • 序言:老撾萬榮一對情侶失蹤城须,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后米苹,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體糕伐,經(jīng)...
    沈念sama閱讀 45,429評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,636評論 3 334
  • 正文 我和宋清朗相戀三年蘸嘶,在試婚紗的時候發(fā)現(xiàn)自己被綠了良瞧。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,785評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡训唱,死狀恐怖褥蚯,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情况增,我是刑警寧澤遵岩,帶...
    沈念sama閱讀 35,492評論 5 345
  • 正文 年R本政府宣布,位于F島的核電站巡通,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏舍哄。R本人自食惡果不足惜宴凉,卻給世界環(huán)境...
    茶點故事閱讀 41,092評論 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望表悬。 院中可真熱鬧弥锄,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,723評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至戒悠,卻和暖如春熬荆,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背绸狐。 一陣腳步聲響...
    開封第一講書人閱讀 32,858評論 1 269
  • 我被黑心中介騙來泰國打工卤恳, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人寒矿。 一個月前我還...
    沈念sama閱讀 47,891評論 2 370
  • 正文 我出身青樓突琳,卻偏偏與公主長得像,于是被迫代替她去往敵國和親符相。 傳聞我的和親對象是個殘疾皇子拆融,可洞房花燭夜當晚...
    茶點故事閱讀 44,713評論 2 354

推薦閱讀更多精彩內(nèi)容