R, Statistics, Data Science.
freedom for the error. s <- min(vh, p) m <- .5 * (abs(vh - p) - 1) N <- .5 * (ve - p - 1) f.approx <- ((2 * N + s + 1) * vs) / ((2 * m + s + 1) * (s - vs)) f.approx ## [1] 4.069718 root.manova$stats[,3][1] ## root$Tree.Number ## 4.069718 The calculated approximate F-statistic matches what was