I am fitting a logistic model using GEE with an interaction term in R:
geeglm(outcome ~ gender*race, family=binomial(link="logit"),data=df, id=cluster,