r - 如何在 R 中对时间序列进行子集化

Question

特别是，我想从 1960 中提取温度测量值的子集。在时间序列中向前 gtemp包裹内astsa :

require(astsa)
gtemp
Time Series:
Start = 1880 
End = 2009 
Frequency = 1 
  [1] -0.28 -0.21 -0.26 -0.27 -0.32 -0.32 -0.29 -0.36 -0.27 -0.17 -0.39 -0.27 -0.32
 [14] -0.33 -0.33 -0.25 -0.14 -0.11 -0.25 -0.15 -0.07 -0.14 -0.24 -0.30 -0.34 -0.24
 [27] -0.19 -0.39 -0.33 -0.35 -0.33 -0.34 -0.32 -0.30 -0.15 -0.10 -0.30 -0.39 -0.33
 [40] -0.20 -0.19 -0.14 -0.26 -0.22 -0.22 -0.17 -0.02 -0.15 -0.12 -0.26 -0.08 -0.02
 [53] -0.08 -0.19 -0.07 -0.12 -0.05  0.07  0.10  0.01  0.04  0.10  0.03  0.09  0.19
 [66]  0.06 -0.05  0.00 -0.04 -0.07 -0.16 -0.04  0.03  0.11 -0.10 -0.10 -0.17  0.08
 [79]  0.08  0.06 -0.01  0.07  0.04  0.08 -0.21 -0.11 -0.03 -0.01 -0.04  0.08  0.03
 [92] -0.10  0.00  0.14 -0.08 -0.05 -0.16  0.12  0.01  0.08  0.18  0.26  0.04  0.26
[105]  0.09  0.05  0.12  0.26  0.31  0.19  0.37  0.35  0.12  0.13  0.23  0.37  0.29
[118]  0.39  0.56  0.32  0.33  0.48  0.56  0.55  0.48  0.62  0.54  0.57  0.43  0.57

个别时间点没有标注年份，所以虽然我可以做 gtemp[3] [1] -0.26 ，我做不到 gtemp[as.date(1960)] ，例如获取 1960 年的值。

我怎样才能带出年份和测量值之间的对应关系，以便以后的子集值？

Answer 1

我们可以利用 window功能

gtemp1 <- window(gtemp, start = 1960)
gtemp1
#Time Series:
#Start = 1960 
#End = 2009 
#Frequency = 1 
#[1] -0.01  0.07  0.04  0.08 -0.21 -0.11 -0.03 -0.01 -0.04  0.08  0.03
#[12]-0.10  0.00  0.14 -0.08 -0.05 -0.16  0.12  0.01  0.08  0.18  0.26  
#[23] 0.04  0.26  0.09  0.05   0.12  0.26  0.31  0.19  0.37  0.35  0.12
#[34] 0.13  0.23  0.37  0.29  0.39  0.56  0.32  0.33  0.48  0.56  0.55
#[45] 0.48  0.62  0.54  0.57  0.43  0.57