[R] Difference between summary.lm() and summary.aov()

Alexander Sirotkin [at Yahoo] alex_s_42 at yahoo.com
Sun Dec 7 13:21:47 CET 2003


John,

What you are saying is that any conclusion I can make
from summary.aov (for instance, to answer a question
if physician is a significant variable) will not be
correct ?


--- John Fox <jfox at mcmaster.ca> wrote:
> Dear Spencer and Alexander,
> 
> In this case, physician is apparently a factor with
> three levels, so 
> summary.aov() gives you a sequential ANOVA,
> equivalent to what you'd get 
> from anova(). There no simple relationship between
> the F-statistic for 
> physician, which has 2 df in the numerator, and the
> two t's. (By the way, I 
> doubt whether a sequential ANOVA is what's wanted
> here.)
> 
> Regards,
>   John
> 
> At 09:17 AM 12/6/2003 -0800, Spencer Graves wrote:
> >      The square of a Student's t with "df" degrees
> of freedom is an F 
> > distribution with 1 and "df" degrees of freedom.
> >      hope this helps.  spencer graves
> >
> >Alexander Sirotkin [at Yahoo] wrote:
> >
> >>I have a simple linear model (fitted with lm())
> with 2
> >>independant
> >>variables : one categorical and one integer.
> >>
> >>When I run summary.lm() on this model, I get a
> >>standard linear
> >>regression summary (in which one categorical
> variable
> >>has to be
> >>converted into many indicator variables) which
> looks
> >>like :
> >>
> >>            Estimate Std. Error t value Pr(>|t|)
> >>(Intercept)  -3595.3     2767.1  -1.299   0.2005
> >>physicianB     802.0     2289.5   0.350   0.7277
> >>physicianC    4906.8     2419.8   2.028   0.0485 *
> >>severity      7554.4      906.3   8.336 1.12e-10
> ***
> >>
> >>and when I run summary.aov() I get similar ANOVA
> table
> >>:
> >>           Df     Sum Sq    Mean Sq F value   
> Pr(>F)
> >>physician    2  294559803  147279901  3.3557  
> 0.04381
> >>*
> >>severity     1 3049694210 3049694210 69.4864
> 1.124e-10
> >>***
> >>Residuals   45 1975007569   43889057
> >>
> >>What is absolutely unclear to me is how F-value
> and
> >>Pr(>F) for the
> >>categorical "physician" variable of the
> summary.aov()
> >>is calculated
> >>from the t-value of the summary.lm() table.
> >>
> >>I looked at the summary.aov() source code but
> still
> >>could not figure
> >>it.
> >>
> >>Thanks a lot.
> >>
> >>__________________________________
> >>

> >>
> >>______________________________________________
> >>R-help at stat.math.ethz.ch mailing list
>
>>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> >>
> >
> >______________________________________________
> >R-help at stat.math.ethz.ch mailing list
>
>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> 
>
-----------------------------------------------------
> 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
>
-----------------------------------------------------
>




More information about the R-help mailing list