anova.r

ANOVA<-function(fit1,fit2){
  temp <- anova(fit2,fit1 )
  fin.aov <- anova(fit1)   
  reg <- temp[2,2:6]
  rownames(reg) <- "Regression"
  reg[1,1:2] <- reg[1,2:3]
  reg[1,3] <- reg[1,2]/reg[1,1]
  colnames(reg) <- colnames(fin.aov)
  res <- fin.aov[tail(nrow(fin.aov),1),]
  tot <- cbind(reg[1,1:2]+res[1,1:2],NA,NA,NA)
  rownames(tot) <- "Total"
  colnames(tot) <- colnames(reg)
  
  rbind(reg,res,tot)
}