rm(list=ls())

require(graphics)

data = read.csv("ILE-2013.csv",header=TRUE)

namevar = c("score","property","corrupt","fiscal","govern","busin","labor",
"monet","trade","invest","finan")

attach(data)

indicadores = data[,4:13]

pairs(indicadores)

hc = hclust(dist(indicadores),"complete")

pdf(file="cluster-ILE.pdf",width=20,height=10)
plot(hc, hang = -1,labels=data[,2],cex=0.4)
abline(h=125,lty=2)
dev.off()

memb = cutree(hc,h=125)

data[(memb==1)&(data[,1]=="Africa"),2]
data[(memb==1)&(data[,1]=="America"),2]
data[(memb==1)&(data[,1]=="Asia"),2]
data[(memb==1)&(data[,1]=="Europa"),2]
data[(memb==1)&(data[,1]=="Oceania"),2]


data[(memb==2)&(data[,1]=="Africa"),2]
data[(memb==2)&(data[,1]=="America"),2]
data[(memb==2)&(data[,1]=="Asia"),2]
data[(memb==2)&(data[,1]=="Europa"),2]
data[(memb==2)&(data[,1]=="Oceania"),2]


data[(memb==3)&(data[,1]=="Africa"),2]
data[(memb==3)&(data[,1]=="America"),2]
data[(memb==3)&(data[,1]=="Asia"),2]
data[(memb==3)&(data[,1]=="Europa"),2]
data[(memb==3)&(data[,1]=="Oceania"),2]

data[(memb==4)&(data[,1]=="Africa"),2]
data[(memb==4)&(data[,1]=="America"),2]
data[(memb==4)&(data[,1]=="Asia"),2]
data[(memb==4)&(data[,1]=="Europa"),2]
data[(memb==4)&(data[,1]=="Oceania"),2]

data[(memb==5)&(data[,1]=="Africa"),2]
data[(memb==5)&(data[,1]=="America"),2]
data[(memb==5)&(data[,1]=="Asia"),2]
data[(memb==5)&(data[,1]=="Europa"),2]
data[(memb==5)&(data[,1]=="Oceania"),2]


# America
hc = hclust(dist(indicadores[data[,1]=="America",]),"complete")
pdf(file="cluster-ILE-america.pdf",width=15,height=10)
plot(hc, hang = -1,labels=data[data[,1]=="America",2])
dev.off()


# Europe
hc = hclust(dist(indicadores[data[,1]=="Europa",]),"complete")
pdf(file="cluster-ILE-europe.pdf",width=15,height=10)
plot(hc, hang = -1,labels=data[data[,1]=="Europa",2])
dev.off()


# Asia
hc = hclust(dist(indicadores[data[,1]=="Asia",]),"complete")
pdf(file="cluster-ILE-asia.pdf",width=15,height=10)
plot(hc, hang = -1,labels=data[data[,1]=="Asia",2])
dev.off()

# Africa
hc = hclust(dist(indicadores[data[,1]=="Africa",]),"complete")
pdf(file="cluster-ILE-africa.pdf",width=15,height=10)
plot(hc, hang = -1,labels=data[data[,1]=="Africa",2])
dev.off()

# Oceania
hc = hclust(dist(indicadores[data[,1]=="Oceania",]),"complete")
pdf(file="cluster-ILE-oceania.pdf",width=15,height=10)
plot(hc, hang = -1,labels=data[data[,1]=="Oceania",2])
dev.off()