Я хочу построить функцию плотности вероятности и кумулятивную функцию распределения для гамма, логнормального, экспоненциального распределения и распределения Парето. Я хочу варьировать параметры и иметь график с 2-3 различными параметрами на одном рисунке. Это должно выглядеть следующим образом (пример распределения Weibull):
Weibull_PDF
library(Cairo)
CairoFonts(regular="DejaVu Sans:style=Regular")
CairoSVG("Weibull PDF.svg")
par(mar=c(3, 3, 1, 1))
x <- seq(0, 2.5, length.out=1000)
plot(x, dweibull(x, .5), type="l", col="blue", xlab="", ylab="", xlim=c(0, 2.5), ylim=c(0, 2.5), xaxs="i", yaxs="i")
lines(x, dweibull(x, 1), type="l", col="red")
lines(x, dweibull(x, 1.5), type="l", col="magenta")
lines(x, dweibull(x, 5), type="l", col="green")
legend("topright", legend=paste("\u03bb = 1, k =", c(.5, 1, 1.5, 5)), lwd=1, col=c("blue", "red", "magenta", "green"))
dev.off()
и
Weibull_CDF
library(Cairo)
CairoFonts(regular="DejaVu Sans:style=Regular")
CairoSVG("Weibull CDF.svg")
par(mar=c(3, 3, 1, 1))
x <- seq(0, 2.5, length.out=1000)
plot(x, pweibull(x, .5), type="l", col="blue", xlab="", ylab="", xlim=c(0, 2.5), ylim=c(0, 1), xaxs="i", yaxs="i")
lines(x, pweibull(x, 1), type="l", col="red")
lines(x, pweibull(x, 1.5), type="l", col="magenta")
lines(x, pweibull(x, 5), type="l", col="green")
legend("bottomright", legend=paste("\u03bb = 1, k =", c(.5, 1, 1.5, 5)), lwd=1, col=c("blue", "red", "magenta", "green"))
dev.off()
Как адаптировать синтаксис, чтобы использовать его для других дистрибутивов?