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 ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
