> print(cbind(SS_TOT_U,DF_TOT_U))
             DF_TOT_U
[1,] 1722811       24
> print(cbind(SS_MODEL,DF_MODEL))
             DF_MODEL
[1,] 1720452        3
> print(cbind(SS_RESIDUAL,DF_RESIDUAL))
              DF_RESIDUAL
[1,] 2359.167          21
> print(cbind(SS_MEAN,DF_MEAN))
             DF_MEAN
[1,] 1622920       1
> print(cbind(SS_REGRESSION,DF_REGRESSION))
             DF_REGRESSION
[1,] 97531.8             2
> 
> # F-test for H0: beta1=beta2=0  (Test Statistic, Critical Value (alpha=0.05,P-value)  
> F_OBS <- (SS_REGRESSION/DF_REGRESSION)/(SS_RESIDUAL/DF_RESIDUAL)  
> F_CRIT <- qf(.95,DF_REGRESSION,DF_RESIDUAL)
> F_P <- 1-pf(F_OBS,DF_REGRESSION,DF_RESIDUAL)
> 
> print(cbind(F_OBS,F_CRIT,F_P))
             F_CRIT  
[1,] 434.087 3.4668 0