# Circular smooths within a GAM-GEE framework

I have a predictor variable which I fit in a GAM as a circular smooth term:

temp <- transform(temp, impact.f = factor(impact))
gam(Dist ~ impact.f + s(TimeAroundHW, by=impact.f, k=k, bs ="cc"),
data=temp, family="Gamma")


I want my standard errors to be robust to residual auto-correlation within blocks of individual as currently there is a high degree of auto-correlation. I also want a population mean response rather than individual effects so I'd like to fit in a GEE framework rather than a mixed-effects model, however I cannot find any examples or possibilities for fitting circular smooths with GEEs. Has anyone done this before?