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 = 'TRUE' > par7 = '1' > par6 = '' > par5 = '1' > par4 = '' > par3 = '0' > par2 = 'periodic' > par1 = '12' > main = 'Seasonal Decomposition by Loess' > par8 <- 'TRUE' > 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 45.04690 213.0277 -22.97458729 Feb 1 46.66551 216.1920 17.84252655 Mar 1 23.95988 219.3562 21.28388152 Apr 1 -14.28222 223.0553 31.92686945 May 1 -31.52655 226.7545 6.17209643 Jun 1 43.58563 231.0483 -33.83391323 Jul 1 18.71923 235.3421 -12.96134674 Aug 1 -11.33991 239.8481 -4.70822824 Sep 1 -29.18637 244.3542 -9.06778610 Oct 1 -46.32479 251.2884 -30.26365459 Nov 1 -23.90257 258.2227 -31.02016323 Dec 1 -21.41472 270.9311 -29.01635762 Jan 2 45.04690 283.6394 -29.18631849 Feb 2 46.66551 298.6275 2.10698707 Mar 2 23.95988 313.6156 0.72453376 Apr 2 -14.28222 328.3168 13.66542691 May 2 -31.52655 343.0180 40.10855910 Jun 2 43.58563 355.7826 -2.76825497 Jul 2 18.71923 368.5473 51.53350709 Aug 2 -11.33991 377.7097 29.23018818 Sep 2 -29.18637 386.8722 5.81419290 Oct 2 -46.32479 389.7023 35.42245544 Nov 2 -23.90257 392.5325 -11.62992218 Dec 2 -21.41472 388.2925 2.12225311 Jan 3 45.04690 384.0524 35.70066192 Feb 3 46.66551 374.7560 57.67853919 Mar 3 23.95988 365.4595 41.88065760 Apr 3 -14.28222 353.9544 26.82777444 May 3 -31.52655 342.4494 15.37713031 Jun 3 43.58563 328.8720 -17.35763179 Jul 3 18.71923 315.2946 -2.41381775 Aug 3 -11.33991 299.7185 -27.07858700 Sep 3 -29.18637 284.1424 -5.95603261 Oct 3 -46.32479 270.3016 -18.47680498 Nov 3 -23.90257 256.4608 3.04178251 Dec 3 -21.41472 245.4772 16.83752636 Jan 4 45.04690 234.4936 -14.64049626 Feb 4 46.66551 226.8428 -19.70827144 Mar 4 23.95988 219.1919 -10.85180548 Apr 4 -14.28222 214.6769 -6.59470454 May 4 -31.52655 210.1619 -1.63536455 Jun 4 43.58563 207.3290 -37.71458510 Jul 4 18.71923 204.4960 -16.01522952 Aug 4 -11.33991 202.5800 -10.64006210 Sep 4 -29.18637 200.6639 17.12242895 Oct 4 -46.32479 199.6069 22.11787759 Nov 4 -23.90257 198.5499 24.35268608 Dec 4 -21.41472 197.8143 3.20044137 Jan 5 45.04690 197.0787 -16.32556983 Feb 5 46.66551 196.0461 -8.71156162 Mar 5 23.95988 195.0134 -18.77331227 Apr 5 -14.28222 193.5048 4.37742288 May 5 -31.52655 191.9962 17.73039706 Jun 5 43.58563 190.5101 -30.89569741 Jul 5 18.71923 189.0240 0.75678425 Aug 5 -11.33991 187.2799 15.85996636 Sep 5 -29.18637 185.5359 16.45047210 Oct 5 -46.32479 183.6327 10.69208851 Nov 5 -23.90257 181.7295 1.57306476 Dec 5 -21.41472 179.3252 -3.41049183 Jan 6 45.04690 176.9209 -8.76781489 Feb 6 46.66551 175.1584 -25.42392574 Mar 6 23.95988 173.3959 -14.55579546 Apr 6 -14.28222 175.0075 15.67474777 May 6 -31.52655 176.6190 8.50753004 Jun 6 43.58563 182.2857 -52.67135777 Jul 6 18.71923 187.9524 -35.67166943 Aug 6 -11.33991 196.8662 -34.32632725 Sep 6 -29.18637 205.7800 -14.69366142 Oct 6 -46.32479 219.0239 -15.49910237 Nov 6 -23.90257 232.2678 -6.66518347 Dec 6 -21.41472 246.4931 11.32163933 Jan 7 45.04690 260.7184 50.33469565 Feb 7 46.66551 274.8367 76.79778840 Mar 7 23.95988 288.9550 90.78512229 Apr 7 -14.28222 301.1842 97.69800719 May 7 -31.52655 313.4134 83.91313113 Jun 7 43.58563 324.5957 -0.08128534 Jul 7 18.71923 335.7779 13.40287434 Aug 7 -11.33991 338.1208 20.21906402 Sep 7 -29.18637 340.4638 32.02257734 Oct 7 -46.32479 336.5056 2.71923150 Nov 7 -23.90257 332.5473 2.85524550 Dec 7 -21.41472 326.5305 -4.21577400 Jan 8 45.04690 320.5137 1.33944001 Feb 8 46.66551 312.9647 -2.73022055 Mar 8 23.95988 305.4158 0.32436003 Apr 8 -14.28222 299.0484 31.43379095 May 8 -31.52655 292.6811 7.84546092 Jun 8 43.58563 287.9439 -42.22955358 Jul 8 18.71923 283.2068 -35.72599194 Aug 8 -11.33991 279.4750 -14.53504430 Sep 8 -29.18637 275.7431 -12.75677302 Oct 8 -46.32479 274.2176 0.50724053 Nov 8 -23.90257 272.6920 4.81061392 Dec 8 -21.41472 273.6871 7.82761044 Jan 9 45.04690 274.6823 -13.12915951 Feb 9 46.66551 275.1917 -12.65716745 Mar 9 23.95988 275.7011 9.83906575 Apr 9 -14.28222 275.3036 9.97859697 May 9 -31.52655 274.9062 36.52036723 Jun 9 43.58563 274.7965 -0.48213984 Jul 9 18.71923 274.6868 4.99392924 Aug 9 -11.33991 274.1622 -16.12228467 Sep 9 -29.18637 273.6375 -17.15117495 Oct 9 -46.32479 272.7325 -17.30771394 Nov 9 -23.90257 271.8275 11.97510692 Dec 9 -21.41472 271.6805 15.73423829 Jan 10 45.04690 271.5335 4.01960319 Feb 10 46.66551 272.2977 -10.46316204 Mar 10 23.95988 273.0618 -14.82168614 Apr 10 -14.28222 275.7855 1.19675329 May 10 -31.52655 278.5091 16.51743177 Jun 10 43.58563 284.6428 -15.12842298 Jul 10 18.71923 290.7765 -25.19570159 Aug 10 -11.33991 300.0547 -36.11483724 Sep 10 -29.18637 309.3330 -29.84664924 Oct 10 -46.32479 323.1799 -30.35510906 Nov 10 -23.90257 337.0268 -0.42420903 Dec 10 -21.41472 355.3924 -0.77763012 Jan 11 45.04690 373.7579 27.59518231 Feb 11 46.66551 384.2975 80.63702205 Mar 11 23.95988 394.8370 96.70310293 Apr 11 -14.28222 403.5940 117.08825452 May 11 -31.52655 412.3509 102.37564516 Jun 11 43.58563 422.8214 55.89301661 Jul 11 18.71923 433.2918 57.78896421 Aug 11 -11.33991 438.9305 33.10940506 Sep 11 -29.18637 444.5692 -9.58283045 Oct 11 -46.32479 436.6054 -15.28064056 Nov 11 -23.90257 428.6417 -26.23909083 Dec 11 -21.41472 418.7276 9.48708076 Jan 12 45.04690 408.8136 13.93948587 Feb 12 46.66551 399.8419 23.29256598 Mar 12 23.95988 390.8702 14.96988722 Apr 12 -14.28222 384.8379 -14.75567926 May 12 -31.52655 378.8056 -14.57900669 Jun 12 43.58563 374.3104 -39.89605580 Jul 12 18.71923 369.8153 -28.03452877 Aug 12 -11.33991 366.7939 -20.75396227 Sep 12 -29.18637 363.7724 -15.08607214 Oct 12 -46.32479 363.7412 5.68354543 Nov 12 -23.90257 363.7101 23.79252285 Dec 12 -21.41472 365.7504 7.76435486 Jan 13 45.04690 367.7907 -0.93757960 Feb 13 46.66551 369.7614 -27.82688863 Mar 13 23.95988 371.7321 20.70804347 Apr 13 -14.28222 374.5743 0.40788831 May 13 -31.52655 377.4166 -7.89002780 Jun 13 43.58563 383.6314 -10.01698020 Jul 13 18.71923 389.8461 -20.16535646 Aug 13 -11.33991 400.0503 -17.61041360 Sep 13 -29.18637 410.2545 -49.56814711 Oct 13 -46.32479 421.8923 -21.86746321 Nov 13 -23.90257 433.5300 -12.92741946 Dec 13 -21.41472 443.0275 25.38725002 Jan 14 45.04690 452.5249 35.92815302 Feb 14 46.66551 459.3181 59.41640998 Mar 14 23.95988 466.1112 52.22890807 Apr 14 -14.28222 468.8356 34.14663909 May 14 -31.52655 471.5599 27.06660914 Jun 14 43.58563 469.0412 18.67312297 Jul 14 18.71923 466.5226 10.85821295 Aug 14 -11.33991 458.3784 -3.03851303 Sep 14 -29.18637 450.2343 -17.64791537 Oct 14 -46.32479 441.6410 -9.01624317 Nov 14 -23.90257 433.0478 -15.04521113 Dec 14 -21.41472 425.6316 -0.11690182 Jan 15 45.04690 418.2155 -1.16235899 Feb 15 46.66551 412.3085 -10.87400572 Mar 15 23.95988 406.4015 1.93858869 Apr 15 -14.28222 402.3135 -1.73123806 May 15 -31.52655 398.2254 28.50117425 Jun 15 43.58563 396.4115 -18.09713224 Jul 15 18.71923 394.5976 -30.41686257 Aug 15 -11.33991 395.0001 0.53982641 Sep 15 -29.18637 395.4025 -20.71616096 Oct 15 -46.32479 397.2191 -27.49434100 Nov 15 -23.90257 399.0357 -2.53316119 Dec 15 -21.41472 401.3295 -3.91478976 Jan 16 45.04690 403.6233 14.02981520 Feb 16 46.66551 405.0644 35.27013588 Mar 16 23.95988 406.5054 13.73469771 Apr 16 -14.28222 406.5716 7.01059790 May 16 -31.52655 406.6378 19.78873713 Jun 16 43.58563 405.6443 6.17007341 Jul 16 18.71923 404.6508 -9.37001416 Aug 16 -11.33991 402.8610 -16.02109936 Sep 16 -29.18637 401.0712 -24.88486092 Oct 16 -46.32479 399.3640 -13.63925399 Nov 16 -23.90257 397.6569 12.04571278 Dec 16 -21.41472 396.1334 4.08128158 Jan 17 45.04690 394.6100 12.14308391 Feb 17 46.66551 392.5832 6.85133656 Mar 17 23.95988 390.5563 7.98383034 Apr 17 -14.28222 387.3512 10.03104968 May 17 -31.52655 384.1460 0.18050807 Jun 17 43.58563 380.1420 21.57236356 Jul 17 18.71923 376.1380 -27.35720481 Aug 17 -11.33991 372.2762 -5.83624995 Sep 17 -29.18637 368.4143 -13.02797146 Oct 17 -46.32479 365.3039 0.82093928 Nov 17 -23.90257 362.1934 -6.49079013 Dec 17 -21.41472 359.7385 2.57622130 Jan 18 45.04690 357.2836 -8.23053376 Feb 18 46.66551 354.4489 16.08557155 Mar 18 23.95988 351.6142 -5.67408200 Apr 18 -14.28222 347.7213 15.76092203 May 18 -31.52655 343.8284 9.09816509 Jun 18 43.58563 338.5880 23.52637804 Jul 18 18.71923 333.3476 -9.16683287 Aug 18 -11.33991 326.6782 1.16167314 Sep 18 -29.18637 320.0089 -6.62249722 Oct 18 -46.32479 313.9435 3.28133079 Nov 18 -23.90257 307.8781 4.82451865 Dec 18 -21.41472 303.8804 -3.66562984 Jan 19 45.04690 299.8826 -20.52954481 Feb 19 46.66551 297.2588 -33.02427834 Mar 19 23.95988 294.6349 -19.59477073 Apr 19 -14.28222 292.7451 -5.46286280 May 19 -31.52655 290.8553 19.97128416 Jun 19 43.58563 289.9591 25.65532251 Jul 19 18.71923 289.0628 -2.78206299 Aug 19 -11.33991 288.3595 5.08037163 Sep 19 -29.18637 287.6562 -8.16987010 Oct 19 -46.32479 286.2974 6.52736480 Nov 19 -23.90257 284.9386 -3.13604045 Dec 19 -21.41472 284.6006 3.31413540 Jan 20 45.04690 284.2626 -13.40945524 Feb 20 46.66551 285.9856 -14.25105974 Mar 20 23.95988 287.7085 -16.26842311 Apr 20 -14.28222 290.5591 -9.87689922 May 20 -31.52655 293.4097 -16.08313629 Jun 20 43.58563 295.4145 23.79985235 Jul 20 18.71923 297.4194 8.76141712 Aug 20 -11.33991 297.7472 7.79269071 Sep 20 -29.18637 298.0751 20.61128794 Oct 20 -46.32479 296.7287 44.79611608 Nov 20 -23.90257 295.3823 18.82030407 Dec 20 -21.41472 293.4646 -0.04983650 Jan 21 45.04690 291.5468 -29.19374354 Feb 21 46.66551 289.6470 -7.61253592 Mar 21 23.95988 287.7472 -18.80708717 Apr 21 -14.28222 286.1795 -22.79727593 May 21 -31.52655 284.6118 -22.68522564 Jun 21 43.58563 283.3391 34.57526519 Jul 21 18.71923 282.0664 20.91433216 Aug 21 -11.33991 280.6498 7.89009781 Sep 21 -29.18637 279.2332 10.65318708 Oct 21 -46.32479 277.4820 19.84279086 Nov 21 -23.90257 275.7308 5.77175449 Dec 21 -21.41472 274.6050 -11.39030944 Jan 22 45.04690 273.4792 -31.02613985 Feb 22 46.66551 274.1809 -28.54639198 Mar 22 23.95988 274.8825 -24.14240298 Apr 22 -14.28222 277.6823 -9.20012080 May 22 -31.52655 280.4822 -18.95559958 Jun 22 43.58563 285.0685 10.34590569 Jul 22 18.71923 289.6548 9.82598710 Aug 22 -11.33991 295.9303 2.40960523 Sep 22 -29.18637 302.2058 22.78054700 Oct 22 -46.32479 310.0922 20.23256221 Nov 22 -23.90257 317.9786 -23.07606274 Dec 22 -21.41472 327.4962 -43.38145198 Jan 23 45.04690 337.0137 -41.46060770 Feb 23 46.66551 348.0733 -15.33878317 Mar 23 23.95988 359.1328 -9.79271751 Apr 23 -14.28222 372.7920 -3.30977131 May 23 -31.52655 386.4511 -16.52458606 Jun 23 43.58563 402.1759 21.13843857 Jul 23 18.71923 417.9007 14.38003935 Aug 23 -11.33991 431.8268 1.51309532 Sep 23 -29.18637 445.7529 12.63347492 Oct 23 -46.32479 456.0804 16.14436730 Nov 23 -23.90257 466.4080 18.19461952 Dec 23 -21.41472 473.6815 11.33318951 Jan 24 45.04690 480.9551 15.39799301 Feb 24 46.66551 486.1547 11.37976059 Mar 24 23.95988 491.3544 2.18576932 Apr 24 -14.28222 494.4561 -10.77388070 May 24 -31.52655 497.5578 -26.63129168 Jun 24 43.58563 498.7937 6.62062752 Jul 24 18.71923 500.0296 14.25112286 Aug 24 -11.33991 500.5068 16.93313698 Sep 24 -29.18637 500.9839 12.20247474 Oct 24 -46.32479 500.6608 2.66394387 Nov 24 -23.90257 500.3378 5.06477285 Dec 24 -21.41472 498.9971 -8.08240263 Jan 25 45.04690 497.6564 1.99665542 Feb 25 46.66551 496.0560 -1.52154851 Mar 25 23.95988 494.4556 3.08448871 Apr 25 -14.28222 492.3640 -8.38177509 May 25 -31.52655 490.2724 -24.34579984 Jun 25 43.58563 486.2324 12.78194772 Jul 25 18.71923 482.1925 16.38827142 Aug 25 -11.33991 476.9028 20.13708050 Sep 25 -29.18637 471.6132 23.37321322 Oct 25 -46.32479 466.0051 27.31967932 Nov 25 -23.90257 460.3971 -9.89449474 Dec 25 -21.41472 454.6869 -21.67213741 Jan 26 45.04690 448.9766 -26.52354656 Feb 26 46.66551 443.9002 -6.06572278 Mar 26 23.95988 438.8238 -11.58365786 Apr 26 -14.28222 435.9096 -4.22739038 May 26 -31.52655 432.9954 -21.56888386 Jun 26 43.58563 432.8766 8.23773761 Jul 26 18.71923 432.7578 3.52293521 Aug 26 -11.33991 434.4805 -2.34056992 Sep 26 -29.18637 436.2031 9.48324858 Oct 26 -46.32479 438.7052 -16.08045218 Nov 26 -23.90257 441.2074 -11.70479310 Dec 26 -21.41472 445.4186 -18.20383192 Jan 27 45.04690 449.6297 6.12336277 Feb 27 46.66551 456.6926 10.64188575 Mar 27 23.95988 463.7555 -12.21535015 Apr 27 -14.28222 474.5876 -30.20533705 May 27 -31.52655 485.4196 -39.49308492 Jun 27 43.58563 499.2879 -4.87351492 Jul 27 18.71923 513.1561 -5.87536879 Aug 27 -11.33991 530.5030 -30.66306819 Sep 27 -29.18637 547.8498 1.53655605 Oct 27 -46.32479 570.6788 -19.95404848 Nov 27 -23.90257 593.5079 -1.10529316 Dec 27 -21.41472 617.6441 14.37061747 Jan 28 45.04690 641.7803 131.17276161 Feb 28 46.66551 660.0835 124.15097190 Mar 28 23.95988 678.3867 133.55342332 Apr 28 -14.28222 695.7278 100.55444870 May 28 -31.52655 713.0688 80.75771312 Jun 28 43.58563 737.7454 75.56899016 Jul 28 18.71923 762.4219 39.75884334 Aug 28 -11.33991 769.9746 10.96528837 Sep 28 -29.18637 777.5273 3.85905703 Oct 28 -46.32479 770.2196 0.50520355 Nov 28 -23.90257 762.9119 -15.90929009 Dec 28 -21.41472 756.6653 -15.75060154 Jan 29 45.04690 750.4188 21.93432053 Feb 29 46.66551 746.5482 10.08625056 Mar 29 23.95988 742.6777 -14.13757827 Apr 29 -14.28222 739.9168 -36.63461892 May 29 -31.52655 737.1560 -75.22942053 Jun 29 43.58563 735.2772 -13.36282427 Jul 29 18.71923 733.3984 5.58234812 Aug 29 -11.33991 732.4315 11.10841508 Sep 29 -29.18637 731.4646 0.32180568 Oct 29 -46.32479 729.9689 -0.34410216 Nov 29 -23.90257 728.4732 4.92934984 Dec 29 -21.41472 724.2188 -0.60406544 Jan 30 45.04690 719.9643 19.78875281 Feb 30 46.66551 714.8795 49.35495753 Mar 30 23.95988 709.7947 21.84540340 Apr 30 -14.28222 703.9097 -32.82750214 May 30 -31.52655 698.0247 -51.39816863 Jun 30 43.58563 690.4919 11.22245935 Jul 30 18.71923 682.9591 -7.57833654 Aug 30 -11.33991 674.7082 12.33166559 Sep 30 -29.18637 666.4574 6.42899136 Oct 30 -46.32479 658.3524 10.07239526 Nov 30 -23.90257 650.2474 8.25515900 Dec 30 -21.41472 642.4329 -33.01822428 Jan 31 45.04690 634.6185 10.03462596 Feb 31 46.66551 628.1310 -0.89655356 Mar 31 23.95988 621.6436 2.29650806 Apr 31 -14.28222 617.5013 -34.41911799 May 31 -31.52655 613.3591 -36.13250501 Jun 31 43.58563 609.1390 -20.12459243 Jul 31 18.71923 604.9189 20.16189629 Aug 31 -11.33991 601.0074 3.43247567 Sep 31 -29.18637 597.0960 11.79037870 Oct 31 -46.32479 593.5190 -1.19417595 Nov 31 -23.90257 589.9419 -3.13937076 Dec 31 -21.41472 586.5613 7.35342353 > m$win s t l 3721 19 13 > m$deg s t l 0 1 1 > m$jump s t l 373 2 2 > m$inner [1] 1 > m$outer [1] 15 > postscript(file="/var/fisher/rcomp/tmp/12ub41355787531.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/2e9121355787531.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/32vae1355787531.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/4tc9c1355787531.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/5zd1l1355787531.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/61k971355787531.tab") > > try(system("convert tmp/12ub41355787531.ps tmp/12ub41355787531.png",intern=TRUE)) character(0) > try(system("convert tmp/2e9121355787531.ps tmp/2e9121355787531.png",intern=TRUE)) character(0) > try(system("convert tmp/32vae1355787531.ps tmp/32vae1355787531.png",intern=TRUE)) character(0) > try(system("convert tmp/4tc9c1355787531.ps tmp/4tc9c1355787531.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.526 0.780 8.311