Performs Shapiro-Wilk for sample sizes between 3 and 5000 (inclusive),
otherwise uses Anderson–Darling. Returns "Normal" if p > 0.05, else "Not Normal".
Usage
shapiro_normality_test(data)
Value
Character scalar: "Normal", "Not Normal", or NA if not applicable.