R中的sapply或其他apply函数

Question

昨天我写了一些代码，这使您困惑。 对于那个很抱歉。 所以我以一种更简单的方式再次写了它。

我的问题是：有没有更简单（或更快速）的方法来实现以下代码？

k <- c(.04, .08, .12, .16, .2);

library(plyr)
valfcn <- function(k, V_next){
  a <- .3;
  b <- .6;

 return_val <- vector()
  for(i in 1:5){
    tmp <- vector()
    for(j in 1:5){
  tmp[j] <- (log(k[i]^a - k[j]) + b*V_next[j]);
}
return_val <- c(return_val,max(tmp[i]))
}
  return_val

}

V0 <- c(rep(0,5))
V1 <- valfcn(k,V0)
V2 <- valfcn(k,V1)

V1
V2

我想使用替代方法，它可能更短但更快，而不是使用for循环方法。

最好！

Answer 1

我相信sapply()根据您的描述不是必需的。 这样的事情可能会满足您的需求：

valfcn <- function(k, V_next){
  a <- .3;
  b <- .6;
  max(log(k^a - k) + b*V_next);
}

在此版本中，对k进行的转换将生成一个向量，然后max()对整个向量进行运算。 不需要使用循环或使用sapply() ，因为max()会处理它。