install.packages("catspec")
library(catspec)
library(survival)
data <- read.table("C:\\ps206hw7.txt", header=T)
data$vote[data$vote==1] <- 0
data$vote[data$vote==3] <- 1
data$vote[data$vote==0] <- 3
data$votef <- as.factor(data$vote)
revdata <- mclgen(data, votef)
revdata$discand <- NA
revdata$discand[revdata$votef==1] <- revdata$dispero[revdata$votef==1]
revdata$discand[revdata$votef==2] <- revdata$disclin[revdata$votef==2]
revdata$discand[revdata$votef==3] <- revdata$disbush[revdata$votef==3]
revdata$age.1[revdata$votef!=2] <- 0
revdata$age.2[revdata$votef!=3] <- 0
revdata$income.1[revdata$votef!=2] <- 0
revdata$income.2[revdata$votef!=3] <- 0
revdata$gender.1[revdata$votef!=2] <- 0
revdata$gender.2[revdata$votef!=3] <- 0
revdata$educ.1[revdata$votef!=2] <- 0
revdata$educ.2[revdata$votef!=3] <- 0
revdata$race.1[revdata$votef!=2] <- 0
revdata$race.2[revdata$votef!=3] <- 0
mod1 <- clogit(depvar~votef+discand +income.1 +
income.2+age.1+age.2+educ.1+educ.2+gender.1+gender.2+strata(id), data = revdata)
summary(mod1)