King Sejong’s National Referendum on Tax Reform (ggplot)
coop711
2017-05-28
Data
원자료는 세종실록, 요약표는 오기수 교수의 논문에서.
| 지역 | 계급 | 관료여부 | 찬반 | 집계 |
|---|---|---|---|---|
| 서울 | 대신 등 | 관료 | 찬성 | 21 |
| 서울 | 대신 등 | 관료 | 반대 | 194 |
| 서울 | 3품이하현직 | 관료 | 찬성 | 259 |
| 서울 | 3품이하현직 | 관료 | 반대 | 393 |
| 서울 | 3품이하전직 | 관료 | 찬성 | 443 |
| 서울 | 3품이하전직 | 관료 | 반대 | 117 |
| 유후사 | 품관촌민 | 품관촌민 | 찬성 | 1123 |
| 유후사 | 품관촌민 | 품관촌민 | 반대 | 71 |
| 경기 | 수령 | 관료 | 찬성 | 29 |
| 경기 | 수령 | 관료 | 반대 | 5 |
| 경기 | 품관촌민 | 품관촌민 | 찬성 | 17076 |
| 경기 | 품관촌민 | 품관촌민 | 반대 | 236 |
| 평안 | 대신 등 | 관료 | 반대 | 1 |
| 평안 | 수령 | 관료 | 찬성 | 6 |
| 평안 | 수령 | 관료 | 반대 | 35 |
| 평안 | 품관촌민 | 품관촌민 | 찬성 | 1326 |
| 평안 | 품관촌민 | 품관촌민 | 반대 | 28474 |
| 황해 | 수령 | 관료 | 찬성 | 17 |
| 황해 | 수령 | 관료 | 반대 | 17 |
| 황해 | 품관촌민 | 품관촌민 | 찬성 | 4454 |
| 황해 | 품관촌민 | 품관촌민 | 반대 | 15601 |
| 충청 | 대신 등 | 관료 | 반대 | 2 |
| 충청 | 수령 | 관료 | 찬성 | 35 |
| 충청 | 수령 | 관료 | 반대 | 26 |
| 충청 | 품관촌민 | 품관촌민 | 찬성 | 6982 |
| 충청 | 품관촌민 | 품관촌민 | 반대 | 14013 |
| 강원 | 수령 | 관료 | 찬성 | 5 |
| 강원 | 수령 | 관료 | 반대 | 10 |
| 강원 | 품관촌민 | 품관촌민 | 찬성 | 939 |
| 강원 | 품관촌민 | 품관촌민 | 반대 | 6888 |
| 함길 | 대신 등 | 관료 | 반대 | 1 |
| 함길 | 수령 | 관료 | 찬성 | 3 |
| 함길 | 수령 | 관료 | 반대 | 14 |
| 함길 | 품관촌민 | 품관촌민 | 찬성 | 75 |
| 함길 | 품관촌민 | 품관촌민 | 반대 | 7387 |
| 경상 | 수령 | 관료 | 찬성 | 55 |
| 경상 | 수령 | 관료 | 반대 | 16 |
| 경상 | 품관촌민 | 품관촌민 | 찬성 | 36262 |
| 경상 | 품관촌민 | 품관촌민 | 반대 | 377 |
| 전라 | 대신 등 | 관료 | 반대 | 2 |
| 전라 | 수령 | 관료 | 찬성 | 42 |
| 전라 | 수령 | 관료 | 반대 | 12 |
| 전라 | 품관촌민 | 품관촌민 | 찬성 | 29505 |
| 전라 | 품관촌민 | 품관촌민 | 반대 | 257 |
전체 찬반
# pander(format(tbl.vote, big.mark = ","))
kable(t(as.matrix(format(tbl.vote,
big.mark = ","))),
caption = "전체 찬반", align = 'c')| 찬성 | 반대 |
|---|---|
| 98,657 | 74,149 |
kable(t(as.table(format(tbl.vote,
big.mark = ","))),
caption = "전체 찬반",
align = 'c')| 찬성 | 반대 |
|---|---|
| 98,657 | 74,149 |
kable(t(as.table(format(prop.table(tbl.vote)*100,
digits = 3,
nsmall = 1))),
caption = "전체 찬반(%)",
align = 'c')| 찬성 | 반대 |
|---|---|
| 57.1 | 42.9 |
원형그래프
source("./theme_kr.R")
# str(theme.kr)
str(tbl.vote)## 'xtabs' int [1:2(1d)] 98657 74149
## - attr(*, "dimnames")=List of 1
## ..$ vote: chr [1:2] "찬성" "반대"
## - attr(*, "call")= language xtabs(formula = counts ~ vote, data = sejong.poll.kr)kable(t(as.matrix(tbl.vote)))| 찬성 | 반대 |
|---|---|
| 98657 | 74149 |
vote.df <- as.data.frame(tbl.vote)
n.vote <- length(levels(vote.df[, 1]))
y.coord <- vote.df$Freq/2 + c(0, cumsum(head(vote.df$Freq, -1)))
pie.label <- paste(levels(vote.df$vote),
format(vote.df$Freq,
big.mark = ","),
sep = "\n")
pie.label## [1] "찬성\n98,657" "반대\n74,149"kable(vote.df,
align = rep("c", 2))| vote | Freq |
|---|---|
| 찬성 | 98657 |
| 반대 | 74149 |
str(vote.df)## 'data.frame': 2 obs. of 2 variables:
## $ vote: Factor w/ 2 levels "찬성","반대": 1 2
## $ Freq: int 98657 74149
p1 <- ggplot(data = vote.df,
mapping = aes(x = "",
y = Freq,
fill = vote))
p2 <- p1 +
geom_bar(width = 1,
stat = "identity",
position = position_stack(reverse = TRUE))
p2
p2.2 <- p1 +
geom_bar(width = 1,
stat = "identity",
position = "dodge")
p2.2
p3.2 <- p2.2 +
theme_bw() + theme.kr
p3.2
p4.2 <- p3.2 +
scale_x_discrete(name = "찬반") +
scale_y_continuous(name = "집계",
breaks = vote.df$Freq,
labels = format(vote.df$Freq, big.mark = ","))
p4.2
p5.2 <- p4.2 +
scale_fill_manual(name = "찬반",
values = rainbow(n.vote)[n.vote:1],
labels = vote.df$vote)
p5.2
p6.2 <- p5.2 +
geom_text(aes(y = vote.df$Freq/2),
label = format(rev(vote.df$Freq), big.mark = ","),
# family = "HCR Dotum LVT",
position = position_dodge(width = 1))
p6.2
p7.2 <- p6.2 +
theme(legend.position = c(0.25, 0.75))
p7.2
p8.2 <- p5.2 +
guides(fill = "none") +
geom_text(aes(y = vote.df$Freq/2),
family = "HCR Dotum LVT",
label = rev(pie.label),
position = position_dodge(width = 1))
p8.2
# ggsave("../pics/sejong_geom_bar_total_ggplot.png", p7.2, dpi = 72)
p3 <- p2 +
theme_bw() + theme.kr
p3 
p4 <- p3 +
scale_x_discrete(name = "찬반") +
scale_y_continuous(name = "집계",
breaks = cumsum(vote.df$Freq),
labels = format(cumsum(vote.df$Freq), big.mark = ","))
p4
p5 <- p4 +
scale_fill_manual(name = "찬반", values = rainbow(2)[2:1], labels = vote.df$vote)
p5
p6 <- p5 +
geom_text(aes(y = y.coord),
label = format(vote.df$Freq, big.mark = ","),
# family = "HCR Dotum LVT",
position = position_identity())
p6
# ggsave("../pics/sejong_geom_bar_total_ggplot_stack.png", p6, dpi = 72)
p7 <- p6 +
theme(legend.position = c(0.2, 0.75))
p7
p8 <- p5 +
guides(fill = "none") +
geom_text(aes(y = y.coord),
label = pie.label,
family = "HCR Dotum LVT",
position = position_identity())
p8
pie.total.1 <- p2 +
coord_polar(theta = "y",
start = 3 * pi /2 ,
direction = -1)
pie.total.1
pie.total.2 <- pie.total.1 +
scale_y_continuous(name = "",
breaks = NULL) +
scale_x_discrete(name = "")
pie.total.2
pie.total.3 <- pie.total.2 +
scale_fill_manual(name = "",
values = rainbow(n.vote)[n.vote:1])
pie.total.3
pie.total.4 <- pie.total.3 +
theme_void(base_family = "HCR Dotum LVT")
pie.total.4
pie.total.5 <- pie.total.4 +
guides(fill = "none")
pie.total.5
pie.total.6 <- pie.total.5 +
geom_text(aes(y = y.coord),
label = pie.label,
family = "HCR Dotum LVT",
position = position_identity())
pie.total.6
pie.total.7 <- pie.total.6 +
ggtitle("전체 찬반") +
theme(plot.margin = unit(c(1, 1, 1.5, 1), "lines"), plot.title = element_text(hjust = 0.5, size = 20))
pie.total.7
ggsave("../pics/sejong_total_pie_ggplot.png", pie.total.7, dpi = 72)
pie.gg()
pie.gg <- function(df, ggtitle = "", font.family = ""){
n <- length(names(df))
y.coord <- df$Freq / 2 + c(0, cumsum(head(df$Freq, -1)))
pie.label <- paste(levels(df$vote), format(df$Freq, big.mark = ","),
sep = "\n")
p1 <- ggplot(df, aes(x = "",
y = Freq,
fill = vote))
p2 <- p1 +
geom_bar(width = 1,
stat = "identity",
position = position_stack(reverse = TRUE))
pie.1 <- p2 +
coord_polar(theta = "y",
start = 3 * pi / 2,
direction = -1)
pie.2 <- pie.1 +
scale_y_continuous(name = "",
breaks = NULL) +
scale_x_discrete(name = "")
pie.3 <- pie.2 +
scale_fill_manual(name = "",
values = rainbow(n)[n:1])
pie.4 <- pie.3 +
theme_void(base_family = font.family)
pie.5 <- pie.4 +
guides(fill = "none")
pie.6 <- pie.5 +
geom_text(aes(y = y.coord),
label = pie.label,
family = font.family,
position = position_identity())
pie.7 <- pie.6 +
ggtitle(ggtitle) +
theme(plot.margin = unit(c(1, 1, 1.5, 1), "lines"),
plot.title = element_text(hjust = 0.5))
return(pie.7)
}
dump("pie.gg",
file = "./pie.gg.R")
pie.gg(vote.df,
ggtitle = "전체 찬반",
font.family = "HCR Dotum LVT")
계급 및 지역별 찬반
계급별 찬반
kable(format(vote.class,
big.mark = ","),
align = "r",
caption = "계급별 찬반")| 대신 등 | 3품이하현직 | 3품이하전직 | 수령 | 품관촌민 | |
|---|---|---|---|---|---|
| 찬성 | 21 | 259 | 443 | 192 | 97,742 |
| 반대 | 200 | 393 | 117 | 135 | 73,304 |
kable(format(prop.table(vote.class,
margin = 2)*100,
digits = 3,
nsmall = 1),
align = "r",
caption = "계급별 찬반(%)")| 대신 등 | 3품이하현직 | 3품이하전직 | 수령 | 품관촌민 | |
|---|---|---|---|---|---|
| 찬성 | 9.5 | 39.7 | 79.1 | 58.7 | 57.1 |
| 반대 | 90.5 | 60.3 | 20.9 | 41.3 | 42.9 |
품관촌민 별도 분석
품관촌민의 수효가 상대적으로 많아서
kable(format(vote.class.2,
big.mark = ","),
align = rep("r", 2),
caption = "품관촌민")| 관료 | 품관촌민 | |
|---|---|---|
| 찬성 | 915 | 97,742 |
| 반대 | 845 | 73,304 |
소계를 교차표의 주변에 계산
kable(format(vote.class.2.am,
big.mark = ","),
caption = "관료와 품관촌민(소계)",
align = rep("r", 3))| 관료 | 품관촌민 | 계 | |
|---|---|---|---|
| 찬성 | 915 | 97,742 | 98,657 |
| 반대 | 845 | 73,304 | 74,149 |
| 계 | 1,760 | 171,046 | 172,806 |
백분율을 계산하여 주변에 집계.
kable(format(prop.table(vote.class.2,
margin = 2)*100,
digits = 3, nsmall = 1),
caption = "관료와 품관촌민(%)",
align = rep("r", 3))| 관료 | 품관촌민 | |
|---|---|---|
| 찬성 | 52.0 | 57.1 |
| 반대 | 48.0 | 42.9 |
data frame 변환
vote.class.2.df <- as.data.frame(vote.class.2)
v.names.class.kr <- c("찬반", "계급", "집계")
kable(format(vote.class.2.df,
big.mark = ","),
align = c('c', 'c', 'r'),
col.names = v.names.class.kr,
caption = "관료와 품관촌민")| 찬반 | 계급 | 집계 |
|---|---|---|
| 찬성 | 관료 | 915 |
| 반대 | 관료 | 845 |
| 찬성 | 품관촌민 | 97,742 |
| 반대 | 품관촌민 | 73,304 |
vote.bureaus.df <- subset(vote.class.2.df,
vote.class.2.df$class.2 == "관료",
select = c("vote", "Freq"))
kable(format(vote.bureaus.df,
big.mark = ","),
align = 'r',
col.names = v.names.class.kr[-2],
caption = "관료의 찬반")| 찬반 | 집계 |
|---|---|
| 찬성 | 915 |
| 반대 | 845 |
vote.commons.df <- subset(vote.class.2.df,
vote.class.2.df$class.2 == "품관촌민",
select = c("vote", "Freq"))
kable(format(vote.commons.df,
big.mark = ","),
align = 'r',
row.names = FALSE,
col.names = v.names.class.kr[-2],
caption = "품관촌민의 찬반")| 찬반 | 집계 |
|---|---|
| 찬성 | 97,742 |
| 반대 | 73,304 |
원형그래프
pie.bureaus <- pie.gg(vote.bureaus.df,
ggtitle = "관료의 찬반",
font.family = "HCR Dotum LVT")
pie.commons <- pie.gg(vote.commons.df,
ggtitle = "품관촌민의 찬반",
font.family = "HCR Dotum LVT")
pies.grid <- grid.arrange(pie.bureaus,
pie.commons,
ncol = 2,
top = "")
pies.grid## TableGrob (2 x 2) "arrange": 3 grobs
## z cells name grob
## 1 1 (2-2,1-1) arrange gtable[layout]
## 2 2 (2-2,2-2) arrange gtable[layout]
## 3 3 (1-1,1-2) arrange text[GRID.text.1265]ggsave("../pics/sejong_bureaus_commons_pie_ggplot.png", pies.grid, width = 8, height = 4, dpi = 72)
지역별 찬반
관료와 품관촌민
kable(vote.region.bureaus,
caption = "관료들의 지역별 찬반")| 서울 | 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|---|
| 찬성 | 723 | 29 | 6 | 17 | 35 | 5 | 3 | 55 | 42 |
| 반대 | 704 | 5 | 36 | 17 | 28 | 10 | 15 | 16 | 14 |
kable(format(prop.table(vote.region.bureaus,
margin = 2) * 100,
digits = 3,
nsmall = 1),
align = rep("r", 9),
caption = "관료들의 지역별 찬반(%)")| 서울 | 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|---|
| 찬성 | 50.7 | 85.3 | 14.3 | 50.0 | 55.6 | 33.3 | 16.7 | 77.5 | 75.0 |
| 반대 | 49.3 | 14.7 | 85.7 | 50.0 | 44.4 | 66.7 | 83.3 | 22.5 | 25.0 |
kable(format(vote.region.commons,
big.mark = ","),
align = rep("r", 9),
caption = "품관촌민들의 지역별 찬반")| 유후사 | 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|---|
| 찬성 | 1,123 | 17,076 | 1,326 | 4,454 | 6,982 | 939 | 75 | 36,262 | 29,505 |
| 반대 | 71 | 236 | 28,474 | 15,601 | 14,013 | 6,888 | 7,387 | 377 | 257 |
kable(format(prop.table(vote.region.commons,
margin = 2) * 100,
digits = 1,
nsmall = 1),
align = rep("r", 9),
caption = "품관촌민들의 지역별 찬반(%)")| 유후사 | 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|---|
| 찬성 | 94.1 | 98.6 | 4.4 | 22.2 | 33.3 | 12.0 | 1.0 | 99.0 | 99.1 |
| 반대 | 5.9 | 1.4 | 95.6 | 77.8 | 66.7 | 88.0 | 99.0 | 1.0 | 0.9 |
서울의 찬반
kable(vote.seoul.class,
caption = "서울의 찬반")| 대신 등 | 3품이하현직 | 3품이하전직 | |
|---|---|---|---|
| 찬성 | 21 | 259 | 443 |
| 반대 | 194 | 393 | 117 |
kable(format(prop.table(vote.seoul.class,
margin = 2) * 100,
digits = 1,
nsmall = 1),
align = rep("r", 3),
caption = "서울의 찬반(%)")| 대신 등 | 3품이하현직 | 3품이하전직 | |
|---|---|---|---|
| 찬성 | 9.8 | 39.7 | 79.1 |
| 반대 | 90.2 | 60.3 | 20.9 |
막대그래프
Stack
vote.seoul.df <- as.data.frame(vote.seoul.class)
n.fill <- length(levels(vote.seoul.df[, 1]))
x.stack <- vote.seoul.df[, 2]
y.stack <- unlist(tapply(vote.seoul.df$Freq,
vote.seoul.df$class,
function(x){x / 2 +
c(0, cumsum(head(x, -1)))}))
y.breaks <- unlist(tapply(vote.seoul.df$Freq,
vote.seoul.df$class,
cumsum))
y.label <- format(y.breaks, big.mark = ",")
b1.seoul <- ggplot(vote.seoul.df,
aes(x = x.stack,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity",
position = position_stack(reverse = TRUE))
b1.seoul
b2.seoul <- b1.seoul +
theme_bw(base_family = "HCR Dotum LVT") +
theme.kr +
scale_x_discrete(name = "계급") +
scale_y_continuous(name = "집계",
breaks = y.breaks,
labels = y.label) +
scale_fill_manual(name = "찬반",
values = rainbow(n.fill)[n.fill:1],
labels = vote.seoul.df$vote)
b2.seoul 
b3.seoul <- b2.seoul +
# geom_text(aes(y = y.stack/2), label = vote.seoul.df$Freq)
geom_text(aes(y = y.stack),
label = vote.seoul.df$Freq,
position = "identity")
b3.seoul
ggsave("../pics/sejong_seoul_barplot_stack_ggplot.png", b3.seoul, width = 6, height = 4.5, dpi = 72)
Dodge
n.fill <- length(levels(vote.seoul.df[, 1]))
b1.seoul.dodge <- ggplot(vote.seoul.df,
aes(x = x.stack,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity", position = "dodge")
b1.seoul.dodge
b2.seoul.dodge <- b1.seoul.dodge +
theme_bw(base_family = "HCR Dotum LVT") +
theme.kr +
scale_x_discrete(name = "계급") +
scale_y_continuous(name = "집계",
breaks = vote.seoul.df$Freq,
labels = vote.seoul.df$Freq) +
scale_fill_manual(name = "찬반",
values = rainbow(n.fill)[n.fill:1],
labels = vote.seoul.df$vote)
b2.seoul.dodge 
# N <- nrow(vote.seoul.df)
# index <- as.vector(matrix(1:N, nrow = 2)[2:1, ])
y.label <- unlist(tapply(vote.seoul.df$Freq,
vote.seoul.df$class,
rev))
b3.seoul.dodge <- b2.seoul.dodge +
geom_text(aes(y = vote.seoul.df$Freq/2),
# label = vote.seoul.df[index, "Freq"],
label = y.label,
position = position_dodge(width = 0.9)) +
ggtitle("서울의 찬반")
b3.seoul.dodge
# ggsave("../pics/sejong_seoul_barplot_dodge_ggplot.png", b3.seoul.dodge, width = 6, height = 4.5, dpi = 72)
Fill
n.fill <- length(levels(vote.seoul.df[, 1]))
y.fill <- unlist(tapply(vote.seoul.df$Freq,
x.stack,
function(x){cumsum(x)/sum(x)}))
p.fill <- unlist(tapply(vote.seoul.df$Freq,
x.stack,
function(x){(x / 2 + c(0, cumsum(head(x, -1))))/sum(x)}))
b1.seoul.fill <- ggplot(vote.seoul.df,
aes(x = x.stack,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity",
position = position_fill(reverse = TRUE))
b1.seoul.fill
b2.seoul.fill <- b1.seoul.fill +
theme_bw(base_family = "HCR Dotum LVT") +
theme.kr +
scale_x_discrete(name = "계급") +
scale_y_continuous(name = "집계(%)",
breaks = y.fill,
labels = format(y.fill * 100,
digits = 2,
nsmall = 1)) +
scale_fill_manual(name = "찬반",
values = rainbow(n.fill)[n.fill:1],
labels = vote.seoul.df$vote)
b2.seoul.fill 
b3.seoul.fill <- b2.seoul.fill +
geom_text(aes(y = p.fill),
label = vote.seoul.df$Freq,
position = "identity") +
ggtitle("서울의 찬반")
b3.seoul.fill
# ggsave("../pics/sejong_seoul_barplot_fill_ggplot.png", b3.seoul.fill, width = 6, height = 4.5, dpi = 72)
barplot.gg
Stack
barplot.gg.stack <- function(df,
base_family = "",
ggtitle = "",
xlab = "",
ylab = "",
fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y <- unlist(tapply(df$Freq,
x,
function(x){x / 2 + c(0, cumsum(head(x, -1)))}))
y.breaks <- unlist(tapply(df$Freq,
x,
cumsum))
y.label <- format(y.breaks, big.mark = ",")
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
diff(y.breaks.sort) < delta
index <- which(diff(y.breaks.sort) > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
y.label <- format(y.breaks, big.mark = ",")
b1 <- ggplot(df,
aes(x = x,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity",
position = position_stack(reverse = TRUE))
b2 <- b1 +
theme_bw(base_family = base_family) +
# theme.kr +
scale_x_discrete(name = xlab) +
scale_y_continuous(name = ylab,
breaks = y.breaks,
labels = y.label) +
scale_fill_manual(name = fill.name,
values = rainbow(n.fill)[n.fill:1],
labels = df$vote,
guide = guide_legend())
b3 <- b2 +
geom_text(aes(y = y),
label = format(df$Freq,
big.mark = ","),
position = "identity") +
ggtitle(ggtitle)
return(b3)
}
Dodge
barplot.gg.dodge <- function(df,
base_family = "",
ggtitle = "",
xlab = "",
ylab = "",
fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y.dodge <- df$Freq
b1 <- ggplot(df,
aes(x = x,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity",
position = "dodge")
b2 <- b1 +
theme_bw(base_family = base_family) +
# theme.kr +
scale_x_discrete(name = xlab) +
scale_y_continuous(name = ylab,
breaks = y.dodge,
labels = format(y.dodge, big.mark = ",")) +
scale_fill_manual(name = fill.name,
values = rainbow(n.fill)[n.fill:1],
labels = df$vote)
# N <- nrow(df)
# index <- as.vector(matrix(1:N, nrow = 2)[2:1, ])
y.label <- unlist(tapply(df$Freq,
df[, 2],
rev))
b3 <- b2 +
geom_text(aes(y = y.dodge/2),
# label = format(df[index, "Freq"], big.mark = ","),
label = y.label,
position = position_dodge(width = 0.9)) +
ggtitle(ggtitle)
return(b3)
}
Fill
barplot.gg.fill <- function(df,
base_family = "",
ggtitle = "",
xlab = "",
ylab = "",
fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y.fill <- unlist(tapply(df$Freq,
x,
function(x){cumsum(x)/sum(x)}))
p.fill <- unlist(tapply(df$Freq,
x,
function(x){(x / 2 + c(0, cumsum(head(x, -1))))/sum(x)}))
b1 <- ggplot(df,
aes(x = x,
y = Freq,
fill = vote)) +
geom_bar(stat = "identity",
position = position_fill(reverse = TRUE))
b2 <- b1 +
theme_bw(base_family = base_family) +
# theme.kr +
scale_x_discrete(name = xlab) +
scale_y_continuous(name = ylab,
breaks = y.fill,
labels = format(y.fill * 100,
digits = 2,
nsmall = 1)) +
scale_fill_manual(name = fill.name,
values = rainbow(n.fill)[n.fill:1],
labels = df$vote,
guide = guide_legend())
b3 <- b2 +
geom_text(aes(y = p.fill),
label = format(df$Freq, big.mark = ","),
position = "identity") +
ggtitle(ggtitle)
return(b3)
}
barplot.gg
barplot.gg <- function(x,
position,
base_family = "",
ggtitle = "",
xlab = "",
ylab = "",
fill.name = ""){
switch(position,
stack = barplot.gg.stack(x,
base_family = base_family,
ggtitle = ggtitle,
xlab = xlab,
ylab = ylab,
fill.name = fill.name),
dodge = barplot.gg.dodge(x,
base_family = base_family,
ggtitle = ggtitle,
xlab = xlab,
ylab = ylab,
fill.name = fill.name),
fill = barplot.gg.fill(x,
base_family = base_family,
ggtitle = ggtitle,
xlab = xlab,
ylab = ylab,
fill.name = fill.name))
}
(bar_seoul_stack <- barplot.gg(vote.seoul.df,
position = "stack",
base_family = "HCR Dotum LVT",
ggtitle = "서울의 찬반",
xlab = "계급",
ylab = "집계",
fill.name = "찬반"))
# ggsave("../pics/sejong_seoul_barplot_stack_ggplotv2.png", bar_seoul_stack, width = 6, height = 4.5, dpi = 72)
(bar_seoul_dodge <- barplot.gg(vote.seoul.df,
position = "dodge",
base_family = "HCR Dotum LVT",
ggtitle = "서울의 찬반",
xlab = "계급",
ylab = "집계",
fill.name = "찬반"))
# ggsave("../pics/sejong_seoul_barplot_dodge_ggplotv2.png", bar_seoul_dodge, width = 6, height = 4.5, dpi = 72)
(bar_seoul_fill <- barplot.gg(vote.seoul.df,
position = "fill",
base_family = "HCR Dotum LVT",
ggtitle = "서울의 찬반",
xlab = "계급",
ylab = "집계(%)",
fill.name = "찬반"))
# ggsave("../pics/sejong_seoul_barplot_fill_ggplotv2.png", bar_seoul_fill, width = 6, height = 4.5, dpi = 72)
dump(c("barplot.gg","barplot.gg.stack", "barplot.gg.dodge", "barplot.gg.fill"),
file = "./barplot.gg.R")
mosaic plot
Data
vote.seoul.class.df <- as.data.frame(vote.seoul.class)
vote.seoul.class.sum <- tapply(vote.seoul.class.df$Freq,
vote.seoul.class.df[, 2],
sum)
vote.seoul.class.p.m <- prop.table(vote.seoul.class.sum)
vote.seoul.class.p <- prop.table(vote.seoul.class)
vote.seoul.class.p.2 <- prop.table(vote.seoul.class,
margin = 2)
vote.seoul.class.p.df <- as.data.frame(vote.seoul.class.p)
vote.seoul.class.p.df$width <- vote.seoul.class.p.m[match(vote.seoul.class.p.df$class,
names(vote.seoul.class.p.m))]
vote.seoul.class.p.2.df <- as.data.frame(vote.seoul.class.p.2)
vote.seoul.class.p.df$height <- vote.seoul.class.p.2.df$Freq
vote.seoul.class.p.df$label.height <- unlist(tapply(vote.seoul.class.p.df$height,
vote.seoul.class.p.df$class,
function(x) x / 2 + c(0, cumsum(head(x, -1)))))
vote.seoul.class.p.df$y.breaks <- unlist(tapply(vote.seoul.class.p.df$height,
vote.seoul.class.p.df$class,
cumsum))
x.center <- vote.seoul.class.p.m / 2 + c(0, cumsum(head(vote.seoul.class.p.m, -1)))
# x.center <- (cumsum(vote.seoul.class.p.m) +
# c(0, head(cumsum(vote.seoul.class.p.m), -1)))/2
vote.seoul.class.p.df$center <- x.center[match(vote.seoul.class.p.df$class,
names(x.center))]
m1 <- ggplot(vote.seoul.class.p.df,
aes(x = center,
y = height,
width = width)) +
geom_bar(aes(fill = vote),
stat = "identity",
col = "white",
size = 2,
position = position_stack(reverse = TRUE))
m1
m2 <- m1 +
theme_bw(base_family = "HCR Dotum LVT")
m2
m3 <- m2 +
geom_text(aes(x = center,
y = 1.05),
label = vote.seoul.class.p.df$class,
family = "HCR Dotum LVT")
m3
m4 <- m3 +
geom_text(aes(x = center,
y = label.height),
label = format(vote.seoul.class.df$Freq,
big.mark = ","),
position = position_identity())
m4
x.breaks <- c(0, ifelse(cumsum(vote.seoul.class.p.m) < 0.1, 0.0, cumsum(vote.seoul.class.p.m)))
x.label <- format(x.breaks * 100,
digits = 3,
nsmall = 1)
y.breaks <- vote.seoul.class.p.df$y.breaks
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
index <- which(diff(y.breaks.sort) > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
y.label <- format(y.breaks * 100,
digits = 2,
nsmall = 1)
m5 <- m4 +
scale_x_continuous(name = "계급(누적 %)",
breaks = x.breaks,
label = x.label) +
scale_y_continuous(name = "찬반(%)",
breaks = y.breaks,
label = y.label) +
scale_fill_manual(name = "찬반",
values = rainbow(2)[2:1],
labels = vote.seoul.df$vote,
guide = guide_legend()) +
ggtitle("서울의 찬반")
m5
mosaic_gg
mosaic_gg <- function(tbl, base_family = "",
ggtitle = "",
xlab = "",
ylab = "",
fill.name = ""){
tbl.df <- as.data.frame(tbl)
N <- length(levels(tbl.df[, 1]))
tbl.sum <- tapply(tbl.df$Freq, tbl.df[, 2], sum)
tbl.p.m <- prop.table(tbl.sum)
tbl.p <- prop.table(tbl)
tbl.p.2 <- prop.table(tbl, margin = 2)
tbl.p.df <- as.data.frame(tbl.p)
tbl.p.df$width <- tbl.p.m[match(tbl.p.df[, 2], names(tbl.p.m))]
tbl.p.df$height <- as.data.frame(tbl.p.2)$Freq
tbl.p.df$label.height <- unlist(tapply(tbl.p.df$height, tbl.p.df[, 2], function(x) x / 2 + c(0, cumsum(head(x, -1)))))
tbl.p.df$y.breaks <- unlist(tapply(tbl.p.df$height, tbl.p.df[, 2], cumsum))
x.center <- tbl.p.m / 2 + c(0, cumsum(head(tbl.p.m, -1)))
# x.center <- (cumsum(tbl.p.m) + c(0, head(cumsum(tbl.p.m), -1)))/2
tbl.p.df$center <- x.center[match(tbl.p.df[, 2], names(x.center))]
m1 <- ggplot(tbl.p.df, aes(x = center, y = height, width = width)) +
geom_bar(aes(fill = tbl.df[1]),
stat = "identity",
col = "white",
size = 1,
position = position_stack(reverse = TRUE))
m2 <- m1 +
theme_bw(base_family = base_family)
m3 <- m2 +
geom_text(aes(x = center, y = 1.05),
label = tbl.p.df[, 2],
family = base_family)
m4 <- m3 +
geom_text(aes(x = center, y = label.height),
label = format(tbl.df$Freq, big.mark = ","),
position = position_identity())
x.breaks <- c(0, ifelse(cumsum(tbl.p.m) < 0.1, 0.0, cumsum(tbl.p.m)))
x.label <- format(x.breaks * 100,
digits = 3,
nsmall = 1)
y.breaks <- tbl.p.df$y.breaks
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
diff(y.breaks.sort) < delta
index <- which(diff(y.breaks.sort) > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
y.label <- format(y.breaks * 100,
digits = 2,
nsmall = 1)
m5 <- m4 +
scale_x_continuous(name = xlab,
breaks = x.breaks,
label = x.label) +
scale_y_continuous(name = ylab,
breaks = y.breaks,
label = y.label) +
scale_fill_manual(name = fill.name,
values = rainbow(N)[N:1],
labels = tbl.df[[1]],
guide = guide_legend()) +
ggtitle(ggtitle) +
theme(plot.margin = unit(c(1, 2, 1, 1), "lines"))
return(m5)
}
dump(list = "mosaic_gg", file = "./mosaic_gg.R")source("./mosaic_gg.R")
m5.seoul <- mosaic_gg(vote.seoul.class,
base_family = "HCR Dotum LVT",
ggtitle = "서울의 찬반",
xlab = "계급",
ylab = "찬반",
fill.name = "찬반")
m5.seoul
ggsave("../pics/sejong_seoul_mosaic_ggplot.png", m5.seoul, width = 6, height = 4.5, dpi = 72)
지역별 관료들의 찬반(서울 제외)
kable(vote.region.bureaus[, -1],
caption = "지역별 관료들의 찬반(서울 제외)")| 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|
| 찬성 | 29 | 6 | 17 | 35 | 5 | 3 | 55 | 42 |
| 반대 | 5 | 36 | 17 | 28 | 10 | 15 | 16 | 14 |
kable(format(prop.table(vote.region.bureaus[, -1],
margin = 2)*100,
digits = 1,
nsmall = 1),
align = "r",
caption = "지역별 관료들의 찬반(서울 제외, %)")| 경기 | 평안 | 황해 | 충청 | 강원 | 함길 | 경상 | 전라 | |
|---|---|---|---|---|---|---|---|---|
| 찬성 | 85.3 | 14.3 | 50.0 | 55.6 | 33.3 | 16.7 | 77.5 | 75.0 |
| 반대 | 14.7 | 85.7 | 50.0 | 44.4 | 66.7 | 83.3 | 22.5 | 25.0 |
막대그래프
vote.region.bureaus.df <- as.data.frame(vote.region.bureaus[, -1])
barplot.gg(vote.region.bureaus.df,
position = "stack",
base_family = "HCR Dotum LVT",
ggtitle = "지방 관료들의 찬반",
xlab = "지역",
ylab = "집계",
fill.name = "찬반")
ggsave("../pics/sejong_bureaus_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)
mosaic plot
m5.bureaus <- mosaic_gg(vote.region.bureaus[, -1],
base_family = "HCR Dotum LVT",
ggtitle = "지방관료의 찬반",
xlab = "지역별 관료들의 비중(누적 %)",
ylab = "찬반(%)",
fill.name = "찬반")
m5.bureaus
ggsave("../pics/sejong_bureaus_mosaic_ggplot.png", m5.bureaus, width = 9, height = 4.5, dpi = 72)
품관촌민들의 지역별 찬반
막대그래프
vote.region.commons.df <- as.data.frame(vote.region.commons)
barplot.gg(vote.region.commons.df,
position = "stack",
base_family = "HCR Dotum LVT",
ggtitle = "품관촌민의 찬반",
xlab = "지역",
ylab = "집계",
fill.name = "찬반")
ggsave("../pics/sejong_commons_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)
mosaic plot
m5.commons <- mosaic_gg(vote.region.commons,
base_family = "HCR Dotum LVT",
ggtitle = "품관촌민의 찬반",
xlab = "지역별 품관촌민의 비중(누적%)",
ylab = "찬반(%)",
fill.name = "찬반")
m5.commons
ggsave("../pics/sejong_commons_mosaic_ggplot.png", m5.commons, width = 12, height = 4.5, dpi = 72)
충청도
관료들의 찬반과 품관촌민의 찬반이 다른 곳.
kable(format(vote.chung.class,
big.mark = ","),
caption = "충청도 사례",
align = "r")| 대신 등 | 수령 | 품관촌민 | |
|---|---|---|---|
| 찬성 | 0 | 35 | 6,982 |
| 반대 | 2 | 26 | 14,013 |
kable(format(prop.table(vote.chung.class,
margin = 2)*100,
digits = 3, nsmall = 1),
caption = "충청도 사례",
align = "r")| 대신 등 | 수령 | 품관촌민 | |
|---|---|---|---|
| 찬성 | 0.0 | 57.4 | 33.3 |
| 반대 | 100.0 | 42.6 | 66.7 |
막대그래프
vote.chung.class.df <- as.data.frame(vote.chung.class)
barplot.gg(vote.chung.class.df,
position = "fill",
base_family = "HCR Dotum LVT",
ggtitle = "충청도의 찬반",
xlab = "계급",
ylab = "집계(%)",
fill.name = "찬반")
ggsave("../pics/sejong_chungcheong_barplot_fill_ggplot.png", width = 6, height = 4.5, dpi = 72)
mosaic plot
m5.chungcheong <- mosaic_gg(vote.chung.class[, -1],
base_family = "HCR Dotum LVT",
ggtitle = "충청도의 찬반",
xlab = "계급(누적%)",
ylab = "찬반(%)",
fill.name = "찬반")
m5.chungcheong
ggsave("../pics/sejong_chungcheong_mosaic_ggplot.png", m5.chungcheong, width = 12, height = 4.5, dpi = 72)자료 정리
save.image(file = "./sejong_ggplot.RData")