I think the .data pronoun doesnt work in that context, so I used a more general tidyeval approach.
However, this meant that I attempt the calculation on all the variables, and some of them required setting success="degree", and one was an ordered factor that I had to throw away the ordering to get any output.
gss_vec <- c("sex", "partyid", "income", "class", "finrela") %>%
set_names()
chi_sq_func <- function(xval){
if(is.ordered(gss[[xval]]))
gss[[xval]] <- factor(gss[[xval]],ordered = FALSE)
eval(expr(specify(x = gss,
formula = college ~ !!sym(xval),
success="degree"))) %>%
hypothesise(null = "independence") %>%
calculate(stat = "Chisq")
}
map_df(.x = gss_vec,
.f = ~chi_sq_func(.x),
.id = "variable")