> P [1] 16841.32 > SSE [1] 1444.654 > SSA [1] 48589.07 > SSB [1] 1141.46 > > SSNA <- P^2/((SSA/B_levels)*(SSB/A_levels)) > SSNA [1] 250.5816 > > SSRES <- SSE-SSNA > SSRES [1] 1194.073 > > MSA <- SSA/(A_levels-1) > MSB <- SSB/(B_levels-1) > MSE <- SSE/((A_levels-1)*(B_levels-1)) > MSNA <- SSNA/1 > MSRES <- SSRES/((A_levels-1)*(B_levels-1)-1) > > F.NA <- MSNA/MSRES > P_F.NA <- 1-pf(F.NA,1,(A_levels-1)*(B_levels-1)-1) > > F.NA [1] 7.34491 > P_F.NA [1] 0.01034526 > >