Guest User

Untitled

a guest
Apr 21st, 2018
92
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.43 KB | None | 0 0
  1. library(ggplot2)
  2. library(doParallel)
  3.  
  4. cluster <- makeCluster(3)
  5.  
  6. calcParticipationNumbers <- function(data) {
  7. data['megjelent_most'] <- NA
  8. data['megjelent_most_pct'] <- NA
  9.  
  10. cities = foreach(c=unique(data$telepules), .combine=rbind) %dopar% {
  11. city_data <- data[data$telepules==c,]
  12. for (k in unique(city_data$id_k)) {
  13. tmp_k <- city_data[city_data$id_k == k,]
  14. tmp_k <- data.table(tmp_k)
  15. tmp_k[,megjelent_most := megjelent - c(0, megjelent[.I -1])]
  16.  
  17. tmp_k$megjelent_most_pct <- tmp_k$megjelent_most / tmp_k$osszes * 100
  18. city_data[city_data$id_k == k,]$megjelent_most <- tmp_k$megjelent_most
  19. city_data[city_data$id_k == k,]$megjelent_most_pct <- tmp_k$megjelent_most_pct
  20. }
  21.  
  22. data.frame(city_data)
  23. }
  24. cities
  25. }
  26.  
  27. drawPlot <- function(city) {
  28. ggplot(test[test$telepules==city & test$id_ido != 8,], aes(x=id_ido, y=megjelent_most_pct)) + geom_point(alpha=.2) + geom_line(aes(group=korzet), alpha=.1) + theme_bw() + geom_smooth(aes(x=as.numeric(id_ido),y=megjelent_most_pct)) + xlab('Ido id') + ylab('Megjelent (%)') + ggtitle(city) + coord_cartesian(ylim=c(0,40))
  29. }
  30.  
  31. tmp <- reszvetel_df[,c('telepules','id_k','osszes','megjelent','id_ido')]
  32. tmp$megjelent <- as.numeric(tmp$megjelent)
  33. tmp$osszes <- as.numeric(tmp$osszes)
  34. tmp$id_ido <- as.factor(tmp$id_ido)
  35.  
  36. tmp['arany'] <- round(tmp$megjelent/tmp$osszes*100, 2)
  37.  
  38. test <- calcParticipationNumbers(tmp)
  39.  
  40. stopCluster(cluster)
  41.  
  42. drawPlot('Sándorfalva')
Add Comment
Please, Sign In to add comment