# Predmet: "Analyza zhlukov a klasifikacia dat"
# Studijny program: "Pravdepodobnost a matematicka statistika"
# Vyucujuci: Radoslav Harman, KAMS FMFI UK Bratislava
#
# Nemetricke mnohorozmerne skalovanie
# Data su z clanku Ekman, G. (1954): Dimensions of color vision, Journal of Psychology, 38, 467-474

color.stimuli <- structure(c(0.86, 0.42, 0.42, 0.18, 0.06, 0.07, 0.04, 0.02, 0.07, 
0.09, 0.12, 0.13, 0.16, 0.5, 0.44, 0.22, 0.09, 0.07, 0.07, 0.02, 
0.04, 0.07, 0.11, 0.13, 0.14, 0.81, 0.47, 0.17, 0.1, 0.08, 0.02, 
0.01, 0.02, 0.01, 0.05, 0.03, 0.54, 0.25, 0.1, 0.09, 0.02, 0.01, 
0, 0.01, 0.02, 0.04, 0.61, 0.31, 0.26, 0.07, 0.02, 0.02, 0.01, 
0.02, 0, 0.62, 0.45, 0.14, 0.08, 0.02, 0.02, 0.02, 0.01, 0.73, 
0.22, 0.14, 0.05, 0.02, 0.02, 0, 0.33, 0.19, 0.04, 0.03, 0.02, 
0.02, 0.58, 0.37, 0.27, 0.2, 0.23, 0.74, 0.5, 0.41, 0.28, 0.76, 
0.62, 0.55, 0.85, 0.68, 0.76), Labels = c("w434", "w445", "w465", 
"w472", "w490", "w504", "w537", "w555", "w584", "w600", "w610", 
"w628", "w651", "w674"), Size = 14L, call = quote(as.dist.default(m = color.stimuli)), class = "dist", Diag = FALSE, Upper = FALSE)

cols<-c("magenta", "blue", "blue", "darkcyan", "darkgreen", "green", 
"lightgreen", "yellowgreen", "yellow", "tan", "orange", "orange", 
"red", "violet") # Toto je len moj velmi priblizny odhad farieb

library(MASS)
col.mds<-isoMDS(1-color.stimuli,y=matrix(runif(28),ncol=2),maxit=100);col.mds$stress
plot(col.mds$points[,1],col.mds$points[,2],col=cols,pch=19,cex=3)