R base graphics gallery

png("fig01.png")

plot(iris[iris$Species == "setosa", ]$Sepal.Length,
     iris[iris$Species == "setosa", ]$Petal.Length,
     col = "red", pch = 1,
     xlab = "Sepal Length", ylab = "Petal Length",
     xlim = range(iris$Sepal.Length), ylim = range(iris$Petal.Length))
points(iris[iris$Species == "versicolor", ]$Sepal.Length,
       iris[iris$Species == "versicolor", ]$Petal.Length,
       col = "green", pch = 2)
points(iris[iris$Species == "virginica", ]$Sepal.Length,
       iris[iris$Species == "virginica", ]$Petal.Length,
       col = "blue", pch = 3)
legend(6.5, 3,
       title = "Species",
       legend = c("Setosa", "Versicolor", "Virginica"),
       col = c("red", "green", "blue"),
       pch = c(1, 2, 3),
       bg = "orange")

dev.off()

png("fig02.png")

x <- seq(0, 2, length.out = 100)
foo <- function(x) x
bar <- function(x) exp(-(x - 1)^2) * sin(x * 8)
points_data <- data.frame(x = seq(0, 2, length = 10), y = rnorm(10))
plot(x, foo(x), type = "l", col = "blue", xlim = c(0, 2), ylim = c(-1.5, 1.5),
     xlab = "x", ylab = "y", main = "Function", cex.main = 1.5)
lines(x, bar(x), col = "red")
points(points_data$x, points_data$y, col = "green")
legend("topright", legend = c("Foo", "Bar", "Baz"), lty = 1, col = c("blue", "red", "green"),
         cex = 0.8, title = expression(bold(italic("Function"))))

dev.off()

png("fig03.png")

data(faithful)

par(bg = "lightgrey")
hist(faithful$waiting,
     main = expression(italic("Waiting Time between Eruptions")),
     breaks = seq(40, 100, by = 5),
     xlab = NULL, ylab = NULL,
     col = "blue", border = "orange", bg = "grey")
axis(side = 1, col.axis = "blue", col = "orange")
title(xlab = "Waiting Time (mins)",
        col.lab = "orange", line = 3)
axis(side = 2, col.axis = "orange", col = "green")
title(ylab = "Frequency",
        col.lab = "green", line = 3)

dev.off()

png("fig04.png")

data(iris)

par(bg = "lightblue", mar = c(5, 4, 4, 2) + 0.1)

boxplot(iris$Sepal.Length ~ iris$Species,
          main = "Distribution of Sepal Length in Iris Dataset",
          col.main = "purple",
          xlab = expression(italic("Species")),
          ylab = "Sepal Length (cm)",
          col = "purple",
          border = "red",
          notch = TRUE,
          boxwex = 0.5,
          staplewex = 0.5,
          outpch = 16,
          outcol = "blue",
          horizontal = FALSE)

dev.off()

png("fig05.png")

nrows <- 5
ncols <- 5
data <- matrix(runif(nrows * ncols), nrow = nrows)
rownames <- as.character(5:9)
colnames <- letters[1:5]

heatmap(data,
          Rowv = TRUE, # include row dendrogram
          Colv = NA,
          xlab = "X Axis Label",
          ylab = "Y Axis Label",
          main = expression(bold(italic("Main Title"))),
          labRow = rownames,
          labCol = colnames,
          col = heat.colors(256))

dev.off()

png("fig06.png")

n <- 100
x <- seq(0, 2, length.out = n)

y_lines <- matrix(NA, nrow = n, ncol = 2)
y_lines[,1] <- sin(x)
y_lines[,2] <- cos(x)
y_points <- matrix(NA, nrow = n, ncol = 2)
y_points[,1] <- rnorm(n = n, mean = sin(x))
y_points[,2] <- rnorm(n = n, mean = cos(x))

matplot(x, y_lines, type = "l", lwd = 3,
        ylim = c(min(y_points), max(y_points)),
        col = c("red", "blue"),
        ylab = "my ylab", xlab = "my xlab")
matpoints(x, y_points, pch = 1, col = c("red", "blue"))
text(x = 1, y = 0,
     labels = "matplot, matpoints and matlines",
     srt = 30, cex = 2.5, font = 4)

dev.off()