printlog
attach(printlog)
plot(week1,week2)

print.lm=lm(week2~week1)
summary(print.lm)
anova(print.lm)


> attach(printlog)
> print.lm=lm(week2~week1)
> summary(print.lm)

Call:
lm(formula = week2 ~ week1)

Residuals:
    Min      1Q  Median      3Q     Max 
-371.55  -64.82  -35.77   42.99  490.27 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 60.92735    9.66554   6.304 1.24e-09 ***
week1        0.44481    0.06538   6.803 7.02e-11 ***
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 116.4 on 260 degrees of freedom
Multiple R-Squared: 0.1511,     Adjusted R-squared: 0.1478 
F-statistic: 46.28 on 1 and 260 DF,  p-value: 7.025e-11 

> anova(print.lm)
Analysis of Variance Table

Response: week2
           Df  Sum Sq Mean Sq F value    Pr(>F)    
week1       1  626519  626519   46.28 7.025e-11 ***
Residuals 260 3519758   13538                      
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 


plot(week1,week2)
lines(week1,fitted(print.lm))

plot(fitted(print.lm),resid(print.lm))
plot(week1,resid(print.lm))

qqnorm(resid(print.lm))