我想感受一下 R 中的函数式编程。为此,我想编写范德蒙德矩阵计算,因为它可能涉及一些构造。
在命令式风格中,这将是:
vandermonde.direct <- function (alpha, n)
{
if (!is.vector(alpha)) stop("argument alpha is not a vector")
if (!is.numeric(alpha)) stop("argument n is not a numeric vector")
m <- length(alpha)
V <- matrix(0, nrow = m, ncol = n)
V[, 1] <- rep(1, m)
j <- 2
while (j <= n) {
V[, j] <- alpha^(j - 1)
j <- j + 1
}
return(V)
}
你会如何在 R 中以函数式风格优雅地编写它?
以下不起作用:
x10 <- runif(10)
n <- 3
Reduce(cbind, aaply(seq_len(n-1),1, function (i) { function (x) {x**i}}), matrix(1,length(x10),1))
因为它告诉我Error: Results must have one or more dimensions.
从i in seq_len(3-1)
函数到函数的函数列表x -> x**i.