Tendencia a grupos!

AndreyMenezes · Feb 28, 2013 · b7c15b5 · b7c15b5
1 parent 713ab2b
commit b7c15b5
Showing 1 changed file with 14 additions and 5 deletions.
diff --git a/US10 - Analise Tendencia A Grupos.R b/US10 - Analise Tendencia A Grupos.R
@@ -1,6 +1,9 @@
-# Codigo que cria as tabelas dos agrupamentos dos dados
+# Codigo que calcula a tendencia de agrupamento dos dados
 # Andrey Menezes - versão 2.0 (Fevereiro 2013)
 
+dados = read.csv("dados/TabelaParaPerfis.csv")
+dados = dados[,2:6]
+
 library("FNN")
 
 hopkins = function(dados, p, objetos, minimos) {
@@ -24,18 +27,20 @@ hopkins = function(dados, p, objetos, minimos) {
 
 hopkins.table = data.frame()
 
+#Calcula a menor distancia de cada linha da tabela para ser utilizado em hopkins
 minimo = knn.dist(dados, k=1, algorithm="kd_tree")
 stats = 0
 
+# Monta a tabela com as estatísticas de hopkins
 col=0
 for(n in 1:10) {
   count=0
   col=col+1
-  for(i in 2:(length(dados)-1)) {
+  for(i in 1:(length(dados)-1)) {
     for(e in (i+1):length(dados)) {
       if(i != e) {
-        objetos = data.frame(runif(15, min(dados[,i]), max(dados[,i])), runif(15, min(dados[,e]), max(dados[,e])))
-        parcial = hopkins(dados[,c(i,e)], 15, objetos, minimo)
+        objetos = data.frame(runif(25, min(dados[,i]), max(dados[,i])), runif(25, min(dados[,e]), max(dados[,e])))
+        parcial = hopkins(dados[,c(i,e)], 25, objetos, minimo)
         if (parcial > stats) {
           stats = parcial
           atrib1 = i
@@ -48,4 +53,8 @@ for(n in 1:10) {
   }
 }
 
-#colnames(hopkins.table) = c()
+colnames(hopkins.table) = c("1_2", "1_3", "1_4", "1_5", "2_3", "2_4", "2_5", "3_4", "3_5", "4_5")
+
+png(filename="BoxSplot Tendencia Agrupamento.png")
+boxplot(hopkins.table)
+dev.off()