[R] Analysing ordinal/nominal data

[Piotr Majdak]

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