.Rhistory

?devtools::install_github
library(QuickUDS)
names(extended_uds)
names(extended_uds[ ,c("arat_pmm":"wahman_teorell_hadenius")])
names(extended_uds[ ,c(12:39)])
names(democracy)
names(democracy[ , c(11:71)])
?opts_chunk
library(QuickUDS)
source('~/.active-rstudio-document', echo=TRUE)
source('~/.active-rstudio-document', echo=TRUE)
data <- extended_uds %>% mutate(adj.z1 = z1 -0.68,
adj.pct025 = pct025 - 0.68,
adj.pct975 = pct975 - 0.68,
index = pnorm(adj.z1),
index.pct025 = pnorm(adj.pct025),
index.pct975 = pnorm(adj.pct975))
countries <- c("United States of America")
data <- data %>% arrange(country_name,year) %>% group_by(country_name) %>% mutate(group = PoliticalDatasets::count_sequence_breaks(year)) %>% ungroup()
ggplot(data %>% filter(country_name %in% countries), aes(x=year,y=index,ymin = index.pct025,ymax = index.pct975, group= group)) +
geom_path() +
geom_ribbon(alpha=0.2) +
theme_bw() +
labs(x = "Year", y = "Extended unified democracy score,\ntransformed to 0-1 probability scale, per year",color="")  +
theme(legend.position="bottom") +
guides(color = guide_legend(ncol=1),fill = guide_legend(nrow=1)) +
facet_wrap(~country_name,ncol=2) +
geom_hline(yintercept = 0.5, color="red") +
geom_vline(xintercept = c(1928,1871,1920,1964), linetype=2)
data <- extended_uds %>% mutate(adj.z1 = z1 -0.68,
adj.pct025 = pct025 - 0.68,
adj.pct975 = pct975 - 0.68,
index = pnorm(adj.z1),
index.pct025 = pnorm(adj.pct025),
index.pct975 = pnorm(adj.pct975))
countries <- c("United States of America")
data <- data %>% arrange(country_name,year) %>% group_by(country_name) %>% mutate(group = PoliticalDatasets::count_sequence_breaks(year)) %>% ungroup()
ggplot(data %>% filter(country_name %in% countries), aes(x=year,y=index,ymin = index.pct025,ymax = index.pct975, group= group)) +
geom_path() +
geom_ribbon(alpha=0.2) +
theme_bw() +
labs(x = "Year", y = "Extended unified democracy score,\ntransformed to 0-1 probability scale, per year",color="")  +
theme(legend.position="bottom") +
guides(color = guide_legend(ncol=1),fill = guide_legend(nrow=1)) +
facet_wrap(~country_name,ncol=2) +
geom_hline(yintercept = 0.5, color="red") +
geom_vline(xintercept = c(1828,1871,1920,1964), linetype=2)
library(ggrepel)
?geom_text_repel
labels <- data.frame(year = c(1828,1871,1920,1954), labels = c("Age of Jackson","Post Civil War era", "Post Women's suffrage era", "Civil rights movement"))
data <- left_join(data,labels)
source('~/.active-rstudio-document', echo=TRUE)
source('~/.active-rstudio-document', echo=TRUE)
pnorm(1)
pnorm(0.68)
source('~/.active-rstudio-document', echo=TRUE)
?row_number
countries <- unique(data$country_name)
periods_of_independence <- extended_uds %>% group_by(country_name,GW_startdate,GW_enddate) %>% summarise(start = lubridate::year(GW_startdate),
end = lubridate::year(GW_enddate))
names(extended_uds)
periods_of_independence <- extended_uds %>% group_by(country_name,group3) %>% summarise(start = min(year), end = max(year))
periods_of_independence
periods_of_independence <- extended_uds %>% group_by(country_name,group3,in_system) %>% summarise(start = min(year), end = max(year))
periods_of_independence
periods_of_independence <- extended_uds %>% group_by(country_name,group3,in_system) %>% summarise(start = min(year), end = max(year)) %>% filter(in_system)
?geom_rect
?democracy
periods_of_independence <- data %>% group_by(country_name,group3,in_system) %>% summarise(start = min(year), end = max(year)) %>% filter(in_system)
source('~/.active-rstudio-document', echo=TRUE)
ggplot() +
geom_path(data %>% filter(country_name %in% countries[113:224]), aes(x=year,y=index,ymin = index.pct025,ymax = index.pct975, group= group)) +
geom_ribbon(data %>% filter(country_name %in% countries[113:224]), aes(x=year,y=index,ymin = index.pct025,ymax = index.pct975, group= group), alpha=0.2)
periods_of_independence <- data %>% group_by(country_name,group3,in_system) %>% summarise(start = min(year), end = max(year)) %>% filter(in_system)
data <- left_join(data,PoliticalDatasets::population.data)
head(data)
data <- data %>% mutate(prop = pop/world.population.spline.interpolated)
source('~/.active-rstudio-document', echo=TRUE)
data <- extended_uds %>% mutate(adj.z1 = z1 -0.68,
adj.pct025 = pct025 - 0.68,
adj.pct975 = pct975 - 0.68,
index = pnorm(adj.z1),
index.pct025 = pnorm(adj.pct025),
index.pct975 = pnorm(adj.pct975))
data <- left_join(data,PoliticalDatasets::population.data)
data <- data %>% mutate(prop = pop/world.population.spline.interpolated)
data <- data %>% mutate(regime = cut(index, 5, include.lowest=TRUE))
ggplot(data = data %>% filter(year > 1815),  aes(x=year, fill=regime,weight=prop)) +
geom_bar() +
theme_bw()+
theme(legend.position="bottom") +
labs(fill="", x="Year", y= "Proportion of the world's population")+
guides(fill=guide_legend(ncol=2,title.position="top"))+
scale_y_continuous(label=scales::percent)
data <- extended_uds %>% mutate(adj.z1 = z1 -0.68,
adj.pct025 = pct025 - 0.68,
adj.pct975 = pct975 - 0.68,
index = pnorm(adj.z1),
index.pct025 = pnorm(adj.pct025),
index.pct975 = pnorm(adj.pct975))
data <- left_join(data,PoliticalDatasets::population.data)
data <- data %>% mutate(prop = pop/world.population.spline.interpolated)
data <- data %>% mutate(regime = cut(index, 5, include.lowest=TRUE))
ggplot(data = data %>% filter(year > 1815),  aes(x=year, fill=regime,weight=prop)) +
geom_bar() +
theme_bw()+
theme(legend.position="bottom") +
labs(fill="", x="Year", y= "Proportion of the world's population")+
guides(fill=guide_legend(ncol=2,title.position="top"))+
scale_y_continuous(label=scales::percent) +
scale_fill_brewer(type="div",palette="RdBu")
data <- extended_uds %>% mutate(adj.z1 = z1 -0.68,
adj.pct025 = pct025 - 0.68,
adj.pct975 = pct975 - 0.68,
index = pnorm(adj.z1),
index.pct025 = pnorm(adj.pct025),
index.pct975 = pnorm(adj.pct975))
data <- left_join(data,PoliticalDatasets::population.data)
data <- data %>% mutate(prop = pop/world.population.spline.interpolated)
data <- data %>% mutate(regime = cut(index, c(0,0.2,0.4,0.6,0.8,1), include.lowest=TRUE))
ggplot(data = data %>% filter(year > 1815),  aes(x=year, fill=regime,weight=prop)) +
geom_bar() +
theme_bw()+
theme(legend.position="bottom") +
labs(fill="", x="Year", y= "Proportion of the world's population")+
guides(fill=guide_legend(ncol=2,title.position="top"))+
scale_y_continuous(label=scales::percent) +
scale_fill_brewer(type="div",palette="RdBu")
data <- left_join(interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, adj.z1.x = adj.z1) %>% select(country_name.x,GWn.x,year,adj.z1.x))
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, adj.z1.x = adj.z1) %>% select(country_name.x,GWn.x,year,adj.z1.x))
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
?rnorm
data <- data %>% group_by(country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,z1.x, se.z1.x),z1.y.norm = rnorm(1000,z1.y, se.z1.y)))
data <- data %>% group_by(country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(1000,.$z1.y, .$se.z1.y)))
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
warnings()
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm), index.y.norm = pnorm(z1.y.norm))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(warname,country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm), index.y.norm = pnorm(z1.y.norm))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
geom_point(data = data, aes(x=index.x.norm,y=index.y.norm), alpha=0.1) +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
library(ggalt)
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
library(ggalt)
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
?geom_circle
?geom_point
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm- 0.68), index.y.norm = pnorm(z1.y.norm - 0.68))
library(ggalt)
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red")
gplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", fill = "Probability density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
coord_cartesian(ylim = c(0,1))
library(ggalt)
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", fill = "Probability density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
coord_cartesian(ylim = c(0,1))
demowars <- data %>% filter(index.x > 0.5, index.y > 0.5)
demowars <- data %>% filter(index.x > 0.5, index.y > 0.5) %>% distinct(warname,index.x,index.y)
demowars
ggplot() +
geom_path(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y, group=warname)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
coord_cartesian(ylim = c(0,1))
ggplot() +
geom_path(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y, group=paste(country.x,country.y))) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
coord_cartesian(ylim = c(0,1))
ggplot() +
geom_path(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y, group=paste(country_name.x,country_name.y))) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
coord_cartesian(ylim = c(0,1))
head(demowars)
demowars <- demowars %>% mutate(label = paste0(warname,": ",country_name.x," vs. ", country_name.y,", ",year))
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(warname,country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm- 0.68), index.y.norm = pnorm(z1.y.norm - 0.68))
library(ggalt)
library(ggrepel)
demowars <- data %>% filter(index.x > 0.5, index.y > 0.5) %>% distinct(warname,index.x,index.y)
demowars <- demowars %>% mutate(label = paste0(warname,": ",country_name.x," vs. ", country_name.y,", ",year))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label))
coord_cartesian(ylim = c(0,1))
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(warname,country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm- 0.68), index.y.norm = pnorm(z1.y.norm - 0.68))
library(ggalt)
library(ggrepel)
demowars <- data %>% filter(index.x > 0.6, index.y > 0.6) %>% distinct(warname,index.x,index.y)
demowars <- demowars %>% mutate(label = paste0(warname,": ",country_name.x," vs. ", country_name.y,", ",year))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label))
coord_cartesian(ylim = c(0,1))
?geom_text_repel
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(warname,country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm- 0.68), index.y.norm = pnorm(z1.y.norm - 0.68))
library(ggalt)
library(ggrepel)
demowars <- data %>% filter(index.x > 0.6, index.y > 0.6) %>% distinct(warname,index.x,index.y)
demowars <- demowars %>% mutate(label = paste0(warname,": ",country_name.x," vs. ", country_name.y,", ",year))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y), aes(x=index.x,y=index.y)) +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label), force = 2)+
coord_cartesian(ylim = c(0,1))
?geom_corss
?geom_crossbar
data <- left_join(PoliticalDatasets::interstate_dyads_yearly, QuickUDS::extended_uds %>% rename(country_name.x = country_name, GWn.x = GWn, z1.x = z1, se.z1.x = se.z1) %>% select(country_name.x,GWn.x,year,z1.x,se.z1.x))
data <- left_join(data, QuickUDS::extended_uds %>% rename(country_name.y = country_name, GWn.y = GWn, z1.y = z1, se.z1.y = se.z1) %>% select(country_name.y,GWn.y,year,z1.y,se.z1.y))
data <- data %>% group_by(warname,country_name.x,GWn.x,year,z1.x,se.z1.x,country_name.y,GWn.y,year,z1.y,se.z1.y) %>% do(data.frame(z1.x.norm = rnorm(500,.$z1.x, .$se.z1.x),z1.y.norm = rnorm(500,.$z1.y, .$se.z1.y)))
data <- data %>% ungroup() %>% mutate(index.x = pnorm(z1.x - 0.68), index.y = pnorm(z1.y - 0.68), index.x.norm = pnorm(z1.x.norm- 0.68), index.y.norm = pnorm(z1.y.norm - 0.68), xmin = pnorm(z1.x - 0.68 - 1.96*se.z1.x), xmax = pnorm(z1.x - 0.68 + 1.96*se.z1.x), ymin = pnorm(z1.y - 0.68 - 1.96*se.z1.y), ymax = pnorm(z1.y - 0.68 + 1.96*se.z1.y))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y),  aes(x=index.x, y = index.y))
geom_rect(data = data %>% distinct(warname,index.x,index.y), aes(xmin=xmin, xmax = xmax, ymin = ymin, ymax = ymax), alpha=0.2, color="blue") +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label), force = 2)+
coord_cartesian(ylim = c(0,1))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y),  aes(x=index.x, y = index.y)) +
geom_rect(data = data %>% distinct(warname,index.x,index.y), aes(xmin=xmin, xmax = xmax, ymin = ymin, ymax = ymax), alpha=0.2, color="blue") +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label), force = 2)+
coord_cartesian(ylim = c(0,1))
ggplot() +
geom_point(data = data %>% distinct(warname,index.x,index.y),  aes(x=index.x, y = index.y)) +
geom_rect(data = data %>% distinct(warname,index.x,index.y), aes(xmin=xmin, xmax = xmax, ymin = ymin, ymax = ymax), alpha=0.2, fill="blue") +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label), force = 2)+
coord_cartesian(ylim = c(0,1))
ggplot() +
geom_rect(data = data %>% distinct(warname,index.x,index.y), aes(xmin=xmin, xmax = xmax, ymin = ymin, ymax = ymax), alpha=0.1, fill="blue") +
stat_bkde2d(data = data, bandwidth=c(0.05,0.05),aes(x=index.x.norm, y=index.y.norm, alpha=..level..), geom="polygon") +
geom_point(data = data %>% distinct(warname,index.x,index.y),  aes(x=index.x, y = index.y)) +
theme_bw() +
theme(legend.position="bottom") +
labs(x="UD score of side A in war\n(0-1 probability scale)",y = "UD score of side B in war\n(0-1 probability scale)", alpha = "Density") +
geom_vline(xintercept=0.5,color="red") +
geom_hline(yintercept=0.5,color="red") +
geom_text_repel(data = demowars, aes(x=index.x,y=index.y, label = label), force = 2)+
coord_cartesian(ylim = c(0,1))
extended_uds %>% mutate(index = pnorm(z1 - 0.68), percentile = percent_rank(index))
extended_uds %>% mutate(index = pnorm(z1 - 0.68), percentile = percent_rank(index)) %>% group_by(percentile) %>% summarise(index)
extended_uds %>% mutate(index = pnorm(z1 - 0.68), percentile = percent_rank(index)) %>% group_by(percentile) %>% summarise(mean(index))
extended_uds %>% mutate(index = pnorm(z1 - 0.68), percentile = percent_rank(index)) %>% group_by(percentile) %>% summarise(m = mean(index)) %>% filter(m > 0.8)
extended_uds %>% mutate(index = pnorm(z1 - 0.68), percentile = percent_rank(index)) %>% group_by(percentile) %>% summarise(m = mean(index)) %>% filter(m > 0.75)
?QuickUDS
library(QuickUDS)
?QuickUDS
??QuickUDS
?democracy