Let k in test_kj be a covariate name

[svteichman]

k can either be numeric, a string that matches a column name of X, or a string that starts one or more column names of X (this last case corresponds to a formula like ~ Group where Group has values in c("A", "B", "C"), and then will test all covariates of the form GroupB, GroupC).