columns <- read.csv("~/Téléchargements/cours/S9/predictive analytics/communities.var", header=FALSE, na.strings="?")
columns=as.vector(t(columns))
communities <- read.csv("~/Téléchargements/cours/S9/predictive analytics/communities.data", header=FALSE, na.strings="?", col.names=columns)
commu_reduced <- communities[,c(4,18,34,35,38,56,97,102,118,120,121,128)]
print(summary(commu_reduced))
descrCor <- cor(na.omit(commu_reduced[,-1]))

summary(pc.communities <- prcomp(na.omit(commu_reduced[,-1]), scale.=TRUE))
print(pc.communities)
plot(pc.communities,type="l")