R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(235.1 + ,280.7 + ,264.6 + ,240.7 + ,201.4 + ,240.8 + ,241.1 + ,223.8 + ,206.1 + ,174.7 + ,203.3 + ,220.5 + ,299.5 + ,347.4 + ,338.3 + ,327.7 + ,351.6 + ,396.6 + ,438.8 + ,395.6 + ,363.5 + ,378.8 + ,357 + ,369 + ,464.8 + ,479.1 + ,431.3 + ,366.5 + ,326.3 + ,355.1 + ,331.6 + ,261.3 + ,249 + ,205.5 + ,235.6 + ,240.9 + ,264.9 + ,253.8 + ,232.3 + ,193.8 + ,177 + ,213.2 + ,207.2 + ,180.6 + ,188.6 + ,175.4 + ,199 + ,179.6 + ,225.8 + ,234 + ,200.2 + ,183.6 + ,178.2 + ,203.2 + ,208.5 + ,191.8 + ,172.8 + ,148 + ,159.4 + ,154.5 + ,213.2 + ,196.4 + ,182.8 + ,176.4 + ,153.6 + ,173.2 + ,171 + ,151.2 + ,161.9 + ,157.2 + ,201.7 + ,236.4 + ,356.1 + ,398.3 + ,403.7 + ,384.6 + ,365.8 + ,368.1 + ,367.9 + ,347 + ,343.3 + ,292.9 + ,311.5 + ,300.9 + ,366.9 + ,356.9 + ,329.7 + ,316.2 + ,269 + ,289.3 + ,266.2 + ,253.6 + ,233.8 + ,228.4 + ,253.6 + ,260.1 + ,306.6 + ,309.2 + ,309.5 + ,271 + ,279.9 + ,317.9 + ,298.4 + ,246.7 + ,227.3 + ,209.1 + ,259.9 + ,266 + ,320.6 + ,308.5 + ,282.2 + ,262.7 + ,263.5 + ,313.1 + ,284.3 + ,252.6 + ,250.3 + ,246.5 + ,312.7 + ,333.2 + ,446.4 + ,511.6 + ,515.5 + ,506.4 + ,483.2 + ,522.3 + ,509.8 + ,460.7 + ,405.8 + ,375 + ,378.5 + ,406.8 + ,467.8 + ,469.8 + ,429.8 + ,355.8 + ,332.7 + ,378 + ,360.5 + ,334.7 + ,319.5 + ,323.1 + ,363.6 + ,352.1 + ,411.9 + ,388.6 + ,416.4 + ,360.7 + ,338 + ,417.2 + ,388.4 + ,371.1 + ,331.5 + ,353.7 + ,396.7 + ,447 + ,533.5 + ,565.4 + ,542.3 + ,488.7 + ,467.1 + ,531.3 + ,496.1 + ,444 + ,403.4 + ,386.3 + ,394.1 + ,404.1 + ,462.1 + ,448.1 + ,432.3 + ,386.3 + ,395.2 + ,421.9 + ,382.9 + ,384.2 + ,345.5 + ,323.4 + ,372.6 + ,376 + ,462.7 + ,487 + ,444.2 + ,399.3 + ,394.9 + ,455.4 + ,414 + ,375.5 + ,347 + ,339.4 + ,385.8 + ,378.8 + ,451.8 + ,446.1 + ,422.5 + ,383.1 + ,352.8 + ,445.3 + ,367.5 + ,355.1 + ,326.2 + ,319.8 + ,331.8 + ,340.9 + ,394.1 + ,417.2 + ,369.9 + ,349.2 + ,321.4 + ,405.7 + ,342.9 + ,316.5 + ,284.2 + ,270.9 + ,288.8 + ,278.8 + ,324.4 + ,310.9 + ,299 + ,273 + ,279.3 + ,359.2 + ,305 + ,282.1 + ,250.3 + ,246.5 + ,257.9 + ,266.5 + ,315.9 + ,318.4 + ,295.4 + ,266.4 + ,245.8 + ,362.8 + ,324.9 + ,294.2 + ,289.5 + ,295.2 + ,290.3 + ,272 + ,307.4 + ,328.7 + ,292.9 + ,249.1 + ,230.4 + ,361.5 + ,321.7 + ,277.2 + ,260.7 + ,251 + ,257.6 + ,241.8 + ,287.5 + ,292.3 + ,274.7 + ,254.2 + ,230 + ,339 + ,318.2 + ,287 + ,295.8 + ,284 + ,271 + ,262.7 + ,340.6 + ,379.4 + ,373.3 + ,355.2 + ,338.4 + ,466.9 + ,451 + ,422 + ,429.2 + ,425.9 + ,460.7 + ,463.6 + ,541.4 + ,544.2 + ,517.5 + ,469.4 + ,439.4 + ,549 + ,533 + ,506.1 + ,484 + ,457 + ,481.5 + ,469.5 + ,544.7 + ,541.2 + ,521.5 + ,469.7 + ,434.4 + ,542.6 + ,517.3 + ,485.7 + ,465.8 + ,447 + ,426.6 + ,411.6 + ,467.5 + ,484.5 + ,451.2 + ,417.4 + ,379.9 + ,484.7 + ,455 + ,420.8 + ,416.5 + ,376.3 + ,405.6 + ,405.8 + ,500.8 + ,514 + ,475.5 + ,430.1 + ,414.4 + ,538 + ,526 + ,488.5 + ,520.2 + ,504.4 + ,568.5 + ,610.6 + ,818 + ,830.9 + ,835.9 + ,782 + ,762.3 + ,856.9 + ,820.9 + ,769.6 + ,752.2 + ,724.4 + ,723.1 + ,719.5 + ,817.4 + ,803.3 + ,752.5 + ,689 + ,630.4 + ,765.5 + ,757.7 + ,732.2 + ,702.6 + ,683.3 + ,709.5 + ,702.2 + ,784.8 + ,810.9 + ,755.6 + ,656.8 + ,615.1 + ,745.3 + ,694.1 + ,675.7 + ,643.7 + ,622.1 + ,634.6 + ,588 + ,689.7 + ,673.9 + ,647.9 + ,568.8 + ,545.7 + ,632.6 + ,643.8 + ,593.1 + ,579.7 + ,546 + ,562.9 + ,572.5) > par8 = 'FALSE' > par7 = '1' > par6 = '' > par5 = '1' > par4 = '' > par3 = '0' > par2 = 'periodic' > par1 = '6' > main = 'Seasonal Decomposition by Loess' > par1 <- as.numeric(par1) #seasonal period > if (par2 != 'periodic') par2 <- as.numeric(par2) #s.window > par3 <- as.numeric(par3) #s.degree > if (par4 == '') par4 <- NULL else par4 <- as.numeric(par4)#t.window > par5 <- as.numeric(par5)#t.degree > if (par6 != '') par6 <- as.numeric(par6)#l.window > par7 <- as.numeric(par7)#l.degree > if (par8 == 'FALSE') par8 <- FALSE else par9 <- TRUE #robust > nx <- length(x) > x <- ts(x,frequency=par1) > if (par6 != '') { + m <- stl(x,s.window=par2, s.degree=par3, t.window=par4, t.degre=par5, l.window=par6, l.degree=par7, robust=par8) + } else { + m <- stl(x,s.window=par2, s.degree=par3, t.window=par4, t.degre=par5, l.degree=par7, robust=par8) + } > m$time.series Time Series: Start = c(1, 1) End = c(62, 6) Frequency = 6 seasonal trend remainder 1.000000 31.1795888 247.3892 -43.46875347 1.166667 19.7193020 244.9113 16.06937110 1.333333 -0.8038918 242.4335 22.97040273 1.500000 -28.8746459 239.3598 30.21484493 1.666667 -28.9970363 236.2861 -5.88907641 1.833333 7.7766795 228.7961 4.22725520 2.000000 31.1795888 221.3060 -11.38560675 2.166667 19.7193020 217.7802 -13.69950822 2.333333 -0.8038918 214.2544 -7.35050264 2.500000 -28.8746459 223.5996 -20.02498621 2.666667 -28.9970363 232.9449 -0.64783333 2.833333 7.7766795 254.9011 -42.17776072 3.000000 31.1795888 276.8573 -8.53688167 3.166667 19.7193020 302.8406 24.84011779 3.333333 -0.8038918 328.8239 10.28002430 3.500000 -28.8746459 348.6641 7.91050637 3.666667 -28.9970363 368.5044 12.09262490 3.833333 7.7766795 375.7917 13.03165797 4.000000 31.1795888 383.0789 24.54149748 4.166667 19.7193020 384.6810 -8.80025546 4.333333 -0.8038918 386.2830 -21.97910136 4.500000 -28.8746459 391.0289 16.64573317 4.666667 -28.9970363 395.7748 -9.77779584 4.833333 7.7766795 403.1307 -41.90738062 5.000000 31.1795888 410.4866 23.13384103 5.166667 19.7193020 406.9015 52.47923403 5.333333 -0.8038918 403.3164 28.78753409 5.500000 -28.8746459 382.8933 12.48131621 5.666667 -28.9970363 362.4703 -7.17326521 5.833333 7.7766795 332.9061 14.41718843 6.000000 31.1795888 303.3420 -2.92155149 6.166667 19.7193020 283.4739 -41.89319342 6.333333 -0.8038918 263.6058 -13.80192831 6.500000 -28.8746459 253.4797 -19.10501417 6.666667 -28.9970363 243.3535 21.24353642 6.833333 7.7766795 240.1110 -6.98763149 7.000000 31.1795888 236.8684 -3.14799296 7.166667 19.7193020 230.9046 3.17614064 7.333333 -0.8038918 224.9407 8.16318129 7.500000 -28.8746459 216.0262 6.64840597 7.666667 -28.9970363 207.1118 -1.11473290 7.833333 7.7766795 200.0527 5.37061126 8.000000 31.1795888 192.9936 -16.97323815 8.166667 19.7193020 191.9872 -31.10646769 8.333333 -0.8038918 190.9807 -1.57679019 8.500000 -28.8746459 193.4971 10.77751892 8.666667 -28.9970363 196.0136 31.98346447 8.833333 7.7766795 199.0876 -27.26429622 9.000000 31.1795888 202.1617 -7.54125048 9.166667 19.7193020 202.2423 12.03836684 9.333333 -0.8038918 202.3230 -1.31910880 9.500000 -28.8746459 200.0271 12.44750316 9.666667 -28.9970363 197.7313 9.46575157 9.833333 7.7766795 192.2570 3.16629480 10.000000 31.1795888 186.7828 -9.46235554 10.166667 19.7193020 181.8628 -9.78206935 10.333333 -0.8038918 176.9428 -3.33887610 10.500000 -28.8746459 175.3588 1.51581098 10.666667 -28.9970363 173.7749 14.62213451 10.833333 7.7766795 175.7100 -28.98672369 11.000000 31.1795888 177.6452 4.37522454 11.166667 19.7193020 179.0013 -2.32063179 11.333333 -0.8038918 180.3575 3.24641893 11.500000 -28.8746459 176.4034 28.87127614 11.666667 -28.9970363 172.4493 10.14776981 11.833333 7.7766795 168.2427 -2.81934807 12.000000 31.1795888 164.0361 -24.21565950 12.166667 19.7193020 170.5881 -39.10735984 12.333333 -0.8038918 177.1400 -14.43615313 12.500000 -28.8746459 204.6972 -18.62258311 12.666667 -28.9970363 232.2544 -1.55737664 12.833333 7.7766795 270.9351 -42.31178212 13.000000 31.1795888 309.6158 15.30461884 13.166667 19.7193020 337.7210 40.85965033 13.333333 -0.8038918 365.8263 38.67758887 13.500000 -28.8746459 370.8581 42.61651977 13.666667 -28.9970363 375.8899 18.90708711 13.833333 7.7766795 365.4270 -5.10370104 14.000000 31.1795888 354.9641 -18.24368275 14.166667 19.7193020 344.4320 -17.15134792 14.333333 -0.8038918 333.9000 10.20389397 14.500000 -28.8746459 330.7161 -8.94144688 14.666667 -28.9970363 327.5322 12.96484872 14.833333 7.7766795 327.6409 -34.51755072 15.000000 31.1795888 327.7496 7.97085628 15.166667 19.7193020 324.5325 12.64821886 15.333333 -0.8038918 321.3154 9.18848849 15.500000 -28.8746459 308.7831 36.29158047 15.666667 -28.9970363 296.2507 1.74630891 15.833333 7.7766795 280.8350 0.68834771 16.000000 31.1795888 265.4192 -30.39880706 16.166667 19.7193020 259.0276 -25.14693147 16.333333 -0.8038918 252.6360 -18.03214884 16.500000 -28.8746459 257.0803 0.19434101 16.666667 -28.9970363 261.5246 21.07246731 16.833333 7.7766795 270.8786 -18.55523528 17.000000 31.1795888 280.2325 -4.81213143 17.166667 19.7193020 286.9311 2.54962496 17.333333 -0.8038918 293.6296 16.67428840 17.500000 -28.8746459 292.3887 7.48592359 17.666667 -28.9970363 291.1478 17.74919523 17.833333 7.7766795 280.3820 29.74135764 18.000000 31.1795888 269.6161 -2.39567351 18.166667 19.7193020 262.7684 -35.78769490 18.333333 -0.8038918 255.9207 -27.81680925 18.500000 -28.8746459 258.4010 -20.42635190 18.666667 -28.9970363 260.8813 28.01574190 18.833333 7.7766795 269.3385 -11.11519279 19.000000 31.1795888 277.7957 11.62467895 19.166667 19.7193020 282.3932 6.38747583 19.333333 -0.8038918 286.9907 -3.98682024 19.500000 -28.8746459 284.2036 7.37105300 19.666667 -28.9970363 281.4165 11.08056270 19.833333 7.7766795 276.3258 28.99748629 20.000000 31.1795888 271.2352 -18.11478369 20.166667 19.7193020 275.6523 -42.77157104 20.333333 -0.8038918 280.0693 -28.96545133 20.500000 -28.8746459 305.2951 -29.92046243 20.666667 -28.9970363 330.5209 11.17616292 20.833333 7.7766795 372.2173 -46.79397380 21.000000 31.1795888 413.9137 1.30669592 21.166667 19.7193020 447.0481 44.83261205 21.333333 -0.8038918 480.1825 36.12143523 21.500000 -28.8746459 490.4011 44.87357185 21.666667 -28.9970363 500.6197 11.57734493 21.833333 7.7766795 485.6622 28.86114678 22.000000 31.1795888 470.7047 7.91575507 22.166667 19.7193020 452.0248 -11.04406433 22.333333 -0.8038918 433.3449 -26.74097667 22.500000 -28.8746459 425.5127 -21.63803423 22.666667 -28.9970363 417.6805 -10.18345534 22.833333 7.7766795 418.0477 -19.02436789 23.000000 31.1795888 418.4149 18.20552600 23.166667 19.7193020 412.7480 37.33273986 23.333333 -0.8038918 407.0810 23.52286076 23.500000 -28.8746459 391.7038 -7.02916912 23.666667 -28.9970363 376.3266 -14.62956256 23.833333 7.7766795 361.6168 8.60652849 24.000000 31.1795888 346.9070 -17.58657402 24.166667 19.7193020 345.6859 -30.70515604 24.333333 -0.8038918 344.4647 -24.16083102 24.500000 -28.8746459 350.4674 1.50725730 24.666667 -28.9970363 356.4701 36.12698208 24.833333 7.7766795 366.3811 -22.05779501 25.000000 31.1795888 376.2922 4.42823435 25.166667 19.7193020 379.9969 -11.11618754 25.333333 -0.8038918 383.7016 33.50229762 25.500000 -28.8746459 382.9120 6.66264012 25.666667 -28.9970363 382.1224 -15.12538092 25.833333 7.7766795 376.1999 33.22340827 26.000000 31.1795888 370.2774 -13.05699610 26.166667 19.7193020 375.2332 -23.85253769 26.333333 -0.8038918 380.1891 -47.88517223 26.500000 -28.8746459 401.6517 -19.07707524 26.666667 -28.9970363 423.1144 2.58265820 26.833333 7.7766795 452.3962 -13.17286657 27.000000 31.1795888 481.6780 20.64241509 27.166667 19.7193020 497.1351 48.54561036 27.333333 -0.8038918 512.5922 30.51171268 27.500000 -28.8746459 507.8041 9.77059276 27.666667 -28.9970363 503.0159 -6.91889070 27.833333 7.7766795 483.8765 39.64678209 28.000000 31.1795888 464.7371 0.18326132 28.166667 19.7193020 448.5564 -24.27565990 28.333333 -0.8038918 432.3756 -28.17167408 28.500000 -28.8746459 424.8161 -9.64148983 28.666667 -28.9970363 417.2567 5.84033087 28.833333 7.7766795 418.9113 -22.58793267 29.000000 31.1795888 420.5658 10.35461022 29.166667 19.7193020 420.9009 7.47980592 29.333333 -0.8038918 421.2360 11.86790867 29.500000 -28.8746459 413.5778 1.59688452 29.666667 -28.9970363 405.9195 18.27749682 29.833333 7.7766795 393.8895 20.23383564 30.000000 31.1795888 381.8594 -30.13901911 30.166667 19.7193020 375.7905 -11.30975851 30.333333 -0.8038918 369.7215 -23.41759086 30.500000 -28.8746459 378.0135 -25.73889608 30.666667 -28.9970363 386.3056 15.29143516 30.833333 7.7766795 401.0570 -32.83363875 31.000000 31.1795888 415.8083 15.71209377 31.166667 19.7193020 424.2686 43.01210584 31.333333 -0.8038918 432.7289 12.27502496 31.500000 -28.8746459 427.5430 0.63166494 31.666667 -28.9970363 422.3571 1.53994138 31.833333 7.7766795 408.4297 39.19359661 32.000000 31.1795888 394.5024 -11.68194172 32.166667 19.7193020 387.3933 -31.61262622 32.333333 -0.8038918 380.2843 -32.48040368 32.500000 -28.8746459 383.3843 -15.10968551 32.666667 -28.9970363 386.4844 28.31266912 32.833333 7.7766795 395.9743 -24.95101642 33.000000 31.1795888 405.4643 15.15610448 33.166667 19.7193020 408.0154 18.36525926 33.333333 -0.8038918 410.5666 12.73732108 33.500000 -28.8746459 402.9074 9.06723706 33.666667 -28.9970363 395.2482 -13.45121050 33.833333 7.7766795 381.5818 55.94151092 34.000000 31.1795888 367.9154 -31.59496122 34.166667 19.7193020 358.8122 -23.43151576 34.333333 -0.8038918 349.7091 -22.70516326 34.500000 -28.8746459 349.9770 -1.30235786 34.666667 -28.9970363 350.2450 10.55208400 34.833333 7.7766795 356.9368 -23.81346781 35.000000 31.1795888 363.6286 -0.70821318 35.166667 19.7193020 366.9331 30.54758766 35.333333 -0.8038918 370.2376 0.46629555 35.500000 -28.8746459 364.5547 13.51994194 35.666667 -28.9970363 358.8718 -8.47477522 35.833333 7.7766795 345.0699 52.85342238 36.000000 31.1795888 331.2680 -19.54757358 36.166667 19.7193020 319.4466 -22.66593228 36.333333 -0.8038918 307.6253 -22.62138393 36.500000 -28.8746459 300.9598 -1.18513240 36.666667 -28.9970363 294.2943 23.50275558 36.833333 7.7766795 294.6227 -23.59937317 37.000000 31.1795888 294.9511 -1.73069548 37.166667 19.7193020 297.8307 -6.64999640 37.333333 -0.8038918 300.7103 -0.90639026 37.500000 -28.8746459 301.2058 0.66881530 37.666667 -28.9970363 301.7014 6.59565732 37.833333 7.7766795 295.7972 55.62609511 38.000000 31.1795888 289.8931 -16.07266066 38.166667 19.7193020 283.0445 -20.66379069 38.333333 -0.8038918 276.1959 -25.09201367 38.500000 -28.8746459 274.9421 0.43257600 38.666667 -28.9970363 273.6882 13.20880214 38.833333 7.7766795 278.7702 -20.04685892 39.000000 31.1795888 283.8521 0.86828646 39.166667 19.7193020 288.5987 10.08199064 39.333333 -0.8038918 293.3453 2.85860188 39.500000 -28.8746459 296.6553 -1.38066805 39.666667 -28.9970363 299.9653 -25.16830152 39.833333 7.7766795 300.4171 54.60622359 40.000000 31.1795888 300.8689 -7.14844487 40.166667 19.7193020 300.7066 -26.22586391 40.333333 -0.8038918 300.5443 -10.24037589 40.500000 -28.8746459 297.7665 26.30815903 40.666667 -28.9970363 294.9887 24.30833040 40.833333 7.7766795 293.6822 -29.45887477 41.000000 31.1795888 292.3757 -16.15527351 41.166667 19.7193020 290.8294 18.15130491 41.333333 -0.8038918 289.2831 4.42079038 41.500000 -28.8746459 290.6687 -12.69400434 41.666667 -28.9970363 292.0542 -32.65716261 41.833333 7.7766795 288.5982 65.12515850 42.000000 31.1795888 285.1421 5.37828605 42.166667 19.7193020 280.9869 -23.50620191 42.333333 -0.8038918 276.8317 -15.32778282 42.500000 -28.8746459 270.9253 8.94929808 42.666667 -28.9970363 265.0190 21.57801544 42.833333 7.7766795 265.8369 -31.81359017 43.000000 31.1795888 266.6548 -10.33438934 43.166667 19.7193020 269.7629 2.81780500 43.333333 -0.8038918 272.8710 2.63290639 43.500000 -28.8746459 278.4692 4.60539982 43.666667 -28.9970363 284.0675 -25.07047030 43.833333 7.7766795 287.2271 43.99617572 44.000000 31.1795888 290.3868 -3.36637182 44.166667 19.7193020 291.4135 -24.13283888 44.333333 -0.8038918 292.4403 4.16360112 44.500000 -28.8746459 295.6671 17.20750385 44.666667 -28.9970363 298.8940 1.10304303 44.833333 7.7766795 311.4186 -56.49523152 45.000000 31.1795888 323.9431 -14.52269962 45.166667 19.7193020 341.8846 17.79614353 45.333333 -0.8038918 359.8260 14.27789374 45.500000 -28.8746459 378.1176 5.95705004 45.666667 -28.9970363 396.4092 -29.01215720 45.833333 7.7766795 407.1676 51.95575660 46.000000 31.1795888 417.9259 1.89447684 46.166667 19.7193020 429.6863 -27.40559441 46.333333 -0.8038918 441.4467 -11.44275860 46.500000 -28.8746459 454.3915 0.38311626 46.666667 -28.9970363 467.3364 22.36062759 46.833333 7.7766795 480.4828 -24.65950602 47.000000 31.1795888 493.6292 16.59116681 47.166667 19.7193020 498.5917 25.88895580 47.333333 -0.8038918 503.5542 14.74965184 47.500000 -28.8746459 504.2054 -5.93074752 47.666667 -28.9970363 504.8565 -36.45951043 47.833333 7.7766795 501.2276 39.99571403 48.000000 31.1795888 497.5987 4.22174493 48.166667 19.7193020 495.9609 -9.58022648 48.333333 -0.8038918 494.3232 -9.51929084 48.500000 -28.8746459 494.1005 -8.22581596 48.666667 -28.9970363 493.8777 16.61929537 48.833333 7.7766795 497.8465 -36.12322252 49.000000 31.1795888 501.8153 11.70506602 49.166667 19.7193020 502.7126 18.76809349 49.333333 -0.8038918 503.6099 18.69402801 49.500000 -28.8746459 501.3347 -2.76004846 49.666667 -28.9970363 499.0595 -35.66248847 49.833333 7.7766795 492.1553 42.66805234 50.000000 31.1795888 485.2510 0.86939958 50.166667 19.7193020 477.1496 -11.16886172 50.333333 -0.8038918 469.0481 -2.44421597 50.500000 -28.8746459 460.1659 15.70876828 50.666667 -28.9970363 451.2836 4.31338899 50.833333 7.7766795 448.0273 -44.20400928 51.000000 31.1795888 444.7710 -8.45060110 51.166667 19.7193020 443.6025 21.17820913 51.333333 -0.8038918 442.4340 9.56992641 51.500000 -28.8746459 441.0238 5.25081214 51.666667 -28.9970363 439.6137 -30.71666567 51.833333 7.7766795 433.3083 43.61504523 52.000000 31.1795888 427.0028 -3.18243742 52.166667 19.7193020 423.5690 -22.48828553 52.333333 -0.8038918 420.1351 -2.83122660 52.500000 -28.8746459 424.4568 -19.28217734 52.666667 -28.9970363 428.7785 5.81850837 52.833333 7.7766795 439.3205 -41.29714203 53.000000 31.1795888 449.8624 19.75801400 53.166667 19.7193020 459.0868 35.19394757 53.333333 -0.8038918 468.3111 7.99278819 53.500000 -28.8746459 474.6787 -15.70402512 53.666667 -28.9970363 481.0462 -37.64920197 53.833333 7.7766795 487.8885 42.33486873 54.000000 31.1795888 494.7307 0.08974587 54.166667 19.7193020 513.7406 -44.95992380 54.333333 -0.8038918 532.7506 -11.74668643 54.500000 -28.8746459 573.0567 -39.78206648 54.666667 -28.9970363 613.3628 -15.86581007 54.833333 7.7766795 664.5789 -61.75553501 55.000000 31.1795888 715.7949 71.02554648 55.166667 19.7193020 753.0790 58.10166338 55.333333 -0.8038918 790.3632 46.34068733 55.500000 -28.8746459 799.1357 11.73894418 55.666667 -28.9970363 807.9082 -16.61116252 55.833333 7.7766795 797.4324 51.69096914 56.000000 31.1795888 786.9565 2.76390724 56.166667 19.7193020 775.2765 -25.39582495 56.333333 -0.8038918 763.5965 -10.59265009 56.500000 -28.8746459 758.6310 -5.35635641 56.666667 -28.9970363 753.6655 -1.56842627 56.833333 7.7766795 752.8677 -41.14436643 57.000000 31.1795888 752.0699 34.15049985 57.166667 19.7193020 746.1047 37.47598514 57.333333 -0.8038918 740.1395 13.16437749 57.500000 -28.8746459 733.4891 -15.61447245 57.666667 -28.9970363 726.8387 -67.44168594 57.833333 7.7766795 721.1671 36.55622074 58.000000 31.1795888 715.4955 11.02493386 58.166667 19.7193020 717.0271 -4.54641527 58.333333 -0.8038918 718.5587 -15.15485736 58.500000 -28.8746459 722.3472 -10.17253335 58.666667 -28.9970363 726.1356 12.36142711 58.833333 7.7766795 731.0228 -36.59946471 59.000000 31.1795888 735.9100 17.71044991 59.166667 19.7193020 730.4896 60.69107374 59.333333 -0.8038918 725.0693 31.33460462 59.500000 -28.8746459 712.9750 -27.30032241 59.666667 -28.9970363 700.8806 -56.78361299 59.833333 7.7766795 686.3632 51.16010537 60.000000 31.1795888 671.8458 -8.92536982 60.166667 19.7193020 663.6778 -7.69713394 60.333333 -0.8038918 655.5099 -11.00599102 60.500000 -28.8746459 648.8819 2.09274641 60.666667 -28.9970363 642.2539 21.34312029 60.833333 7.7766795 639.2547 -59.03137322 61.000000 31.1795888 636.2555 22.26493970 61.166667 19.7193020 630.0350 24.14574461 61.333333 -0.8038918 623.8144 24.88945657 61.500000 -28.8746459 617.5012 -19.82659636 61.666667 -28.9970363 611.1880 -36.49101284 61.833333 7.7766795 602.6245 22.19882129 62.000000 31.1795888 594.0609 18.55946186 62.166667 19.7193020 589.8617 -16.48104296 62.333333 -0.8038918 585.6625 -5.15864074 62.500000 -28.8746459 580.7335 -5.85889347 62.666667 -28.9970363 575.8045 16.09249025 62.833333 7.7766795 570.9230 -6.19971590 > m$win s t l 3721 11 7 > m$deg s t l 0 1 1 > m$jump s t l 373 2 1 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/fisher/rcomp/tmp/19gi51353953346.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(m,main=main) > dev.off() null device 1 > mylagmax <- nx/2 > postscript(file="/var/fisher/rcomp/tmp/2g1ky1353953346.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > op <- par(mfrow = c(2,2)) > acf(as.numeric(x),lag.max = mylagmax,main='Observed') > acf(as.numeric(m$time.series[,'trend']),na.action=na.pass,lag.max = mylagmax,main='Trend') > acf(as.numeric(m$time.series[,'seasonal']),na.action=na.pass,lag.max = mylagmax,main='Seasonal') > acf(as.numeric(m$time.series[,'remainder']),na.action=na.pass,lag.max = mylagmax,main='Remainder') > par(op) > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/3buot1353953346.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > op <- par(mfrow = c(2,2)) > spectrum(as.numeric(x),main='Observed') > spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'trend']),'trend']),main='Trend') > spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'seasonal']),'seasonal']),main='Seasonal') > spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'remainder']),'remainder']),main='Remainder') > par(op) > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/4qk241353953346.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > op <- par(mfrow = c(2,2)) > cpgram(as.numeric(x),main='Observed') > cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'trend']),'trend']),main='Trend') > cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'seasonal']),'seasonal']),main='Seasonal') > cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'remainder']),'remainder']),main='Remainder') > par(op) > dev.off() null device 1 > > #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Seasonal Decomposition by Loess - Parameters',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Component',header=TRUE) > a<-table.element(a,'Window',header=TRUE) > a<-table.element(a,'Degree',header=TRUE) > a<-table.element(a,'Jump',header=TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Seasonal',header=TRUE) > a<-table.element(a,m$win['s']) > a<-table.element(a,m$deg['s']) > a<-table.element(a,m$jump['s']) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Trend',header=TRUE) > a<-table.element(a,m$win['t']) > a<-table.element(a,m$deg['t']) > a<-table.element(a,m$jump['t']) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Low-pass',header=TRUE) > a<-table.element(a,m$win['l']) > a<-table.element(a,m$deg['l']) > a<-table.element(a,m$jump['l']) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/fisher/rcomp/tmp/5xj7j1353953346.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Seasonal Decomposition by Loess - Time Series Components',6,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'t',header=TRUE) > a<-table.element(a,'Observed',header=TRUE) > a<-table.element(a,'Fitted',header=TRUE) > a<-table.element(a,'Seasonal',header=TRUE) > a<-table.element(a,'Trend',header=TRUE) > a<-table.element(a,'Remainder',header=TRUE) > a<-table.row.end(a) > for (i in 1:nx) { + a<-table.row.start(a) + a<-table.element(a,i,header=TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]+m$time.series[i,'remainder']) + a<-table.element(a,m$time.series[i,'seasonal']) + a<-table.element(a,m$time.series[i,'trend']) + a<-table.element(a,m$time.series[i,'remainder']) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/fisher/rcomp/tmp/6y8uq1353953346.tab") > > try(system("convert tmp/19gi51353953346.ps tmp/19gi51353953346.png",intern=TRUE)) character(0) > try(system("convert tmp/2g1ky1353953346.ps tmp/2g1ky1353953346.png",intern=TRUE)) character(0) > try(system("convert tmp/3buot1353953346.ps tmp/3buot1353953346.png",intern=TRUE)) character(0) > try(system("convert tmp/4qk241353953346.ps tmp/4qk241353953346.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.033 0.646 7.659