linear test for trend from a gee model

I am trying to run a trend test using a binary outcome and an ordinal exposure (+ other covariates in R)

I haven't been able to find much online about how to do this. I have a colleague who did it in SAS but I've tried to redo it in R and no luck.