#################################################
# TERRORISM
#################################################

terror <- read.csv('/Users/lesliejohns/Dropbox/EIRT/Data/data-terror.csv', header=T, na.strings='.')
terror <- subset(terror, terror$year<2016)
data.us <- subset(terror, terror$us==1)
data.both <- subset(terror, terror$both==1)
data.foreign <- subset(terror, terror$foreign==1)

years <- seq(1997, 2015, length=19)
count.us <- rep(0, length=length(years))
count.both <- rep(0, length=length(years))
count.foreign <- rep(0, length=length(years))
for(i in 1:length(years)) {
      for(j in 1:length(data.us$year)) {
      if(years[i]==data.us$year[j]) {count.us[i] <- count.us[i]+1}
      }}
for(i in 1:length(years)) {
      for(j in 1:length(data.both$year)) {
      if(years[i]==data.both$year[j]) {count.both[i] <- count.both[i]+1}
      }}
for(i in 1:length(years)) {
      for(j in 1:length(data.foreign$year)) {
      if(years[i]==data.foreign$year[j]) {count.foreign[i] <- count.foreign[i]+1}
      }}
counts <- rbind(count.us, count.both, count.foreign)
colors = c("dim grey", "dark grey", "light grey")

pdf("/Users/lesliejohns/Dropbox/EIRT/graphics/fig-terror.pdf")
barplot(counts, space=0, col=colors, xaxt='n', yaxt='n')
axis(side=1, at=c(4,9,14,19), labels=c('2000', '2005', '2010', '2015'), tick=FALSE)
axis(side=2, at=c(0, 3, 6, 9, 12, 15, 18), labels=c(0, 3, 6, 9, 12, 15, 18), lwd=0, lwd.ticks=1, line=-1)
segments(0,0,0,17)
segments(0,17,19,17)
segments(19,17, 19, 0)
mtext(side=1, text="Year", line=3)
mtext(side=2, text="Number of New Cases", line=2)
legend(c(0,16), c("US only", "US & foreign", "Foreign only"), fill=colors, bty='n')
total <- count.us + count.both + count.foreign
lines(lowess(total), lwd=2)
dev.off()

#################################################
# ALIEN TORT STATUTE
#################################################

ATS <- read.csv('/Users/lesliejohns/Dropbox/EIRT/Data/data-ats.csv', header=T, na.strings='.')
ATS <- subset(ATS, ATS$year<2016)

years <- seq(1978, 2015, length=38)
count <- rep(0, length=length(years))
for(i in 1:length(years)) {
      for(j in 1:length(ATS$year)) {
      if(years[i]==ATS$year[j]) {count[i] <- count[i]+1}
      }}

pdf("/Users/lesliejohns/Dropbox/EIRT/graphics/fig-ats.pdf")
barplot(count, space=0, col='light grey', xaxt='n', yaxt='n')
axis(side=1, at=c(3,13,23,33), labels=c('1980', '1990', '2000', '2010'), tick=FALSE)
axis(side=2, at=c(0, 5, 10, 15, 20), labels=c(0, 5, 10, 15, 20), lwd=0, lwd.ticks=1, line=-1)
segments(0,0,0,21)
segments(0,21,38,21)
segments(38,21,38,0)
lines(lowess(count), lwd=2)
dev.off()

#################################################
# FDI ARBITRATION ENFORCEMENT
#################################################

fdi <- read.csv('/Users/Leslie/Dropbox/EIRT/Data/data-fdi.csv', header=T, na.strings='.')
fdi <- subset(fdi, fdi$Year <2016)

years <- seq(1975, 2015, length=41)
count <- rep(0, length=length(years))
for(i in 1:length(years)) {
      for(j in 1:length(fdi$ID)) {
      if(years[i]==fdi$Year[j]) {count[i] <- count[i]+1}
      }}

pdf("/Users/lesliejohns/Dropbox/EIRT/graphics/fig-fdi.pdf")
barplot(count, space=0, col='light grey', xaxt='n', yaxt='n')
axis(side=1, at=c(0,10,20,30,40), labels=c('1975', '1985', '1995', '2005', '2015'), tick=FALSE)
axis(side=2, at=c(0, 1, 2, 3, 4, 5, 6, 7, 8), labels=c(0, 1, 2, 3, 4, 5, 6, 7, 8), lwd=0, lwd.ticks=1, line=-1)
segments(0,0,0,7)
segments(0,7,40,7)
mtext(side=1, text="Year", line=3)
mtext(side=2, text="Number of New Cases", line=2)
lines(lowess(count), lwd=2)
dev.off()