[R] Any R function for self-controlled case series method /effect
Jari Haukka
jari.haukka at ktl.fi
Tue Mar 20 12:58:17 CET 2007
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
>>
>> ______________________________________________
>> R-help at stat.math.ethz.ch 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.
>>
>>
>
More information about the R-help
mailing list