Então, olhei para esta pergunta postada antes, que foi usada para somar a cada 2 valores em cada linha de uma matriz. Aqui está o link:soma colunas específicas entre linhas. Também dei uma olhada em outra pergunta aqui:R Soma todas as k colunas na matriz o que é mais parecido com o meu. Não foi possível obter a solução neste caso para funcionar. Aqui está o código com o qual estou trabalhando ...

y <- matrix(1:27, nrow = 3)
y

m1 <- as.matrix(y)
n <- 3
dim(m1) <- c(nrow(m1)/n, ncol(m1), n)
res <- matrix(rowSums(apply(m1, 1, I)), ncol=n)
identical(res[1,],rowSums(y[1:3,]))


sapply(split.default(y, 0:(length(y)-1) %/% 3), rowSums)

Acabei de receber uma mensagem de erro ao aplicar isso. A saída desejada é uma matriz com os seguintes valores:

      [,1] [,2] [,3]
[1,]   12   39   66
[2,]   15   42   69
[3,]   18   45   72