[R] Any R function for self-controlled case series method /effect absorption?

[Jari Haukka]

Hello,

Has anyone written R functions for applying self-controlled case series 
methods (http://statistics.open.ac.uk/sccs/).

In fact only thing needed is to modify glm function to allow absorption 
of effect. Eg. in Poisson model individual effect is used as factor, but 
it is considered as nuisance term where parameter estimates are not needed.

Could anyone point how absorbing individual effect could be carried out 
in glm.

There is already code for  Stata 
(http://statistics.open.ac.uk/sccs/stata/aglm.ado), Genstat 
(http://statistics.open.ac.uk/sccs/genstat%5Csccs.gen), Glim 
(http://statistics.open.ac.uk/sccs/glim/SCCS.MAC) , and SAS 
(http://statistics.open.ac.uk/sccs/sas/macro/poisreg.sas).


Regards,
Jari Haukka