I was able to reproduce the reprex, but I'm going crazy trying to apply it to my data. Something is off and I can't figure it out:
colnames(x=logit)[colnames(x=logit)=="W_Male"] <- 'Male Gender'
colnames(x=logit)[colnames(x=logit)=="W_AGE_60"] <- 'Age>60'
colnames(x=logit)[colnames(x=logit)=="W_FACILITY_ACADEMIC"] <- 'Academic Facility'
colnames(x=logit)[colnames(x=logit)=="W_RACE_WHITE"] <- 'White Race'
colnames(x=logit)[colnames(x=logit)=="W_RURAL_250"] <- 'Rural'
colnames(x=logit)[colnames(x=logit)=="W_DIST_60"] <- 'Distance >60 Miles'
colnames(x=logit)[colnames(x=logit)=="W_EOE_Max"] <- 'Extraocular Extension'
colnames(x=logit)[colnames(x=logit)=="Z_year2010"] <- 'Diagnosed after 2010'
colnames(x=logit)[colnames(x=logit)=="W_INSURANCE_UN"] <- 'Uninsured'
colnames(x=logit)[colnames(x=logit)=="W_HS_13"] <- 'Uneducated Zip Code'
colnames(x=logit)[colnames(x=logit)=="W_MED_INC_48"] <- 'Zip code income'
colnames(x=logit)[colnames(x=logit)=="W_CharlsonScore"] <- 'Charlson/Deyo score'
colnames(x=logit)[colnames(x=logit)=="W_SITE_CILIARY"] <- 'Ciliary body involvement'
colnames(x=logit)[colnames(x=logit)=="Z_SZ_class"] <- 'T stage'
forest_model(model = glm(formula = (T_TXT_NGP ~'Male Gender'+'Age>60'+'Academic Facility'+'White Race'+'Rural'+'Distance >60 Miles'+'Extraocular Extension'+'Diagnosed after 2010'+'Uninsured'+'Uneducated Zip Code'+'Zip code income'+'Charlson/Deyo score'+'Ciliary body involvement'+'T stage'), data = logit, family = "binomial"))
Error in terms.formula(formula, data = data) :
invalid model formula in ExtractVars