ElectionData/BuildDataset.R at master · nunompmoniz/ElectionData · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
votes <- read.csv("~/Desktop/Open Data/ElectionData/SourceData/votes.csv")
overall <- read.csv("~/Desktop/Open Data/ElectionData/SourceData/overall_results.csv")

votes$time <- as.POSIXct(votes$time)
overall$time <- as.POSIXct(overall$time)

times <- unique(overall$time)

overall <- overall[overall$time %in% times[16:length(times)],]

times <- unique(overall$time)

final <- votes[votes$time == times[length(times)],c("District","Party","Mandates")]

res <- c()

for(t in 1:length(times)) {

  distr <- overall[overall$time==times[t],]
  distr.names <- distr$territoryName

  timedif <- as.numeric(round(difftime(times[t],times[1],units="mins")))

  for(i in 1:length(distr.names)) {

    info.distr <- distr[distr$territoryName==distr.names[i],c(1,3,5:15,18:24)]

    votes.distr <- votes[votes$time==unique(distr$time) & votes$District==distr.names[i],3:7]

    if(nrow(votes.distr)>0) {

      maxElected <- sum(final[final$District==distr.names[i],]$Mandates)

      hondt.m <- matrix(nrow=maxElected,ncol=nrow(votes.distr))
      colnames(hondt.m) <- votes.distr$Party; rownames(hondt.m) <- 1:nrow(hondt.m)

      for(m in 1:maxElected) {
        hondt.m[m,] <- votes.distr$Votes/m
      }

      threshold <- sort(hondt.m,decreasing = TRUE)[maxElected]

      votes.distr["Hondt"] <- apply(hondt.m,2,FUN=function(x) length(x[x>=threshold]))

      ###

      parties <- votes.distr$Party

      for(p in 1:length(parties)) {

        final.mandates <- final[final$District==distr.names[i] & final$Party==parties[p],]$Mandates

        res.row <- cbind(TimeElapsed=timedif,info.distr,votes.distr[votes.distr$Party==parties[p],],FinalMandates=final.mandates)

        res <- rbind(res,res.row)

      }

    }

  }

}

write.csv(res,file="ElectionData.csv",row.names=FALSE)