在R中绘制双轴图表：基于生产总量W数据集的技术问询

Got it, let's break down how to build a dual-axis chart in R using your sens dataset. First, I’m assuming your data frame includes columns for W (total production), t, and h (the two variables influencing W). We’ll use ggplot2—the go-to package for flexible, polished visualizations—since it makes dual-axis setups straightforward once you get the hang of scaling.

Step 1: Load Required Packages #

First, make sure you have these packages installed (if not, run install.packages(c("ggplot2", "scales")) first):

library(ggplot2)
library(scales)

Step 2: Verify Your Data Structure #

Before plotting, double-check your data has all the columns you need. Run this to inspect the first few rows and column names:

# Check data structure
str(sens)
# View first 6 rows
head(sens)

This confirms you have t, h, and W columns to work with.

Step 3: Build the Dual-Axis Chart #

Dual-axis charts work by plotting two datasets on the same x-axis, then scaling one of the y-variables to match the range of the other (so they fit on the same plot). Below is an example where we use t as the x-axis, left y-axis for W (production), and right y-axis for h:

# Create the base plot
dual_axis_plot <- ggplot(sens, aes(x = t)) +
  # First layer: W vs t (solid line + points)
  geom_line(aes(y = W, color = "Total Production (W)"), linewidth = 1) +
  geom_point(aes(y = W, color = "Total Production (W)"), size = 2) +
  
  # Second layer: h vs t (dashed line + points, scaled to W's range)
  geom_line(aes(y = h * (max(sens$W)/max(sens$h)), color = "Variable h"), 
            linewidth = 1, linetype = "dashed") +
  geom_point(aes(y = h * (max(sens$W)/max(sens$h)), color = "Variable h"), size = 2) +
  
  # Configure y-axes: left for W, right for scaled h
  scale_y_continuous(
    name = "Total Production (W)",
    sec.axis = sec_axis(~ . / (max(sens$W)/max(sens$h)), 
                        name = "Variable h")
  ) +
  
  # Customize colors and labels
  scale_color_manual(values = c("Total Production (W)" = "#2980b9", "Variable h" = "#e74c3c")) +
  labs(
    title = "Production (W) vs Variables t and h",
    x = "Variable t",
    color = "Metric"
  ) +
  
  # Clean up the theme
  theme_minimal() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 14, face = "bold"),
    axis.title = element_text(size = 12),
    legend.position = "top",
    legend.title = element_text(face = "bold")
  )

# Display the plot
print(dual_axis_plot)

Key Notes & Adjustments #

• Swap Variables: If you want h as the x-axis instead of t, just replace x = t with x = h, and adjust the second layer to use t instead of h.
• Scaling: The scaling factor max(sens$W)/max(sens$h) ensures h fits within the range of W. If your variables have different ranges, this keeps the plot readable.
• Caution: Dual-axis charts can be misleading if not used carefully—readers might misinterpret the relationship between variables. Only use this if you truly need to compare two metrics with vastly different scales.

If your dataset is missing t or h columns (or has different column names), just adjust the code to match your actual column names!

内容的提问来源于stack exchange，提问作者anup