如何在Stata中转换日期并绘制日度数据年度时间序列图？

一、解决substr()类型不匹配错误 #

• 错误原因：substr()是字符串处理函数，你的date是数字型变量，直接传入会触发类型不匹配。
• 两种解决方法：
  1. 先转字符串再提取年份：

     gen year_str = substr(string(date, "%8.0f"), 1, 4)
     destring year_str, replace  // 转成数值型年份
     

     注：%8.0f确保数字转字符串时保留8位，避免短位数日期的格式问题。
  2. 更推荐方案：转成标准日期格式后提取年份（更准确，后续绘图也能直接复用）：

     gen date_real = date(string(date, "%8.0f"), "YMD")  // 转成Stata可识别的日期格式
     format date_real %td  // 设置日期显示格式
     gen year = year(date_real)  // 一键提取年份
     

二、仅用年份绘制日度数据的合理性分析 #

完全不合理。日度数据每年约有300个观测值，若将X轴设为年份，同一年的所有数据点会重叠在同一X位置，导致图形无法展示日度波动，最终呈现杂乱无章的效果（正如你当前的图）。正确思路是基于日期格式变量绘图，让工具自动按年份标记刻度，兼顾日度细节和年份维度的可读性。

三、便捷绘图解决方法 #

以下是Stata和R语言的最优绘图方案：

Stata 实现 #

// 1. 先将数字型日期转为标准日期变量
gen date_real = date(string(date, "%8.0f"), "YMD")
format date_real %td

// 2. 绘制优化后的时间序列图
line spread date_real, ///
xlabel(#10, format(%tdCY))  // 自动生成约10个刻度，显示为年份（如2011、2012）
title("期限利差日度时间序列") ///
ytitle("期限利差(%)") ///
xtitle("年份") ///
legend(off) ///
graphregion(color(white))  // 优化背景整洁度

R（ggplot2）实现 #

library(ggplot2)

// 1. 转换日期格式
df$date_real <- as.Date(as.character(df$date), format = "%Y%m%d")

// 2. 绘制优化后的时间序列图
ggplot(df, aes(x = date_real, y = spread)) +
  geom_line(color = "#2c3e50") +
  scale_x_date(
    date_labels = "%Y",  // 刻度显示为年份
    date_breaks = "1 year"  // 每年设置一个刻度
  ) +
  labs(
    title = "期限利差日度时间序列",
    x = "年份",
    y = "期限利差(%)"
  ) +
  theme_minimal()

这样绘制的图既保留了日度数据的波动细节，X轴又以年份为刻度标记，清晰易读，不会出现数据重叠问题。

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