# NOT RUN {
# Import data.
fname <- system.file("extdata", "finndiane.txt", package = "Numero")
dataset <- read.delim(file = fname)
# Prepare training data.
trvars <- c("CHOL", "HDL2C", "TG", "CREAT", "uALB")
trdata <- scale.default(dataset[,trvars])
# K-means clustering.
km <- nroKmeans(data = trdata)
# Self-organizing map.
sm <- nroKohonen(seeds = km)
sm <- nroTrain(map = sm, data = trdata)
# Assign data points into districts.
matches <- nroMatch(centroids = sm, data = trdata)
# District averages for all variables.
planes <- nroAggregate(topology = sm, districts = matches, data = dataset)
# District colors for cholesterol.
chol <- nroColorize(values = planes$CHOL)
print(head(chol))
# District colors for all variables.
colrs <- nroColorize(values = planes)
print(head(colrs))
# }
Run the code above in your browser using DataLab