Methods and Practices for Plotting Multiple Curves in the Same Graph in R

Introduction

In data visualization and statistical analysis, it is often necessary to plot multiple curves in the same graph for comparative analysis. R, as a powerful statistical analysis tool, provides multiple methods to achieve this functionality. This article starts from the base plotting system and delves into the technical details of plotting multiple curves in the same graph.

Problem Background and Core Challenges

A common issue faced by users when plotting in R is that consecutive calls to the plot() function each create a new graphics window, preventing multiple curves from being displayed in the same graph. This contrasts sharply with MATLAB's hold on functionality and requires different approaches to achieve the same effect.

Solutions in the Base Plotting System

Using the lines() Function for Curve Superposition

In R's base plotting system, the most direct and effective method is to use the plot() function to create the initial graph, then use the lines() or points() functions to add subsequent curves. This approach does not create new graphics windows but instead superimposes directly on the existing graph.

# Generate sample data
x <- seq(-2, 2, 0.05)
y1 <- pnorm(x)  # Standard normal cumulative distribution function
y2 <- pnorm(x, 1, 1)  # Normal cumulative distribution function with mean=1, sd=1

# Create base graph
plot(x, y1, type = "l", col = "red", lwd = 2, 
     xlab = "X Values", ylab = "Cumulative Probability", 
     main = "Two Normal Cumulative Distribution Function Curves")

# Add second curve
lines(x, y2, col = "green", lwd = 2)

# Add legend
legend("topleft", legend = c("Standard Normal", "N(1,1) Distribution"), 
       col = c("red", "green"), lty = 1, lwd = 2)

The advantage of this method lies in its simplicity and intuitiveness. The lines() function is specifically designed to add lines to existing graphs, maintaining graph consistency and integrity.

Using the points() Function to Add Scatter Points

In addition to lines, the points() function can be used to add scatter points to existing graphs:

# Create base scatter plot
plot(x, y1, col = "red", pch = 16, cex = 0.8,
     xlab = "X Values", ylab = "Cumulative Probability")

# Add second set of scatter points
points(x, y2, col = "green", pch = 17, cex = 0.8)

Alternative Approach Using the par() Function

Another method involves using the par(new = TRUE) parameter, which allows superimposing new graphs on existing ones:

# Plot first curve
plot(x, y1, type = "l", col = "red", 
     xlab = "", ylab = "", axes = FALSE)

# Set superposition parameters
par(new = TRUE)

# Plot second curve
plot(x, y2, type = "l", col = "green", 
     xlab = "X Values", ylab = "Cumulative Probability")

It is important to note that this method requires manual handling of axis labels and ticks to avoid overlap and confusion.

Comparative Analysis with Other Tools

Comparison with Matplotlib

In Python's Matplotlib library, implementing multiple curve superposition is more straightforward, requiring only consecutive calls to the plot() function:

import matplotlib.pyplot as plt
import numpy as np

x = np.arange(0, 2*np.pi, 0.05)
y1 = np.sin(x)
y2 = np.cos(x)

plt.plot(x, y1, color='red', label='sin')
plt.plot(x, y2, color='green', label='cos')
plt.legend()
plt.show()

This design in Matplotlib better aligns with user intuition, but R's lines() method offers greater flexibility in controlling graph details.

Comparison with Tableau

In Tableau, multiple curve superposition can be achieved through dual-axis charts or measure value boxes. This approach is more common in business intelligence applications but offers less programming flexibility compared to R.

Advanced Applications and Best Practices

Multi-Curve Color Management

When plotting multiple curves, proper color management is crucial:

# Define color vector
colors <- c("red", "green", "blue", "orange", "purple")

# Create base graph
plot(x, y1, type = "l", col = colors[1], lwd = 2,
     xlab = "X Values", ylab = "Function Values", 
     ylim = range(c(y1, y2, y3)))

# Add multiple curves
lines(x, y2, col = colors[2], lwd = 2)
lines(x, y3, col = colors[3], lwd = 2)

Dynamic Curve Generation and Plotting

In practical applications, dynamic generation and plotting of multiple curves is often required:

# Generate multiple normal distribution curves
mu_values <- c(0, 0.5, 1, 1.5)
sigma_values <- c(0.5, 1, 1.5, 2)

# Create base graph
plot(x, dnorm(x, mu_values[1], sigma_values[1]), 
     type = "l", col = 1, lwd = 2,
     xlab = "X Values", ylab = "Probability Density")

# Loop to add other curves
for(i in 2:length(mu_values)) {
    lines(x, dnorm(x, mu_values[i], sigma_values[i]), 
          col = i, lwd = 2)
}

Common Issues and Solutions

Axis Range Adjustment

When curve value ranges differ significantly, manual setting of appropriate axis ranges is necessary:

# Calculate value range for all curves
y_range <- range(c(y1, y2, y3))

plot(x, y1, type = "l", col = "red", 
     ylim = y_range, xlab = "X Values", ylab = "Function Values")
lines(x, y2, col = "green")
lines(x, y3, col = "blue")

Legend Optimization

Clear legends are crucial for the readability of multi-curve graphs:

legend("topright", 
       legend = c("Curve 1", "Curve 2", "Curve 3"),
       col = c("red", "green", "blue"),
       lty = 1, lwd = 2, 
       bty = "n",  # No border
       cex = 0.8)  # Font size

Performance Optimization Recommendations

When handling large datasets, consider the following optimization measures:

# Reduce data point density
x_sparse <- seq(-2, 2, 0.1)
y1_sparse <- pnorm(x_sparse)
y2_sparse <- pnorm(x_sparse, 1, 1)

plot(x_sparse, y1_sparse, type = "l", col = "red")
lines(x_sparse, y2_sparse, col = "green")

Conclusion

In R, using the lines() and points() functions to plot multiple curves in the same graph is the most reliable and efficient method. This approach not only maintains code simplicity but also provides sufficient flexibility to customize various aspects of the graph. Compared to the par(new = TRUE) method, the lines() method is more stable and easier to maintain. In practical applications, combining proper color management, axis settings, and legend design can create professional-level multi-curve comparison graphs.