# # Attempt to Run simple Probit Copying Ernesto's Code # library(arm) setwd("C:/calvo") # # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # Read in data from STATA output # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # rcx.file <- "c:/calvo/h106_data.txt" # # X[,1] = 1 if bush vote >= 50% # X[,2] = 1 if Gore vote >= 50% # X[,3] = Bush Percentage in CD # X[,4] = Gore Percentage in CD # X[,5] = Black Percentage in CD # X[,6] = 1 if Southern State (11 states of Confederacy + OK + KY # X[,7] = Hispanic Percentage in CD # X[,8] = Median Family Income (in thousands) in CD # X[,9] = Percent Owner-Occupied Housing # X[,10] = DW-NOMINATE 1st Dimension # X[,11] = DW-NOMINATE 2nd Dimension # # Standard fields and their widths # rcx.fields <- c("ybush","ygore","bushvote","gorevote","black","south","hispanic","income","owner","dwnom1n","dwnom2n") rcx.fieldWidths <- c(9,11,10,10,11,10,11,11,11,11,11) # # Input File # TT <- read.fwf(file=rcx.file,widths=rcx.fieldWidths,as.is=TRUE,col.names=rcx.fields) dim(TT) nrow <- length(TT[,1]) #ncol <- length(TT[1,]) #data = read.dta("hdmg106_2009_fixed.dta") #attach(data) V <- cbind(TT[,1],TT[,5],TT[,6],TT[,7],TT[,8],TT[,9],TT[,10],TT[,11]) N = nrow K = ncol(V) data.data = list(N=N,K=K,V=V) data.inits = function() {list(beta=runif(K,-2,2), delta=runif(N,-1,1))} data.parameters = c("beta") wide.sim = bugs(data.data, data.inits, data.parameters,"h106_probit_RtoWINBUGS_model.txt", n.chains=4, n.thin=1, n.burnin=15000,n.iter=40000, debug=T) detach(data)