Eu tenho 2 matrizes de parâmetros diferentes: M1 e M3 com as mesmas dimensões. Eu gostaria de fazer uma coluna mais sábia em R.

M1<- matrix( c(2,3, 1, 4, 3, 3, 1,1, 5, 7), nrow=5, ncol=2) 
M3<- matrix( c(1, 3, 1,5, 7,3, 1, 3, 3, 4), nrow=5, ncol=2)

Vou querer fazer um teste de causalidade de granger para determinar se M2 granger causa M1. Minhas matrizes reais contêm mais colunas e linhas, mas este é apenas um exemplo. O código original entre dois vetores está abaixo:

library(lmtest)
data(ChickEgg)
grangertest(chicken ~ egg, order = 3, data = ChickEgg)

Como escrevo isso para uma análise sábia de colunas, de modo que uma matriz com 2 linhas ("F [2]" e "Pr (> F) [2]") e duas colunas seja retornada como resultado, por favor?