[R] dropterm() in MASS
John Fox
jfox at mcmaster.ca
Mon Apr 29 06:36:00 CEST 2002
Dear Patrick,
At 11:51 AM 4/29/2002 +1200, Patrick Connolly wrote:
>To compare two different models, I've compared the result of using
>dropterm() on both.
>
>
>Single term deletions
>
>Model:
>growth ~ days + I(days^0.5)
> Df Sum of Sq RSS AIC
><none> 2.8750 -0.2290
>days 1 4.8594 7.7344 4.6984
>I(days^0.5) 1 0.0234 2.8984 -2.1722
>
>AND
>
>Single term deletions
>
>Model:
>growth ~ days + I(days^2)
> Df Sum of Sq RSS AIC
><none> 2.8750 -0.2290
>days 1 5.8338 8.7088 5.5290
>I(days^2) 1 0.0234 2.8984 -2.1722
>
>I understood that the AIC column was the change in AIC when that term
>was dropped from the model.
>
>QUESTION:
>
>Why would (should) the effect of dropping quite different terms be
>exactly the same? It's consistent with all the data of this type I've
>tried it on.
Notice that the residual sums of squares for the two models are identical,
so the peculiar result here pertains not just to dropping a term but also
to the general fit of the model. What does the variable days look like? Are
there just three different values of this variable?
I hope that this helps,
John
-----------------------------------------------------
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario, Canada L8S 4M4
email: jfox at mcmaster.ca
phone: 905-525-9140x23604
web: www.socsci.mcmaster.ca/jfox
-----------------------------------------------------
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list