# Leitura de dados - - - - - - - - - - - - - - -
base <- read.table("Siderurgia.csv",header=TRUE,sep=";",dec=",")
# Ajuste da base - - - - - - - - 
base.q = base[,-1]
row.names(base.q) = base[,1]
base.z = scale(base.q)
# Cálculos de distâncias - - - - - - - - - - - -
d = dist(base.z,method="euclidean")
d2 = d^2
# Análise de cluster hierárquico - - - - - - - - 
clust = hclust(d,method="ward")
plot(clust)
grupos = cutree(clust,k=4)
# Resultados - - - - - - - - - - - - - - - - - -
clust[]
grupos[]
# Adiciona o grupo à base de dados "base"
base$grupo = grupos
base
# Boxplot- - - - - - - - - - - - - - - - - - - - 
attach(base) 
boxplot(Fatur~grupo)
boxplot(Rent~grupo)
boxplot(Endiv~grupo)
boxplot(Empreg~grupo)
detach(base)