R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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(52.61 + ,65.04 + ,67.54 + ,63.58 + ,57.35 + ,54.93 + ,54.30 + ,58.89 + ,65.95 + ,82.65 + ,100.08 + ,100.68 + ,97.53 + ,92.29 + ,85.08 + ,91.61 + ,93.61 + ,90.40 + ,99.31 + ,107.71 + ,106.18 + ,98.80 + ,99.58 + ,98.85 + ,92.69 + ,91.82 + ,92.63 + ,98.41 + ,94.56 + ,85.78 + ,84.59 + ,83.49 + ,84.68 + ,80.12 + ,84.37 + ,85.94 + ,87.07 + ,84.52 + ,83.13 + ,75.95 + ,70.12 + ,78.10 + ,83.06 + ,87.92 + ,90.21 + ,89.95 + ,97.08 + ,102.08 + ,100.64 + ,97.73 + ,97.61 + ,100.32 + ,102.04 + ,107.80 + ,111.51 + ,110.18 + ,110.08 + ,117.40 + ,119.82 + ,118.79 + ,113.18 + ,122.76 + ,120.43 + ,129.16 + ,132.48 + ,135.68 + ,141.49 + ,122.40 + ,137.06 + ,144.84 + ,154.64 + ,148.04 + ,152.76 + ,172.00 + ,169.03 + ,179.68 + ,190.38 + ,233.23 + ,231.45 + ,244.87 + ,299.12 + ,385.01 + ,381.48 + ,321.56 + ,317.27 + ,323.09 + ,392.72 + ,372.37 + ,386.52 + ,412.83 + ,404.91 + ,406.73 + ,392.41 + ,363.31 + ,357.95 + ,375.10 + ,369.74 + ,386.14 + ,353.40 + ,346.87 + ,362.53 + ,349.87 + ,347.03 + ,332.94 + ,327.48 + ,327.92 + ,308.91 + ,285.71 + ,318.81 + ,284.76 + ,301.04 + ,315.16 + ,388.34 + ,383.37 + ,416.77 + ,423.24 + ,429.90 + ,486.07 + ,394.41 + ,410.93 + ,430.88 + ,447.29 + ,431.65 + ,456.53 + ,452.93 + ,440.90 + ,416.46 + ,451.49 + ,432.00 + ,436.19 + ,428.55 + ,421.40 + ,425.18 + ,437.24 + ,431.92 + ,412.65 + ,419.37 + ,436.40 + ,421.37 + ,423.66 + ,402.45 + ,402.82 + ,400.46 + ,425.73 + ,417.93 + ,403.43 + ,404.96 + ,393.64 + ,399.98 + ,375.93 + ,366.57 + ,353.90 + ,347.51 + ,364.10 + ,328.64 + ,348.01 + ,329.63 + ,350.96 + ,336.16 + ,332.15 + ,349.46 + ,383.64 + ,369.82 + ,345.50 + ,337.80 + ,334.76 + ,338.02 + ,346.74 + ,371.84 + ,375.90 + ,373.31 + ,391.91 + ,374.28 + ,384.69 + ,372.16 + ,371.97 + ,351.76 + ,352.89 + ,330.48 + ,347.70 + ,345.58 + ,360.76 + ,364.40 + ,374.62 + ,369.07 + ,341.80 + ,337.87 + ,336.58 + ,332.66 + ,335.74 + ,321.64 + ,329.38 + ,321.84 + ,324.56 + ,330.90 + ,310.91 + ,318.07 + ,312.36 + ,315.19 + ,332.89 + ,310.67 + ,321.26 + ,316.15 + ,283.87 + ,280.65 + ,280.21 + ,265.93 + ,267.80 + ,278.03 + ,291.86 + ,262.61 + ,264.80 + ,265.67 + ,251.05 + ,256.11 + ,279.75 + ,282.52 + ,288.89 + ,308.46 + ,292.89 + ,280.79 + ,273.61 + ,276.67 + ,277.92 + ,250.28 + ,264.70 + ,268.95 + ,261.69 + ,257.99 + ,251.28 + ,243.14 + ,246.81 + ,224.50 + ,241.25 + ,254.97 + ,261.39 + ,266.67 + ,264.28 + ,270.45 + ,274.97 + ,281.13 + ,300.65 + ,321.12 + ,354.79 + ,318.97 + ,298.71 + ,318.85 + ,327.89 + ,348.19 + ,335.18 + ,332.98 + ,331.04 + ,317.52 + ,325.31 + ,317.59 + ,313.37 + ,313.00 + ,314.77 + ,298.37 + ,311.10 + ,308.79 + ,297.30 + ,293.58 + ,291.35 + ,291.51 + ,289.94 + ,287.07 + ,280.74 + ,294.95 + ,288.98 + ,285.63 + ,294.55 + ,290.67 + ,314.78 + ,306.50 + ,304.48 + ,308.65 + ,307.01 + ,298.59 + ,293.51 + ,294.90 + ,296.14 + ,294.25 + ,291.75 + ,290.49 + ,288.68 + ,310.07 + ,297.45 + ,300.81 + ,301.56 + ,296.89 + ,305.23 + ,298.45 + ,298.75 + ,273.02 + ,266.62 + ,266.06 + ,284.48 + ,275.71 + ,284.19 + ,284.81 + ,267.29 + ,272.95 + ,262.35 + ,246.34 + ,251.03 + ,247.54 + ,254.80 + ,245.08 + ,251.30 + ,261.48 + ,258.85 + ,270.89 + ,257.55 + ,253.08 + ,238.81 + ,241.22 + ,280.75 + ,284.56 + ,289.35 + ,289.56 + ,289.55 + ,305.00 + ,289.22 + ,301.82 + ,293.56 + ,300.59 + ,298.67 + ,311.55 + ,310.08 + ,312.06 + ,309.13 + ,292.31 + ,284.41 + ,290.02 + ,291.52 + ,296.81 + ,315.60 + ,319.63 + ,303.89 + ,300.53 + ,321.84 + ,309.48 + ,307.68 + ,310.53 + ,327.91 + ,343.18 + ,345.48 + ,342.03 + ,349.57 + ,322.50 + ,310.74 + ,318.96 + ,327.53 + ,320.00 + ,320.72 + ,330.86 + ,342.34 + ,322.37 + ,306.86 + ,301.75 + ,307.27 + ,301.30 + ,315.18 + ,342.11 + ,333.18 + ,332.26 + ,332.32 + ,330.00 + ,321.78 + ,318.59 + ,344.78 + ,324.09 + ,322.03 + ,325.32 + ,325.10 + ,335.10 + ,334.66 + ,334.54 + ,341.15 + ,320.47 + ,323.85 + ,328.06 + ,328.93 + ,337.50 + ,335.65 + ,361.05 + ,353.19 + ,352.28 + ,392.53 + ,393.03 + ,420.42 + ,434.91 + ,468.38 + ,466.35 + ,480.93 + ,511.25 + ,508.39 + ,479.80 + ,495.63 + ,487.09 + ,473.06 + ,473.03 + ,487.87 + ,479.28 + ,500.60 + ,502.82 + ,497.13 + ,496.06 + ,489.80 + ,481.66 + ,486.17 + ,492.94 + ,522.45 + ,545.71 + ,533.77 + ,570.26 + ,623.56 + ,639.94 + ,589.13 + ,559.45 + ,569.96 + ,590.43 + ,588.37 + ,565.80 + ,629.69 + ,576.28 + ,641.89 + ,625.70 + ,717.52 + ,749.58 + ,690.29 + ,666.55 + ,689.18 + ,666.24 + ,662.32 + ,665.83 + ,681.23 + ,704.87 + ,783.13 + ,757.97 + ,775.93 + ,812.08 + ,824.40 + ,886.89 + ,984.07 + ,1015.59 + ,897.30 + ,980.37 + ,957.37 + ,968.96 + ,1062.80 + ,1047.67 + ,967.91 + ,1021.58 + ,1014.02 + ,1034.98 + ,1068.80 + ,1038.38 + ,1133.26 + ,1259.55 + ,1207.42 + ,1234.59 + ,1297.03) > par8 = 'FALSE' > par7 = '1' > par6 = '' > par5 = '1' > par4 = '' > par3 = '1' > par2 = 'periodic' > par1 = '12' > main = 'Seasonal Decomposition by Loess' > #'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 -2.7764725 53.30156 2.084908656 Feb 1 1.0828964 55.83058 8.126522171 Mar 1 -1.8182369 58.35960 10.998637850 Apr 1 -4.9768164 61.13599 7.420822759 May 1 0.1797171 63.91239 -6.742105305 Jun 1 1.8639596 66.80133 -13.735286286 Jul 1 -2.3418062 69.69027 -13.048458944 Aug 1 4.8502087 72.42217 -18.382375804 Sep 1 -1.2633427 75.15407 -7.940726352 Oct 1 4.9867812 78.07773 -0.414515963 Nov 1 2.7726846 81.00140 16.305914877 Dec 1 -2.5595726 84.35359 18.885986151 Jan 2 -2.7764725 87.70577 12.600700129 Feb 2 1.0828964 90.57250 0.634603391 Mar 2 -1.8182369 93.43923 -6.540991183 Apr 2 -4.9768164 94.81385 1.772962206 May 2 0.1797171 96.18848 -2.758197378 Jun 2 1.8639596 96.47171 -7.935668333 Jul 2 -2.3418062 96.75494 4.896869036 Aug 2 4.8502087 96.96508 5.894708203 Sep 2 -1.2633427 97.17523 10.268113682 Oct 2 4.9867812 97.32790 -3.514679756 Nov 2 2.7726846 97.48057 -0.673252743 Dec 2 -2.5595726 96.87531 4.534260684 Jan 3 -2.7764725 96.27006 -0.803583183 Feb 3 1.0828964 94.85346 -4.116353461 Mar 3 -1.8182369 93.43686 1.011378424 Apr 3 -4.9768164 91.85523 11.531585171 May 3 0.1797171 90.27360 4.106678945 Jun 3 1.8639596 89.06776 -5.151719049 Jul 3 -2.3418062 87.86191 -0.930108720 Aug 3 4.8502087 86.87673 -8.236936766 Sep 3 -1.2633427 85.89154 0.051801499 Oct 3 4.9867812 84.82097 -9.687754443 Nov 3 2.7726846 83.75041 -2.153089934 Dec 3 -2.5595726 83.05267 5.446906145 Jan 4 -2.7764725 82.35493 7.491544929 Feb 4 1.0828964 82.36672 1.070388457 Mar 4 -1.8182369 82.37850 2.569734150 Apr 4 -4.9768164 82.91902 -1.992204918 May 4 0.1797171 83.45954 -13.519256960 Jun 4 1.8639596 84.50474 -8.268699225 Jul 4 -2.3418062 85.54994 -0.148133167 Aug 4 4.8502087 87.13952 -4.069731181 Sep 4 -1.2633427 88.72911 2.744237116 Oct 4 4.9867812 90.86330 -5.900080836 Nov 4 2.7726846 92.99749 1.309821663 Dec 4 -2.5595726 95.30509 9.334485882 Jan 5 -2.7764725 97.61268 5.803792807 Feb 5 1.0828964 99.62059 -2.973490944 Mar 5 -1.8182369 101.62851 -2.200272529 Apr 5 -4.9768164 103.39003 1.906788765 May 5 0.1797171 105.15155 -3.291262913 Jun 5 1.8639596 106.77630 -0.840259289 Jul 5 -2.3418062 108.40105 5.450752658 Aug 5 4.8502087 110.21132 -4.881526917 Sep 5 -1.2633427 112.02158 -0.678240180 Oct 5 4.9867812 114.18394 -1.770722079 Nov 5 2.7726846 116.34630 0.701016473 Dec 5 -2.5595726 118.65965 2.689923175 Jan 6 -2.7764725 120.97300 -5.016527419 Feb 6 1.0828964 123.10218 -1.425077765 Mar 6 -1.8182369 125.23136 -2.983125947 Apr 6 -4.9768164 127.50192 6.634896279 May 6 0.1797171 129.77248 2.527805532 Jun 6 1.8639596 132.53780 1.278236078 Jul 6 -2.3418062 135.30313 8.528674949 Aug 6 4.8502087 138.81872 -21.268923897 Sep 6 -1.2633427 142.33430 -4.010956431 Oct 6 4.9867812 146.97432 -7.121105827 Nov 6 2.7726846 151.61435 0.252965229 Dec 6 -2.5595726 158.57543 -7.975856470 Jan 7 -2.7764725 165.53651 -10.000035463 Feb 7 1.0828964 176.37718 -5.460080212 Mar 7 -1.8182369 187.21786 -16.369622797 Apr 7 -4.9768164 202.90772 -18.250901046 May 7 0.1797171 218.59758 -28.397292267 Jun 7 1.8639596 235.61344 -4.247397824 Jul 7 -2.3418062 252.62930 -18.837495057 Aug 7 4.8502087 269.12283 -29.103036901 Sep 7 -1.2633427 285.61636 14.766987567 Oct 7 4.9867812 301.70589 78.317329160 Nov 7 2.7726846 317.79542 60.911891205 Dec 7 -2.5595726 332.31399 -8.194418893 Jan 8 -2.7764725 346.83256 -26.786086286 Feb 8 1.0828964 356.72872 -34.721611766 Mar 8 -1.8182369 366.62487 27.913364918 Apr 8 -4.9768164 370.99534 6.351479398 May 8 0.1797171 375.36580 10.974480905 Jun 8 1.8639596 377.85874 33.107298552 Jul 8 -2.3418062 380.35168 26.900124523 Aug 8 4.8502087 381.32236 20.557431865 Sep 8 -1.2633427 382.29304 11.380305518 Oct 8 4.9867812 379.95910 -21.635879808 Nov 8 2.7726846 377.62516 -22.447844683 Dec 8 -2.5595726 373.45942 4.200151453 Jan 9 -2.7764725 369.29368 3.222790295 Feb 9 1.0828964 364.74504 20.312058973 Mar 9 -1.8182369 360.19641 -4.978170185 Apr 9 -4.9768164 355.51447 -3.667655041 May 9 0.1797171 350.83254 11.517747131 Jun 9 1.8639596 345.18446 2.821579726 Jul 9 -2.3418062 339.53639 9.835420645 Aug 9 4.8502087 333.67440 -5.584606164 Sep 9 -1.2633427 327.81241 0.930933339 Oct 9 4.9867812 325.03163 -2.098411639 Nov 9 2.7726846 322.25085 -16.113536165 Dec 9 -2.5595726 324.58700 -36.317429234 Jan 10 -2.7764725 326.92315 -5.336679597 Feb 10 1.0828964 334.63725 -50.960151320 Mar 10 -1.8182369 342.35136 -39.493120880 Apr 10 -4.9768164 353.31262 -33.175799848 May 10 0.1797171 364.27387 23.886408210 Jun 10 1.8639596 376.05439 5.451651454 Jul 10 -2.3418062 387.83490 31.276903021 Aug 10 4.8502087 398.84185 19.547942128 Sep 10 -1.2633427 409.84880 21.314547545 Oct 10 4.9867812 418.20921 62.874012444 Nov 10 2.7726846 426.56962 -34.932302206 Dec 10 -2.5595726 430.52155 -17.031979557 Jan 11 -2.7764725 434.47349 -0.817014203 Feb 11 1.0828964 435.61441 10.592690050 Mar 11 -1.8182369 436.75534 -3.287103532 Apr 11 -4.9768164 436.89905 24.607768010 May 11 0.1797171 437.04276 15.707526580 Jun 11 1.8639596 436.98130 2.054745102 Jul 11 -2.3418062 436.91983 -18.118028052 Aug 11 4.8502087 435.80372 10.836072228 Sep 11 -1.2633427 434.68760 -1.424261180 Oct 11 4.9867812 433.03383 -1.830615406 Nov 11 2.7726846 431.38006 -5.602749181 Dec 11 -2.5595726 430.08699 -6.127417715 Jan 12 -2.7764725 428.79392 -0.837443543 Feb 12 1.0828964 427.41110 8.746001107 Mar 12 -1.8182369 426.02829 7.709947921 Apr 12 -4.9768164 424.08814 -6.461328263 May 12 0.1797171 422.14800 -2.957717421 Jun 12 1.8639596 420.43719 14.098853893 Jul 12 -2.3418062 418.72637 4.985433531 Aug 12 4.8502087 416.98447 1.825320274 Sep 12 -1.2633427 415.24257 -11.529226671 Oct 12 4.9867812 413.07861 -15.245386916 Nov 12 2.7726846 410.91464 -13.227326710 Dec 12 -2.5595726 407.62397 20.665604611 Jan 13 -2.7764725 404.33329 16.373178638 Feb 13 1.0828964 399.90794 2.439163087 Mar 13 -1.8182369 395.48259 11.295649701 Apr 13 -4.9768164 390.02927 8.587545679 May 13 0.1797171 384.57595 15.224328684 Jun 13 1.8639596 378.15734 -4.091299706 Jul 13 -2.3418062 371.73873 -2.826919773 Aug 13 4.8502087 365.63090 -16.581113004 Sep 13 -1.2633427 359.52308 -10.749739923 Oct 13 4.9867812 355.19652 3.916700509 Nov 13 2.7726846 350.86995 -25.002638608 Dec 13 -2.5595726 349.38332 1.186253898 Jan 14 -2.7764725 347.89668 -15.490210892 Feb 14 1.0828964 347.85847 2.018635859 Mar 14 -1.8182369 347.82025 -9.842015225 Apr 14 -4.9768164 347.70679 -10.579971922 May 14 0.1797171 347.59332 1.686958409 Jun 14 1.8639596 348.27395 33.502086871 Jul 14 -2.3418062 348.95458 23.207223657 Aug 14 4.8502087 350.97220 -10.322405730 Sep 14 -1.2633427 352.98981 -13.926468806 Oct 14 4.9867812 355.50492 -25.731696865 Nov 14 2.7726846 358.02002 -22.772704473 Dec 14 -2.5595726 360.22243 -10.922857732 Jan 15 -2.7764725 362.42484 12.191631714 Feb 15 1.0828964 364.54321 10.273890311 Mar 15 -1.8182369 366.66159 8.466651073 Apr 15 -4.9768164 367.55475 29.332070164 May 15 0.1797171 368.44791 5.652376281 Jun 15 1.8639596 367.20632 15.619717375 Jul 15 -2.3418062 365.96474 8.537066793 Aug 15 4.8502087 363.86421 3.255582483 Sep 15 -1.2633427 361.76368 -8.740335516 Oct 15 4.9867812 360.12607 -12.222855083 Nov 15 2.7726846 358.48847 -30.781154199 Dec 15 -2.5595726 356.88848 -6.628908938 Jan 16 -2.7764725 355.28849 -6.932020973 Feb 16 1.0828964 353.71417 5.962933268 Mar 16 -1.8182369 352.13985 14.078389673 Apr 16 -4.9768164 350.61100 28.985812352 May 16 0.1797171 349.08216 19.808122059 Jun 16 1.8639596 346.82899 -6.892952629 Jul 16 -2.3418062 344.57583 -4.364018993 Aug 16 4.8502087 341.22873 -9.498940535 Sep 16 -1.2633427 337.88164 -3.958295766 Oct 16 4.9867812 334.31151 -3.558294184 Nov 16 2.7726846 330.74139 -11.874072150 Dec 16 -2.5595726 328.19137 3.748206253 Jan 17 -2.7764725 325.64135 -1.024872639 Feb 17 1.0828964 324.17583 -0.698730173 Mar 17 -1.8182369 322.71032 10.007914457 Apr 17 -4.9768164 321.35009 -5.463277242 May 17 0.1797171 319.98986 -2.099581913 Jun 17 1.8639596 317.58969 -7.093650077 Jul 17 -2.3418062 315.18952 2.342290082 Aug 17 4.8502087 311.48617 16.553618459 Sep 17 -1.2633427 307.78283 4.150513147 Oct 17 4.9867812 303.72914 12.544074540 Nov 17 2.7726846 299.67546 13.701856384 Dec 17 -2.5595726 295.64652 -9.216943574 Jan 18 -2.7764725 291.61757 -8.191100826 Feb 18 1.0828964 287.13780 -8.010693097 Mar 18 -1.8182369 282.65802 -14.909783202 Apr 18 -4.9768164 278.52552 -5.748699815 May 18 0.1797171 274.39301 3.457270599 Jun 18 1.8639596 272.47124 17.524795969 Jul 18 -2.3418062 270.54948 -5.597670338 Aug 18 4.8502087 271.21492 -11.265133060 Sep 18 -1.2633427 271.88037 -4.947029470 Oct 18 4.9867812 273.42526 -27.362045390 Nov 18 2.7726846 274.97016 -21.632840859 Dec 18 -2.5595726 276.19009 6.119487229 Jan 19 -2.7764725 277.41001 7.886458021 Feb 19 1.0828964 278.26339 9.543711865 Mar 19 -1.8182369 279.11677 31.161467873 Apr 19 -4.9768164 279.03807 18.828748212 May 19 0.1797171 278.95937 1.650915577 Jun 19 1.8639596 277.18473 -5.438694067 Jul 19 -2.3418062 275.41010 3.601704612 Aug 19 4.8502087 272.02716 1.042631573 Sep 19 -1.2633427 268.64422 -17.100875154 Oct 19 4.9867812 264.76327 -5.050049990 Nov 19 2.7726846 260.88232 5.294995625 Dec 19 -2.5595726 257.99022 6.259349885 Jan 20 -2.7764725 255.09813 5.668346850 Feb 20 1.0828964 253.78046 -3.583356913 Mar 20 -1.8182369 252.46280 -7.504558512 Apr 20 -4.9768164 252.53616 -0.749344939 May 20 0.1797171 252.60953 -28.289244339 Jun 20 1.8639596 254.10878 -14.722737702 Jul 20 -2.3418062 255.60803 1.703777258 Aug 20 4.8502087 260.03010 -3.490311045 Sep 20 -1.2633427 264.45218 3.481166964 Oct 20 4.9867812 271.43172 -12.138497129 Nov 20 2.7726846 278.41126 -10.733940770 Dec 20 -2.5595726 285.30194 -7.772371518 Jan 21 -2.7764725 292.19263 -8.286159560 Feb 21 1.0828964 298.35361 1.213493280 Mar 21 -1.8182369 304.51459 18.423648285 Apr 21 -4.9768164 310.21051 49.556302917 May 21 0.1797171 315.90644 2.883844576 Jun 21 1.8639596 319.85632 -23.010282750 Jul 21 -2.3418062 323.80621 -2.614401752 Aug 21 4.8502087 325.06469 -2.024897831 Sep 21 -1.2633427 326.32317 23.130172402 Oct 21 4.9867812 325.89875 4.294468230 Nov 21 2.7726846 325.47433 4.732984511 Dec 21 -2.5595726 324.53752 9.062047756 Jan 22 -2.7764725 323.60072 -3.304246294 Feb 22 1.0828964 321.54139 2.685713360 Mar 22 -1.8182369 319.48206 -0.073824822 Apr 22 -4.9768164 316.48405 1.862770366 May 22 0.1797171 313.48603 -0.665747419 Jun 22 1.8639596 310.50811 2.397934898 Jul 22 -2.3418062 307.53018 -6.818374462 Aug 22 4.8502087 304.76175 1.488042223 Sep 22 -1.2633427 301.99332 8.060025219 Oct 22 4.9867812 299.64709 -7.333875784 Nov 22 2.7726846 297.30087 -6.493556336 Dec 22 -2.5595726 295.49721 -1.587640563 Jan 23 -2.7764725 293.69355 0.592917915 Feb 23 1.0828964 292.66471 -3.807606780 Mar 23 -1.8182369 291.63587 -2.747629311 Apr 23 -4.9768164 291.80406 -6.087240232 May 23 0.1797171 291.97225 2.798035874 Jun 23 1.8639596 293.19231 -6.076271794 Jul 23 -2.3418062 294.41238 -6.440571138 Aug 23 4.8502087 295.87307 -6.173279364 Sep 23 -1.2633427 297.33376 -5.400421278 Oct 23 4.9867812 298.46583 11.327386391 Nov 23 2.7726846 299.59790 4.129414512 Dec 23 -2.5595726 300.13010 6.909471436 Jan 24 -2.7764725 300.66230 10.764171064 Feb 24 1.0828964 300.35647 5.570631780 Mar 24 -1.8182369 300.05064 0.357594660 Apr 24 -4.9768164 299.09046 -0.603640011 May 24 0.1797171 298.13027 -3.409987654 Jun 24 1.8639596 297.36006 -3.084022948 Jul 24 -2.3418062 296.58986 0.001950082 Aug 24 4.8502087 296.45118 -9.551388047 Sep 24 -1.2633427 296.31250 -4.559159864 Oct 24 4.9867812 296.74001 -13.046795308 Nov 24 2.7726846 297.16753 10.129789698 Dec 24 -2.5595726 297.15983 2.849745661 Jan 25 -2.7764725 297.15213 6.434344329 Feb 25 1.0828964 295.89420 4.582899989 Mar 25 -1.8182369 294.63628 4.071957813 Apr 25 -4.9768164 292.56418 17.642636152 May 25 0.1797171 290.49208 7.778201517 Jun 25 1.8639596 288.25370 8.632343481 Jul 25 -2.3418062 286.01531 -10.653506231 Aug 25 4.8502087 283.63157 -21.861778677 Sep 25 -1.2633427 281.24783 -13.924484811 Oct 25 4.9867812 278.45006 1.043159084 Nov 25 2.7726846 275.65229 -2.714976570 Dec 25 -2.5595726 272.98796 13.761615728 Jan 26 -2.7764725 270.32362 17.262850730 Feb 26 1.0828964 268.02616 -1.819053343 Mar 26 -1.8182369 265.72869 9.039544747 Apr 26 -4.9768164 263.31813 4.008682587 May 26 0.1797171 260.90758 -14.747292546 Jun 26 1.8639596 258.98439 -9.818348226 Jul 26 -2.3418062 257.06120 -7.179395583 Aug 26 4.8502087 256.00468 -6.054890070 Sep 26 -1.2633427 254.94816 -8.604818245 Oct 26 4.9867812 254.64087 -8.327652136 Nov 26 2.7726846 254.33358 4.373734425 Dec 26 -2.5595726 255.59494 5.814630514 Jan 27 -2.7764725 256.85630 16.810169309 Feb 27 1.0828964 259.52313 -3.056029047 Mar 27 -1.8182369 262.18996 -7.291725239 Apr 27 -4.9768164 265.40732 -21.620506726 May 27 0.1797171 268.62468 -27.584401186 Jun 27 1.8639596 271.97965 6.906388026 Jul 27 -2.3418062 275.33462 11.567185561 Aug 27 4.8502087 279.24628 5.253508174 Sep 27 -1.2633427 283.15795 7.665397098 Oct 27 4.9867812 287.31150 -2.748285328 Nov 27 2.7726846 291.46506 10.762252696 Dec 27 -2.5595726 294.61045 -2.830875788 Jan 28 -2.7764725 297.75583 6.840638432 Feb 28 1.0828964 299.30007 -6.822965926 Mar 28 -1.8182369 300.84430 1.563931880 Apr 28 -4.9768164 300.98347 2.663348754 May 28 0.1797171 301.12263 10.247652656 Jun 28 1.8639596 300.75479 7.461254688 Jul 28 -2.3418062 300.38694 14.014865043 Aug 28 4.8502087 300.59097 3.688820352 Sep 28 -1.2633427 300.79500 -7.221658027 Oct 28 4.9867812 301.16187 -21.738647776 Nov 28 2.7726846 301.52873 -14.281417073 Dec 28 -2.5595726 302.04504 -7.965469877 Jan 29 -2.7764725 302.56135 -2.974879976 Feb 29 1.0828964 303.86255 10.654551194 Mar 29 -1.8182369 305.16375 16.284484528 Apr 29 -4.9768164 307.85253 1.014283369 May 29 0.1797171 310.54131 -10.191030763 Jun 29 1.8639596 313.87836 6.097679811 Jul 29 -2.3418062 317.21541 -5.393601292 Aug 29 4.8502087 319.96280 -17.133005527 Sep 29 -1.2633427 322.71019 -10.916843451 Oct 29 4.9867812 324.42974 -1.506523268 Nov 29 2.7726846 326.14930 14.258017367 Dec 29 -2.5595726 327.26385 20.775724915 Jan 30 -2.7764725 328.37840 16.428075168 Feb 30 1.0828964 328.94619 19.540915659 Mar 30 -1.8182369 329.51398 -5.195741686 Apr 30 -4.9768164 328.95341 -13.236592667 May 30 0.1797171 328.39284 -9.612556621 Jun 30 1.8639596 326.34263 -0.676591246 Jul 30 -2.3418062 324.29242 -1.950617548 Aug 30 4.8502087 322.19052 -6.320724223 Sep 30 -1.2633427 320.08861 12.034735413 Oct 30 4.9867812 319.47367 17.879549783 Nov 30 2.7726846 318.85873 0.738584604 Dec 30 -2.5595726 319.29508 -9.875505202 Jan 31 -2.7764725 319.73142 -15.204952303 Feb 31 1.0828964 320.26562 -14.078520565 Mar 31 -1.8182369 320.79982 -17.681586663 Apr 31 -4.9768164 321.35943 -1.202611157 May 31 0.1797171 321.91903 20.011251377 Jun 31 1.8639596 323.36724 7.948801423 Jul 31 -2.3418062 324.81545 9.786359792 Aug 31 4.8502087 326.18256 1.287234862 Sep 31 -1.2633427 327.54967 3.713676243 Oct 31 4.9867812 328.01591 -11.222695668 Nov 31 2.7726846 328.48216 -12.664847127 Dec 31 -2.5595726 328.61597 18.723603760 Jan 32 -2.7764725 328.74978 -1.883302647 Feb 32 1.0828964 329.08517 -8.138063980 Mar 32 -1.8182369 329.42056 -2.282323149 Apr 32 -4.9768164 329.52146 0.555361122 May 32 0.1797171 329.62235 5.297932420 Jun 32 1.8639596 329.76363 3.032414389 Jul 32 -2.3418062 329.90490 6.976904682 Aug 32 4.8502087 330.95162 5.348167221 Sep 32 -1.2633427 331.99835 -10.265003929 Oct 32 4.9867812 334.12900 -15.265781778 Nov 32 2.7726846 336.25965 -10.972339175 Dec 32 -2.5595726 340.26239 -8.772812764 Jan 33 -2.7764725 344.26512 -3.988643648 Feb 33 1.0828964 351.49820 -16.931094936 Mar 33 -1.8182369 358.73128 4.136955940 Apr 33 -4.9768164 369.42035 -11.253532995 May 33 0.1797171 380.10942 -28.009134903 Jun 33 1.8639596 393.16077 -2.494729769 Jul 33 -2.3418062 406.21212 -10.840316312 Aug 33 4.8502087 419.76183 -4.192040874 Sep 33 -1.2633427 433.31154 2.861800875 Oct 33 4.9867812 445.18260 18.210616345 Nov 33 2.7726846 457.05366 6.523652267 Dec 33 -2.5595726 465.50208 17.987494537 Jan 34 -2.7764725 473.95049 40.075979514 Feb 34 1.0828964 478.72532 28.581782715 Mar 34 -1.8182369 483.50015 -1.881911919 Apr 34 -4.9768164 485.63879 14.968022703 May 34 0.1797171 487.77744 -0.867155648 Jun 34 1.8639596 488.31667 -17.120628246 Jul 34 -2.3418062 488.85590 -13.484092520 Aug 34 4.8502087 488.68522 -5.665432362 Sep 34 -1.2633427 488.51455 -7.971205892 Oct 34 4.9867812 489.18632 6.426900696 Nov 34 2.7726846 489.85809 10.189227736 Dec 34 -2.5595726 492.57842 7.111149302 Jan 35 -2.7764725 495.29876 3.537713574 Feb 35 1.0828964 500.53268 -11.815572545 Mar 35 -1.8182369 505.76659 -22.288356499 Apr 35 -4.9768164 514.08234 -22.935526506 May 35 0.1797171 522.39809 -29.637809485 Jun 35 1.8639596 531.78134 -11.195304148 Jul 35 -2.3418062 541.16460 6.887209512 Aug 35 4.8502087 550.05973 -21.139934889 Sep 35 -1.2633427 558.95486 12.568487022 Oct 35 4.9867812 566.64707 51.926151185 Nov 35 2.7726846 574.33928 62.828035800 Dec 35 -2.5595726 580.24601 11.443562111 Jan 36 -2.7764725 586.15274 -23.926268874 Feb 36 1.0828964 590.86768 -21.990581082 Mar 36 -1.8182369 595.58263 -3.334391125 Apr 36 -4.9768164 602.42934 -9.082528396 May 36 0.1797171 609.27606 -43.655778640 Jun 36 1.8639596 618.89021 8.935830781 Jul 36 -2.3418062 628.50436 -49.882551475 Aug 36 4.8502087 638.49618 -1.456385787 Sep 36 -1.2633427 648.48800 -21.524653788 Oct 36 4.9867812 656.41411 56.119111508 Nov 36 2.7726846 664.34022 82.467097255 Dec 36 -2.5595726 671.45479 21.394781833 Jan 37 -2.7764725 678.56936 -9.242890883 Feb 37 1.0828964 685.66633 2.430774498 Mar 37 -1.8182369 692.76329 -24.705057956 Apr 37 -4.9768164 700.41563 -33.118810188 May 37 0.1797171 708.06796 -42.417675392 Jun 37 1.8639596 720.97931 -41.613266355 Jul 37 -2.3418062 733.89066 -26.678848994 Aug 37 4.8502087 755.83602 22.443773857 Sep 37 -1.2633427 777.78138 -18.548036981 Oct 37 4.9867812 803.44121 -32.497994051 Nov 37 2.7726846 829.10105 -19.793730670 Dec 37 -2.5595726 853.73298 -26.773412316 Jan 38 -2.7764725 878.36492 11.301548743 Feb 38 1.0828964 901.09466 81.892444311 Mar 38 -1.8182369 923.82439 93.583842043 Apr 38 -4.9768164 942.77038 -40.493560983 May 38 0.1797171 961.71636 18.473923018 Jun 38 1.8639596 974.85018 -19.344142159 Jul 38 -2.3418062 987.98401 -16.682199012 Aug 38 4.8502087 997.53421 60.415576722 Sep 38 -1.2633427 1007.08442 41.848918768 Oct 38 4.9867812 1020.95581 -58.032591021 Nov 38 2.7726846 1034.82720 -16.019880360 Dec 38 -2.5595726 1058.03462 -41.455045432 Jan 39 -2.7764725 1081.24204 -43.485567799 Feb 39 1.0828964 1104.78615 -37.069050875 Mar 39 -1.8182369 1128.33027 -88.132031787 Apr 39 -4.9768164 1152.53872 -14.301907518 May 39 0.1797171 1176.74718 82.623103778 Jun 39 1.8639596 1202.20298 3.353062598 Jul 39 -2.3418062 1227.65878 9.273029741 Aug 39 4.8502087 1254.06322 38.116573704 > m$win s t l 4641 19 13 > m$deg s t l 0 1 1 > m$jump s t l 465 2 2 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/wessaorg/rcomp/tmp/1ipwg1324472358.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/2rb0x1324472358.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/3qcr51324472358.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/48fm61324472358.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/51mwk1324472358.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/6izcw1324472358.tab") > > try(system("convert tmp/1ipwg1324472358.ps tmp/1ipwg1324472358.png",intern=TRUE)) character(0) > try(system("convert tmp/2rb0x1324472358.ps tmp/2rb0x1324472358.png",intern=TRUE)) character(0) > try(system("convert tmp/3qcr51324472358.ps tmp/3qcr51324472358.png",intern=TRUE)) character(0) > try(system("convert tmp/48fm61324472358.ps tmp/48fm61324472358.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 6.826 0.303 7.130