我想知道如何同时使用sapply
和lapply
以便我可以避免编写名为GG
的函数,如下所示?
GG = function(x, y) dnorm(250, mean = x, sd = y)*dnorm(265, mean = x, sd = y) *
dnorm(259, mean = x, sd = y)
P.S。我知道如果我上面的函数中只有x
变化,以下情况可能有效:
function(x) sapply(lapply(x, dnorm, x = c(250, 265, 259), 10), prod)
但在我的情况下x
和y
都有所不同。
答案 0 :(得分:3)
我们可以将Map
与Reduce
base R
一起使用。使用Map
的原因是函数可以应用于传递给它的对象的相应元素。在这里,dnorm
是获取' x'的每个对应元素的函数。并且' y'作为mean
和sd
参数,它有一个常数向量" x" (c(250, 265, 259)
)。 Map
的输出为list
,我们Reduce
通过乘以(list
)
*
的相应元素添加到单个元素中
GG1 <- function(x, y) Reduce(`*`, Map(dnorm, x = c(250, 265, 259),
mean = list(x), sd = list(y)))
identical(GG(24, 12), GG1(24, 12))
#[1] TRUE
identical(GG(32, 15), GG1(32, 15))
#[1] TRUE
根据OP的评论,
x <- seq(10, 40, length= 30)
y <- x
z <- outer(x, y, GG1)
persp(x, y, z , theta = 0, phi = 20, expand = 0.5, col = 'pink')