Мощность матрицы в R

Это неправильный ответ, но он может быть хорошим местом для обсуждения и понимания внутренней работы R. Этот вид ошибок уже подкрался раньше в другом пакете I использовал.

Во-первых, обратите внимание, что простое присвоение матрицы сначала новой переменной не помогает:

> A <- B <-matrix(c(1,3,0,2,8,4,1,1,1),nrow=3)
> r1 <- A %^% 5
> A
     [,1] [,2] [,3]
[1,]  691 1926  312
[2,] 2331 6502 1056
[3,] 1116 3108  505
> B
     [,1] [,2] [,3]
[1,]  691 1926  312
[2,] 2331 6502 1056
[3,] 1116 3108  505

Я предполагаю, что R пытается быть умным, передавая по ссылке вместо значений. Чтобы заставить это работать, вам нужно что-то сделать, чтобы отличить A от B:

`%m%` <- function(x, k) {
    tmp <- x*1
    res <- tmp%^%k
    res
}
> B <-matrix(c(1,3,0,2,8,4,1,1,1),nrow=3)
> r2 <- B %m% 5
> B
     [,1] [,2] [,3]
[1,]    1    2    1
[2,]    3    8    1
[3,]    0    4    1

Каков явный способ сделать это?

Наконец, в коде C для пакета есть следующий комментарий:

• NB: x будет изменен! При необходимости вызывающий должен сделать копию

Но я не понимаю, почему R позволяет коду C / Fortran иметь побочные эффекты в глобальной среде.

Хотя исходный код не отображается в пакете, поскольку он упакован в файл .dll , я считаю, что в пакете используется алгоритм алгоритм быстрого возведения в степень , который вы можете изучить, взяв вместо этого функцию с именем matpowfast .

Вам нужны две переменные:

1. result , чтобы сохранить результат,
2. mat в качестве промежуточной переменной.

Чтобы вычислить A ^ 6 , поскольку 6 = 110 (двоичная запись), в итоге result = A ^ 6 и mat = A ^ 4 . То же самое для A ^ 5 .

Вы можете легко проверить, соответствует ли mat = A ^ 8 , когда вы попытаетесь вычислить A ^ n для любого 8 . Если да, то у вас есть объяснение.