Hi! I want to plot the probability density function and the cumulative distribution function for the gamma, lognormal, exponential and Pareto distribution. I want to vary the parameters and and have a plot with 2-3 different parameters in the same figure. It should look like this (example weibull distribution):
http://en.wikipedia.org/wiki/File:Weibull_PDF.svg 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() and http://en.wikipedia.org/wiki/File:Weibull_CDF.svg 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() How do I adapt the syntax to use it for the other distributions? -- View this message in context: http://r.789695.n4.nabble.com/Plotting-probability-density-and-cumulative-distribution-function-tp4303198p4303198.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
