nes$income5.n = as.numeric(nes$income5)
nes$dem_age_group_r.n = as.numeric(nes$dem_age_group_r)

liberalsftgr.w = svyglm(ftgr_liberals ~ pid_3 + gender + income5.n + dem_age_group_r.n, design = nesD, na.action="na.omit")
nesD= svydesign(id=~1, data=nes, weights=~wt)
summary(lm(nes$ftgr_liberals ~ pid_3 + gender + income5.n + dem_age_group_r.n, data = nes))
summary(liberalsftgr.w)
printC(liberalsftgr.w)