# Produces four plots of different density functions in a 2x2 window. # The exponential, gamma, beta, and normal densities are shown. # Set up a 2x2 plot. par(mfrow = c(2, 2)) # Setting the parameters 'xaxs' and 'yaxs' to 'i' (for "internal") lets us # specify the plot windows exactly. By default, plot windows are enlarged # by about 6%. Using the "internal" setting prevents this enlargement. par(xaxs = "i", yaxs = "i") # Set some other style parameters. line.width = 2 label.size = 1.3 # Plot the exponential density. x <- seq(from = 0, to = 8, by = 0.01) exponential.pdf <- dexp(x) plot(x, exponential.pdf, ylim = c(0, 1), type = "l", lwd = line.width, main = "", xlab = "", ylab = "", bty = "n", cex.lab = label.size) # Plot the Gamma density. x <- seq(from = 0, to = 10, by = 0.01) gamma.pdf <- dgamma(x, 2.5, 1.25) plot(x, gamma.pdf, ylim = c(0, 0.4), type = "l", lwd = line.width, main = "", xlab = "", ylab = "", bty = "n", cex.lab = label.size) # Plot the Beta density. x <- seq(from = 0, to = 1, by = 0.01) beta.pdf = dbeta(x, 4, 4) plot(x, beta.pdf, ylim = c(0, 2.5), type = "l", lwd = line.width, main = "", xlab = "", ylab = "", bty = "n", cex.lab = label.size) # Plot the Normal density. x <- seq(from = -4, to = 4, by = 0.01) normal.pdf <- dnorm(x) plot(x, normal.pdf, ylim = c(0, 0.4), type = "l", lwd = 2, main = "", xlab = "", ylab = "", bty = "n", cex.lab = 1.3) # Print the image to a ps. dev.print(device = postscript, "3.3.eps", horizontal = TRUE)