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?


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.