0
$\begingroup$

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?

$\endgroup$

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.