Question

我有以下数据框：

> coc_comp_model[1:3,]
    Relationship Output  Input |r-Value|   Y-Intercept   Gradient
1   DG-r ~ DG-cl   DG-r  DG-cl 0.8271167  0.0027217513 12.9901380
2 CA3-r ~ CA3-cl  CA3-r CA3-cl 0.7461309  0.0350767684 27.6107963
3 CA2-r ~ CA2-cl  CA2-r CA2-cl 0.9732584 -0.0040992226 35.8299582

我想为数据框的每一行创建简单的函数。这是我尝试过的：

for(i in 1:nrow(coc_comp_model)) {
coc_glm_f[i] <- function(x)
x*coc_comp_model$Gradient[i] + coc_comp_model$Y-Intercept[i]
}

也尝试制作一个函数向量，它也可以工作。

感谢您阅读本文/帮助。

Answer 1

这样的事情：

myfunc<-function(datrow, x){
  x*as.numeric(datrow[6]) + as.numeric(datrow[5] )
}

然后您可以使用apply在每一行上调用它，根据需要更改x：

apply(hzdata, 1, myfunc, x = 0.5)

注意：使用dput()分享您的数据要比子集中的粘贴容易得多。

Answer 2

没有函数向量这样的东西。 R中有6种原子向量类型：raw，logical，integer，double，complex和character，还有异构列表类型，最后还有鲜为人知的表达式类型，它基本上是解析树的向量（如从调用substitute()函数开始。这些都是R中的矢量类型。

printAndType <- function(x) { print(x); typeof(x); };
printAndType(as.raw(1:3));
## [1] 01 02 03
## [1] "raw"
printAndType(c(T,F));
## [1]  TRUE FALSE
## [1] "logical"
printAndType(1:3);
## [1] 1 2 3
## [1] "integer"
printAndType(as.double(1:3));
## [1] 1 2 3
## [1] "double"
printAndType(c(1i,2i,3i));
## [1] 0+1i 0+2i 0+3i
## [1] "complex"
printAndType(letters[1:3]);
## [1] "a" "b" "c"
## [1] "character"
printAndType(list(c(T,F),1:3,letters[1:3]));
## [[1]]
## [1]  TRUE FALSE
##
## [[2]]
## [1] 1 2 3
##
## [[3]]
## [1] "a" "b" "c"
##
## [1] "list"
printAndType(expression(a+1,sum(1,2+3*4),if (T) 1 else 2));
## expression(a + 1, sum(1, 2 + 3 * 4), if (T) 1 else 2)
## [1] "expression"

如果要在单个对象中存储多个函数，则必须使用列表，并且必须使用左值中的双括号索引运算符来分配它：

fl <- list();
for (i in 1:3) fl[[i]] <- (function(i) { force(i); function(a) a+i; })(i);
fl;
## [[1]]
## function (a)
## a + i
## <environment: 0x600da11a0>
##
## [[2]]
## function (a)
## a + i
## <environment: 0x600da1ab0>
##
## [[3]]
## function (a)
## a + i
## <environment: 0x600da23f8>
sapply(fl,function(f) environment(f)$i);
## [1] 1 2 3
sapply(fl,function(f) f(3));
## [1] 4 5 6

在上面的代码中，我还演示了关闭循环变量的正确方法。这需要创建一个临时函数评估环境来保存i的副本，然后返回的函数将关闭该评估环境，以便它可以访问特定于迭代的i。这适用于支持动态函数和闭包的其他语言，例如JavaScript。在R中还有一个额外的要求，即强制承诺通过force()来解决，否则，对于每个生成的函数，只有在特定生成函数的第一次评估之前，承诺才会得到解决。那时锁定该特定生成函数的promise目标（在这种情况下为全局i变量）的当前值。还应该提到的是，这是一个非常浪费的设计，为每次迭代生成一个临时函数并对其进行评估，从而生成一个带有循环变量副本的新评估环境。

如果您想使用此设计，那么您的代码将变为：

coc_glm_f <- list();
for (i in 1:nrow(coc_comp_model)) {
    coc_glm_f[[i]] <- (function(i) { force(i); function(x) x*coc_comp_model$Gradient[i] + coc_comp_model$`Y-Intercept`[i]; })(i);
};

但是，为data.frame的每一行创建一个单独的函数可能没有意义。如果您希望x参数采用标量值（我的意思是单元素向量），那么您可以按如下方式定义函数：

coc_glm_f <- function(x) x*coc_comp_model$Gradient + coc_comp_model$`Y-Intercept`;

此函数是矢量化的，这意味着您可以传递x的向量，其中x的每个元素都对应一行coc_comp_model。例如：

coc_comp_model <- data.frame(Relationship=c('DG-r ~ DG-cl','CA3-r ~ CA3-cl','CA2-r ~ CA2-cl'),Output=c('DG-r','CA3-r','CA2-r'),Input=c('DG-cl','CA3-cl','CA2-cl'),`|r-Value|`=c(0.8271167,0.7461309,0.9732584),`Y-Intercept`=c(0.0027217513,0.0350767684,-0.0040992226),Gradient=c(12.9901380,27.6107963,35.8299582),check.names=F);
coc_glm_f(seq_len(nrow(coc_comp_model)));
## [1]  12.99286  55.25667 107.48578

R - 在for循环中命名多个函数

2 个答案: