Интересно, как извлечь часть Multivariate Tests: Site из вывода fm1 в следующем MWE. library(car) fm1 <- summary(Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery))) fm1 Type II MANOVA Tests: Sum of squares and products for error:...

Мне нужно вычислить сумму квадратов матрицы перекрестного произведения (на самом деле след этой матрицы) в многомерной линейной модели с Y (n x q) и X (n x p). Стандартный код R для этого: require(MASS) require(car) # Example data q <- 10 n...