[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