Singular value decomposition of datasauRus

datasauRus_svd.R

library(datasauRus)
library(ggplot2)
library(dplyr)

X = dplyr::filter(datasaurus_dozen,dataset=="dino") %>% 
  dplyr::select(-dataset) %>% 
  as.matrix()
res_svd <- svd(as.matrix(X))

df_svd <- bind_rows(data.frame(res_svd$u, matrix="U"),
                    data.frame(res_svd$u%*%diag(res_svd$d), matrix="UD"),
                    data.frame(res_svd$u%*%res_svd$v, matrix="UV"),
                    data.frame(res_svd$u%*%diag(res_svd$d)%*%res_svd$v, matrix="UDV"))


ggplot(data = df_svd)+
  geom_point(aes(x = X1, y = X2))+
  facet_wrap(~matrix,scales = "free")+
  theme_classic(16)
ggsave("dino_svd.png")