A = [1,2,3,4,5,6]; B = [2,3,4,5,6,7,8]; C = [3,4,5,6,7,8,9,10]; Xa = mean(A); Va = var(A); Xb = mean(B); Vb = var(B); Xc = mean(C); Vc = var(C); Na = length(A); Nb = length(B); Nc = length(C); k = 3; N = Na + Nb + Nc G = (Na*Xa + Nb*Xb + Nc*Xc) / N MSSb = (Na*(Xa-G)^2 + Nb*(Xb-G)^2 + Nc*(Xc-G)^2) / (k-1) MSSw = ((Na-1)*Va + (Nb-1)*Vb + (Nc-1)*Vc) / (N-k) F = MSSb / MSSw