Dear Vito, I found found package "gnm" with function "gnm". The "eliminate" argument does exactly what I needed. Conditional Poisson models are handy to fit by gnm.
Regards, Jari vito muggeo wrote: > Dear Jari > The problem is to build the dataset to apply the conditional logit > model. However, as far as I know, no R function exists. > > BTW if you are dealing with time series of pollution and health, the > following two papers might be of interest of you: > It appears that the time series approach could be preferred. > > > Heather J. Whitaker, Mounia N. Hocine, C. Paddy Farrington > On case-crossover methods for environmental time series data > Environmetrics > Volume 18, Issue 2, Date: March 2007, Pages: 157-171 > > Lu, Zeger. On the equivalence of case-crossover and time series > methods in environmental epidemiology Biostatistics, Early view > > > Jari Haukka wrote: >> 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. >> >> > ______________________________________________ [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.
