From j.c.l.ooms at uu.nl  Tue Sep 30 15:28:14 2008
From: j.c.l.ooms at uu.nl (Jeroen Ooms)
Date: Tue, 30 Sep 2008 06:28:14 -0700 (PDT)
Subject: [Rd]  weird behavior of drop1() for polr models (MASS)
Message-ID: <19742120.post@talk.nabble.com>


I would like to do a SS type III analysis on a proportional odds logistic
regression model. I use drop1(), but dropterm() shows the same behaviour. It
works as expected for regular main effects models, however when the model
includes an interaction effect it seems to have problems with matching the
parameters to the predictor terms. An example:

library("MASS");
options(contrasts = c("contr.treatment", "contr.poly"));

house.plr1 <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data =
housing);
drop1(house.plr1,attributes(house.plr1$terms)$term.labels,test="Chisq");

house.plr2 <- polr(Sat ~ Infl * Type + Cont, weights = Freq, data =
housing); 
drop1(house.plr2,attributes(house.plr2$terms)$term.labels,test="Chisq");

Notice that model 2 has a * instead of a + between predictors Infl and Type.
In model 1, estimated parameters are nicely attributed to the right term,
however in house.plr2, only 2 of the 4 terms are evaluated.

I am using R version 2.7.2 (2008-08-25) i386-pc-mingw32, and MASS_7.2-44.

-- 
View this message in context: http://www.nabble.com/weird-behavior-of-drop1%28%29-for-polr-models-%28MASS%29-tp19742120p19742120.html
Sent from the R devel mailing list archive at Nabble.com.



[Rd] weird behavior of drop1() for polr models (MASS)
Jeroen Ooms