Блог

Я хочу создавать попарные произведения переменных чисел матриц/векторов в базе R. У меня есть только это уродливое решение (уродливое <<- ), но интуитивно я думаю, что существует более приятный — возможно, рекурсивный — способ или, возможно, даже функция. Мне нужна парная версия prod .

 f1 <- function(...) {
  input <- list(...)
  output <- input[[1]]
  sapply(2:length(input), function(m) output <<- output*input[[m]])
  return(output)
}

m1 <- matrix(1:6, ncol = 2)
m2 <- matrix(6:1, ncol = 2)
m3 <- 1/matrix(6:1, ncol = 2)

all(f1(m1,m2,m3) == m1*m2*m3) #[1] TRUE
 

Воспользуйся Reduce :

 f1 <- function(...) {
  Reduce(`*`, list(...))
}

all(f1(m1,m2,m3) == m1*m2*m3)
#[1] TRUE