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 = '6' > par1 = '12' > main = 'Seasonal Decomposition by Loess' > par8 <- 'FALSE' > par7 <- '1' > par6 <- '' > par5 <- '1' > par4 <- '' > par3 <- '0' > par2 <- '6' > par1 <- '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > 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 seasonal trend remainder Jan 1 37.57226955 189.0271 8.500657325 Feb 1 60.93370676 196.5908 23.175506819 Mar 1 35.27362497 204.1545 25.171875315 Apr 1 2.23399731 211.7182 26.747789688 May 1 -15.82861959 220.0157 -2.787080728 Jun 1 19.07419831 228.3132 -6.587385943 Jul 1 22.42640368 236.6107 -17.937078634 Aug 1 -15.14368342 245.4800 -6.536364043 Sep 1 -30.81302798 254.3494 -17.436392007 Oct 1 -49.54761251 263.2188 -38.971179988 Nov 1 -34.48913860 272.3762 -34.587060455 Dec 1 -28.34370049 281.5336 -32.689905123 Jan 2 37.88861516 290.6910 -29.079627334 Feb 2 59.44531189 303.1207 -15.166035730 Mar 2 33.81114956 315.5504 -11.061585076 Apr 2 1.20467808 327.9801 -1.484825270 May 2 -16.22107291 339.4620 28.359117932 Jun 2 18.60192652 350.9438 27.054310715 Jul 2 21.69965931 362.4256 54.674770144 Aug 2 -15.08983951 369.0420 41.647801838 Sep 2 -29.42483407 375.6585 17.266329275 Oct 2 -47.98121356 382.2750 44.506241634 Nov 2 -32.15079408 380.6932 8.457615402 Dec 2 -26.80538641 379.1114 16.694000974 Jan 3 38.00057987 377.5296 49.269827936 Feb 3 56.26905633 369.0976 53.733387415 Mar 3 30.57173210 360.6655 40.062747587 Apr 3 -1.16110043 352.2335 15.427616062 May 3 -16.32940907 340.0510 2.578434321 Jun 3 18.31319988 327.8685 8.918334990 Jul 3 21.69390404 315.6860 -5.779859543 Aug 3 -14.29345506 302.2236 -26.630114628 Sep 3 -26.48414216 288.7612 -13.277041714 Oct 3 -43.65252695 275.2988 -26.146271106 Nov 3 -27.03660231 263.9695 -1.332857019 Dec 3 -22.57809797 252.6401 10.837977373 Jan 4 41.80887738 241.3108 -18.219659242 Feb 4 50.71388702 234.1893 -31.103191840 Mar 4 24.56019755 227.0678 -19.328025328 Apr 4 -5.88008951 219.9464 -20.266261220 May 4 -19.03555623 215.4045 -19.368979688 Jun 4 14.84632871 210.8627 -12.509049808 Jul 4 15.69062230 206.3209 -14.811528575 Aug 4 -17.58272862 203.7590 -5.576311093 Sep 4 -22.47877439 201.1972 9.881601237 Oct 4 -39.70528645 198.6353 16.469979867 Nov 4 -17.67395489 198.0657 18.608259826 Dec 4 -16.10472576 197.4961 -1.791357778 Jan 5 42.28927837 196.9265 -13.415750375 Feb 5 44.72425296 196.1362 -6.860499393 Mar 5 23.35883516 195.3460 -18.504856015 Apr 5 0.16643245 194.5558 -11.122227726 May 5 -16.81705853 192.4001 2.616913108 Jun 5 8.26928084 190.2445 4.686223592 Jul 5 5.80770400 188.0888 14.603450287 Aug 5 -20.68492581 185.8964 26.588561125 Sep 5 -23.62317433 183.7039 12.719290681 Oct 5 -45.48818691 181.5114 11.976784295 Nov 5 -21.55861607 180.1495 0.809095125 Dec 5 -21.29679339 178.7876 -2.990845895 Jan 6 45.37272991 177.4258 -9.598487527 Feb 6 49.20294955 178.5435 -31.346443937 Mar 6 30.98728352 179.6612 -27.848514683 Apr 6 12.64620381 180.7790 -17.025171737 May 6 -8.91273320 187.8951 -25.382324252 Jun 6 6.99096124 195.0111 -28.802108217 Jul 6 2.37904003 202.1272 -33.506276528 Aug 6 -19.32980138 214.8229 -44.293103574 Sep 6 -25.71983106 227.5186 -39.898742342 Oct 6 -50.64315804 240.2142 -32.371083820 Nov 6 -27.64174772 254.6180 -25.276217558 Dec 6 -25.26208144 269.0217 -7.359607252 Jan 7 48.03381899 283.4254 24.640768899 Feb 7 52.17438344 295.5725 50.553122437 Mar 7 37.72710928 307.7196 58.253314582 Apr 7 19.47649222 319.8667 45.256849633 May 7 -5.62534202 326.3213 45.104016410 Jun 7 9.48249145 332.7760 25.841515484 Jul 7 -0.56798326 339.2307 29.237322729 Aug 7 -23.69644581 338.2624 32.434066788 Sep 7 -32.91102567 337.2941 38.916928159 Oct 7 -56.26333017 336.3258 12.837514163 Nov 7 -28.66340315 330.6809 9.482490879 Dec 7 -22.85384148 325.0360 -1.282167049 Jan 8 48.21342882 319.3911 -0.704533611 Feb 8 52.13650576 313.1975 -8.434019274 Mar 8 39.99583349 307.0039 -17.299755721 Apr 8 17.00054997 300.8103 -1.610880917 May 8 -2.86164809 295.6092 -23.747565172 Jun 8 18.69517090 290.4081 -19.803266475 Jul 8 2.38239993 285.2070 -21.389377827 Aug 8 -26.33499713 282.3829 -2.447884751 Sep 8 -38.43534606 279.5588 -7.323439799 Oct 8 -59.72103017 276.7347 11.386340321 Nov 8 -24.91158920 276.1448 2.366803611 Dec 8 -19.95570250 275.5549 4.500821170 Jan 9 46.33844657 274.9650 -14.703423638 Feb 9 50.71414189 275.2704 -16.784511452 Mar 9 36.56365747 275.5758 -2.639419523 Apr 9 12.73047205 275.8812 -17.611626592 May 9 -2.55007733 275.5284 6.921649322 Jun 9 28.62494362 275.1757 14.099354909 Jul 9 6.64982963 274.8230 16.927195429 Aug 9 -28.83586953 274.2265 1.309418257 Sep 9 -48.21518822 273.6299 1.885260620 Oct 9 -66.25406887 273.0334 2.320664937 Nov 9 -28.60516567 273.2510 15.254179102 Dec 9 -20.78003321 273.4686 13.311464000 Jan 10 45.92872797 273.6862 0.985120173 Feb 10 56.43646804 276.3482 -24.284624040 Mar 10 40.63777366 279.0102 -37.447933815 Apr 10 11.51267603 281.6722 -30.484840336 May 10 -1.55766337 289.7598 -24.702133201 Jun 10 36.46989142 297.8474 -21.217320260 Jul 10 13.38976893 305.9351 -35.024830048 Aug 10 -27.87737825 319.6486 -39.171255127 Sep 10 -53.56876944 333.3622 -29.493436191 Oct 10 -69.86962676 347.0758 -30.706151122 Nov 10 -32.76640076 362.4992 -17.032824728 Dec 10 -20.98118938 377.9227 -23.741483719 Jan 11 49.37464620 393.3461 3.679233101 Feb 11 60.91225776 405.1264 45.561326222 Mar 11 45.32312214 416.9067 53.270166532 Apr 11 6.75244368 428.6870 70.960549683 May 11 -9.85740941 434.0039 59.053484106 Jun 11 34.53445947 439.3208 48.444696562 Jul 11 11.35031418 444.6378 53.811923196 Aug 11 -27.37913303 441.8155 46.263627449 Sep 11 -56.51791900 438.9932 23.324670462 Oct 11 -68.03058940 436.1710 6.859597910 Nov 11 -33.16178944 428.1844 -16.522642543 Dec 11 -19.97155814 420.1979 6.573685664 Jan 12 53.12933207 412.2113 2.459354965 Feb 12 63.77480244 404.6520 1.373183178 Mar 12 52.82232669 397.0927 -20.115042494 Apr 12 5.03914429 389.5334 -38.772561511 May 12 -17.72176982 384.0755 -33.653713367 Jun 12 34.77422462 378.6175 -35.391773778 Jul 12 10.53029230 373.1596 -23.189907418 Aug 12 -23.91314628 370.5614 -11.948244234 Sep 12 -58.42533021 367.9632 9.962164301 Oct 12 -63.27748966 365.3649 21.012548354 Nov 12 -33.45316517 366.1982 30.854993792 Dec 12 -14.64281557 367.0314 -0.288585875 Jan 13 56.38181494 367.8646 -12.346446464 Feb 13 62.71255453 371.5685 -45.681019450 Mar 13 52.53236681 375.2723 -11.404665124 Apr 13 -2.46622523 378.9761 -15.809906488 May 13 -22.28100826 385.4301 -25.149045813 Jun 13 34.71509971 391.8840 -9.399076140 Jul 13 7.38007848 398.3379 -17.317977268 Aug 13 -22.42124604 407.1318 -13.610550169 Sep 13 -57.54171778 415.9257 -26.883975841 Oct 13 -59.98422929 424.7196 -11.035361746 Nov 13 -29.95439868 432.9243 -6.269937137 Dec 13 -12.69136884 441.1291 18.562288229 Jan 14 57.15849747 449.3338 27.007677132 Feb 14 59.35968858 453.6753 52.365028956 Mar 14 47.49314960 458.0167 36.790110856 Apr 14 -5.95679635 462.3582 32.298599735 May 14 -18.63309520 461.4031 24.329956040 Jun 14 37.27032350 460.4481 33.581594797 Jul 14 3.62086049 459.4930 32.986115266 Aug 14 -22.33851177 453.9979 12.340584946 Sep 14 -58.90770270 448.5028 13.804873306 Oct 14 -64.70490107 443.0077 7.997169098 Nov 14 -30.53677706 436.0956 -11.458775540 Dec 14 -16.15685792 429.1834 -8.926515309 Jan 15 57.78110682 422.2712 -17.952300679 Feb 15 65.27930238 416.8497 -34.029028797 Mar 15 44.27779407 411.4283 -23.406053036 Apr 15 -3.60971803 406.0068 -16.097073493 May 15 -12.49271293 403.5896 4.103161383 Jun 15 44.93413036 401.1723 -24.206441922 Jul 15 4.04169367 398.7551 -19.896765254 Aug 15 -20.58142839 398.9120 5.869402200 Sep 15 -55.05919138 399.0690 1.490210589 Oct 15 -64.93038586 399.2259 -10.895549524 Nov 15 -29.71696785 400.1488 2.168128626 Dec 15 -20.33445873 401.0717 -4.737284327 Jan 16 53.81690952 401.9946 6.888443590 Feb 16 63.15640675 402.6541 21.189536831 Mar 16 35.63252675 403.3135 5.254007314 Apr 16 -4.67765574 403.9729 0.004780272 May 16 -14.34407449 403.8625 5.381550237 Jun 16 50.46560209 403.7522 1.182224872 Jul 16 0.84735866 403.6418 9.510819509 Aug 16 -18.67530367 402.2892 -8.113939110 Sep 16 -49.25223934 400.9367 -4.684424394 Oct 16 -59.55873254 399.5841 -0.625352141 Nov 16 -25.53345700 397.3113 14.022135436 Dec 16 -21.89128431 395.0386 5.652725853 Jan 17 47.74967473 392.7658 11.284529929 Feb 17 58.02740298 390.1625 -2.089860422 Mar 17 27.89264016 387.5591 7.048240291 Apr 17 -4.91502953 384.9558 3.059247883 May 17 -17.80854873 382.1058 -11.497246460 Jun 17 59.08152881 379.2558 6.962662467 Jul 17 3.24412395 376.4058 -12.149946220 Aug 17 -16.57911569 373.2477 -1.568628483 Sep 17 -44.98130884 370.0897 1.091642764 Oct 17 -52.25117016 366.9316 5.119582171 Nov 17 -23.73040457 363.5339 -8.003526301 Dec 17 -21.75897865 360.1363 2.522704900 Jan 18 40.05706092 356.7386 -2.695677557 Feb 18 46.92130854 352.7846 17.494074427 Mar 18 20.12816116 348.8306 0.941221419 Apr 18 -7.48066886 344.8766 11.804051048 May 18 -21.63194996 340.1225 2.909470766 Jun 18 67.03911535 335.3683 3.292544072 Jul 18 9.73706643 330.6142 2.548731612 Aug 18 -11.34034059 325.8697 1.970620480 Sep 18 -37.02201146 321.1252 0.096773203 Oct 18 -41.17490890 316.3808 -4.305847513 Nov 18 -23.55788596 312.1295 0.228370095 Dec 18 -22.11247795 307.8783 -6.965797370 Jan 19 30.44120395 303.6270 -9.668238726 Feb 19 36.79319668 300.3537 -26.246853454 Mar 19 11.51315873 297.0803 -9.593437491 Apr 19 -14.91765731 293.8069 -5.889243448 May 19 -28.42844503 292.0567 15.671727910 Jun 19 70.49892287 290.3065 -1.605456361 Jul 19 19.35651849 288.5563 -2.912868340 Aug 19 -5.69303265 288.1192 -0.326133794 Sep 19 -26.99398928 287.6820 -10.387993766 Oct 19 -29.75691623 287.2448 -10.987883410 Nov 19 -19.71439499 287.6907 -10.076341822 Dec 19 -23.58397485 288.1367 1.947300849 Jan 20 22.84962095 288.5826 4.467767876 Feb 20 29.65444799 289.7096 -0.964036895 Mar 20 5.62328389 290.8366 -1.059850535 Apr 20 -23.44895938 291.9635 -2.114584996 May 20 -36.78377260 292.9905 -10.406736956 Jun 20 71.19219335 294.0175 -2.409668088 Jul 20 28.79058575 295.0444 1.064974336 Aug 20 -2.68745959 295.1147 1.772776162 Sep 20 -18.39155631 295.1849 12.706629364 Oct 20 -22.34799278 295.2552 22.292822321 Nov 20 -18.97768244 294.5960 14.681708634 Dec 20 -28.26222800 293.9368 6.325450855 Jan 21 18.28751670 293.2776 -4.165097180 Feb 21 27.84831661 291.7283 9.123410406 Mar 21 4.40130651 290.1790 -1.680271993 Apr 21 -26.30009655 288.6297 -13.229561438 May 21 -45.82581908 286.3961 -10.170282854 Jun 21 68.94161736 284.1625 8.395836751 Jul 21 33.78208635 281.9290 5.988923811 Aug 21 -2.52581890 280.4202 -0.694399623 Sep 21 -11.14141746 278.9114 -7.070029739 Oct 21 -18.68272292 277.4027 -7.719952960 Nov 21 -19.93027713 277.3458 0.184510689 Dec 21 -33.33997654 277.2889 -2.148880473 Jan 22 18.35284619 277.2319 -8.084793770 Feb 22 31.95235272 278.4547 -18.107044950 Mar 22 7.46605411 279.6774 -12.443490986 Apr 22 -25.77150672 280.9002 -0.928674811 May 22 -51.12397530 284.3240 -3.200042193 Jun 22 65.13132478 287.7479 -13.879178232 Jul 22 35.21482056 291.1717 -8.186509974 Aug 22 -3.08090487 297.8578 -7.776888687 Sep 22 -10.35017778 304.5439 1.606280074 Oct 22 -23.92623027 311.2300 -3.303771578 Nov 22 -19.99702827 320.7187 -29.721691214 Dec 22 -33.61760108 330.2074 -33.889836033 Jan 23 26.58024543 339.6962 -25.676400177 Feb 23 36.58379104 351.4644 -8.648156821 Mar 23 13.05524032 363.2326 -2.987817130 Apr 23 -24.14276818 375.0008 4.341980342 May 23 -53.04387550 388.0382 3.405724307 Jun 23 59.07763993 401.0755 6.746845535 Jul 23 33.80159142 414.1129 3.085530691 Aug 23 -1.66105203 425.7799 -2.118824814 Sep 23 -11.27170213 437.4469 3.024826333 Oct 23 -29.63019011 449.1139 6.416315362 Nov 23 -22.47352580 457.4884 25.685138141 Dec 23 -34.30575149 465.8629 32.042850914 Jan 24 33.07178314 474.2374 34.090803357 Feb 24 41.14754177 479.5421 23.510407668 Mar 24 17.81415127 484.8467 14.839161104 Apr 24 -23.24732298 490.1513 2.495998298 May 24 -54.67061015 492.5802 1.490377929 Jun 24 55.68314844 495.0091 -1.692288188 Jul 24 32.28070526 497.4380 3.281247465 Aug 24 -0.03578026 497.6120 8.523743017 Sep 24 -11.69944387 497.7860 -2.086583335 Oct 24 -32.73463124 497.9600 -8.225385929 Nov 24 -22.46387307 497.3897 6.574161678 Dec 24 -31.84327320 496.8194 4.523867582 Jan 25 39.30505557 496.2491 9.145844576 Feb 25 46.72363888 494.4230 0.053380210 Mar 25 19.36608110 492.5969 9.537056932 Apr 25 -24.74749929 490.7707 3.676756264 May 25 -57.34281460 486.8795 4.863290158 Jun 25 51.96934944 482.9883 7.642344706 Jul 25 28.24987227 479.0971 9.953040465 Aug 25 -5.47227534 474.3704 16.801828147 Sep 25 -16.39284024 469.6438 12.549033112 Oct 25 -43.27663316 464.9172 25.359466092 Nov 25 -28.22300521 460.5879 -5.764889013 Dec 25 -33.85474310 456.2586 -10.803878280 Jan 26 51.56552987 451.9293 -35.994878400 Feb 26 57.21011523 448.5578 -21.267899112 Mar 26 27.26191983 445.1862 -21.248139057 Apr 26 -20.93945656 441.8147 -3.475198009 May 26 -53.95440544 440.1812 -6.326795935 Jun 26 51.70115250 438.5477 -5.548900675 Jul 26 24.11176274 436.9143 -6.026057711 Aug 26 -16.57303030 438.1755 -0.802466292 Sep 26 -24.05922663 439.4367 1.122528421 Oct 26 -55.38643641 440.6979 -9.011463421 Nov 26 -39.33778707 445.3383 -0.400524905 Dec 26 -39.85848959 449.9787 -4.320234519 Jan 27 66.27832688 454.6191 -20.097463133 Feb 27 69.68380751 464.4693 -20.153123707 Mar 27 38.03064154 474.3195 -36.850137686 Apr 27 -15.28329439 484.1697 -38.786381700 May 27 -50.60685028 501.7523 -36.745483858 Jun 27 53.63337052 519.3350 -34.968362709 Jul 27 23.31726264 536.9177 -34.234912879 Aug 27 -21.75219409 560.2993 -50.047058872 Sep 27 -29.83507402 583.6809 -33.645781675 Oct 27 -62.30050682 607.0625 -40.361951595 Nov 27 -42.38707803 630.5635 -19.676463457 Dec 27 -39.64433035 654.0646 -3.820294205 Jan 28 76.97667065 677.5657 63.457621725 Feb 28 78.71161777 696.0790 56.109357396 Mar 28 42.33671703 714.5923 78.970940921 Apr 28 -20.30090618 733.1057 69.195246923 May 28 -58.82682491 743.0357 78.091150702 Jun 28 52.33867452 752.9657 51.595636318 Jul 28 22.99575127 762.8957 35.008544614 Aug 28 -18.06846139 763.2905 24.377960313 Sep 28 -31.07769470 763.6853 19.592396654 Oct 28 -59.15475279 764.0801 19.474657767 Nov 28 -36.04442679 759.5909 -0.446477377 Dec 28 -36.64662922 755.1017 1.044915902 Jan 29 73.88167853 750.6125 -7.094200992 Feb 29 76.78305098 746.3455 -19.828549205 Mar 29 40.71228887 742.0785 -30.290762859 Apr 29 -29.29381434 737.8114 -19.517635413 May 29 -65.08370978 734.8830 -39.399246870 Jun 29 48.01445350 731.9545 -14.468917055 Jul 29 24.18934416 729.0260 4.484685387 Aug 29 -12.18590309 726.7447 17.641173317 Sep 29 -27.75682606 724.4635 5.893336962 Oct 29 -52.77579215 722.1822 13.893543732 Nov 29 -31.16790997 719.6520 21.015928839 Dec 29 -35.95599530 717.1217 21.034281452 Jan 30 72.08715818 714.5914 -1.878604741 Feb 30 74.33339301 710.0264 26.540189064 Mar 30 38.95912967 705.4614 11.179481034 Apr 30 -32.89277856 700.8964 -11.203582105 May 30 -67.59625657 693.6566 -10.960386685 Jun 30 45.45123983 686.4169 13.431834330 Jul 30 25.13546138 679.1772 -10.212669802 Aug 30 -10.18451964 671.3080 14.576524125 Sep 30 -26.41519880 663.4388 6.676416201 Oct 30 -50.63603760 655.5696 17.166467915 Nov 30 -29.43099331 649.1209 14.910115030 Dec 30 -34.48395666 642.6722 -20.188230217 Jan 31 71.18087475 636.2235 -17.704370225 Feb 31 73.56612665 631.1595 -30.825676119 Mar 31 38.11039188 626.0956 -16.305995354 Apr 31 -34.91239997 621.0317 -17.319257503 May 31 -69.28883558 616.0708 -1.081965486 Jun 31 44.42345613 611.1099 -22.933400782 Jul 31 25.46767011 606.1491 12.183241641 Aug 31 -8.75187660 601.4812 0.370684380 Sep 31 -25.52877834 596.8133 8.415482146 Oct 31 -49.11524389 592.1454 2.969843733 Nov 31 -28.05855974 587.7810 3.177555621 Dec 31 -33.90791427 583.4166 22.991306186 > m$win s t l 6 25 13 > m$deg s t l 0 1 1 > m$jump s t l 1 3 2 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/wessaorg/rcomp/tmp/1d1ek1352541591.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/wessaorg/rcomp/tmp/2vhfw1352541591.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/wessaorg/rcomp/tmp/3dwey1352541591.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/wessaorg/rcomp/tmp/4bp2v1352541591.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/5qjkx1352541591.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/wessaorg/rcomp/tmp/6fbce1352541591.tab") > > try(system("convert tmp/1d1ek1352541591.ps tmp/1d1ek1352541591.png",intern=TRUE)) character(0) > try(system("convert tmp/2vhfw1352541591.ps tmp/2vhfw1352541591.png",intern=TRUE)) character(0) > try(system("convert tmp/3dwey1352541591.ps tmp/3dwey1352541591.png",intern=TRUE)) character(0) > try(system("convert tmp/4bp2v1352541591.ps tmp/4bp2v1352541591.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.927 0.630 8.544