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)