[R-sig-ME] Conditional logistic regression vs lmer

Federico Calboli f.calboli at imperial.ac.uk
Tue Feb 23 18:32:05 CET 2010


Hi,

I am doing a conditional logistic regression with clogit() of library survival, but I have a number of random factors I'd like to add to my model and I was thinking of lmer...

clg.tigeR = clogit(c.c ~ tIgE.Resp + strata(match), data.ready)

and all is good. I also have two more variables, both random which I'd like to use, but I am not 100% sure what's the way to go with lmer. I suspect I would need to specify a different intercept for each strata, but seems to baffle me at the moment. The code below might do what I want, but I have no idea if that's correct (i.e. does it match the clogit one + two more random effects?):

lmer.tigeR = lmer(c.c ~ tIgE.Resp + (tIgE.Resp|match) + (1|Run) + (1|Box), family = binomial, data.ready)

lmer.tigeR is computed without problems btw.

Any illumination would be gratefully accepted.

Best

Federico


--
Federico C. F. Calboli
Department of Epidemiology and Biostatistics
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG

Tel +44 (0)20 75941602   Fax +44 (0)20 75943193

f.calboli [.a.t] imperial.ac.uk
f.calboli [.a.t] gmail.com




More information about the R-sig-mixed-models mailing list