Re: [R] Analysing ordinal/nominal data

2005-06-17 Thread Prof Brian Ripley
On Thu, 16 Jun 2005, Piotr Majdak wrote: I'm looking for a solution to analyse data, which consists of dichotomous responses (yes/no) for 2 multinomial ordinal variables. Please explain how you get a binary response for a `multinomial ordinal variables'? If you intend these variables to be

Re: [R] Analysing ordinal/nominal data

2005-06-17 Thread Piotr Majdak
Prof Brian Ripley wrote: On Thu, 16 Jun 2005, Piotr Majdak wrote: I'm looking for a solution to analyse data, which consists of dichotomous responses (yes/no) for 2 multinomial ordinal variables. Please explain how you get a binary response for a `multinomial ordinal variables'? If

Re: [R] Analysing ordinal/nominal data

2005-06-17 Thread Prof Brian Ripley
As I suggested before, a binomial logistic model is appropriate here, not a Poisson log-linear one. (They are equivalent, but the binomial version is easier to interpret and less wasteful to fit.) You have still not defined v' and w', nor the scores (are they estimated or not). But the model

Re: [R] Analysing ordinal/nominal data

2005-06-17 Thread Piotr Majdak
Prof Brian Ripley wrote: You have still not defined v' and w', nor the scores (are they estimated or not). But the model I suggested is such a model with scores 1,2,... Sorry for that, here it is: scores v and w: integer scores, reflecting the ordering of columns/rows. Agresti suggests to

Re: [R] Analysing ordinal/nominal data

2005-06-17 Thread Piotr Majdak
Hi Brian (and the list of course!), I still have problems analysing data in R, because I don't know how to tell glm() use row-effect model, please. The models are well defined by Agresti, but can't get the link from the theory to the implementations in R. Different names, definitions and no

[R] Analysing ordinal/nominal data

2005-06-16 Thread 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