R语言ggstatsplot画回归拟合图

suppressMessages(library(ggstatsplot)) #不显示加载包附带的信息
library(ggplot2)

#将mtcars数据集的行名转换成mtcars数据集的一个列，将数据集存储为mtcars_new。
mtcars_new <- mtcars %>%
  tibble::rownames_to_column(., var = "car") %>%  #将mtcars的行名存储为‘car’
  tibble::as_tibble(x = .)  #将dataframe转换为tibble格式

mtcars_new #查看数据
# A tibble: 32 x 12
   car                 mpg   cyl  disp    hp  drat    wt  qsec    vs    am  gear  carb
                          
 1 Mazda RX4          21       6  160    110  3.9   2.62  16.5     0     1     4     4
 2 Mazda RX4 Wag      21       6  160    110  3.9   2.88  17.0     0     1     4     4
 3 Datsun 710         22.8     4  108     93  3.85  2.32  18.6     1     1     4     1
 4 Hornet 4 Drive     21.4     6  258    110  3.08  3.22  19.4     1     0     3     1
 5 Hornet Sportabout  18.7     8  360    175  3.15  3.44  17.0     0     0     3     2
 6 Valiant            18.1     6  225    105  2.76  3.46  20.2     1     0     3     1
 7 Duster 360         14.3     8  360    245  3.21  3.57  15.8     0     0     3     4
 8 Merc 240D          24.4     4  147.    62  3.69  3.19  20       1     0     4     2
 9 Merc 230           22.8     4  141.    95  3.92  3.15  22.9     1     0     4     2
10 Merc 280           19.2     6  168.   123  3.92  3.44  18.3     1     0     4     4
# ... with 22 more rows

ggstatsplot::ggscatterstats(
  data = mtcars_new,
  x = wt,
  y = mpg,
  label.var = car, #标签名称为car变量
  label.expression = wt < 4 & mpg < 20, #显示出来的标签为wt < 4 和mpg < 24的car， 其他的car不显示。
  axes.range.restrict = TRUE,
  centrality.para = "median", #集中参数为中位数
)

image.png

#将边际图形设置为密度曲线
ggstatsplot::ggscatterstats(
  data = mtcars_new,
  x = wt,
  y = mpg,
  marginal.type = 'density',
  label.var = car, #标签名称为car变量
  label.expression = wt < 4 & mpg < 20, #显示出来的标签为wt < 4 和mpg < 24的car， 其他的car不显示。
  centrality.para = "mean", #集中参数为均数
)

image.png