> Bear <- lm(AGE~WEIGHT) # lm is for "linear model" and does linear regressions.  Here, y=AGE and x=WEIGHT.
> summary(Bear)

Call:
lm(formula = AGE ~ WEIGHT)

Residuals:
    Min      1Q  Median      3Q     Max 
-29.245 -11.919  -4.586   1.613  80.995 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  5.59389    5.57349   1.004     0.32    
WEIGHT       0.20736    0.02544   8.152 7.33e-11 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 22.56 on 52 degrees of freedom
Multiple R-Squared: 0.561,      Adjusted R-squared: 0.5526 
F-statistic: 66.46 on 1 and 52 DF,  p-value: 7.325e-11 




> anova(Bear) # can pick out estimate of variance (but sd picked out of summary)
Analysis of Variance Table

Response: AGE
          Df Sum Sq Mean Sq F value    Pr(>F)    
WEIGHT     1  33810   33810  66.457 7.325e-11 ***
Residuals 52  26455     509                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1