shapiro.test(..) cannot deal with more than 5000 data points
This is a safety limitation. Please read this: Perform a Shapiro-Wilk Normality Test
Other tests of normality do not have this limitation such as the Kolmogorov-Smirnov test:
ks.test(x=rnorm(10^4),y='pnorm',alternative='two.sided')