Як відтінку графа за допомогою кривої () в R - r, сюжеті, графіку, області

Я будую стандартний нормальний розподіл.

curve(dnorm(x), from=-4, to=4,
main = "The Standard Normal Distibution",
ylab = "Probability Density",
xlab = "X")

З педагогічних міркувань я хочу затінити площу нижче певного квантилу на мій вибір. Як я можу це зробити?

Відповіді:

Якщо ви хочете використовувати curve і базовий сюжет, тоді ви можете написати невелику функцію самостійно polygon:

colorArea <- function(from, to, density, ..., col="blue", dens=NULL){
y_seq <- seq(from, to, length.out=500)
d <- c(0, density(y_seq, ...), 0)
polygon(c(from, y_seq, to), d, col=col, density=dens)
}

Маленький приклад випливає:

curve(dnorm(x), from=-4, to=4,
main = "The Standard Normal Distibution",
ylab = "Probability Density",
xlab = "X")

colorArea(from=-4, to=qnorm(0.025), dnorm)
colorArea(from=qnorm(0.975), to=4, dnorm, mean=0, sd=1, col=2, dens=20)