Предположим, что A - некоторая квадратная матрица. Как легко измерить эту матрицу в R?
Я уже пробовал два пути: Trial 1 с хакером for-loop и Trial 2 немного элегантнее, но он все еще далек от простоты A k.
Испытание 1
set.seed(10)
t(matrix(rnorm(16),ncol=4,nrow=4)) -> a
for(i in 1:2){a <- a %*% a}
Испытание 2
a <- t(matrix(c(0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0),nrow=4))
i <- diag(4)
(function(n) {if (n<=1) a else (i+a) %*% Recall(n-1)})(10)