US10

AndreyMenezes · Feb 27, 2013 · fccd6c0 · fccd6c0
1 parent e6b9769
commit fccd6c0
Show file tree

Hide file tree

Showing 3 changed files with 57 additions and 7 deletions.
diff --git a/US10 - Analise Tendencia A Grupos.R b/US10 - Analise Tendencia A Grupos.R
@@ -0,0 +1,50 @@
+# Codigo que cria as tabelas dos agrupamentos dos dados
+# Andrey Menezes - versão 2.0 (Fevereiro 2013)
+
+library("FNN")
+
+hopkins = function(dados, p, objetos, minimos) {
+
+  amostra = sample(length(dados[,1]), size=p, replace = FALSE)
+
+  soma_w = 0
+  for(i in 1:length(amostra)) {
+    soma_w = soma_w + minimos[i]
+  }
+
+  soma_u = 0
+  for(i in 1:length(objetos[,1])) {
+    dist = knnx.dist(as.matrix(dados),as.matrix(objetos),2,"kd_tree")
+    soma_u = sum(dist[,2])
+  }
+  if(soma_u > 0) {
+    return(soma_u/(soma_u + soma_w))
+  } else {return(0)}
+}
+
+hopkins.table = data.frame()
+
+minimo = knn.dist(dados, k=1, algorithm="kd_tree")
+stats = 0
+
+col=0
+for(n in 1:10) {
+  count=0
+  col=col+1
+  for(i in 2:(length(dados)-1)) {
+    for(e in (i+1):length(dados)) {
+      if(i != e) {
+        objetos = data.frame(runif(60, min(dados[,i]), max(dados[,i])), runif(60, min(dados[,e]), max(dados[,e])))
+        parcial = hopkins(dados[,c(i,e)], 60, objetos, minimo)
+        if (parcial > stats) {
+          stats = parcial
+          atrib1 = i
+          atrib2 = e
+        }
+        count = count+1
+        hopkins.table[col,count] = parcial
+      }
+    }
+  }
+}
+colnames(hopkins.table) = c("A_Q", "A_C", "A_TA", "A_AS", "A_QS", "A_CS", "Q_C", "Q_TA", "Q_AS", "Q_QS", "Q_CS", "C_TA", "C_AS", "C_QS", "C_CS", "TA_AS", "TA_QS", "TA_CS", "AS_QS", "AS_CS", "QS_CS")
diff --git a/US2 - Automatização.R b/US2 - Automatização.R
@@ -1,5 +1,5 @@
 # Codigo que cria as tabelas dos agrupamentos dos dados
-# Andrey Menezes - versão 2.0 (Fevereiro 2013)
+# Andrey Menezes - versão 2.0 (Fevereiro 2013)
 
 #dados usados
 data_exerc = read.csv("../dados/exercicios-20112.csv")
@@ -9,11 +9,11 @@ prova.um = read.csv("../dados/Prova1.csv")
 prova.dois = read.csv("../dados/Prova2.csv")
 prova.tres = read.csv("../dados/Prova3.csv")
 
-#Script de automatização UC3 e UC4
+#Script de automatiza??o UC3 e UC4
 source("../Processamento/processamento.R")
 
-#Script de automatização UC1
+#Script de automatiza??o UC1
 source("../Processamento/UC1-A Priori.R")
 
-#Script de automatização UC2
+#Script de automatiza??o UC2
 source("../Processamento/UC2.R")
diff --git a/US7 - Analise NumeroSubmissões-Tempo.R b/US7 - Analise NumeroSubmissões-Tempo.R
@@ -1,5 +1,5 @@
-# Codigo para analise de número de submissões pelo tempo
-# Andrey Menezes - versão 2.0 (Fevereiro 2013)
+# Codigo para analise de número de submissões pelo tempo
+# Andrey Menezes - versão 2.0 (Fevereiro 2013)
 
 library(plyr)
 
@@ -11,4 +11,4 @@ colnames(dados) = c("matricula","data")
 submission = count(dados, "data")
 
 a = qplot(submission$data, submission$freq, geom="histogram", xlab="Tempo", ylab="Numero de submissoes")
-a + geom_vline(xintercept = "9/26/2011", colour="green", linetype = "longdash")                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
+a + geom_vline(xintercept = "9/26/2011", colour="green", linetype = "longdash")