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
"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
______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.