# # mimic_keyvotes.r # # ANALYZING SPATIAL MODELS OF CHOICE AND JUDGMENT WITH R # Dave Armstrong, Ryan Bakker, Royce Carroll, Christopher Hare, Keith T. Poole, and Howard Rosenthal # http://voteview.com/asmcjr.asp # # # Remove all objects just to be safe. rm(list=ls(all=TRUE)) # library(foreign) library(rjags) library(coda) library(wnominate) # H111 <- readKH("ftp://voteview.com/dtaord/hou111kh.ord") # keyvotes <- c(49,69,103,334,351,475,678,680,717,882,885,1150,1301,1302,1320,1376,1398,1415,1457,1496,1501,1537,1541,1544,1608,1630) # H111$votes <- H111$votes[,keyvotes] H111$m <- length(keyvotes) # result <- wnominate(H111, minvotes=15, dims=1, polarity=2) wnom1 <- result$legislators$coord1D wnom1[is.na(wnom1)] <- 0 # rollcalls <- H111$votes rollcalls[rollcalls==0 | rollcalls==7 | rollcalls==8 | rollcalls==9] <- NA rollcalls[rollcalls==1] <- 1 rollcalls[rollcalls==6] <- 0 # rownames(rollcalls) <- colnames(rollcalls) <- NULL rollcalls <- as.matrix(rollcalls) # N <- nrow(rollcalls) Q <- ncol(rollcalls) y <- rollcalls # inits <- function() {list (ideo=c(-1,1,wnom1[3:N]), alpha=rnorm(Q,0,1), beta=rnorm(Q,0,1))} # MIMIC.sim <- jags.model( url("ftp://voteview.com/wf1/mimic_keyvotes.bug"), data = list('y' = y, 'Q' = Q, 'N' = N), inits = inits, n.chains = 1, n.adapt = 10) # ideo <- coda.samples(MIMIC.sim, 'ideo', 10, thin=1) # Normalize "ideo" across iterations ideo2 <- t(apply(ideo[[1]], 1, scale)) #