Hi,
These are materials for logistic regression course II: Conditional logistic regression, 2010. This year I use R version 2.12.0 and the regression outputs are not different from the 2009 version where R version 2.10.0 was used. However, the different in glm() output at the 3 or later digits after the decimal point is very usual when you change R version. This is a technical issue, not a bug, of the maximum likelihood estimation method. None of the results is wrong. It means that you should not pay attention too much on later digits after the second decimal points.
If you are eager to know why the maximum likelihood estimation method is that tricky, please read some textbooks specific to this matter.
Conditional logistic regression requires library survival that already exists on your R library folder. Just call library(survival) or require(survival) and you are ready to use clogit function.
agechd.dta
cca-match.dta
The module can be downloaded here.
logistic1012-2.pdf
You can also download ICE modules from
r-ice.project.net
The required modules for this course are ice and epid.
Finally, you can follow my R script file below. Don't forget to change working directory to yours.
conditional.R
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment