On 4/2/2012 3:34 AM, Christofer Bogaso wrote:
Dear all, can somebody give me some pointer how I can fit a
"linear-by-linear association model" (i.e. loglinear model for the
ordinal variables) in R? A brief description can be found here
'https://onlinecourses.science.psu.edu/stat504/node/141'.
Thanks for your help
See my short course, http://www.datavis.ca/courses/VCD/
in particular, http://www.datavis.ca/courses/VCD/R/mental-glm.R
# fit linear x linear (uniform) association. Use integer scores for
rows/cols
Cscore <- as.numeric(Mental$ses)
Rscore <- as.numeric(Mental$mental)
linlin <- glm(Freq ~ mental + ses + Rscore:Cscore,
family = poisson, data = Mental)
--
Michael Friendly Email: friendly AT yorku DOT ca
Professor, Psychology Dept.
York University Voice: 416 736-5115 x66249 Fax: 416 736-5814
4700 Keele Street Web: http://www.datavis.ca
Toronto, ONT M3J 1P3 CANADA
______________________________________________
R-help@r-project.org 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.
