R,未知数 vectors/matrices 的成对乘积
R, pair-wise product of unknown number of vectors/matrices
我想在基数 R 中生成可变数量 matrices/vectors 的成对乘积。
我只有这个丑陋的解决方案(丑陋的是 <<-
),但直觉上认为存在更好的 - 可能是递归的 - 方式,甚至可能是一个函数。我需要 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
我想在基数 R 中生成可变数量 matrices/vectors 的成对乘积。
我只有这个丑陋的解决方案(丑陋的是 <<-
),但直觉上认为存在更好的 - 可能是递归的 - 方式,甚至可能是一个函数。我需要 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