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) > par8 = 'FALSE' > par7 = '1' > par6 = '' > par5 = '1' > par4 = '' > par3 = '0' > par2 = 'periodic' > par1 = '12' > main = 'Seasonal Decomposition by Loess' > par8 <- 'FALSE' > par7 <- '1' > par6 <- '' > par5 <- '1' > par4 <- '' > par3 <- '0' > par2 <- 'periodic' > 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 46.126734 209.0927 -20.11944274 Feb 1 54.944778 212.7586 12.99664915 Mar 1 31.598316 216.4244 16.57724743 Apr 1 -8.339539 220.5017 28.53782168 May 1 -29.893525 224.5790 6.71452690 Jun 1 41.489279 229.1346 -29.82391133 Jul 1 16.349030 233.6903 -8.93929714 Aug 1 -15.202202 238.4357 0.56650141 Sep 1 -33.490109 243.1811 -3.59102519 Oct 1 -49.202316 250.2550 -26.35272718 Nov 1 -27.837860 257.3290 -26.19109260 Dec 1 -26.542588 270.4341 -23.39153801 Jan 2 46.126734 283.5393 -30.16603370 Feb 2 54.944778 299.1676 -6.71238379 Mar 2 31.598316 314.7959 -8.09422749 Apr 2 -8.339539 330.4575 5.58200582 May 2 -29.893525 346.1192 35.37437011 Jun 2 41.489279 359.6833 -4.57258986 Jul 2 16.349030 373.2475 49.20350259 Aug 2 -15.202202 382.7240 28.07816686 Sep 2 -33.490109 392.2006 4.78950598 Oct 2 -49.202316 395.2793 32.72301221 Nov 2 -27.837860 398.3580 -13.52014499 Dec 2 -26.542588 394.2365 1.30609665 Jan 3 46.126734 390.1150 28.55828802 Feb 3 54.944778 380.5655 43.58973304 Mar 3 31.598316 371.0160 28.68568447 Apr 3 -8.339539 358.6552 16.18438016 May 3 -29.893525 346.2943 9.89920685 Jun 3 41.489279 331.8400 -18.22927085 Jul 3 16.349030 317.3857 -2.13469613 Aug 3 -15.202202 301.3570 -24.85476871 Sep 3 -33.490109 285.3283 -2.83816643 Oct 3 -49.202316 270.8836 -16.18124879 Nov 3 -27.837860 256.4389 6.99900542 Dec 3 -26.542588 245.0606 22.38200980 Jan 4 46.126734 233.6823 -14.90903609 Feb 4 54.944778 225.8005 -26.94524843 Mar 4 31.598316 217.9186 -17.21695437 Apr 4 -8.339539 213.4437 -11.30418352 May 4 -29.893525 208.9688 -2.07528170 Jun 4 41.489279 206.4705 -34.75978511 Jul 4 16.349030 203.9722 -13.12123612 Aug 4 -15.202202 202.3462 -6.54398449 Sep 4 -33.490109 200.7202 21.36994200 Oct 4 -49.202316 199.6364 24.96589645 Nov 4 -27.837860 198.5527 28.28518747 Dec 4 -26.542588 197.6105 8.53205342 Jan 5 46.126734 196.6684 -16.99513090 Feb 5 54.944778 195.4459 -16.39068080 Mar 5 31.598316 194.2234 -25.62172429 Apr 5 -8.339539 192.7055 -0.76592645 May 5 -29.893525 191.1875 16.90600238 Jun 5 41.489279 189.9326 -28.22190370 Jul 5 16.349030 188.6777 3.47324263 Aug 5 -15.202202 187.0605 19.94173199 Sep 5 -33.490109 185.4432 20.84689621 Oct 5 -49.202316 183.3135 13.88886477 Nov 5 -27.837860 181.1837 6.05416989 Dec 5 -26.542588 178.2286 2.81398107 Jan 6 46.126734 175.2735 -8.20025803 Feb 6 54.944778 172.7688 -31.31357528 Mar 6 31.598316 170.2641 -19.06238613 Apr 6 -8.339539 171.1680 13.57157365 May 6 -29.893525 172.0719 11.42166442 Jun 6 41.489279 179.4224 -47.71171282 Jul 6 16.349030 186.7730 -32.12203764 Aug 6 -15.202202 201.0319 -34.62968075 Sep 6 -33.490109 215.2908 -19.90064900 Oct 6 -49.202316 233.9219 -27.51959703 Nov 6 -27.837860 252.5531 -23.01520850 Dec 6 -26.542588 271.4826 -8.53997955 Jan 7 46.126734 290.4121 19.56119912 Feb 7 54.944778 306.8457 36.50954348 Mar 7 31.598316 323.2793 48.82239425 Apr 7 -8.339539 334.4191 58.52045021 May 7 -29.893525 345.5589 50.13463717 Jun 7 41.489279 349.4676 -22.85687192 Jul 7 16.349030 353.3763 -1.82532860 Aug 7 -15.202202 350.4052 11.79698397 Sep 7 -33.490109 347.4341 29.35597140 Oct 7 -49.202316 341.3191 0.78318722 Nov 7 -27.837860 335.2041 4.13373962 Dec 7 -26.542588 327.7027 -0.26006679 Jan 8 46.126734 320.2012 0.57207652 Feb 8 54.944778 312.0954 -10.14021037 Mar 8 31.598316 303.9897 -5.88799085 Apr 8 -8.339539 297.2512 27.28837738 May 8 -29.893525 290.5126 8.38087659 Jun 8 41.489279 285.9038 -38.09309060 Jul 8 16.349030 281.2950 -31.44400537 Aug 8 -15.202202 277.9576 -9.15538893 Sep 8 -33.490109 274.6202 -7.33009763 Oct 8 -49.202316 273.5609 4.04141306 Nov 8 -27.837860 272.5016 8.93626033 Dec 8 -26.542588 273.6282 13.01437886 Jan 9 46.126734 274.7548 -14.28155289 Feb 9 54.944778 275.3926 -21.13739852 Mar 9 31.598316 276.0304 1.87126226 Apr 9 -8.339539 275.7407 3.59887505 May 9 -29.893525 275.4509 34.34261882 Jun 9 41.489279 275.5695 0.84125173 Jul 9 16.349030 275.6880 6.36293705 Aug 9 -15.202202 275.2198 -13.31755071 Sep 9 -33.490109 274.7515 -13.96136362 Oct 9 -49.202316 273.5519 -15.24961004 Nov 9 -27.837860 272.3524 15.38548010 Dec 9 -26.542588 271.7534 20.78917746 Jan 10 46.126734 271.1544 3.31882455 Feb 10 54.944778 271.7381 -18.18286585 Mar 10 31.598316 272.3217 -21.72004986 Apr 10 -8.339539 274.7846 -3.74503082 May 10 -29.893525 277.2474 16.14611921 Jun 10 41.489279 284.2742 -12.66348711 Jul 10 16.349030 291.3010 -23.35004100 Aug 10 -15.202202 305.4530 -37.65080616 Sep 10 -33.490109 319.6050 -35.81489647 Oct 10 -49.202316 338.9805 -43.27821301 Nov 10 -27.837860 358.3561 -17.81819299 Dec 10 -26.542588 378.9361 -19.19354991 Jan 11 46.126734 399.5162 0.75704290 Feb 11 54.944778 417.2321 39.42310013 Mar 11 31.598316 434.9480 48.95366377 Apr 11 -8.339539 445.6620 69.07750876 May 11 -29.893525 456.3760 56.71748473 Jun 11 41.489279 459.0554 21.75531493 Jul 11 16.349030 461.7348 31.71619754 Aug 11 -15.202202 456.9927 18.90952147 Sep 11 -33.490109 452.2506 -12.96047974 Oct 11 -49.202316 442.1374 -17.93508151 Nov 11 -27.837860 432.0242 -25.68634671 Dec 11 -26.542588 420.4073 12.93531073 Jan 12 46.126734 408.7903 12.88291789 Feb 12 54.944778 399.0920 15.76324747 Mar 12 31.598316 389.3936 8.80808344 Apr 12 -8.339539 383.1830 -19.04342403 May 12 -29.893525 376.9723 -14.37880053 Jun 12 41.489279 372.7172 -36.20650383 Jul 12 16.349030 368.4621 -24.31115472 Aug 12 -15.202202 365.8122 -15.91004352 Sep 12 -33.490109 363.1624 -10.17225745 Oct 12 -49.202316 363.2625 9.03979559 Nov 12 -27.837860 363.3627 28.07518520 Dec 12 -26.542588 365.2688 13.37379047 Jan 13 46.126734 367.1749 -1.40165452 Feb 13 54.944778 368.9588 -35.30358238 Mar 13 31.598316 370.7427 14.05899617 Apr 13 -8.339539 373.3711 -4.33157458 May 13 -29.893525 375.9995 -8.10601435 Jun 13 41.489279 382.7120 -7.00130920 Jul 13 16.349030 389.4245 -17.37355163 Aug 13 -15.202202 400.5490 -14.24679237 Sep 13 -33.490109 411.6735 -46.68335826 Oct 13 -49.202316 423.8385 -20.93614944 Nov 13 -27.837860 436.0035 -11.46560405 Dec 13 -26.542588 447.0178 26.52476735 Jan 14 46.126734 458.0322 29.34108848 Feb 14 54.944778 465.4926 44.96261602 Mar 14 31.598316 472.9530 37.74864997 Apr 14 -8.339539 474.6837 22.35586188 May 14 -29.893525 476.4143 20.57920477 Jun 14 41.489279 471.9832 17.82749632 Jul 14 16.349030 467.5521 12.19884029 Aug 14 -15.202202 459.2964 -0.09422082 Sep 14 -33.490109 451.0407 -14.15060707 Oct 14 -49.202316 442.3659 -6.86355315 Nov 14 -27.837860 433.6910 -11.75316267 Dec 14 -26.542588 425.9409 4.70168231 Jan 15 46.126734 418.1908 -2.21752299 Feb 15 54.944778 412.0542 -18.89898192 Mar 15 31.598316 405.9176 -5.21593444 Apr 15 -8.339539 401.7489 -7.10931315 May 15 -29.893525 397.5801 27.51343911 Jun 15 41.489279 395.9635 -15.55278039 Jul 15 16.349030 394.3469 -27.79594748 Aug 15 -15.202202 395.0895 4.31271957 Sep 15 -33.490109 395.8320 -16.84193853 Oct 15 -49.202316 397.6986 -25.09631384 Nov 15 -27.837860 399.5652 0.87264742 Dec 15 -26.542588 401.6947 0.84785277 Jan 16 46.126734 403.8243 12.74900786 Feb 16 54.944778 405.1815 26.87370119 Mar 16 31.598316 406.5388 6.06290093 Apr 16 -8.339539 406.6519 0.98764740 May 16 -29.893525 406.7650 18.02852485 Jun 16 41.489279 405.8682 8.04248312 Jul 16 16.349030 404.9715 -7.32050620 Aug 16 -15.202202 403.2118 -12.50959120 Sep 16 -33.490109 401.4521 -20.96200136 Oct 16 -49.202316 399.5608 -10.95845487 Nov 16 -27.837860 397.6694 15.96842818 Dec 16 -26.542588 395.8147 9.52787364 Jan 17 46.126734 393.9600 11.71326883 Feb 17 54.944778 391.8158 -0.66059104 Mar 17 31.598316 389.6716 1.23005550 Apr 17 -8.339539 386.6205 4.81902399 May 17 -29.893525 383.5694 -0.87587653 Jun 17 41.489279 379.8298 23.98091174 Jul 17 16.349030 376.0902 -24.93924757 Aug 17 -15.202202 372.4366 -2.13442637 Sep 17 -33.490109 368.7830 -9.09293032 Oct 17 -49.202316 365.6265 3.37585530 Nov 17 -27.837860 362.4699 -2.83202251 Dec 17 -26.542588 359.8229 7.61973468 Jan 18 46.126734 357.1758 -9.20255839 Feb 18 54.944778 354.3239 7.93127604 Mar 18 31.598316 351.4721 -13.17038312 Apr 18 -8.339539 347.6969 9.84264115 May 18 -29.893525 343.9217 7.37179640 Jun 18 41.489279 338.8416 25.36908823 Jul 18 16.349030 333.7615 -7.21056752 Aug 18 -15.202202 327.1189 4.58325589 Sep 18 -33.490109 320.4764 -2.78624584 Oct 18 -49.202316 314.1426 5.95968794 Nov 18 -27.837860 307.8089 8.82895830 Dec 18 -26.542588 303.3976 1.94501620 Jan 19 46.126734 298.9862 -20.71297617 Feb 19 54.944778 296.2399 -40.28466302 Mar 19 31.598316 293.4935 -26.09184347 Apr 19 -8.339539 291.8614 -10.52181610 May 19 -29.893525 290.2292 18.96434226 Jun 19 41.489279 289.7608 27.94988948 Jul 19 16.349030 289.2925 -0.64151090 Aug 19 -15.202202 288.8280 8.47421358 Sep 19 -33.490109 288.3635 -4.57338709 Oct 19 -49.202316 286.9166 8.78571564 Nov 19 -27.837860 285.4697 0.26815493 Dec 19 -26.542588 284.8119 8.23064644 Jan 20 46.126734 284.1542 -14.38091233 Feb 20 54.944778 285.8173 -22.36206769 Mar 20 31.598316 287.4804 -23.67871665 Apr 20 -8.339539 290.6498 -15.91021690 May 20 -29.893525 293.8191 -18.12558617 Jun 20 41.489279 296.3949 24.91583821 Jul 20 16.349030 298.9707 9.58031500 Aug 20 -15.202202 299.6434 9.75875853 Sep 20 -33.490109 300.3162 22.67387692 Oct 20 -49.202316 298.8030 45.59934856 Nov 20 -27.837860 297.2897 20.84815677 Dec 20 -26.542588 295.0713 3.47127847 Jan 21 46.126734 292.8529 -31.57965010 Feb 21 54.944778 290.6407 -16.88543483 Mar 21 31.598316 288.4284 -27.12671315 Apr 21 -8.339539 286.6784 -29.23890347 May 21 -29.893525 284.9285 -24.63496280 Jun 21 41.489279 283.7282 36.28255930 Jul 21 16.349030 282.5278 22.82313381 Aug 21 -15.202202 281.0951 11.30714065 Sep 21 -33.490109 279.6623 14.52782236 Oct 21 -49.202316 277.5662 22.63613321 Nov 21 -27.837860 275.4701 9.96778062 Dec 21 -26.542588 273.8124 -5.46983093 Jan 22 46.126734 272.1548 -30.78149275 Feb 22 54.944778 272.6931 -35.33790159 Mar 22 31.598316 273.2315 -30.12980401 Apr 22 -8.339539 276.2529 -13.71334531 May 22 -29.893525 279.2743 -19.38075562 Jun 22 41.489279 283.8170 13.69368625 Jul 22 16.349030 288.3598 13.49118053 Aug 22 -15.202202 294.4892 7.71297300 Sep 22 -33.490109 300.6187 28.67144033 Oct 22 -49.202316 308.0241 25.17821783 Nov 22 -27.837860 315.4295 -16.59166810 Dec 22 -26.542588 324.5740 -35.33145376 Jan 23 46.126734 333.7186 -39.24528968 Feb 23 54.944778 345.1328 -20.67755395 Mar 23 31.598316 356.5470 -14.84531182 Apr 23 -8.339539 370.9326 -7.39310795 May 23 -29.893525 385.3183 -17.02477309 Jun 23 41.489279 401.5684 23.84235730 Jul 23 16.349030 417.8184 16.83254011 Aug 23 -15.202202 432.1531 5.04905416 Sep 23 -33.490109 446.4879 16.20224307 Oct 23 -49.202316 456.5716 18.53067181 Nov 23 -27.837860 466.6554 21.88243713 Dec 23 -26.542588 473.5796 16.56298609 Jan 24 46.126734 480.5038 14.76948479 Feb 24 54.944778 485.5358 3.71943132 Mar 24 31.598316 490.5678 -4.66611574 Apr 24 -8.339539 493.7508 -16.01121822 May 24 -29.893525 496.9337 -27.64018972 Jun 24 41.489279 498.4682 9.04251757 Jul 24 16.349030 500.0027 16.64827728 Aug 24 -15.202202 500.6722 20.62998386 Sep 24 -33.490109 501.3417 16.14836530 Oct 24 -49.202316 500.8623 5.33997603 Nov 24 -27.837860 500.3829 8.95492333 Dec 24 -26.542588 498.8286 -2.78596524 Jan 25 46.126734 497.2742 1.29909592 Feb 25 54.944778 495.5363 -9.28105617 Mar 25 31.598316 493.7984 -3.89670186 Apr 25 -8.339539 491.8430 -13.80348686 May 25 -29.893525 489.8877 -25.59414087 Jun 25 41.489279 486.2379 14.87280730 Jul 25 16.349030 482.5882 18.36280788 Aug 25 -15.202202 477.4985 23.40366320 Sep 25 -33.490109 472.4089 26.88119338 Oct 25 -49.202316 466.5612 29.64110882 Nov 25 -27.837860 460.7135 -6.27563918 Dec 25 -26.542588 454.6536 -16.51097329 Jan 26 46.126734 448.5936 -27.22035768 Feb 26 54.944778 443.3873 -13.83211720 Mar 26 31.598316 438.1811 -18.57937031 Apr 26 -8.339539 435.2929 -9.55338442 May 26 -29.893525 432.4048 -22.61126754 Jun 26 41.489279 432.5311 10.67964531 Jul 26 16.349030 432.6574 5.99361056 Aug 26 -15.202202 434.5511 1.45109668 Sep 26 -33.490109 436.4449 13.54525766 Oct 26 -49.202316 438.7072 -13.20485318 Nov 26 -27.837860 440.9695 -7.53162744 Dec 26 -26.542588 444.6018 -12.25923274 Jan 27 46.126734 448.2342 6.43911168 Feb 27 54.944778 454.8052 4.25002125 Mar 27 31.598316 461.3762 -17.47456278 Apr 27 -8.339539 472.2963 -33.85679741 May 27 -29.893525 483.2164 -38.92290105 Jun 27 41.489279 501.3156 -4.80489381 Jul 27 16.349030 519.4148 -9.76383416 Aug 27 -15.202202 545.6511 -41.94887324 Sep 27 -33.490109 571.8873 -18.19723746 Oct 27 -49.202316 602.1517 -48.54940134 Nov 27 -27.837860 632.4161 -36.07822865 Dec 27 -26.542588 661.2341 -24.09146926 Jan 28 46.126734 690.0520 81.82123986 Feb 28 54.944778 714.1426 61.81263298 Mar 28 31.598316 738.2332 66.06853250 Apr 28 -8.339539 754.3146 36.02492895 May 28 -29.893525 770.3961 21.79745639 Jun 28 41.489279 776.2548 39.15594409 Jul 28 16.349030 782.1135 22.43748421 Aug 28 -15.202202 779.2891 5.51313636 Sep 28 -33.490109 776.4646 9.22546337 Oct 28 -49.202316 768.9246 4.67771118 Nov 28 -27.837860 761.3846 -10.44670445 Dec 28 -26.542588 753.4157 -7.37309997 Jan 29 46.126734 745.4468 25.82645424 Feb 29 54.944778 740.1011 8.25410725 Mar 29 31.598316 734.7554 -13.85373334 Apr 29 -8.339539 731.9095 -34.57000999 May 29 -29.893525 729.0637 -68.77015566 Jun 29 41.489279 727.9314 -3.92071318 Jul 29 16.349030 726.7992 14.55178172 Aug 29 -15.202202 727.1543 20.24785260 Sep 29 -33.490109 727.5095 8.58059834 Oct 29 -49.202316 726.8860 5.61627019 Nov 29 -27.837860 726.2626 11.07527861 Dec 29 -26.542588 722.8410 5.90158053 Jan 30 46.126734 719.4194 19.25383219 Feb 30 54.944778 714.0567 41.89850646 Mar 30 31.598316 708.6940 15.30768714 Apr 30 -8.339539 702.4843 -37.34471622 May 30 -29.893525 696.2745 -51.28098860 Jun 30 41.489279 688.4492 15.36148434 Jul 30 16.349030 680.6240 -2.87299030 Aug 30 -15.202202 672.0621 18.84008910 Sep 30 -33.490109 663.5003 13.68984335 Oct 30 -49.202316 654.2711 17.03125370 Nov 30 -27.837860 645.0419 17.39600063 Dec 30 -26.542588 635.6919 -21.14935107 Jan 31 46.126734 626.3420 17.23124697 Feb 31 54.944778 616.7476 2.20758542 Mar 31 31.598316 607.1533 9.14843026 Apr 31 -8.339539 597.2322 -20.09270415 May 31 -29.893525 587.3112 -11.71770758 > m$win s t l 3651 19 13 > m$deg s t l 0 1 1 > m$jump s t l 366 2 2 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/wessaorg/rcomp/tmp/1pfwq1352745677.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/2yyqt1352745677.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/3rkw81352745677.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/4ylq01352745677.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/560p91352745677.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/6hcx21352745678.tab") > > try(system("convert tmp/1pfwq1352745677.ps tmp/1pfwq1352745677.png",intern=TRUE)) character(0) > try(system("convert tmp/2yyqt1352745677.ps tmp/2yyqt1352745677.png",intern=TRUE)) character(0) > try(system("convert tmp/3rkw81352745677.ps tmp/3rkw81352745677.png",intern=TRUE)) character(0) > try(system("convert tmp/4ylq01352745677.ps tmp/4ylq01352745677.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 6.953 0.445 7.437