映射基础知识？ - 如何从两个向量和一个函数创建一个矩阵

Question

我正在尝试创建一个 data.frame 来创建图表。我有一个函数和两个向量，我想将它们用作两个输入。这有点简化，但基本上我所拥有的是：

relGPA <- seq(-1.5,1.5,.2)
avgGPA <- c(-2,0,2)

f <- function(relGPA, avgGPA) 1/(1+exp(sum(relGPA*pred.model$coef[1],avgGPA*pred.model$coef[2])))

我想要的只是一个 data.frame，其中 3 列用于 avgGPA 值，16 行用于 relGPA 值以及单元格中的结果值。

对于这是多么基本的事情，我深表歉意，但我向你保证，我已经尽力在没有你的帮助的情况下实现这一点。我已经尝试遵循 sapply 和 mapply 手册页上的示例，但我对 R 有点太陌生，无法了解我正在尝试做什么。

谢谢！

Answer 1

无法使用所提供的信息进行测试，但这应该可行：

expGPA  <- outer(relGPA, avgGPA, FUN=f) # See below for way to make this "work"

当您想要生成组合时，另一个有用的函数是 expand.grid ，这将为您提供“长形式”：

expGPA2 <-expand.grid(relGPA, avgGPA)
expGPA2$fn <- apply(expGPA2, 1, f)

长形式是什么网格和 ggplot 将期望作为更高级别绘图的输入格式。

编辑：可能需要构造一个更具体的方法来将列引用传递给函数，如 djhurio 所指出的和 Sam Swift 使用 Vectorize 策略（已解决）。在 apply 的情况下，sum 函数可以像上面描述的那样开箱即用，但除法运算符则不然，所以这里有一个可以推广的进一步示例具有多个参数的更复杂的函数。程序员所需要的只是“apply()”-ed”函数中适当参数的列号，因为（不幸的是）列名称没有传递到 x 参数：

> expGPA2$fn <- apply(expGPA2, 1, function(x) x[1]/x[2])
> str(expGPA2)
'data.frame':   48 obs. of  3 variables:
 $ Var1: num  -1.5 -1.3 -1.1 -0.9 -0.7 ...
 $ Var2: num  -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 ...
 $ fn  : num  0.75 0.65 0.55 0.45 0.35 ...
 - attr(*, "out.attrs")=List of 2
  ..$ dim     : int  16 3
  ..$ dimnames:List of 2
  .. ..$ Var1: chr  "Var1=-1.5" "Var1=-1.3" "Var1=-1.1" "Var1=-0.9" ...
  .. ..$ Var2: chr  "Var2=-2" "Var2= 0" "Var2= 2"

Edit2 : (2013-01-05) 一年后看到这个，我意识到 SamSwift 的函数可以通过使其主体使用“+”而不是 sum 进行矢量化：

 1/(1+exp( relGPA*pred.model$coef[1] + avgGPA*pred.model$coef[2]) # all vectorized fns