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)