Максимум и минимум можно получить с помощью функции optimize
.
optimize(rho, c(-6,6))
$minimum
[1] -1.000002
$objective
[1] -0.5
optimize(rho, c(-6,6), maximum=TRUE)
$maximum
[1] 1.000002
$objective
[1] 0.5
Вы можете добавить его к графику, используя geom_point
.
MIN = optimize(rho, c(-6,6))
MAX = optimize(rho, c(-6,6), maximum=TRUE)
Points = data.frame(x=c(MIN[[1]][1], MAX[[1]][1]),
y=c(MIN[[2]][1], MAX[[2]][1]))
P1 = ggplot(data.frame(x=c(-6,6)), aes(x=x)) + stat_function(fun=rho, geom="line") + ylab(expression(rho)) + xlab(expression(theta)) +
ggtitle("Autocorrelation at Lag 1 of MA(1)") +
geom_hline(yintercept=0.2, color = "blue", linetype = "dashed") +
scale_color_manual(values = c(NA, "red"))
P1 + geom_point(data=Points, size=2, aes(x,y) )