組間比較不只有箱線圖,還有這些

收集了幾張漂亮的組間比較圖们妥,調(diào)整了一些細(xì)節(jié)猜扮,分享給大家。使用的數(shù)據(jù)時(shí)R語言內(nèi)置數(shù)據(jù)iris监婶,以下代碼都可以直接運(yùn)行旅赢。

1.ggstatplot

難點(diǎn)是這個(gè)R包安裝時(shí)通常會遇到一些依賴包安裝不成功的問題,需要多折騰幾下压储,安好了使用起來就非常輕松咯

library(ggstatsplot)
ggbetweenstats(iris,x = "Species",y = "Sepal.Length")

2.經(jīng)典箱線圖疊加點(diǎn)圖

library(ggplot2)
library(ggpubr)
library(paletteer)
ggplot(iris,aes(x = Species,y=Sepal.Length))+
  geom_boxplot(aes(fill = Species))+
  geom_jitter(shape = 21,size = 2,color = "black",aes(fill = Species),stroke = 1.5)+
  scale_fill_paletteer_d("basetheme::minimal")+
  stat_compare_means(method = "wilcox.test",comparisons = list(c("setosa","virginica"),c("setosa","versicolor"),c("virginica","versicolor")))+
  theme_bw()+
  theme(legend.position = c("top"),panel.grid = element_blank())
ggplot(iris,aes(x = Species,y=Sepal.Length))+
  geom_boxplot(aes(fill = Species))+
  geom_dotplot(binaxis = "y",binwidth = 0.12,stackdir = "center",stroke = 1.5,aes(fill = Species))+
  scale_fill_paletteer_d("basetheme::minimal")+
  stat_compare_means(method = "wilcox.test",comparisons = list(c("setosa","virginica"),c("setosa","versicolor"),c("virginica","versicolor")))+
  theme_bw()+
  theme(legend.position = c("top"),panel.grid = element_blank())

jitter太跳脫了鲜漩,dotplot又比較呆板,有個(gè)折中的圖:蜜蜂圖

3.箱線圖疊加蜜蜂圖

library(ggbeeswarm)
ggplot(iris,aes(x = Species,y = Sepal.Length,fill = Species))+
  geom_boxplot()+
  geom_beeswarm(size = 3,cex = 3,shape = 21,stroke = 1.5)+
  theme_bw()+
  theme(legend.position = c("top"),panel.grid = element_blank())+
  stat_compare_means(method = "wilcox.test",comparisons = list(c("setosa","virginica"),c("setosa","versicolor"),c("virginica","versicolor")))+
  scale_fill_paletteer_d("basetheme::minimal")

4.微笑版密蜂圖加分位數(shù)線

有些糾結(jié)到底應(yīng)該是疊加分位數(shù)線集惋,還是疊加誤差棒孕似,發(fā)現(xiàn)其實(shí)兩個(gè)都說的過去,干脆都畫一下咯刮刑。

ggplot(iris,aes(x = Species,y = Sepal.Length,fill = Species))+
  geom_quasirandom(method = "smiley",size = 3,width = 0.25,shape = 21,stroke = 1.5)+
  theme_bw()+
  theme(legend.position = c("top"),panel.grid = element_blank())+
  stat_compare_means(method = "wilcox.test",comparisons = list(c("setosa","virginica"),c("setosa","versicolor"),c("virginica","versicolor")))+
  stat_summary(fun = median, fun.min = median, fun.max = median, 
               geom = 'crossbar', width = 0.3, size = 0.4,color = 'black') +
  stat_summary(fun.data = function(x) median_hilow(x, 0.5), 
               geom = 'errorbar', width = 0.25, size = 1,color = 'black')+
  scale_fill_paletteer_d("basetheme::minimal")

5.微笑版蜜蜂圖疊加誤差棒

library(dplyr)
df2 <- group_by(iris,Species)%>%summarise(sd = sd(Sepal.Length),
                                          Sepal.Length=mean(Sepal.Length))
head(df2)
## # A tibble: 3 x 3
##   Species       sd Sepal.Length
##   <fct>      <dbl>        <dbl>
## 1 setosa     0.352         5.01
## 2 versicolor 0.516         5.94
## 3 virginica  0.636         6.59
ggplot(iris,aes(x = Species,y = Sepal.Length,fill = Species))+
  geom_quasirandom(method = "smiley",size = 3,width = 0.25,shape = 21,stroke = 1.5)+
  theme_bw()+
  theme(legend.position = c("top"),panel.grid = element_blank())+
  stat_compare_means(method = "wilcox.test",comparisons = list(c("setosa","virginica"),c("setosa","versicolor"),c("virginica","versicolor")))+
  stat_summary(fun = median, fun.min = median, fun.max = median, 
               geom = 'crossbar', width = 0.3, size = 0.4,color = 'black')+
  geom_errorbar(dat = df2,aes(ymin=Sepal.Length-sd, ymax=Sepal.Length+sd), width=.2)+
  scale_fill_paletteer_d("basetheme::minimal")

參考代碼:http://www.sthda.com/english/wiki/ggplot2-error-bars-quick-start-guide-r-software-and-data-visualization
https://paulvanderlaken.com/2019/01/25/visualization-innovation-waffleplots-and-swarmplots-aka-beeplots/
https://mp.weixin.qq.com/s/8LwTRKTlOR0CsQDUc15sBA

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末喉祭,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子雷绢,更是在濱河造成了極大的恐慌泛烙,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,919評論 6 502
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件翘紊,死亡現(xiàn)場離奇詭異蔽氨,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,567評論 3 392
  • 文/潘曉璐 我一進(jìn)店門鹉究,熙熙樓的掌柜王于貴愁眉苦臉地迎上來宇立,“玉大人,你說我怎么就攤上這事自赔÷栲冢” “怎么了?”我有些...
    開封第一講書人閱讀 163,316評論 0 353
  • 文/不壞的土叔 我叫張陵绍妨,是天一觀的道長润脸。 經(jīng)常有香客問我,道長他去,這世上最難降的妖魔是什么毙驯? 我笑而不...
    開封第一講書人閱讀 58,294評論 1 292
  • 正文 為了忘掉前任,我火速辦了婚禮孤页,結(jié)果婚禮上尔苦,老公的妹妹穿的比我還像新娘。我一直安慰自己行施,他們只是感情好允坚,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,318評論 6 390
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著蛾号,像睡著了一般稠项。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上鲜结,一...
    開封第一講書人閱讀 51,245評論 1 299
  • 那天展运,我揣著相機(jī)與錄音,去河邊找鬼精刷。 笑死拗胜,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的怒允。 我是一名探鬼主播埂软,決...
    沈念sama閱讀 40,120評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼纫事!你這毒婦竟也來了勘畔?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 38,964評論 0 275
  • 序言:老撾萬榮一對情侶失蹤丽惶,失蹤者是張志新(化名)和其女友劉穎炫七,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體钾唬,經(jīng)...
    沈念sama閱讀 45,376評論 1 313
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡万哪,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,592評論 2 333
  • 正文 我和宋清朗相戀三年侠驯,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片壤圃。...
    茶點(diǎn)故事閱讀 39,764評論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡陵霉,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出伍绳,到底是詐尸還是另有隱情,我是刑警寧澤乍桂,帶...
    沈念sama閱讀 35,460評論 5 344
  • 正文 年R本政府宣布冲杀,位于F島的核電站,受9級特大地震影響睹酌,放射性物質(zhì)發(fā)生泄漏权谁。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,070評論 3 327
  • 文/蒙蒙 一憋沿、第九天 我趴在偏房一處隱蔽的房頂上張望旺芽。 院中可真熱鬧,春花似錦辐啄、人聲如沸采章。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,697評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽悯舟。三九已至,卻和暖如春砸民,著一層夾襖步出監(jiān)牢的瞬間抵怎,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,846評論 1 269
  • 我被黑心中介騙來泰國打工岭参, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留反惕,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,819評論 2 370
  • 正文 我出身青樓演侯,卻偏偏與公主長得像姿染,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個(gè)殘疾皇子蚌本,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,665評論 2 354

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

  • 給直方圖和線圖添加誤差棒 準(zhǔn)備數(shù)據(jù) 這里使用ToothGrowth 數(shù)據(jù)集盔粹。它描述了維他命C對Guinea豬牙齒的...
    王詩翔閱讀 7,410評論 2 19
  • 話說“一圖勝千言”,在各類數(shù)據(jù)分析報(bào)告中經(jīng)常會看見各種各樣的圖形程癌,例如折線圖舷嗡、條形圖、箱線圖嵌莉、點(diǎn)圖等进萄。 其中折線圖...
    生信補(bǔ)給站閱讀 9,116評論 0 5
  • 本系列文章的目的:我將首先簡要介紹一下ggplot2的基本思路,然后使用前面介紹的TCGA數(shù)據(jù)模仿articles...
    ZZZZZZ_XX閱讀 3,935評論 0 14
  • 我是黑夜里大雨紛飛的人啊 1 “又到一年六月,有人笑有人哭中鼠,有人歡樂有人憂愁可婶,有人驚喜有人失落,有的覺得收獲滿滿有...
    陌忘宇閱讀 8,536評論 28 53
  • 信任包括信任自己和信任他人 很多時(shí)候援雇,很多事情矛渴,失敗、遺憾惫搏、錯過具温,源于不自信,不信任他人 覺得自己做不成筐赔,別人做不...
    吳氵晃閱讀 6,187評論 4 8