library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paset("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep="")
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(max(var)-min(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
(out[2]-out[1])/(bounds[i,2]-bounds[i,1])
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep="")
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(max(var)-min(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
(out[2]-out[1])/(bounds[i,2]-bounds[i,1])
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(max(var)-min(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
(out[2]-out[1])/(bounds[i,2]-bounds[i,1])
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(max(var)-min(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(max(var),min(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015")
vars     <-c("internet_users100","pidentity")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","Adams_et_al_AJPS_2006")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015","rep_adams_2006.dta")
vars     <-c("internet_users100","pidentity","vshift")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006")
vars     <-c("internet_users100","pidentity","vshift")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA)),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta")
path <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006")
vars     <-c("internet_users100","pidentity","vshift")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006")
vars     <-c("internet_users100","pidentity","vshift")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)")
vars     <-c("internet_users100","pidentity","vshift","racresent1")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)")
vars     <-c("internet_users100","pidentity","vshift","racresent1")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta",
"rep_bodea_io_2015a.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)",
"Bodea_IO_2015")
vars     <-c("internet_users100","pidentity","vshift","racresent1",
"polity2_cen")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Data/Included/Malesky_APSR_2012")
library(foreign)
library(coda)
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta",
"rep_bodea_io_2015a.dta","rep_truex_2014a.dta","rep_truex_2014a.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)",
"Bodea_IO_2015","Truex_APSR_2014","Truex_APSR_2014")
vars     <-c("internet_users100","pidentity","vshift","racresent1",
"polity2_cen","so_portion","rev2007")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:length(datasets)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",path[i],sep=""))
d <- read.dta(paste(datasets[i]))
var <- d[,paste(vars[i])]
if(is.na(bounds[i,1])){
bounds[i,1:2] <- c(min(var),max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Code")
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv("moderatorlist.csv")
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv("moderatorlist.csv",header=F)
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F)
dd
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F,sep=",")
dd
dd[1]
dd$V1[1]
paste(dd$V1[1])
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F,sep=",")
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta",
"rep_bodea_io_2015a.dta","rep_truex_2014a.dta","rep_truex_2014a.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)",
"Bodea_IO_2015","Truex_APSR_2014","Truex_APSR_2014")
vars     <-c("internet_users100","pidentity","vshift","racresent1",
"polity2_cen","so_portion","rev2007")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999)){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
dd
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var))
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var)
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(bounds[i,2]-bounds[i,1])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F,sep=",")
datasets <-c("rep_malesky_2012.dta","rep_huddy_2015a.dta","rep_adams_2006.dta","rep_banks_2012.dta",
"rep_bodea_io_2015a.dta","rep_truex_2014a.dta","rep_truex_2014a.dta")
path     <-c("Malesky_APSR_2012","Huddy_APSR_2015","Adams_et_al_AJPS_2006","Banks_Valentino_AJPS_2012/dataverse_files (6)",
"Bodea_IO_2015","Truex_APSR_2014","Truex_APSR_2014")
vars     <-c("internet_users100","pidentity","vshift","racresent1",
"polity2_cen","so_portion","rev2007")
bounds   <- matrix(c(
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA),
c(NA,NA)
),
ncol=2)
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var)
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(dd$V5[i]-dd$V4[i])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F,sep=",")
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var)
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(dd$V5[i]-dd$V4[i])
cat(" \n ",datasets[i],": ",stat,"\n")
}
setwd("~/Dropbox/Interaction Paper/Code")
library(foreign)
library(coda)
dd <- read.csv2("moderatorlist.csv",header=F,sep=",")
for(i in 1:nrow(dd)){
setwd(paste("~/Dropbox/Interaction Paper/Data/Included/",dd$V1[i],sep=""))
d <- read.dta(paste(dd$V2[i]))
var <- d[,paste(dd$V3[i])]
if(dd$V4[i]==999){
dd$V4[i] <- min(var)
dd$V5[i] <- max(var)
}
out <- HPDinterval(as.mcmc(var), prob=0.8)
stat <- (out[2]-out[1])/(dd$V5[i]-dd$V4[i])
cat(" \n ",datasets[i],": ",stat,"\n")
}
dim(d)
dim(d)
data = d
nbins = 3
cuts.X
nbins
groupX
formula0
formula1
Z
names(data.aug)
formula0
FE
Z
formula0
formula1
v[1:6,1:6]
v
