# Shapiro Wilk normality test for one variable
iris %>% shapiro_test(Sepal.Length)
# Shapiro Wilk normality test for two variables
iris %>% shapiro_test(Sepal.Length, Petal.Width)
# Multivariate normality test
mshapiro_test(iris[, 1:3])
Run the code above in your browser using DataLab