In [8]:
data.frame(PC1 = pca$x[,1], PC2 = pca$x[,2],
           label=factor(mnist$train$label)) %>%
  sample_n(1000) %>%
  ggplot(aes(PC1, PC2, fill=label)) +
    geom_point(cex=3, pch=21)