Re: [R] Analysing ordinal/nominal data

From: Piotr Majdak <piotr_at_majdak.com>
Date: Fri 17 Jun 2005 - 20:54:15 EST

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 use 1,2,3,... as you did. v' and w': mean of v and w, respectively

> 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.)

After your suggestions I've read the logit-Chapter in Agresti (1984) and tried to fit a logit model to my data:

log(m_ij2/m_ij1) = intercept + beta^PR_i(u_i-u') +

                    beta^ENV_j(v_j-v')

with:
mijk: expected frequencies for cell with indicies i,j,k u_i: scores of PR with u={0,1,2,3}
v_j: scores of ENV with v={0,1,2}
beta^PR: association parameter for PR with responses beta^ENV: association parameter for ENV with responses

I wrote in R:

count=c(250,274,285,241,279,299,247,246,255,280,355,362,

         230,207,195,239,200,181,233,235,224,200,125,118)

PR=as.integer(gl(4,3,12))-1
ENV=as.integer(gl(3,1,12))-1
countM=matrix(count,12,2)
fit=glm(countM ~ as.integer(PR)+as.integer(ENV),
	family=binomial(link=logit))

summary(fit)

R gives me:

Coefficients:

                 Estimate Std. Error z value Pr(>|z|)
(Intercept)      0.18435    0.07737   2.383 0.017179 *
as.integer(PR) -0.10538 0.03085 -3.416 0.000635 *** as.integer(ENV) -0.07075 0.04748 -1.490 0.136191

Looking at my data, I know that for PR=0 I have a very weak dependence on ENV and for PR=3 very strong dependence on ENV. How can I get this interpretation from the summary above? I'm new in R: is fit$linearpredictors what I'm looking for?

-- 
Piotr Majdak
Institut für Schallforschung
Österreichische Akademie der Wissenschaften
Reichsratsstr. 17
A-1010 Wien
Tel.: +43-1-4277-29511
Fax: +43-1-4277-9296
E-Mail: piotr@majdak.com
WWW: http://www.kfs.oeaw.ac.at

______________________________________________
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
Received on Fri Jun 17 21:26:16 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:47 EST