[R] Decomposing irregular time series
Jue Lin-Ye
jl.iccp at gmail.com
Wed May 6 15:04:41 CEST 2015
Dear fellow members of the r-help list,
May someone be so kind as to point out how an irregular time series
should be decomposed into season and trend?
My code, the related error and my data are enclosed herein.
Best regards and thanks in advance!
--
Jue
> x<-read.table(file="Data3.txt",header=T)
> for(i5 in 1:nrow(
x)){
+ time3<-ISOdatetime(
x[i5,"AAi"],X_season1[i5,"MMi"],x[i5,"DDi"],x[i5,"HHi"],0,0)
+ }
> zoo(
x[,"dur_storm"],time3)->ZOO1
> stl(ZOO1)
Error in stl(ZOO1) : series is not periodic or has less than two periods
In addition: Warning message:
In as.ts.zoo(x) : ‘x’ does not have an underlying regularity
-------------- next part --------------
"AAi" "MMi" "DDi" "HHi" "dur_storm"
2051 1 14 16 107
2051 1 22 0 167
2051 1 31 12 87
2051 2 16 17 1027
2051 8 17 18 67
2051 11 14 9 977
2051 11 21 15 357
2051 11 24 11 587
2051 12 8 21 107
2051 12 22 18 217
2052 1 20 22 197
2052 2 27 16 127
2052 3 4 21 217
2052 3 26 11 97
2052 7 25 18 97
2052 8 24 14 137
2052 10 14 7 127
2052 10 26 5 157
2052 11 20 23 527
2052 12 9 12 307
2052 12 28 17 107
2053 1 2 7 377
2053 2 6 6 667
2053 2 23 15 457
2053 3 9 5 347
2053 4 15 18 247
2053 4 30 22 127
2053 10 9 8 97
2053 10 31 23 437
2053 11 17 23 447
2053 12 31 2 177
2054 1 26 8 97
2054 1 31 12 147
2054 2 5 23 67
2054 2 13 14 137
2054 2 28 8 67
2054 4 14 21 177
2054 4 16 17 247
2054 10 13 21 207
2054 10 18 21 87
2054 11 1 14 337
2054 11 4 9 157
2054 11 12 20 97
2054 11 23 5 597
2054 12 3 1 147
2055 1 9 0 327
2055 1 27 10 287
2055 2 15 18 127
2055 3 3 17 127
2055 4 1 21 107
2055 4 16 3 117
2055 4 18 16 797
2055 5 16 10 67
2055 7 7 15 97
2055 9 10 13 307
2055 12 23 13 397
2055 12 29 0 127
2056 1 17 0 107
2056 1 20 10 107
2056 3 18 7 577
2056 3 24 7 457
2056 4 5 11 167
2056 4 18 20 477
2056 5 18 16 67
2056 11 25 10 237
2056 12 15 3 197
2056 12 27 3 167
2057 1 2 21 87
2057 1 15 8 1197
2057 2 3 4 107
2057 2 14 7 497
2057 2 23 23 307
2057 2 27 20 787
2057 3 6 17 67
2057 3 20 16 227
2057 5 4 14 147
2057 5 18 1 217
2057 10 21 12 567
2057 11 7 7 487
2057 12 10 1 107
2057 12 26 9 87
2058 1 10 1 817
2058 1 29 6 117
2058 1 31 8 137
2058 2 14 13 97
2058 2 15 18 77
2058 2 19 4 247
2058 3 21 14 67
2058 3 23 12 107
2058 4 7 2 187
2058 6 22 20 67
2058 9 29 5 67
2058 10 4 14 137
2058 11 7 6 567
2058 11 27 8 77
2058 12 29 15 67
2059 2 8 18 397
2059 2 16 2 287
2059 2 19 19 307
2059 4 15 0 447
2059 8 8 8 177
2059 10 21 15 77
2059 10 28 12 487
2059 11 7 10 167
2059 11 14 1 297
2059 11 22 23 67
2059 12 20 3 857
2060 1 17 4 377
2060 2 10 9 187
2060 2 21 8 507
2060 3 1 8 347
2060 3 14 8 327
2060 4 7 5 297
2060 4 17 20 507
2060 4 22 3 127
2060 10 1 16 217
2060 11 22 21 137
2060 12 17 8 117
2060 12 26 13 157
2060 12 31 16 87
2061 1 3 21 577
2061 1 21 10 187
2061 4 8 21 237
2061 4 13 21 327
2061 4 16 19 617
2061 5 24 17 397
2061 11 13 8 357
2061 11 17 8 117
2062 4 6 20 167
2062 4 11 10 377
2062 7 16 4 87
2062 8 29 18 447
2062 10 16 19 107
2062 11 1 18 127
2062 12 20 13 667
2062 12 27 8 267
2063 2 8 23 357
2063 3 13 17 257
2063 3 21 18 227
2063 3 26 16 167
2063 4 2 1 247
2063 4 9 16 277
2063 10 30 17 97
2063 11 11 7 87
2064 1 20 19 77
2064 2 8 14 337
2064 2 21 5 157
2064 3 16 21 227
2064 8 28 5 147
2064 11 3 2 137
2064 12 17 7 167
2065 1 5 19 97
2065 1 21 17 297
2065 1 31 8 397
2065 2 13 20 127
2065 2 24 7 337
2065 3 15 16 197
2065 3 27 21 207
2065 4 9 6 217
2065 4 26 13 97
2065 5 12 17 97
2065 9 7 15 87
2065 11 4 14 237
2065 11 18 1 157
2065 11 19 12 487
2065 12 18 10 77
2066 1 3 17 217
2066 2 22 5 297
2066 3 1 8 407
2066 3 22 17 77
2066 4 3 4 147
2066 5 4 23 397
2066 6 5 10 87
2066 10 9 14 167
2066 10 16 17 67
2066 11 3 18 807
2066 12 16 15 77
2067 1 17 2 107
2067 1 22 2 117
2067 1 28 12 147
2067 2 13 2 67
2067 2 16 2 297
2067 2 25 9 187
2067 3 27 11 97
2067 4 16 14 97
2067 4 26 5 217
2067 4 28 11 237
2067 5 5 7 197
2067 10 2 22 357
2067 10 18 13 327
2067 11 23 17 197
2067 11 30 10 317
2067 12 8 20 287
2067 12 12 10 97
2067 12 14 7 287
2068 1 4 22 217
2068 2 3 14 237
2068 2 13 15 287
2068 2 27 7 117
2068 3 7 4 397
2068 10 16 17 97
2068 10 18 14 207
2068 10 20 19 257
2068 10 29 4 207
2068 11 28 2 167
2068 11 30 5 227
2068 12 10 4 117
2069 1 4 8 257
2069 1 20 20 87
2069 4 6 18 107
2069 4 12 21 207
2069 5 2 8 217
2069 5 30 15 127
2069 6 30 6 77
2069 9 17 3 127
2069 10 4 17 67
2069 11 7 8 197
2069 11 12 5 477
2069 11 26 4 127
2069 11 30 7 257
2070 1 14 18 117
2070 1 29 11 717
2070 3 11 1 317
2070 3 21 2 67
2070 3 31 23 477
2070 4 13 11 297
2070 10 24 2 477
2070 12 15 17 67
2071 1 20 1 417
2071 1 22 12 117
2071 2 20 1 547
2071 2 25 16 357
2071 3 14 14 217
2071 3 24 10 137
2071 3 28 19 377
2071 4 2 3 167
2071 6 6 15 127
2071 10 2 0 127
2072 1 2 8 267
2072 2 12 0 347
2072 2 27 14 137
2072 5 7 9 117
2072 8 5 20 77
2072 10 9 7 97
2072 10 11 13 137
2072 10 17 4 397
2072 10 22 14 127
2072 10 28 1 227
2072 11 21 13 107
2072 11 26 3 147
2072 12 14 23 197
2072 12 19 3 217
2072 12 31 20 87
2073 2 13 15 117
2073 2 18 1 237
2073 3 2 0 277
2073 3 4 7 117
2073 3 5 14 97
2073 3 13 22 87
2073 3 19 8 127
2073 3 20 10 827
2073 3 28 16 77
2073 3 30 5 267
2073 10 7 20 87
2073 10 10 7 177
2073 10 13 3 167
2073 10 20 19 167
2073 12 5 8 97
2073 12 8 16 107
2074 1 8 4 67
2074 1 14 1 207
2074 1 24 21 227
2074 2 2 21 157
2074 3 23 8 317
2074 3 30 18 347
2074 8 22 20 197
2074 9 30 19 207
2074 11 10 0 617
2074 12 2 20 77
2074 12 17 7 197
2074 12 21 11 137
2074 12 24 9 267
2075 1 4 5 137
2075 1 10 15 207
2075 1 20 10 127
2075 1 23 21 97
2075 1 29 18 457
2075 4 12 20 447
2075 5 3 4 187
2075 5 8 1 157
2075 5 23 14 117
2075 9 30 19 127
2075 10 21 18 187
2075 11 5 10 117
2075 11 21 16 97
2075 11 27 16 377
2075 12 28 5 67
2076 1 21 9 347
2076 1 25 6 667
2076 3 5 3 237
2076 3 7 17 337
2076 3 25 1 167
2076 4 2 4 127
2076 9 17 20 67
2076 11 1 16 427
2076 12 31 18 187
2077 1 16 14 137
2077 2 1 20 187
2077 2 25 6 497
2077 3 22 15 267
2077 3 30 1 537
2077 4 24 10 197
2077 4 30 15 167
2077 6 6 16 197
2077 9 2 12 147
2077 9 26 14 367
2077 10 21 14 207
2077 12 9 6 317
2078 1 3 2 357
2078 1 5 12 187
2078 2 10 16 107
2078 2 15 10 307
2078 2 23 4 177
2078 3 22 8 177
2078 4 4 11 97
2078 9 7 18 207
2078 9 14 18 137
2078 11 6 18 257
2078 11 14 7 197
2078 12 20 15 857
2079 1 12 17 177
2079 1 29 3 207
2079 3 15 14 127
2079 3 19 4 247
2079 4 13 14 317
2079 6 17 3 97
2079 10 15 22 77
2079 10 19 20 137
2079 11 25 11 267
2079 12 20 9 97
2080 1 10 14 177
2080 1 20 4 147
2080 2 29 16 177
2080 3 27 16 97
2080 4 5 14 117
2080 4 28 12 67
2080 10 2 9 147
2080 10 27 7 117
2080 10 29 16 257
2080 12 23 13 87
2081 1 7 11 137
2081 3 3 11 197
2081 3 20 5 157
2081 3 25 11 167
2081 3 26 17 607
2081 4 2 0 627
2081 4 29 6 77
2081 9 2 17 377
2081 9 24 12 157
2081 10 4 9 357
2081 10 6 13 117
2081 10 8 17 397
2081 10 12 19 217
2081 10 15 7 77
2081 11 11 13 137
2081 11 13 19 307
2081 12 2 2 137
2081 12 30 10 487
2082 2 13 2 97
2082 2 26 17 67
2082 5 4 17 157
2082 5 12 19 67
2082 10 19 17 67
2082 10 27 16 177
2082 11 7 23 137
2082 12 16 8 67
2082 12 22 14 67
2082 12 23 22 707
2083 1 16 14 457
2083 1 28 12 87
2083 2 16 14 87
2083 3 2 5 477
2083 3 8 22 257
2083 3 11 17 87
2083 4 10 23 227
2083 5 12 2 77
2083 9 18 2 167
2083 10 22 9 317
2083 11 10 4 87
2083 11 18 21 327
2083 12 7 15 267
2084 1 16 11 757
2084 1 21 23 77
2084 1 26 9 127
2084 2 5 8 197
2084 3 10 17 77
2084 3 27 17 207
2084 4 1 17 77
2084 4 15 13 137
2084 5 6 11 257
2084 5 16 2 237
2084 10 20 22 87
2084 10 26 19 407
2084 11 14 9 77
2084 11 23 14 127
2084 12 3 6 137
2084 12 16 19 77
2084 12 31 4 297
2085 1 19 3 337
2085 2 21 22 317
2085 3 14 18 317
2085 4 7 18 97
2085 8 21 10 177
2085 11 23 18 417
2085 12 29 17 267
2086 1 22 2 427
2086 5 8 12 127
2086 10 12 2 417
2086 11 3 7 117
2086 11 8 8 387
2086 12 4 14 767
2086 12 19 20 87
2087 1 8 19 107
2087 1 13 14 657
2087 2 12 1 227
2087 2 14 14 327
2087 3 2 13 407
2087 3 8 17 327
2087 3 14 1 147
2087 3 17 14 207
2087 3 22 3 107
2087 3 24 23 237
2087 4 26 15 317
2087 6 25 16 327
2087 8 29 16 87
2087 10 30 23 267
2087 11 14 4 267
2088 1 6 3 117
2088 2 23 22 567
2088 3 5 5 237
2088 5 7 9 147
2088 5 12 8 77
2088 10 15 15 447
2088 10 19 9 87
2088 12 13 9 67
2088 12 19 2 287
2089 1 2 1 207
2089 1 6 20 127
2089 1 13 13 597
2089 1 19 14 617
2089 1 28 13 447
2089 2 17 23 237
2089 2 25 3 107
2089 3 7 4 67
2089 3 14 8 207
2089 4 9 4 427
2089 9 23 17 67
2089 9 26 8 147
2089 10 1 7 107
2089 10 31 3 367
2089 11 28 0 147
2089 11 30 10 507
2090 1 5 21 217
2090 1 13 22 77
2090 2 15 23 127
2090 2 21 6 97
2090 3 9 19 217
2090 3 11 14 67
2090 3 16 8 127
2090 3 22 12 467
2090 3 28 12 137
2090 3 30 12 147
2090 4 19 22 147
2090 12 3 9 387
2090 12 15 14 87
2090 12 25 1 377
2090 12 28 0 827
2091 1 7 22 467
2091 1 12 2 447
2091 1 25 12 667
2091 2 20 19 147
2091 3 17 20 197
2091 3 22 17 167
2091 4 6 19 197
2091 4 18 19 427
2091 6 14 21 97
2091 6 29 11 177
2091 8 26 7 87
2091 9 3 19 107
2091 9 17 17 377
2091 10 14 6 237
2091 12 10 5 227
2091 12 16 0 247
2092 1 21 22 467
2092 1 29 0 177
2092 2 11 15 177
2092 2 27 13 607
2092 5 5 17 87
2092 8 15 5 207
2092 9 5 11 187
2092 9 19 8 107
2092 9 22 14 407
2092 10 8 18 137
2092 10 10 10 177
2092 10 12 14 107
2092 10 31 7 107
2092 11 7 18 67
2092 11 29 1 187
2092 12 3 5 157
2092 12 10 13 87
2092 12 11 17 307
2093 1 14 13 67
2093 3 4 0 277
2093 3 14 7 697
2093 5 20 13 77
2093 8 19 18 107
2093 10 6 18 187
2093 12 7 16 417
2093 12 14 11 97
2094 1 14 19 257
2094 1 17 17 447
2094 2 27 21 267
2094 3 5 10 167
2094 3 7 22 97
2094 3 10 8 177
2094 3 12 19 117
2094 4 21 1 207
2094 5 13 7 197
2094 9 14 16 107
2094 10 5 11 87
2094 11 17 9 457
2094 12 11 13 487
2094 12 26 3 87
2094 12 27 19 167
2095 2 27 20 297
2095 3 8 11 197
2095 3 13 3 97
2095 4 2 8 287
2095 4 18 11 77
2095 10 11 8 547
2095 10 19 0 247
2095 11 19 9 167
2095 11 20 21 137
2096 1 31 11 247
2096 2 24 14 87
2096 3 6 20 157
2096 3 27 16 127
2096 4 30 0 137
2096 10 5 20 227
2096 10 15 4 1227
2096 11 3 20 457
2096 11 12 7 117
2096 12 6 2 137
2096 12 28 14 77
2097 1 15 4 317
2097 1 25 18 157
2097 2 21 18 117
2097 4 7 11 907
2097 4 14 13 147
2097 4 26 13 187
2097 9 26 7 77
2097 10 20 17 617
2097 10 25 21 207
2097 12 2 20 197
2097 12 14 8 147
2098 1 13 9 107
2098 1 16 7 327
2098 1 19 8 417
2098 2 12 3 107
2098 3 4 9 1157
2098 3 16 13 87
2098 4 15 4 587
2098 5 10 19 67
2098 5 11 16 67
2098 9 9 18 347
2098 9 29 5 347
2098 11 16 16 87
2098 11 20 12 687
2098 12 13 16 137
2099 1 16 16 157
2099 2 28 12 287
2099 3 11 8 267
2099 3 15 17 67
2099 3 29 1 497
2099 5 17 10 107
2099 5 23 2 507
2099 8 1 11 157
2099 8 30 0 107
2099 9 10 9 187
2099 11 12 4 247
2099 11 23 23 277
2099 12 9 12 277
2099 12 22 15 167
2099 12 26 17 747
2100 3 12 14 107
2100 4 4 10 397
2100 4 18 13 1067
2100 10 20 6 157
2100 12 30 3 297
More information about the R-help
mailing list