A tibble data frame with one row for each value of the by variable,
or one row overall if there is no by variable. For the y variable whose
normality is being tested, for each subset of values corresponding to the
values of they by variable, or for all values if there is no by
variable, return the Shapiro-Wilk statistic, and the corresponding p-value.