如何计算多个因素和水平的平均值？

Question

我有一个包含多列的数据集（Orig.data）：第一列=主题，随后的列是因子（Factor1，Factor2，...），每个因子具有2个级别（A，B分别为Factor1、1例如，＆2为Factor2）；然后将结果放在不同的列（因变量）中。

每种情况均包括两项试验。因此，前两行代表主题1的试验1＆2，因子1-A，因子2-1，第3和第4行是Sub1因子1-A的试验，因子2-2，第5和第6行是因子1-B，因子2-1的Sub1。 ， 等等。

我想计算前两行，后两行的平均值，依此类推。使用下面的代码，可以从结构上获得所需的输出，但是计算完全不可用：

MySummaryFunc <- function(x) {c(mean = mean(x, na.rm = T), len = sum(!is.na(x)))}

data_avg = summaryBy(Depend.Var1 ~ Subject + Factor1 + Factor2, 
                     data = Orig.data, 
                     FUN = MySummaryFunc, keep.names = TRUE)

有人可以帮助我解决问题吗？