library(gtsummary)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(ggplot2)  
mtcars %>%
  select(mpg, cyl, gear) %>%
  tbl_summary(by = cyl) %>%
  add_p()
Characteristic 4
N = 111		 6
N = 71		 8
N = 141		 p-value2
mpg 26.0 (22.8, 30.4) 19.7 (18.1, 21.0) 15.2 (14.3, 16.4) <0.001
gear


<0.001
    3 1 (9.1%) 2 (29%) 12 (86%)
    4 8 (73%) 4 (57%) 0 (0%)
    5 2 (18%) 1 (14%) 2 (14%)
1 Median (Q1, Q3); n (%)
2 Kruskal-Wallis rank sum test; Fisher’s exact test 
mtcars %>%
  tbl_summary()
Characteristic N = 321
mpg 19.2 (15.4, 22.8)
cyl
    4 11 (34%)
    6 7 (22%)
    8 14 (44%)
disp 196 (121, 334)
hp 123 (96, 180)
drat 3.70 (3.08, 3.92)
wt 3.33 (2.54, 3.65)
qsec 17.71 (16.89, 18.90)
vs 14 (44%)
am 13 (41%)
gear
    3 15 (47%)
    4 12 (38%)
    5 5 (16%)
carb
    1 7 (22%)
    2 10 (31%)
    3 3 (9.4%)
    4 10 (31%)
    6 1 (3.1%)
    8 1 (3.1%)
1 Median (Q1, Q3); n (%)