Adds state secciones nyt maps

code/plot2018pr.r

@@ -926,8 +926,7 @@ summary(pr.se$camloneg[pr.se$camloneg<0])
 summary(pr.se$cprineg[pr.se$cprineg<0])
 summary(pr.se$cpanneg[pr.se$cpanneg<0])
 
-
-edon <- 1
+edon <- 31
 edo <- edos[edon]
 
 #################
@@ -972,10 +971,12 @@ ed.map <- tmp
 plot(mu.map, lwd = .5, border = "gray")
 axis(1); axis(2)
 
-# amlo v pan
+###################################
+# amlo v pri secciones por estado #
+###################################
 library(graphics)
-#png("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPri.png", width = 20, height = 20, units = "cm", res = 196)
-#pdf("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPri.pdf", width = 10, height = 7)
+#png(paste("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPri-", edo, ".png", sep = ""), width = 20, height = 20, units = "cm", res = 196)
+#pdf(paste("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPri-", edo, ".pdf", sep = ""), width = 10, height = 7)
 par(mar = c(0,0,0,0))
 #plot(shave(mu.map, p = .95), lwd = .5, border = "gray")
 plot(se.map, lwd = .25, border = "lightgray")
@@ -986,30 +987,95 @@ plot(ed.map, lwd = .5, lty = 3, border = "gray", add = TRUE)
 # add arrows
 # start-end of arrows
 #cx <- 150000; cy <- cx/3
-xlength <- (tmp.ranges[2] - tmp.ranges[1]) / 8
-for (i in 1:nrow(se.map@data)){
-    # start-end of arrows
-    #i <- 1 # debug
-    cx <- xlength*se.map$camlopos; cy <- cx/3
-    arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
-           coordinates(se.map)[,1]-cx, coordinates(se.map)[,2]+cy,
-           length = .025, angle = 10,
-           col = brown, lwd = (.1+cx/xlength))
-    cx <- xlength*se.map$cpanpos; cy <- cx/3
-    arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
-           coordinates(se.map)[,1]+cx, coordinates(se.map)[,2]+cy,
-           length = .025, angle = 10,
-           col = colpan, lwd = (.1+cx/xlength))
-}
+xlength <- (tmp.ranges[2] - tmp.ranges[1]) / 10
+#
+# start-end of arrows
+cx <- xlength*se.map$camlopos; cy <- cx/3
+arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
+       coordinates(se.map)[,1]-cx, coordinates(se.map)[,2]+cy,
+       length = .025, angle = 10,
+       col = brown, lwd = (.1+cx/(xlength*.4)))
+cx <- xlength*se.map$cpripos; cy <- cx/3
+arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
+       coordinates(se.map)[,1]+cx, coordinates(se.map)[,2]+cy,
+       length = .025, angle = 10,
+       col = colpri, lwd = (.1+cx/(xlength*.4)))
 # legend
 #xl <-  -11440000; yl <- 2555000
-xl <-  tmp.ranges[2]-(tmp.ranges[2]-tmp.ranges[1])*.1; yl <- tmp.ranges[4]-(tmp.ranges[4]-tmp.ranges[3])*.1
-arrows(xl,yl,xl-xlength*.6,yl+xlength*.6/3,length = .025, angle = 10, col = brown, lwd = 2)
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.20
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.15
+arrows(xl,yl,xl-xlength*.5,yl+xlength*.5/3,length = .025, angle = 10, col = brown, lwd = 2)
 text(xl, yl, pos = 4, labels = "AMLO creció", cex = .75)
-arrows(xl+310000,yl-150000,xl+310000+150000,yl-150000+50000,length = .025, angle = 10, col = colpan, lwd = .75)
-text(xl+310000,yl-150000, pos = 2, labels = "Anaya creció", cex = .75)
-text(xl+180000, yl+280000, labels = "Cambio desde", font = 2)
-text(xl+180000, yl+180000, labels = "2012 en secciones", font = 2)
-text(-13000000, 1550000, labels = "Preparado por Eric Magar con datos del INE (ericmagar.com)", col = "lightgray", pos = 4, cex = .65)
-#dev.off()
+#
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.07
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.18
+arrows(xl,yl,xl+xlength*.5,yl+xlength*.5/3,length = .025, angle = 10, col = colpri, lwd = 2)
+text(xl,yl, pos = 2, labels = "Meade creció", cex = .75)
+#
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.135
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.07
+text(xl, yl, labels = "Cambio desde", font = 2)
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.135
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.10
+text(xl, yl, labels = "2012 en secciones", font = 2)
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.9
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.95
+#
+text(xl, yl, labels = "Preparado por Eric Magar con datos del INE (ericmagar.com)", col = "lightgray", pos = 4, cex = .65)
+dev.off()
+
+
+###################################
+# amlo v pan secciones por estado #
+###################################
+library(graphics)
+#png(paste("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPan-", edo, ".png", sep = ""), width = 20, height = 20, units = "cm", res = 196)
+#pdf(paste("/home/eric/Desktop/MXelsCalendGovt/elecReturns/graph/nytPlusPan-", edo, ".pdf", sep = ""), width = 10, height = 7)
+par(mar = c(0,0,0,0))
+#plot(shave(mu.map, p = .95), lwd = .5, border = "gray")
+plot(se.map, lwd = .25, border = "lightgray")
+tmp.ranges <- par("usr") # keep calculated xy ranges to compute arrow length
+plot(mu.map, lwd = .5,          border = "gray", add = TRUE)
+plot(mu.map, lwd = .5, lty = 3, border = "white", add = TRUE)
+plot(ed.map, lwd = .5, lty = 3, border = "gray", add = TRUE)
+# add arrows
+# start-end of arrows
+#cx <- 150000; cy <- cx/3
+xlength <- (tmp.ranges[2] - tmp.ranges[1]) / 10
+#
+# start-end of arrows
+cx <- xlength*se.map$camlopos; cy <- cx/3
+arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
+       coordinates(se.map)[,1]-cx, coordinates(se.map)[,2]+cy,
+       length = .025, angle = 10,
+       col = brown, lwd = (.1+cx/(xlength*.4)))
+cx <- xlength*se.map$cpanpos; cy <- cx/3
+arrows(coordinates(se.map)[,1],    coordinates(se.map)[,2],
+       coordinates(se.map)[,1]+cx, coordinates(se.map)[,2]+cy,
+       length = .025, angle = 10,
+       col = colpan, lwd = (.1+cx/(xlength*.4)))
+# legend
+#xl <-  -11440000; yl <- 2555000
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.2
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.15
+arrows(xl,yl,xl-xlength*.5,yl+xlength*.5/3,length = .025, angle = 10, col = brown, lwd = 2)
+text(xl, yl, pos = 4, labels = "AMLO creció", cex = .75)
+#
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.07
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.18
+arrows(xl,yl,xl+xlength*.5,yl+xlength*.5/3,length = .025, angle = 10, col = colpan, lwd = 2)
+#arrows(xl+310000,yl-150000,xl+310000+150000,yl-150000+50000,length = .025, angle = 10, col = colpan, lwd = .75)
+text(xl,yl, pos = 2, labels = "Anaya creció", cex = .75)
+#
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.135
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.07
+text(xl, yl, labels = "Cambio desde", font = 2)
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.135
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.10
+text(xl, yl, labels = "2012 en secciones", font = 2)
+xl <- tmp.ranges[2] - (tmp.ranges[2] - tmp.ranges[1])*.9
+yl <- tmp.ranges[4] - (tmp.ranges[4] - tmp.ranges[3])*.95
+text(xl, yl, labels = "Preparado por Eric Magar con datos del INE (ericmagar.com)", col = "lightgray", pos = 4, cex = .65)
+dev.off()
+