R語言利用ggplot2繪制QQ圖和箱線圖詳解

更新時間：2022年06月09日 14:12:32 作者：zoujiahui_2018

這篇文章主要為大家介紹了R語言如何利用ggplot2繪制QQ圖和箱線圖，文中的示例代碼講解詳細(xì)，對我們學(xué)習(xí)R語言有一定的幫助，需要的可以參考一下

繪制qq圖

在ggplot2中繪制qq圖需要兩步，geom_qq()將繪制樣本分位點(diǎn)，geom_qq_line()將繪制標(biāo)準(zhǔn)正態(tài)線

函數(shù)介紹

geom_qq()

geom_qq(
  mapping = NULL,
  data = NULL,
  geom = "point",
  position = "identity",
  ...,
  distribution = stats::qnorm,
  dparams = list(),
  na.rm = FALSE,
  show.legend = NA,
  inherit.aes = TRUE
)

geom_qq_line(
  mapping = NULL,
  data = NULL,
  geom = "path",
  position = "identity",
  ...,
  distribution = stats::qnorm,
  dparams = list(),
  line.p = c(0.25, 0.75),
  fullrange = FALSE,
  na.rm = FALSE,
  show.legend = NA,
  inherit.aes = TRUE
)

參數(shù)介紹

**aes()**中的映射參數(shù)必須包含sample，可選參數(shù)有g(shù)roup,x,y distribution

Distribution function to use, if x not specified
dparams Additional parameters passed on to distribution function.
line.p Vector of quantiles to use when fitting the Q-Q line, defaults defaults to c(.25, .75).
fullrange Should the q-q line span the full range of the plot, or just the data

注意事項(xiàng)

**aes()**中的映射參數(shù)必須包含sample

例子

Using to explore the distribution of a variable

ggplot(mtcars, aes(sample = mpg)) +
  stat_qq() +
  stat_qq_line()
ggplot(mtcars, aes(sample = mpg, colour = factor(cyl))) +
  stat_qq() +
  stat_qq_line()

繪制boxplot

函數(shù)介紹

geom_boxplot(
  mapping = NULL,
  data = NULL,
  stat = "boxplot",
  position = "dodge2",
  ...,
  outlier.colour = NULL,
  outlier.color = NULL,
  outlier.fill = NULL,
  outlier.shape = 19,
  outlier.size = 1.5,
  outlier.stroke = 0.5,
  outlier.alpha = NULL,
  notch = FALSE,
  notchwidth = 0.5,
  varwidth = FALSE,
  na.rm = FALSE,
  orientation = NA,
  show.legend = NA,
  inherit.aes = TRUE
)

參數(shù)介紹

aes()可接收的參數(shù)有：

x or y，利用x將會是橫向箱線圖，y的是縱向
lower or xlower
upper or xupper
middle or xmiddle
ymin or xmin
ymax or xmax
alpha
colour
fill
group
linetype
shape
size
weight

notch If FALSE (default) make a standard box plot. If TRUE, make a notched box plot. Notches are used to compare groups; if the notches
of two boxes do not overlap, this suggests that the medians are
significantly different.
notchwidth For a notched box plot, width of the notch relative to the body (defaults to notchwidth = 0.5).
varwidth If FALSE (default) make a standard box plot. If TRUE, boxes are drawn with widths proportional to the square-roots of the
number of observations in the groups (possibly weighted, using the
weight aesthetic).

例子

p <- ggplot(mpg, aes(x=class, y=hwy))
p + geom_boxplot()

ggplot(mpg, aes(x=hwy, y=class)) + geom_boxplot()

p <- ggplot(mpg, aes(x=class, y=hwy))
p + geom_boxplot(notch = TRUE,varwidth = TRUE,fill = "white", colour = "#3366FF")

ggplot(diamonds, aes(carat, price)) +
  geom_boxplot(aes(group = cut_width(carat, 0.25)))

p <- ggplot(mpg, aes(x=class, y=hwy))
p + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)

利用分位點(diǎn)繪制箱線圖

y <- rnorm(100)
df <- data.frame(
  x = 1,
  y0 = min(y),
  y25 = quantile(y, 0.25),
  y50 = median(y),
  y75 = quantile(y, 0.75),
  y100 = max(y)
)
ggplot(df, aes(x)) +
  geom_boxplot(
    aes(ymin = y0, lower = y25, middle = y50, upper = y75, ymax = y100),
    stat = "identity"
  )

將QQ圖和箱線圖進(jìn)行融合

函數(shù)介紹

該函數(shù)是來自于qqboxplot包，因此使用前需要安裝

geom_qqboxplot(
  mapping = NULL,
  data = NULL,
  stat = "qqboxplot",
  position = "dodge2",
  ...,
  outlier.colour = NULL,
  outlier.color = NULL,
  outlier.fill = NULL,
  outlier.shape = 19,
  outlier.size = 1.5,
  outlier.stroke = 0.5,
  outlier.alpha = NULL,
  notch = FALSE,
  notchwidth = 0.5,
  varwidth = FALSE,
  na.rm = FALSE,
  show.legend = NA,
  inherit.aes = TRUE
)

參數(shù)介紹

大部分參數(shù)和geom_qq()和geom_boxplot()中的參數(shù)含義相同

reference_dist 表示參數(shù)比較的標(biāo)準(zhǔn)分布名稱，如果有參數(shù)需要有dparams

compdata 用于比較的標(biāo)準(zhǔn)樣本數(shù)據(jù)，是個向量

注意事項(xiàng)

aes()函數(shù)中的y不可缺

例子

library(dplyr)
library(ggplot2)
library(qqboxplot)

simulated_data=tibble(y=c(rnorm(1000, mean=2), rt(1000, 16), rt(500, 4), 
                          rt(1000, 8), rt(1000, 32)),
                      group=c(rep("normal, mean=2", 1000), 
                              rep("t distribution, df=16", 1000), 
                              rep("t distribution, df=4", 500), 
                              rep("t distribution, df=8", 1000), 
                              rep("t distribution, df=32", 1000)))
p <- ggplot2::ggplot(simulated_data, ggplot2::aes(factor(group,
                                                         levels=c("normal, mean=2", "t distribution, df=32", "t distribution, df=16",
                                                                  "t distribution, df=8", "t distribution, df=4")), y=y))
p + geom_qqboxplot()
p + geom_qqboxplot(reference_dist = "norm")


p + geom_qqboxplot(compdata = comparison_dataset)