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 = '13' > par5 = '1' > par4 = '20' > par3 = '0' > par2 = 'periodic' > par1 = '12' > main = 'Seasonal Decomposition by Loess' > par8 <- 'FALSE' > par7 <- '1' > par6 <- '13' > par5 <- '1' > par4 <- '20' > 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.944015 203.2626 -14.10665193 Feb 1 54.644196 208.1086 17.94718624 Mar 1 31.179871 212.9546 20.46553134 Apr 1 -8.882152 218.3216 31.26055679 May 1 -30.560306 223.6886 8.27171311 Jun 1 41.095776 229.6285 -29.92424812 Jul 1 17.255063 235.5684 -11.72341472 Aug 1 -14.806234 241.9258 -3.31954703 Sep 1 -32.828831 248.2832 -9.35437912 Oct 1 -49.056164 254.6898 -30.93361887 Nov 1 -27.770598 261.0964 -30.02575709 Dec 1 -26.214634 273.6025 -26.88788721 Jan 2 45.944015 286.1087 -32.55270188 Feb 2 54.644196 300.6612 -7.90542518 Mar 2 31.179871 315.2138 -8.09364155 Apr 2 -8.882152 329.4898 7.09238198 May 2 -30.560306 343.7658 38.39453639 Jun 2 41.095776 356.6935 -1.18929449 Jul 2 17.255063 369.6213 51.92366926 Aug 2 -14.806234 378.8084 31.59778649 Sep 2 -32.828831 387.9956 8.33320394 Oct 2 -49.056164 390.9329 36.92321490 Nov 2 -27.770598 393.8703 -9.09967262 Dec 2 -26.214634 389.9814 5.23327614 Jan 3 45.944015 386.0924 32.76354036 Feb 3 54.644196 377.1303 47.32550294 Mar 3 31.179871 368.1682 31.95197246 Apr 3 -8.882152 356.4815 18.90063498 May 3 -30.560306 344.7949 12.06542838 Jun 3 41.095776 330.9887 -16.98447850 Jul 3 17.255063 317.1825 -2.83759074 Aug 3 -14.806234 301.9967 -25.89042850 Sep 3 -32.828831 286.8108 -4.98196604 Oct 3 -49.056164 272.4747 -17.91851669 Nov 3 -27.770598 258.1386 5.23203420 Dec 3 -26.214634 246.8836 20.23098853 Jan 4 45.944015 235.6287 -16.67274167 Feb 4 54.644196 228.1506 -28.99477896 Mar 4 31.179871 220.6724 -19.55230931 Apr 4 -8.882152 216.0390 -13.35688900 May 4 -30.560306 211.4056 -3.84533782 Jun 4 41.095776 208.0516 -35.94733020 Jul 4 17.255063 204.6975 -14.75252794 Aug 4 -14.806234 202.4308 -7.02453930 Sep 4 -32.828831 200.1641 21.26474955 Oct 4 -49.056164 198.9884 25.46780426 Nov 4 -27.770598 197.8126 28.95796050 Dec 4 -26.214634 197.1469 8.66769231 Jan 5 45.944015 196.4812 -16.62526041 Feb 5 54.644196 195.7848 -16.42898146 Mar 5 31.179871 195.0883 -26.06819557 Apr 5 -8.882152 193.5431 -1.06099161 May 5 -30.560306 191.9980 16.76234323 Jun 5 41.095776 189.9114 -27.80720745 Jul 5 17.255063 187.8249 3.42003650 Aug 5 -14.806234 186.0539 20.55231128 Sep 5 -32.828831 184.2829 21.34588627 Oct 5 -49.056164 182.3768 14.67933784 Nov 5 -27.770598 180.4707 6.69989093 Dec 5 -26.214634 178.0501 2.66451980 Jan 6 45.944015 175.6295 -8.37353588 Feb 6 54.644196 174.1896 -32.43381782 Mar 6 31.179871 172.7497 -21.12959282 Apr 6 -8.882152 174.5328 10.74934386 May 6 -30.560306 176.3159 7.84441142 Jun 6 41.095776 184.1009 -51.99666414 Jul 6 17.255063 191.8859 -38.14094506 Aug 6 -14.806234 205.8324 -39.82615606 Sep 6 -32.828831 219.7789 -25.05006684 Oct 6 -49.056164 236.7648 -30.50866145 Nov 6 -27.770598 253.7508 -24.28015454 Dec 6 -26.214634 270.6925 -8.07784058 Jan 7 45.944015 287.6342 22.52178883 Feb 7 54.644196 302.9161 40.73972842 Mar 7 31.179871 318.1980 54.32217495 Apr 7 -8.882152 329.0312 64.45092362 May 7 -30.560306 339.8645 56.49580316 Jun 7 41.095776 344.4189 -17.41466606 Jul 7 17.255063 348.9733 1.67165936 Aug 7 -14.806234 347.4143 14.39198381 Sep 7 -32.828831 345.8552 30.27360847 Oct 7 -49.056164 340.1207 1.83549841 Nov 7 -27.770598 334.3861 4.88448988 Dec 7 -26.214634 326.9102 0.20448000 Jan 8 45.944015 319.4342 1.52178558 Feb 8 54.644196 312.1399 -9.88406763 Mar 8 31.179871 304.8455 -6.32541390 Apr 8 -8.882152 298.5368 26.54539462 May 8 -30.560306 292.2280 7.33233402 Jun 8 41.095776 287.4468 -39.24260249 Jul 8 17.255063 282.6657 -33.72074438 Aug 8 -14.806234 279.4742 -11.06794343 Sep 8 -32.828831 276.2827 -9.65384227 Oct 8 -49.056164 275.0820 2.37415959 Nov 8 -27.770598 273.8813 7.48926298 Dec 8 -26.214634 274.1186 12.19604256 Jan 9 45.944015 274.3558 -13.69986240 Feb 9 54.644196 274.8412 -20.28538318 Mar 9 31.179871 275.3265 2.99360297 Apr 9 -8.882152 275.4066 4.47554486 May 9 -30.560306 275.4867 34.97361762 Jun 9 41.095776 275.3552 1.44902631 Jul 9 17.255063 275.2237 5.92122964 Aug 9 -14.806234 274.7948 -13.28855866 Sep 9 -32.828831 274.3659 -14.23704674 Oct 9 -49.056164 273.6583 -15.50209339 Nov 9 -27.770598 272.9506 14.71996149 Dec 9 -26.214634 272.3109 19.90369972 Jan 10 45.944015 271.6712 2.98475340 Feb 10 54.644196 272.4364 -18.58057626 Mar 10 31.179871 273.2015 -22.18139899 Apr 10 -8.882152 276.6815 -5.09938675 May 10 -30.560306 280.1615 13.89875636 Jun 10 41.095776 288.2638 -16.25961722 Jul 10 17.255063 296.3661 -29.32119618 Aug 10 -14.806234 310.4815 -43.07527595 Sep 10 -32.828831 324.5969 -41.46805551 Oct 10 -49.056164 342.7066 -47.15042071 Nov 10 -27.770598 360.8163 -20.34568438 Dec 10 -26.214634 379.1000 -19.68538649 Jan 11 45.944015 397.3838 3.07222686 Feb 11 54.644196 412.8446 44.11119807 Mar 11 31.179871 428.3055 56.01467622 Apr 11 -8.882152 438.6090 76.67312637 May 11 -30.560306 448.9126 64.84770741 Jun 11 41.095776 452.7207 28.48354457 Jul 11 17.255063 456.5288 36.01617637 Aug 11 -14.806234 453.2383 22.26796041 Sep 11 -32.828831 449.9478 -11.31895533 Oct 11 -49.056164 440.8087 -16.75254569 Nov 11 -27.770598 431.6696 -25.39903453 Dec 11 -26.214634 420.5671 12.44752610 Jan 12 45.944015 409.4646 12.39140220 Feb 12 54.644196 400.2933 14.86249472 Mar 12 31.179871 391.1220 7.49809418 Apr 12 -8.882152 385.1214 -20.43929468 May 12 -30.560306 379.1209 -15.86055267 Jun 12 41.095776 374.8901 -37.98583541 Jul 12 17.255063 370.6593 -27.41432352 Aug 12 -14.806234 367.7569 -18.25069799 Sep 12 -32.828831 364.8546 -12.52577224 Oct 12 -49.056164 364.2980 7.85817449 Nov 12 -27.770598 363.7414 27.62922274 Dec 12 -26.214634 365.2377 13.07693222 Jan 13 45.944015 366.7340 -0.77804284 Feb 13 54.644196 368.9808 -35.02501529 Mar 13 31.179871 371.2276 13.99251920 Apr 13 -8.882152 375.0052 -5.42306778 May 13 -30.560306 378.7828 -10.22252389 Jun 13 41.095776 385.9076 -9.80333831 Jul 13 17.255063 393.0323 -21.88735809 Aug 13 -14.806234 403.3491 -17.44285876 Sep 13 -32.828831 413.6659 -49.33705922 Oct 13 -49.056164 424.6735 -21.91736521 Nov 13 -27.770598 435.6812 -11.21056968 Dec 13 -26.214634 445.1583 28.05632574 Jan 14 45.944015 454.6354 32.92053662 Feb 14 54.644196 461.0961 49.65970898 Mar 14 31.179871 467.5567 43.56338828 Apr 14 -8.882152 469.6324 27.94973205 May 14 -30.560306 471.7081 25.95220670 Jun 14 41.095776 468.5174 21.68682762 Jul 14 17.255063 465.3267 13.51824317 Aug 14 -14.806234 458.2488 0.55744760 Sep 14 -32.828831 451.1709 -14.94204777 Oct 14 -49.056164 442.9643 -7.60816202 Nov 14 -27.770598 434.7578 -12.88717475 Dec 14 -26.214634 427.0905 3.22414887 Jan 15 45.944015 419.4232 -3.26721206 Feb 15 54.644196 413.1969 -19.74105830 Mar 15 31.179871 406.9705 -5.85039761 Apr 15 -8.882152 403.0147 -7.83254934 May 15 -30.560306 399.0589 26.70142980 Jun 15 41.095776 397.6530 -16.84881605 Jul 15 17.255063 396.2472 -30.60226728 Aug 15 -14.806234 396.7139 2.29232042 Sep 15 -32.828831 397.1806 -18.85179167 Oct 15 -49.056164 398.6394 -26.18327875 Nov 15 -27.770598 400.0983 0.27233570 Dec 15 -26.214634 401.4530 0.76163529 Jan 16 45.944015 402.8077 13.94825034 Feb 16 54.644196 403.8294 28.52638432 Mar 16 31.179871 404.8511 8.16902524 Apr 16 -8.882152 405.3438 2.83838016 May 16 -30.560306 405.8364 19.62386596 Jun 16 41.095776 405.4103 8.89389354 Jul 16 17.255063 404.9842 -8.23928424 Aug 16 -14.806234 403.4035 -13.09728125 Sep 16 -32.828831 401.8228 -21.99397804 Oct 16 -49.056164 399.7988 -11.34262311 Nov 16 -27.770598 397.7748 15.79583333 Dec 16 -26.214634 395.5580 9.45664004 Jan 17 45.944015 393.3412 12.51476219 Feb 17 54.644196 390.9980 0.45780941 Mar 17 31.179871 388.6548 2.66536357 Apr 17 -8.882152 385.9351 6.04703883 May 17 -30.560306 383.2155 0.14484496 Jun 17 41.095776 379.7741 24.43009785 Jul 17 17.255063 376.3328 -26.08785464 Aug 17 -14.806234 372.7571 -2.85084233 Sep 17 -32.828831 369.1814 -10.15252981 Oct 17 -49.056164 366.0945 2.76170167 Nov 17 -27.770598 363.0076 -3.43696532 Dec 17 -26.214634 360.0360 7.07859743 Jan 18 45.944015 357.0645 -8.90852436 Feb 18 54.644196 353.9007 8.65507981 Mar 18 31.179871 350.7369 -12.01680909 Apr 18 -8.882152 346.8809 11.20124744 May 18 -30.560306 343.0249 8.93543484 Jun 18 41.095776 337.7049 26.89933120 Jul 18 17.255063 332.3849 -6.73997781 Aug 18 -14.806234 326.2263 5.07998330 Sep 18 -32.828831 320.0676 -3.03875538 Oct 18 -49.056164 314.6588 5.29731470 Nov 18 -27.770598 309.2501 7.32048630 Dec 18 -26.214634 305.0784 -0.06372138 Jan 19 45.944015 300.9066 -22.45061360 Feb 19 54.644196 297.9048 -41.64903824 Mar 19 31.179871 294.9031 -27.08295594 Apr 19 -8.882152 292.9262 -11.04406959 May 19 -30.560306 290.9494 18.91094763 Jun 19 41.095776 289.7157 28.38856492 Jul 19 17.255063 288.4820 -0.73702316 Aug 19 -14.806234 287.7500 9.15623636 Sep 19 -32.828831 287.0180 -3.88920391 Oct 19 -49.056164 286.5505 9.00563501 Nov 19 -27.770598 286.0830 -0.41242455 Dec 19 -26.214634 286.1175 6.59712954 Jan 20 45.944015 286.1520 -16.19600091 Feb 20 54.644196 287.7698 -24.01404141 Mar 20 31.179871 289.3877 -25.16757497 Apr 20 -8.882152 291.7669 -16.48472434 May 20 -30.560306 294.1460 -17.78574283 Jun 20 41.095776 295.5810 26.12324622 Jul 20 17.255063 297.0159 10.62902990 Aug 20 -14.806234 297.1383 11.86794650 Sep 20 -32.828831 297.2607 25.06816331 Oct 20 -49.056164 296.6558 47.60035091 Nov 20 -27.770598 296.0510 22.01964003 Dec 20 -26.214634 294.9087 3.30592072 Jan 21 45.944015 293.7665 -32.31048314 Feb 21 54.644196 292.0149 -17.95905055 Mar 21 31.179871 290.2632 -28.54311102 Apr 21 -8.882152 288.0855 -30.10334370 May 21 -30.560306 285.9078 -24.94744551 Jun 21 41.095776 283.5474 36.85682349 Jul 21 17.255063 281.1871 23.25788712 Aug 21 -14.806234 279.4682 12.53800500 Sep 21 -32.828831 277.7494 15.77942310 Oct 21 -49.056164 276.6248 23.43132478 Nov 21 -27.770598 275.5003 9.87032799 Dec 21 -26.214634 275.1003 -7.08570214 Jan 22 45.944015 274.7004 -33.14441680 Feb 22 54.644196 275.5102 -37.85434963 Mar 22 31.179871 276.3199 -32.79977553 Apr 22 -8.882152 278.4875 -15.40533371 May 22 -30.560306 280.6551 -20.09476101 Jun 22 41.095776 284.3524 13.55187001 Jul 22 17.255063 288.0496 12.89529567 Aug 22 -14.806234 294.1440 7.66224354 Sep 22 -32.828831 300.2383 28.39049162 Oct 22 -49.056164 308.5279 24.52824156 Nov 22 -27.770598 316.8175 -18.04690699 Dec 22 -26.214634 326.6989 -37.78426606 Jan 23 45.944015 336.5803 -41.92430967 Feb 23 54.644196 348.0471 -23.29124791 Mar 23 31.179871 359.5138 -17.39367921 Apr 23 -8.882152 372.9797 -8.89752584 May 23 -30.560306 386.4455 -17.48524159 Jun 23 41.095776 401.2098 24.59443787 Jul 23 17.255063 415.9740 17.77091195 Aug 23 -14.806234 429.5227 7.28354946 Sep 23 -32.828831 443.0713 18.95748718 Oct 23 -49.056164 453.4180 21.53814504 Nov 23 -27.770598 463.7647 24.70590442 Dec 23 -26.214634 471.4727 18.34192467 Jan 24 45.944015 479.1807 16.27526038 Feb 24 54.644196 484.5998 4.95595673 Mar 24 31.179871 490.0190 -3.69883999 Apr 24 -8.882152 493.2665 -14.98431278 May 24 -30.560306 496.5140 -26.55365469 Jun 24 41.095776 497.8587 10.04552178 Jul 24 17.255063 499.2034 16.54149289 Aug 24 -14.806234 499.4803 21.42597049 Sep 24 -32.828831 499.7571 17.07174830 Oct 24 -49.056164 499.4537 6.60244144 Nov 24 -27.770598 499.1504 10.12023610 Dec 24 -26.214634 498.3845 -2.66982021 Jan 25 45.944015 497.6185 1.13743895 Feb 25 54.644196 496.1011 -9.54525625 Mar 25 31.179871 494.5836 -4.26344451 Apr 25 -8.882152 491.8711 -13.28896660 May 25 -30.560306 489.1587 -24.19835781 Jun 25 41.095776 484.9892 16.51502515 Jul 25 17.255063 480.8197 19.22520274 Aug 25 -14.806234 475.7502 24.75598494 Sep 25 -32.828831 470.6808 27.94806735 Oct 25 -49.056164 465.4835 30.57263211 Nov 25 -27.770598 460.2863 -5.91570161 Dec 25 -26.214634 455.2459 -17.43130091 Jan 26 45.944015 450.2056 -28.64958475 Feb 26 54.644196 445.4600 -15.60420093 Mar 26 31.179871 440.7144 -20.69431018 Apr 26 -8.882152 437.5065 -11.22430672 May 26 -30.560306 434.2985 -23.83817239 Jun 26 41.095776 433.8457 9.75850771 Jul 26 17.255063 433.3930 4.35198245 Aug 26 -14.806234 434.8898 0.71648069 Sep 26 -32.828831 436.3866 12.94227915 Oct 26 -49.056164 439.0036 -13.64748193 Nov 26 -27.770598 441.6207 -8.25014148 Dec 26 -26.214634 445.9569 -13.94230114 Jan 27 45.944015 450.2931 4.56285465 Feb 27 54.644196 457.3454 2.01036148 Mar 27 31.179871 464.3978 -20.07762476 Apr 27 -8.882152 476.1246 -37.14243654 May 27 -30.560306 487.8514 -42.89111744 Jun 27 41.095776 506.8071 -9.90288142 Jul 27 17.255063 525.7628 -17.01785077 Aug 27 -14.806234 550.9982 -47.69196898 Sep 27 -32.828831 576.2336 -23.20478697 Oct 27 -49.056164 604.2780 -50.82185327 Nov 27 -27.770598 632.3224 -36.05181805 Dec 27 -26.214634 659.1118 -22.29712183 Jan 28 45.944015 685.9011 86.15488985 Feb 28 54.644196 708.2842 67.97164876 Mar 28 31.179871 730.6672 74.05291461 Apr 28 -8.882152 746.4538 44.42840125 May 28 -30.560306 762.2403 30.62001877 Jun 28 41.095776 769.5796 46.22461398 Jul 28 17.255063 776.9189 26.72600382 Aug 28 -14.806234 775.4066 8.99961370 Sep 28 -32.828831 773.8943 11.13452380 Oct 28 -49.056164 767.2624 6.19377337 Nov 28 -27.770598 760.6305 -9.75987554 Dec 28 -26.214634 753.8771 -8.16241718 Jan 29 45.944015 747.1236 24.33235664 Feb 29 54.644196 742.1390 6.51683577 Mar 29 31.179871 737.1543 -15.83417815 Apr 29 -8.882152 734.1368 -36.25460236 May 29 -30.560306 731.1192 -70.15889570 Jun 29 41.095776 729.6997 -5.29547417 Jul 29 17.255063 728.2802 12.16474198 Aug 29 -14.806234 727.2876 19.71860289 Sep 29 -32.828831 726.2951 9.13376400 Oct 29 -49.056164 724.7363 7.61987392 Nov 29 -27.770598 723.1775 14.09308537 Dec 29 -26.214634 720.4537 7.96095127 Jan 30 45.944015 717.7299 21.12613264 Feb 30 54.644196 713.1052 43.15064503 Mar 30 31.179871 708.4805 15.93966435 Apr 30 -8.882152 702.2466 -36.56446111 May 30 -30.560306 696.0128 -50.35245569 Jun 30 41.095776 688.3475 15.85668803 Jul 30 17.255063 680.6823 -3.83737361 Aug 30 -14.806234 672.0855 18.42073458 Sep 30 -32.828831 663.4887 13.04014298 Oct 30 -49.056164 655.3012 15.85495719 Nov 30 -27.770598 647.1137 15.25687293 Dec 30 -26.214634 640.2845 -26.06988011 Jan 31 45.944015 633.4553 10.30068230 Feb 31 54.644196 627.9982 -8.74236268 Mar 31 31.179871 622.5410 -5.82090072 Apr 31 -8.882152 618.1226 -40.44047003 May 31 -30.560306 613.7042 -37.44390846 Jun 31 41.095776 609.3486 -17.84437424 Jul 31 17.255063 604.9930 21.55195461 Aug 31 -14.806234 601.0091 6.89711438 Sep 31 -32.828831 597.0253 15.50357436 Oct 31 -49.056164 593.4606 1.59552327 Nov 31 -27.770598 589.8960 0.77457371 Dec 31 -26.214634 586.6065 12.10812638 > m$win s t l 3721 20 13 > m$deg s t l 0 1 1 > m$jump s t l 373 2 2 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/wessaorg/rcomp/tmp/1j7l71352547230.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/2if4k1352547230.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/371ff1352547230.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/4du2e1352547230.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/5xysh1352547230.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/68pp01352547230.tab") > > try(system("convert tmp/1j7l71352547230.ps tmp/1j7l71352547230.png",intern=TRUE)) character(0) > try(system("convert tmp/2if4k1352547230.ps tmp/2if4k1352547230.png",intern=TRUE)) character(0) > try(system("convert tmp/371ff1352547230.ps tmp/371ff1352547230.png",intern=TRUE)) character(0) > try(system("convert tmp/4du2e1352547230.ps tmp/4du2e1352547230.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.999 0.481 8.463