Technical Methods for Plotting Multiple Curves with Consistent Scales in R

The Challenge of Scale Consistency in Multi-Curve Plotting

In data visualization practice, it is often necessary to plot multiple data series on the same graph for comparative analysis. However, when different data series have significantly different numerical ranges, directly using R's plot function for overlay plotting leads to scale mismatch issues. For example, consider the following two data series: y1 <- c(100, 200, 300, 400, 500) and y2 <- c(1, 2, 3, 4, 5), which share the same x-values x <- c(1, 2, 3, 4, 5).

Limitations of Traditional Approaches

Beginners often attempt to overlay multiple graphs using the par(new = TRUE) parameter as shown below:

# First graph
plot(x, y1)

# Second graph
par(new = TRUE)
plot(x, y2, axes = FALSE, xlab = "", ylab = "")

Although this method overlays both graphs in the same coordinate system, since the second graph hides its own axes, the two data series actually use different y-axis scales. Specifically, the first graph establishes the y-axis based on the range of y1 (100-500), while the second graph, although not displaying axes, is still plotted based on the range of y2 (1-5), causing y2 data points to be compressed at the bottom of the graph and preventing effective comparison with y1.

Correct Methods for Scale Unification

To address this issue, the key is to ensure that all graphs use the same y-axis range. This can be achieved by setting an appropriate ylim parameter during the first call to the plot function, which should encompass the range of all data series. Here is the corrected code:

# First graph, setting y-axis range to include all data
plot(x, y1, ylim = range(c(y1, y2)))

# Second graph, using the same y-axis range
par(new = TRUE)
plot(x, y2, ylim = range(c(y1, y2)), axes = FALSE, xlab = "", ylab = "")

Here, range(c(y1, y2)) calculates the combined range of y1 and y2 (1-500), ensuring both graphs share the same y-axis scale. Through this approach, y2 data points will be correctly distributed in the graph, forming meaningful comparisons with y1.

Concise Solution Using matplot

In addition to the above method, R provides the more concise matplot function specifically designed for handling multi-curve plotting. This function is optimized for matrix data and automatically manages scale unification for multiple data series. Here are usage examples:

# Basic matplot plotting
matplot(x, cbind(y1, y2))

# matplot plotting with customized point shapes
matplot(x, cbind(y1, y2), pch = 1)

The matplot function takes x and a matrix composed of y1 and y2 cbind(y1, y2) as input, automatically calculating an appropriate y-axis range to accommodate all data series. By default, different series are distinguished by colors, while the pch = 1 parameter sets all data points to hollow circles. This method not only provides cleaner code but also avoids the complexity of manually managing graph overlays.

Technical Summary

The core of achieving consistent scales in multi-curve plotting lies in ensuring all data series are presented within the same coordinate system range. Key steps include: 1) Using the range() function to calculate the combined range of all data series; 2) Explicitly setting this range through the ylim parameter during the first plot call; 3) Maintaining the same ylim settings for subsequent overlay graphs. For more complex scenarios, similar principles may need to be applied to x-axis range unification.

In practical applications, the choice between plot overlays and matplot depends on specific requirements: the former offers greater flexibility allowing fine-grained control over each series; the latter is more suitable for quickly plotting standardized multi-series graphs. Regardless of the method chosen, understanding the principles of scale unification is fundamental to achieving effective data visualization.