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
<chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
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", #集中参数为中位数
)
#将边际图形设置为密度曲线
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", #集中参数为均数
)