所以我有一个函数列表.我想创建一个for循环,它返回(显然作为一个函数)它们的总和.
为了在for循环中创建函数列表,我正在使用此代码
##CODE
f=dnorm
h=function(x){log(f(x))}
S=c(-3,-2,-1,0,1,2,3)
K=matrix(rep(1:length(S),2),ncol=2)
for(i in 1:length(S)){
K[i,]=c(S[i],h(S[i]))
}
funcs=list()
## LOOP TO DEFINE THE LINES
for(i in 1:6){
## Make function name
funcName <- paste( 'hl', i,i+1, sep = '' )
## Make function
func1 = paste('function(x){ (K[',i,'+1,2]-K[',i,',2])/(K[',i,'+1,1]-K[',i,',1])*x+
K[',i,'+1,2]-((K[',i,'+1,2]-K[',i,',2])/(K[',i,'+1,1]-K[',i,',1]))*K[',i,'+1,1]}',sep
= '')
funcs[[funcName]] = eval(parse(text=func1))
}
Run Code Online (Sandbox Code Playgroud)
它创建了6个函数的列表.我怎样才能得到他们的金额?我尝试使用apply命令但是我的语法不正确或者它们不起作用.
PS我实际上是在尝试为ars命令编写一个代码.
正如尼克指出的那样,"功能的总和"没有意义.我疯狂地猜测你想在某个时刻(在S?)评估函数,然后取这些值的总和.这应该可以解决问题.
rowSums(sapply(funcs, function(f) f(S)))
Run Code Online (Sandbox Code Playgroud)
您的大部分代码都可以更加干净地编写,并以矢量化方式编写.
f <- dnorm
h <- function(x) log(f(x))
S <- -3:3
K <- cbind(S, h(S)) #No need to define this twice; no need to use a loop
i <- seq_len(6)
funcNames <- paste('hl', i, i+1, sep = '') #paste is vectorised
#You can avoid using `paste`/`eval`/`parse` with this function to create the functions
#Can possibly be done even more cleanly by using local
makeFunc <- function(i)
{
evalq(substitute(
function(x)
{
(K[i + 1, 2] - K[i, 2]) / (K[i + 1, 1] - K[i, 1]) * x +
K[i + 1, 2] -
((K[i + 1, 2] - K[i, 2]) / (K[i + 1, 1] - K[i, 1])) * K[i + 1, 1]
},
list(i = i)
))
}
funcs <- lapply(i, makeFunc)
names(funcs) <- funcNames
rowSums(sapply(funcs, function(f) f(S)))
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
1526 次 |
| 最近记录: |