# Using libraries
library(mva)                # Loads the multivariate analysis library
library(help=mva)

search()                    # The search path

hist(choc$Fat)

choc <- read.table("/home/units/r.data/chocolates.dat", header=T, row.names=1)
choc
plot(choc)

choc[1,]
choc[1:3,]
choc[,1]

choc.hc <- hclust( dist(choc) )
plot(choc.hc)

choc.km <- kmeans(choc, 3)
names(choc.km)

choc.km$cluster
rownames(choc)[choc.km$cluster==1]
rownames(choc)[choc.km$cluster==2]
rownames(choc)[choc.km$cluster==3]