我需要创建一个函数(x,k),x =长度为n的向量,k =整数。 该函数需要给我一个维度矩阵[nxk],并且每列中的条目需要被提升到该列数的幂(即第一列中的条目是x,第2列中的条目是x ^ 2等)。 我很难搞清楚如何构建一个按列进行此类操作的函数。
非常感谢你。
答案 0 :(得分:1)
这样的事情可能会利用outer
,因为将函数应用于两个向量而返回matrix
。
matpower <- function(x,k) outer(x,seq_len(k),`^`)
matpower(1:4,4)
# [,1] [,2] [,3] [,4]
#[1,] 1 1 1 1
#[2,] 2 4 8 16
#[3,] 3 9 27 81
#[4,] 4 16 64 256