NCI voting

library(haven)
library(tidyverse)
library(dplyr)

left_join(NCI_17_18, state_disprov, by = "STATENUM")

glimpse(NCI_17_18)

glimpse(state_disprov)

NCI_dis <- left_join(NCI_17_18, state_disprov, by = "STATENUM")
glimpse(NCI_dis)

NCI.select <- NCI_dis %>%
  select(.,
         STATENUM,
         STATENAM,
         DIS_PROV,
         VOTE15,
         AGE,
         RACE15_AMERIN,
         RACE15_ASIAN,
         RACE15_BLACK,
         RACE15_DK,
         RACE15_HAWAII,
         RACE15_WHITE,
         RACE15_LATINO,
         RACE15_OTHER,
         LGLSTAT15,
         LEVELMR08,
         EXPRESS15,
         MOBILITY08,
         GENDER)
glimpse(NCI.select)

age.cat <- rep(NA,25671)
age.cat[NCI.select$AGE <= 29] <- 1
age.cat[NCI.select$AGE >= 30 & NCI.select$AGE <=44] <- 2
age.cat[NCI.select$AGE >=45 & NCI.select$AGE <=59] <- 3
age.cat[NCI.select$AGE >=60] <- 4
NCI.select$age.cat <- as.factor(age.cat)

glimpse(NCI.select)

table(NCI.select$VOTE15)

```{r}
NCI.select %>%
mutate(vote.fac = as_factor(VOTE15))

guardian_dic <- recode_factor(NCI.select$LGLSTAT15, "1"="0", "2"="1", "3"="1", "4"="1")
guardian_dic

NCI.select$guardian_dic <- as.factor(guardian_dic)

glimpse(NCI.select)

guardian_dic <- recode_factor(NCI.select$LGLSTAT15, "1"="0", "2"="1", "3"="1", "4"="1")
guardian_dic

```{r}
comm_dic <- recode_factor(NCI.select$EXPRESS15, "1"="0", "2"="1", "3"="1", "4"="1", "5"="1")
comm_dic

NCI.select$comm_dic <- as.factor(comm_dic)

mobility_dic <- recode_factor(NCI.select$MOBILITY08, "1"="0", "2"="1", "3"="1")
mobility_dic

NCI.select$mobility_dic <- as.factor(mobility_dic)

table(NCI.select$mobility_dic)

library(lme4)

library(lmerTest)
library(Hmisc)

describe(NCI.select)

nci.clean <- NCI.select %>%
  
  filter(!(LEVELMR08 %in%
             c("5", "6", "99")))%>%
  mutate(., 
         vote.fac = as_factor(VOTE15),
         )

model.null <- glmer(vote.fac ~ (1|STATENUM), family = binomial, data = nci.clean)
summary(model.null)

model.2 <- glmer(vote.fac ~ age.cat + LEVELMR08 + comm_dic + mobility_dic + GENDER + RACE15_AMERIN + RACE15_ASIAN + RACE15_BLACK + RACE15_HAWAII + RACE15_WHITE + RACE15_LATINO + RACE15_OTHER + guardian_dic + (1|STATENUM), family = binomial, data = nci.clean)
summary(model.2)