Question

我想生成一个矩阵，其值是对具有两个参数的函数的评估。

到目前为止，以下是我想要的，但我怀疑在R中有一种更好的方法可以做到这一点，而不使用for循环。

x <- c(0.01, 0.1,1,10,100)
n <- length(x)
M <- matrix(NA, nrow=n, ncol=n)
for (i in 1:n){
 for (j in 1:n){
  M[i,j]<-x[i]^x[j]
 }
}

Answer 1

使用outer功能：

> outer(x,x,"^")
          [,1]      [,2]  [,3]  [,4]   [,5]
[1,] 0.9549926 0.6309573 1e-02 1e-20 1e-200
[2,] 0.9772372 0.7943282 1e-01 1e-10 1e-100
[3,] 1.0000000 1.0000000 1e+00 1e+00  1e+00
[4,] 1.0232930 1.2589254 1e+01 1e+10 1e+100
[5,] 1.0471285 1.5848932 1e+02 1e+20 1e+200
> identical(M,outer(x,x,"^"))
[1] TRUE

您可以使用函数名称或运算符的引用名称。另请注意，我测试outer的答案与循环中的答案相同。一直在测试！

Answer 2

使用outer功能：

outer(x,x,"^")

          [,1]      [,2]  [,3]  [,4]   [,5]
[1,] 0.9549926 0.6309573 1e-02 1e-20 1e-200
[2,] 0.9772372 0.7943282 1e-01 1e-10 1e-100
[3,] 1.0000000 1.0000000 1e+00 1e+00  1e+00
[4,] 1.0232930 1.2589254 1e+01 1e+10 1e+100
[5,] 1.0471285 1.5848932 1e+02 1e+20 1e+200

在没有for循环的情况下在R中填充函数求值矩阵

2 个答案: