У меня есть 2 матрицы разных параметров: M1 и M3 с одинаковыми размерами. Я хотел бы сделать мудрый столбец grangertest в 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)

Я хочу сделать тест причинности Грейнджер, чтобы определить, вызывает ли М2 грейнджер M1. Мои настоящие матрицы содержат больше столбцов и строк, но это только пример. Оригинальный код между двумя векторами ниже:

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

Как мне написать это для анализа по столбцам, чтобы в качестве результата возвращалась матрица с 2 строками ("F [2]" и "Pr (> F) [2]") и двумя столбцами?