r - 使用 R 更改时间序列图中的 X 轴值

Question

我是 R 的新用户，我需要有关时间序列图的帮助。我创建了一个时间序列图，但无法弄清楚如何更改我的 x 轴值以对应于我的样本日期。我的数据如下:

 Year  Month  Level
 2009   8     350
 2009   9     210
 2009   10    173
 2009   11    166
 2009   12    153
 2010   1     141
 2010   2     129
 2010   3     124
 2010   4     103
 2010   5     69
 2010   6     51
 2010   7     49
 2010   8     51
 2010   9     51

假设此数据保存为名称“data.csv”

data = read.table("data.csv", sep = ",", header = T)
data.ts = ts(data, frequency = 1)
plot(dat.mission.ts[, 3], ylab = "level", main = "main", axes = T)

我也试过将 start = c(2009, 8) 输入到 ts 函数中，但我仍然得到错误的值

当我绘制此图时，我的 x 轴与 2009 年 8 月至 2010 年 9 月不相关。它将按年增加或仅按小数增加。我在网上和通过 ?帮助 R，但找不到重新标记我的轴值的方法。任何帮助将不胜感激。

Answer 1

使用 base 编码，您可以通过几个步骤完成此操作。如本 SO answer 中所述，如果您同时使用 as.Date 和 paste 函数并合并一天(即第一天月份；“1”)。出于此答案的目的，我将简单地将您提供的数据称为 df:

df$date<-with(df,as.Date(paste(Year,Month,'1',sep='-'),format='%Y-%m-%d'))
df
   Year Month Level       date
1  2009     8   350 2009-08-01
2  2009     9   210 2009-09-01
3  2009    10   173 2009-10-01
4  2009    11   166 2009-11-01
5  2009    12   153 2009-12-01
6  2010     1   141 2010-01-01
7  2010     2   129 2010-02-01
8  2010     3   124 2010-03-01
9  2010     4   103 2010-04-01
10 2010     5    69 2010-05-01
11 2010     6    51 2010-06-01
12 2010     7    49 2010-07-01
13 2010     8    51 2010-08-01
14 2010     9    51 2010-09-01

然后，您可以使用基本的 plot、axis 和 mtext 函数来控制数据和坐标轴的可视化方式。例如:

xmin<-min(df$date,na.rm=T);xmax<-max(df$date,na.rm=T) #ESTABLISH X-VALUES (MIN & MAX)
ymin<-min(df$Level,na.rm=T);ymax<-max(df$Level,na.rm=T) #ESTABLISH Y-VALUES (MIN & MAX)

xseq<-seq.Date(xmin,xmax,by='1 month') #CREATE DATE SEQUENCE THAT INCREASES BY MONTH FROM DATE MINIMUM TO MAXIMUM
yseq<-round(seq(0,ymax,by=50),0) # CREATE SEQUENCE FROM 0-350 BY 50

par(mar=c(1,1,0,0),oma=c(6,5,3,2)) #CONTROLS YOUR IMAGE MARGINS
plot(Level~date,data=df,type='b',ylim=c(0,ymax),axes=F,xlab='',ylab='');box() #PLOT LEVEL AS A FUNCTION OF DATE, REMOVE AXES FOR FUTURE CUSTOMIZATION 
axis.Date(side=1,at=xseq,format='%Y-%m',labels=T,las=3) #ADD X-AXIS LABELS WITH "YEAR-MONTH" FORMAT
axis(side=2,at=yseq,las=2) #ADD Y-AXIS LABELS
mtext('Date (Year-Month)',side=1,line=5) #X-AXIS LABEL
mtext('Level',side=2,line=4) #Y-AXIS LABEL