sdc <- read.csv("SDC-bisque-scores.txt",sep="\t")
pca <- prcomp(t(data.matrix(sdc)),scale.=T)
plot(x=pca$x[,1],y=pca$x[,2])

heatmap(data.matrix(sdc))
d <- dist(t(data.matrix(sdc)))
h <- hclust(d,method = "ward.D")
plot(h)
clusters <- cutree(h,3)
clusters

plot(x=pca$x[,1],y=pca$x[,2],pch=19,col=clusters)

library(Rtsne)
tsne <- Rtsne(t(data.matrix(sdc)), perplexity = 7)
plot(x=tsne$Y[,1],y=tsne$Y[,2],pch=19,col=clusters)

tsne <- Rtsne(pca$x[,1:10], perplexity = 7)
plot(x=tsne$Y[,1],y=tsne$Y[,2],pch=19,col=clusters)