在R中绘制参数均值

Question

我有一个名为r的大型实际1-d数据集。 我想绘图：

 mean(log(1+a*r)) vs a, with a > -1 . 

我该怎么做？

Answer 1

您可以使用sapply评估一组值的功能。

x <- seq(-1,1,length=20)
y <- sapply(x, function(a) mean(log(1+a*r)) )
plot(x,y, type="l")

或者，您可以定义一个函数来计算此数量 并使用curve绘制它。 该函数必须进行矢量化。

r <- runif(100)
f <- function(a) mean(log(1+a*r))
f <- Vectorize(f)
curve(f, xlim=c(-1,1), las=1)