R, Statistics, Data Science.
sum(b) E[j,i] <- sum(b) } } Verify the results by comparing the computed matrices to the output of summary(manova()). root.manova$SS[1] ## $`root$Tree.Number` ## [,1] [,2] [,3] [,4] ## [1,] 0.07356042 0.5373852 0.3322646 0.208470 ## [2,] 0.53738521 4.1996619 2.3553885 1.637108 ## [3,] 0.33226458