Как подмножество временных рядов в 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 году.

Как я могу принести o ut соответствие между годом и измерениями, чтобы более поздние значения подмножества?

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