Rのパラメトリック平均をプロット

Question

私はrと呼ばれる1dの大きなデータセットを持っています。 私はプロットをしたいと思います：

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)