[R] Difficult to compare models
Ronaldo Reis Jr.
chrysopa at insecta.ufv.br
Mon Apr 21 20:49:45 CEST 2003
Hi,
I have the follow situation in a ancova analysis:
y (count data)
x1 (continuous data)
x2 (qualitative data)
x2 have 10 levels (a,b,c,d,e,f,g,h,i,j)
I make the model
m1 <- glm(y~x1*x2, family=poisson)
then I want to remove the levels h, i and j (make these levels with
intercept = 0 and slope = 0)
A make a new datafile without these levels. This new table is smallest than
the original table.
I refit the model:
m2 <- glm(y~x1*x2, family=poisson)
Then I want to compare m1 with m2, I need to know if the removing of the 3
levels are significant. But howto compare this if the total GL and the total
deviance are different.
Existe another mean to make this? Maybe using offset? How?
Thanks for all
Ronaldo
--
You can observe a lot just by watching. -- Yogi Berra
--
| // | \\ [*****************************][*******************]
|| ( õ õ ) [Ronaldo Reis Júnior ][PentiumIII-600 ]
| V [UFV/DBA-Entomologia ][HD: 30 + 10 Gb ]
|| / \ [36571-000 Viçosa - MG ][RAM: 128 Mb ]
| /(.''`.)\ [Fone: 31-3899-2532 ][Video: SiS620-8Mb ]
||/(: :' :)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar]
|/ (`. `'` ) \[ICQ#: 5692561 ][Kernel: 2.4.18 ]
|| ( `- ) [*****************************][*******************]
||| _/ \_Powered by GNU/Debian W/Sarge D+ || Lxuser#: 205366
More information about the R-help
mailing list