График параметрического среднего в R

Question

У меня большой реальный 1-й набор данных, который называется r.Я хотел бы сюжет:

 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)