Тестирование на нормальность и отображение комбинированных результатов для нескольких столбцов в кадре данных (R)

Question

Я пытаюсь объединить два графика (ggqplot и ggdensity) и включить значение p из выходных данных теста shapiro для нескольких столбцов кадра данных.

Ввод - dataframe (df) список интересующих столбцов - "varlist"

Я извлек имена столбцов, с которыми я хочу извлечь (пример) varlist <- names(df[c(10:20)] Я хочу L oop? или используйте другой метод для получения вывода значений qqplot, ggdensity и P для теста shapiro, отображаемых вместе.

library(ggpubr)

ggqplot(df$i)
ggdensity(df$i)
shapiro.test(df$i)

Это было бы полезно для других, желающих быстро получить все эти показатели на своих данных , Я посмотрел, но не могу найти, как собрать это воедино. Я думаю, что функция, которая делает каждый граф и объединяет информацию, а затем передать через него я из varlist?

Я извлек эту функцию из предыдущего ответа, который быстро получает результаты теста Шапиро - это может быть полезно?

myfun <- function(x, group) {
  data.frame(x, group) %>%
    group_by(group) %>%
    summarise(
      statistic = ifelse(sd(x)!=0,shapiro.test(x)$statistic,NA), 
      p.value = ifelse(sd(x)!=0,shapiro.test(x)$p.value,NA)
    )
}