Группировка графика Лайкерта с помощью sjPlot

Я хотел бы создать Likert-Plots с помощью удивительного sjPlot-Package (немного похожего, но другого пакета: Построение похожих на шкалу вопросов и их группировка), где относительная частота каждого элемента сгруппирована по фактору (например, полу). Есть ли способ сделать это без этого обходного пути (который ему так не нравится)?

library(tidyverse)
library(sjPlot)

plot_df <- tibble(VAR_1 = sample(c(1:4), 100, replace = TRUE, prob = c(0.1, 0.3, 0.4, 0.2)),
                 VAR_2 = sample(c(1:4), 100, replace = TRUE, prob = c(0.4, 0.25, 0.25, 0.2)),
                 sex = sample(c("male", "female"), 100, replace = TRUE, prob = c(0.49, 0.51)))


plot_df <- plot_df %>% 
 transmute(VAR_1_male = case_when(sex == "male" ~ VAR_1,
                                  TRUE  ~ NA_integer_),
           VAR_1_female = case_when(sex == "female" ~ VAR_1,
                                    TRUE ~ NA_integer_),
           VAR_2_male = case_when(sex == "male" ~ VAR_2,
                                  TRUE  ~ NA_integer_),
           VAR_2_female = case_when(sex == "female" ~ VAR_2,
                                    TRUE ~ NA_integer_))

plot_df %>% plot_likert()