Tengo 2 matrices de diferentes parámetros: M1 y M3 con las mismas dimensiones. Me gustaría hacer una grangertest de columna en 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)

Querré hacer una prueba de causalidad de granger para determinar si M2 granger causa M1. Mis matrices reales contienen más columnas y filas, pero esto es solo un ejemplo. El código original entre dos vectores está abajo:

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

¿Cómo escribo esto para un análisis inteligente de columnas de modo que una matriz con 2 filas ("F [2]" y "Pr (> F) [2]") y dos columnas se devuelva como resultados, por favor?