如何绘制1880-2017年月度气温时间序列折线图？解决x/y长度不匹配报错

Hey there! Let's troubleshoot that frustrating error you're getting when trying to switch from a scatter plot to a line plot for your 1880-2017 monthly temperature data. Here are the most common causes and fixes:

1. First, Verify Your Data Lengths Match #

The error directly tells you that your x-axis (time) and y-axis (temperature) vectors don't have the same number of observations. Start by checking this with simple code:

# Replace with your actual column names
length(your_data$date_column)
length(your_data$temperature_column)

If these numbers don't match, you've found the root issue. This usually happens because:

• Missing values in one column got dropped accidentally (e.g., you filtered temperatures but didn't apply the same filter to dates)
• Your date column has invalid entries that failed to parse into a time format, leaving behind NA values that get excluded silently

2. Clean and Standardize Your Time Column #

Most of the time, this error pops up because your date data isn't properly formatted as a time object. If your dates are stored as strings (like "1880-01"), convert them to a time format first to avoid parsing issues:

Option 1: Use zoo's as.yearmon (great for monthly data) #

library(zoo)
# Convert string dates to year-month objects
your_data$date <- as.yearmon(your_data$date_column, format = "%Y-%m")
# Check for NA values from bad date formats
sum(is.na(your_data$date))

Option 2: Base R Date object (add a dummy day like the 1st) #

your_data$date <- as.Date(paste0(your_data$date_column, "-01"), format = "%Y-%m-%d")
sum(is.na(your_data$date))

If there are NA values here, go back to your raw data and fix any malformed date entries.

3. Handle Missing Values Consistently #

If you have missing temperatures, make sure you filter out the corresponding date rows too. Don't just drop NA temps without touching dates—this will create a length mismatch:

# Remove rows with NA in either date or temperature
clean_data <- na.omit(your_data)
# Double-check lengths now
length(clean_data$date) == length(clean_data$temperature)

4. Plot Explicitly with X and Y Defined #

Avoid vague plot calls like plot(your_data, type="l")—this can confuse R when your data frame has multiple columns. Instead, explicitly specify your x and y axes:

plot(x = clean_data$date, 
     y = clean_data$temperature,
     type = "l",
     xlab = "Year",
     ylab = "Monthly Temperature (°C)",
     main = "1880-2017 Monthly Temperature Trend")

Bonus: Use a Time Series Object for Easier Plotting #

For monthly data, converting your temperature column to a ts object simplifies plotting entirely—no need to worry about date columns at all:

# Create a time series starting at January 1880, 12 observations per year
temp_time_series <- ts(clean_data$temperature, start = c(1880, 1), frequency = 12)
# Plot directly
plot(temp_time_series, type = "l",
     xlab = "Year",
     ylab = "Monthly Temperature (°C)",
     main = "1880-2017 Monthly Temperature Trend")

Give these steps a try—you should be able to get that line plot working in no time!

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