Hi,
These are materials for logistic regression course I, 2011. This year I use R version 2.14.0 and the regression outputs are slightly different from those produced by the 2010 version where R version 2.12.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.
anc.dta
agechd.dta
cca.dta
lowbwt.dta
The revised module can be downloaded here.
logistic1112-1.pdf
You can also download ICE modules from
www.r-ice-project.net
The required modules for this course are ice and epid which can be installed into R by typing the following two lines on your R console.
install.packages("ice",repos="http://r-ice-project.net")
install.packages("epid",repos="http://r-ice-project.net")
Finally, you can try following my R script file. Please change the working directory to yours, otherwise the script will not run correctly.
exercises.R
Thursday, December 15, 2011
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment