Hi! I'm looking for a solution to analyse data, which consists of dichotomous responses (yes/no) for 2 multinomial ordinal variables. I was trying glm() and got hierarhical models treating all variables as nominal, but I can't figure out how to tell glm() to use a model for ordinal data like this:
log(Mij) = intercept + X + Y + Z + beta*(x-x')*(y-y') where beta is a regression factor for interaction between X and Y. Do you know a trick to code it in R or point me to some documentation? thanks a lot, Piotr Majdak ______________________________________________ R-help@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