我有一个函数列表
funs <- list(fn1 = function(x) x^2,
fn2 = function(x) x^3,
fn3 = function(x) sin(x),
fn4 = function(x) x+1)
#in reality these are all f = splinefun()
我有一个数据框:
mydata <- data.frame(x1 = c(1, 2, 3, 2),
x2 = c(3, 2, 1, 0),
x3 = c(1, 2, 2, 3),
x4 = c(1, 2, 1, 2))
#actually a 500x15 dataframe of 500 samples from 15 parameters
对于每个i行,我想评估函数j在每个j列并对结果求和:
unlist(funs)
attach(mydata)
a <- rep(NA,4)
for (i in 1:4) {
a[i] <- sum(fn1(x1[i]), fn2(x2[i]), fn3(x3[i]), fn4(x4[i]))
}
我怎样才能有效地做到这一点?现在是实施的适当时机吗?plyr
功能?如果是这样,怎么办?
额外问题:为什么a[4]
NA
?
现在是使用以下函数的合适时机吗plyr
,如果可以,我该怎么做?