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(235.1
+ ,280.7
+ ,264.6
+ ,240.7
+ ,201.4
+ ,240.8
+ ,241.1
+ ,223.8
+ ,206.1
+ ,174.7
+ ,203.3
+ ,220.5
+ ,299.5
+ ,347.4
+ ,338.3
+ ,327.7
+ ,351.6
+ ,396.6
+ ,438.8
+ ,395.6
+ ,363.5
+ ,378.8
+ ,357
+ ,369
+ ,464.8
+ ,479.1
+ ,431.3
+ ,366.5
+ ,326.3
+ ,355.1
+ ,331.6
+ ,261.3
+ ,249
+ ,205.5
+ ,235.6
+ ,240.9
+ ,264.9
+ ,253.8
+ ,232.3
+ ,193.8
+ ,177
+ ,213.2
+ ,207.2
+ ,180.6
+ ,188.6
+ ,175.4
+ ,199
+ ,179.6
+ ,225.8
+ ,234
+ ,200.2
+ ,183.6
+ ,178.2
+ ,203.2
+ ,208.5
+ ,191.8
+ ,172.8
+ ,148
+ ,159.4
+ ,154.5
+ ,213.2
+ ,196.4
+ ,182.8
+ ,176.4
+ ,153.6
+ ,173.2
+ ,171
+ ,151.2
+ ,161.9
+ ,157.2
+ ,201.7
+ ,236.4
+ ,356.1
+ ,398.3
+ ,403.7
+ ,384.6
+ ,365.8
+ ,368.1
+ ,367.9
+ ,347
+ ,343.3
+ ,292.9
+ ,311.5
+ ,300.9
+ ,366.9
+ ,356.9
+ ,329.7
+ ,316.2
+ ,269
+ ,289.3
+ ,266.2
+ ,253.6
+ ,233.8
+ ,228.4
+ ,253.6
+ ,260.1
+ ,306.6
+ ,309.2
+ ,309.5
+ ,271
+ ,279.9
+ ,317.9
+ ,298.4
+ ,246.7
+ ,227.3
+ ,209.1
+ ,259.9
+ ,266
+ ,320.6
+ ,308.5
+ ,282.2
+ ,262.7
+ ,263.5
+ ,313.1
+ ,284.3
+ ,252.6
+ ,250.3
+ ,246.5
+ ,312.7
+ ,333.2
+ ,446.4
+ ,511.6
+ ,515.5
+ ,506.4
+ ,483.2
+ ,522.3
+ ,509.8
+ ,460.7
+ ,405.8
+ ,375
+ ,378.5
+ ,406.8
+ ,467.8
+ ,469.8
+ ,429.8
+ ,355.8
+ ,332.7
+ ,378
+ ,360.5
+ ,334.7
+ ,319.5
+ ,323.1
+ ,363.6
+ ,352.1
+ ,411.9
+ ,388.6
+ ,416.4
+ ,360.7
+ ,338
+ ,417.2
+ ,388.4
+ ,371.1
+ ,331.5
+ ,353.7
+ ,396.7
+ ,447
+ ,533.5
+ ,565.4
+ ,542.3
+ ,488.7
+ ,467.1
+ ,531.3
+ ,496.1
+ ,444
+ ,403.4
+ ,386.3
+ ,394.1
+ ,404.1
+ ,462.1
+ ,448.1
+ ,432.3
+ ,386.3
+ ,395.2
+ ,421.9
+ ,382.9
+ ,384.2
+ ,345.5
+ ,323.4
+ ,372.6
+ ,376
+ ,462.7
+ ,487
+ ,444.2
+ ,399.3
+ ,394.9
+ ,455.4
+ ,414
+ ,375.5
+ ,347
+ ,339.4
+ ,385.8
+ ,378.8
+ ,451.8
+ ,446.1
+ ,422.5
+ ,383.1
+ ,352.8
+ ,445.3
+ ,367.5
+ ,355.1
+ ,326.2
+ ,319.8
+ ,331.8
+ ,340.9
+ ,394.1
+ ,417.2
+ ,369.9
+ ,349.2
+ ,321.4
+ ,405.7
+ ,342.9
+ ,316.5
+ ,284.2
+ ,270.9
+ ,288.8
+ ,278.8
+ ,324.4
+ ,310.9
+ ,299
+ ,273
+ ,279.3
+ ,359.2
+ ,305
+ ,282.1
+ ,250.3
+ ,246.5
+ ,257.9
+ ,266.5
+ ,315.9
+ ,318.4
+ ,295.4
+ ,266.4
+ ,245.8
+ ,362.8
+ ,324.9
+ ,294.2
+ ,289.5
+ ,295.2
+ ,290.3
+ ,272
+ ,307.4
+ ,328.7
+ ,292.9
+ ,249.1
+ ,230.4
+ ,361.5
+ ,321.7
+ ,277.2
+ ,260.7
+ ,251
+ ,257.6
+ ,241.8
+ ,287.5
+ ,292.3
+ ,274.7
+ ,254.2
+ ,230
+ ,339
+ ,318.2
+ ,287
+ ,295.8
+ ,284
+ ,271
+ ,262.7
+ ,340.6
+ ,379.4
+ ,373.3
+ ,355.2
+ ,338.4
+ ,466.9
+ ,451
+ ,422
+ ,429.2
+ ,425.9
+ ,460.7
+ ,463.6
+ ,541.4
+ ,544.2
+ ,517.5
+ ,469.4
+ ,439.4
+ ,549
+ ,533
+ ,506.1
+ ,484
+ ,457
+ ,481.5
+ ,469.5
+ ,544.7
+ ,541.2
+ ,521.5
+ ,469.7
+ ,434.4
+ ,542.6
+ ,517.3
+ ,485.7
+ ,465.8
+ ,447
+ ,426.6
+ ,411.6
+ ,467.5
+ ,484.5
+ ,451.2
+ ,417.4
+ ,379.9
+ ,484.7
+ ,455
+ ,420.8
+ ,416.5
+ ,376.3
+ ,405.6
+ ,405.8
+ ,500.8
+ ,514
+ ,475.5
+ ,430.1
+ ,414.4
+ ,538
+ ,526
+ ,488.5
+ ,520.2
+ ,504.4
+ ,568.5
+ ,610.6
+ ,818
+ ,830.9
+ ,835.9
+ ,782
+ ,762.3
+ ,856.9
+ ,820.9
+ ,769.6
+ ,752.2
+ ,724.4
+ ,723.1
+ ,719.5
+ ,817.4
+ ,803.3
+ ,752.5
+ ,689
+ ,630.4
+ ,765.5
+ ,757.7
+ ,732.2
+ ,702.6
+ ,683.3
+ ,709.5
+ ,702.2
+ ,784.8
+ ,810.9
+ ,755.6
+ ,656.8
+ ,615.1
+ ,745.3
+ ,694.1
+ ,675.7
+ ,643.7
+ ,622.1
+ ,634.6
+ ,588
+ ,689.7
+ ,673.9
+ ,647.9
+ ,568.8
+ ,545.7
+ ,632.6
+ ,643.8
+ ,593.1
+ ,579.7
+ ,546
+ ,562.9
+ ,572.5)
> par8 = 'FALSE'
> par7 = '1'
> par6 = ''
> par5 = '0'
> par4 = ''
> par3 = '0'
> 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  45.540853 231.0872 -41.52802069
Feb  1  54.276511 231.6971  -5.27364381
Mar  1  30.847662 232.3071   1.44524000
Apr  1  -9.102209 233.1958  16.60642563
May  1 -30.668211 234.0845  -2.01625787
Jun  1  41.196236 235.6256 -36.02188061
Jul  1  17.563888 237.1668 -13.63070872
Aug  1 -14.418589 240.2135  -1.99488181
Sep  1 -32.362366 243.2601  -4.79775469
Oct  1 -48.791405 250.2753 -26.78386344
Nov  1 -27.707545 257.2904 -26.28287067
Dec  1 -26.374826 270.4400 -23.56512744
Jan  2  45.540853 283.5895 -29.63034368
Feb  2  54.276511 299.2361  -6.11264020
Mar  2  30.847662 314.8828  -7.43042978
Apr  2  -9.102209 330.5174   6.28480379
May  2 -30.668211 346.1520  36.11616823
Jun  2  41.196236 359.6769  -4.27312309
Jul  2  17.563888 373.2017  48.03438023
Aug  2 -14.418589 382.6585  27.36008955
Sep  2 -32.362366 392.1153   3.74709907
Oct  2 -48.791405 395.2174  32.37403725
Nov  2 -27.707545 398.3195 -13.61192306
Dec  2 -26.374826 394.2423   1.13250723
Jan  3  45.540853 390.1652  29.09397803
Feb  3  54.276511 380.6340  44.18947664
Mar  3  30.847662 371.1029  29.34948219
Apr  3  -9.102209 358.7150  16.88717814
May  3 -30.668211 346.3272  10.64100497
Jun  3  41.196236 331.8336 -17.92980408
Jul  3  17.563888 317.3399  -3.30381849
Aug  3 -14.418589 301.2914 -25.57284603
Sep  3 -32.362366 285.2429  -3.88057335
Oct  3 -48.791405 270.8216 -16.53022376
Nov  3 -27.707545 256.4003   6.90722735
Dec  3 -26.374826 245.0664  22.20842038
Jan  4  45.540853 233.7325 -14.37334607
Feb  4  54.276511 225.8690 -26.34550482
Mar  4  30.847662 218.0055 -16.55315664
Apr  4  -9.102209 213.5036 -10.60138554
May  4 -30.668211 209.0017  -1.33348357
Jun  4  41.196236 206.4641 -34.46031834
Jul  4  17.563888 203.9265 -14.29035848
Aug  4 -14.418589 202.2807  -7.26206181
Sep  4 -32.362366 200.6348  20.32753507
Oct  4 -48.791405 199.5745  24.61692147
Nov  4 -27.707545 198.5141  28.19340940
Dec  4 -26.374826 197.6164   8.35846400
Jan  5  45.540853 196.7186 -16.45944088
Feb  5  54.276511 195.5144 -15.79093719
Mar  5  30.847662 194.3103 -24.95792656
Apr  5  -9.102209 192.7653  -0.06312846
May  5 -30.668211 191.2204  17.64780052
Jun  5  41.196236 189.9262 -27.92243692
Jul  5  17.563888 188.6320   2.30412027
Aug  5 -14.418589 186.9949  19.22365466
Sep  5 -32.362366 185.3579  19.80448927
Oct  5 -48.791405 183.2515  13.53988978
Nov  5 -27.707545 181.1452   5.96239182
Dec  5 -26.374826 178.2344   2.64039164
Jan  6  45.540853 175.3237  -7.66456801
Feb  6  54.276511 172.8373 -30.71383166
Mar  6  30.847662 170.3509 -18.39858838
Apr  6  -9.102209 171.2278  14.27437165
May  6 -30.668211 172.1047  12.16346256
Jun  6  41.196236 179.4160 -47.41224604
Jul  6  17.563888 186.7273 -33.29116001
Aug  6 -14.418589 200.9663 -35.34775808
Sep  6 -32.362366 215.2054 -20.94305594
Oct  6 -48.791405 233.8600 -27.86857202
Nov  6 -27.707545 252.5145 -23.10698658
Dec  6 -26.374826 271.4884  -8.71356898
Jan  7  45.540853 290.4623  20.09688915
Feb  7  54.276511 306.9142  37.10928710
Mar  7  30.847662 323.3661  49.48619199
Apr  7  -9.102209 334.4790  59.22324822
May  7 -30.668211 345.5918  50.87643532
Jun  7  41.196236 349.4612 -22.55740514
Jul  7  17.563888 353.3306  -2.99445097
Aug  7 -14.418589 350.3397  11.07890663
Sep  7 -32.362366 347.3488  28.31356444
Oct  7 -48.791405 341.2572   0.43421223
Nov  7 -27.707545 335.1656   4.04196153
Dec  7 -26.374826 327.7085  -0.43365622
Jan  8  45.540853 320.2514   1.10776655
Feb  8  54.276511 312.1640  -9.54046674
Mar  8  30.847662 304.0765  -5.22419310
Apr  8  -9.102209 297.3110  27.99117539
May  8 -30.668211 290.5455   9.12267474
Jun  8  41.196236 285.8974 -37.79362381
Jul  8  17.563888 281.2492 -32.61312774
Aug  8 -14.418589 277.8921  -9.87346627
Sep  8 -32.362366 274.5349  -8.37250459
Oct  8 -48.791405 273.4990   3.69243806
Nov  8 -27.707545 272.4631   8.84448224
Dec  8 -26.374826 273.6340  12.84078943
Jan  9  45.540853 274.8050 -13.74586286
Feb  9  54.276511 275.4611 -20.53765489
Mar  9  30.847662 276.1173   2.53506002
Apr  9  -9.102209 275.8005   4.30167307
May  9 -30.668211 275.4838  35.08441699
Jun  9  41.196236 275.5630   1.14071851
Jul  9  17.563888 275.6423   5.19381467
Aug  9 -14.418589 275.1542 -14.03562807
Sep  9 -32.362366 274.6661 -15.00377059
Oct  9 -48.791405 273.4900 -15.59858506
Nov  9 -27.707545 272.3138  15.29370201
Dec  9 -26.374826 271.7592  20.61558804
Jan 10  45.540853 271.2046   3.85451459
Feb 10  54.276511 271.8066 -17.58312222
Mar 10  30.847662 272.4086 -21.05625209
Apr 10  -9.102209 274.8444  -3.04223279
May 10 -30.668211 277.2803  16.88791738
Jun 10  41.196236 284.2678 -12.36402032
Jul 10  17.563888 291.2553 -24.51916339
Aug 10 -14.418589 305.3875 -38.36888353
Sep 10 -32.362366 319.5197 -36.85730345
Oct 10 -48.791405 338.9186 -43.62718803
Nov 10 -27.707545 358.3175 -17.90997108
Dec 10 -26.374826 378.9420 -19.36713933
Jan 11  45.540853 399.5664   1.29273294
Feb 11  54.276511 417.3006  40.02284377
Mar 11  30.847662 435.0349  49.61746155
Apr 11  -9.102209 445.7219  69.78030679
May 11 -30.668211 456.4089  57.45928291
Jun 11  41.196236 459.0490  22.05478172
Jul 11  17.563888 461.6890  30.54707516
Aug 11 -14.418589 456.9271  18.19144410
Sep 11 -32.362366 452.1653 -14.00288674
Oct 11 -48.791405 442.0755 -18.28405653
Nov 11 -27.707545 431.9857 -25.77812481
Dec 11 -26.374826 420.4131  12.76172130
Jan 12  45.540853 408.8405  13.41860793
Feb 12  54.276511 399.1605  16.36299111
Mar 12  30.847662 389.4805   9.47188123
Apr 12  -9.102209 383.2428 -18.34062600
May 12 -30.668211 377.0052 -13.63700234
Jun 12  41.196236 372.7108 -35.90703704
Jul 12  17.563888 368.4164 -25.48027711
Aug 12 -14.418589 365.7467 -16.62812089
Sep 12 -32.362366 363.0770 -11.21466445
Oct 12 -48.791405 363.2006   8.69082056
Nov 12 -27.707545 363.3241  27.98340709
Dec 12 -26.374826 365.2746  13.20020105
Jan 13  45.540853 367.2251  -0.86596448
Feb 13  54.276511 369.0273 -34.70383873
Mar 13  30.847662 370.8295  14.72279396
Apr 13  -9.102209 373.4310  -3.62877654
May 13 -30.668211 376.0324  -7.36421616
Jun 13  41.196236 382.7056  -6.70184241
Jul 13  17.563888 389.3788 -18.54267402
Aug 13 -14.418589 400.4835 -14.96486975
Sep 13 -32.362366 411.5881 -47.72576527
Oct 13 -48.791405 423.7765 -21.28512448
Nov 13 -27.707545 435.9649 -11.55738216
Dec 13 -26.374826 447.0236  26.35117792
Jan 14  45.540853 458.0824  29.87677852
Feb 14  54.276511 465.5611  45.56235968
Mar 14  30.847662 473.0399  38.41244777
Apr 14  -9.102209 474.7435  23.05865993
May 14 -30.668211 476.4472  21.32100297
Jun 14  41.196236 471.9768  18.12696312
Jul 14  17.563888 467.5064  11.02971790
Aug 14 -14.418589 459.2309  -0.81229820
Sep 14 -32.362366 450.9554 -15.19301409
Oct 14 -48.791405 442.3039  -7.21252820
Nov 14 -27.707545 433.6525 -11.84494078
Dec 14 -26.374826 425.9467   4.52809288
Jan 15  45.540853 418.2410  -1.68183295
Feb 15  54.276511 412.1227 -18.29923826
Mar 15  30.847662 406.0045  -4.55213663
Apr 15  -9.102209 401.8087  -6.40651509
May 15 -30.668211 397.6130  28.25523732
Jun 15  41.196236 395.9571 -15.25331359
Jul 15  17.563888 394.3012 -28.96506988
Aug 15 -14.418589 395.0239   3.59464217
Sep 15 -32.362366 395.7467 -17.88434556
Oct 15 -48.791405 397.6367 -25.44528889
Nov 15 -27.707545 399.5267   0.78086930
Dec 15 -26.374826 401.7006   0.67426334
Jan 16  45.540853 403.8744  13.28469790
Feb 16  54.276511 405.2500  27.47344485
Mar 16  30.847662 406.6256   6.72669874
Apr 16  -9.102209 406.7118   1.69044546
May 16 -30.668211 406.7979  18.77032306
Jun 16  41.196236 405.8618   8.34194992
Jul 16  17.563888 404.9257  -8.48962859
Aug 16 -14.418589 403.1463 -13.22766860
Sep 16 -32.362366 401.3668 -22.00440840
Oct 16 -48.791405 399.4988 -11.30742993
Nov 16 -27.707545 397.6309  15.87665006
Dec 16 -26.374826 395.8205   9.35428421
Jan 17  45.540853 394.0102  12.24895887
Feb 17  54.276511 391.8843  -0.06084737
Mar 17  30.847662 389.7585   1.89385332
Apr 17  -9.102209 386.6804   5.52182207
May 17 -30.668211 383.6023  -0.13407831
Jun 17  41.196236 379.8234  24.28037854
Jul 17  17.563888 376.0445 -26.10836997
Aug 17 -14.418589 372.3711  -2.85250378
Sep 17 -32.362366 368.6977 -10.13533736
Oct 17 -48.791405 365.5645   3.02688024
Nov 17 -27.707545 362.4313  -2.92380063
Dec 17 -26.374826 359.8287   7.44614525
Jan 18  45.540853 357.2260  -8.66686835
Feb 18  54.276511 354.3925   8.53101971
Mar 18  30.847662 351.5589 -12.50658529
Apr 18  -9.102209 347.7568  10.54543923
May 18 -30.668211 343.9546   8.11359462
Jun 18  41.196236 338.8352  25.66855503
Jul 18  17.563888 333.7158  -8.37968992
Aug 18 -14.418589 327.0534   3.86517848
Sep 18 -32.362366 320.3910  -3.82865290
Oct 18 -48.791405 314.0807   5.61071288
Nov 18 -27.707545 307.7704   8.73718017
Dec 18 -26.374826 303.4034   1.77142677
Jan 19  45.540853 299.0364 -20.17728612
Feb 19  54.276511 296.3084 -39.68491935
Mar 19  30.847662 293.5804 -25.42804564
Apr 19  -9.102209 291.9212  -9.81901801
May 19 -30.668211 290.2621  19.70614049
Jun 19  41.196236 289.7544  28.24935628
Jul 19  17.563888 289.2467  -1.81063331
Aug 19 -14.418589 288.7625   7.75613616
Sep 19 -32.362366 288.2782  -5.61579415
Oct 19 -48.791405 286.8547   8.43674056
Nov 19 -27.707545 285.4312   0.17637680
Dec 19 -26.374826 284.8178   8.05705700
Jan 20  45.540853 284.2044 -13.84522228
Feb 20  54.276511 285.8858 -21.76232401
Mar 20  30.847662 287.5673 -23.01491881
Apr 20  -9.102209 290.7096 -15.20741881
May 20 -30.668211 293.8520 -17.38378793
Jun 20  41.196236 296.3885  25.21530501
Jul 20  17.563888 298.9249   8.41119259
Aug 20 -14.418589 299.5779   9.04068111
Sep 20 -32.362366 300.2309  21.63146985
Oct 20 -48.791405 298.7410  45.25037348
Nov 20 -27.707545 297.2512  20.75637864
Dec 20 -26.374826 295.0771   3.29768903
Jan 21  45.540853 292.9031 -31.04396005
Feb 21  54.276511 290.7092 -16.28569115
Mar 21  30.847662 288.5153 -26.46291531
Apr 21  -9.102209 286.7383 -28.53610537
May 21 -30.668211 284.9614 -23.89316455
Jun 21  41.196236 283.7217  36.58202610
Jul 21  17.563888 282.4821  21.65401139
Aug 21 -14.418589 281.0295  10.58906323
Sep 21 -32.362366 279.5770  13.48541528
Oct 21 -48.791405 277.5042  22.28715812
Nov 21 -27.707545 275.4315   9.87600248
Dec 21 -26.374826 273.8182  -5.64342037
Jan 22  45.540853 272.2049 -30.24580270
Feb 22  54.276511 272.7616 -34.73815790
Mar 22  30.847662 273.3183 -29.46600616
Apr 22  -9.102209 276.3128 -13.01054720
May 22 -30.668211 279.3072 -18.63895737
Jun 22  41.196236 283.8106  13.99315305
Jul 22  17.563888 288.3141  12.32205811
Aug 22 -14.418589 294.4237   6.99489556
Sep 22 -32.362366 300.5333  27.62903323
Oct 22 -48.791405 307.9622  24.82924273
Nov 22 -27.707545 315.3910 -16.68344624
Dec 22 -26.374826 324.5799 -35.50504319
Jan 23  45.540853 333.7687 -38.70959963
Feb 23  54.276511 345.2013 -20.07781026
Mar 23  30.847662 356.6339 -14.18151396
Apr 23  -9.102209 370.9925  -6.69030983
May 23 -30.668211 385.3512 -16.28297483
Jun 23  41.196236 401.5619  24.14182411
Jul 23  17.563888 417.7727  15.66341769
Aug 23 -14.418589 432.0876   4.33097672
Sep 23 -32.362366 446.4025  15.15983597
Oct 23 -48.791405 456.5097  18.18169671
Nov 23 -27.707545 466.6169  21.79065898
Dec 23 -26.374826 473.5854  16.38939665
Jan 24  45.540853 480.5540  15.30517485
Feb 24  54.276511 485.6043   4.31917502
Mar 24  30.847662 490.6547  -4.00231787
Apr 24  -9.102209 493.8106 -15.30842010
May 24 -30.668211 496.9666 -26.89839146
Jun 24  41.196236 498.4618   9.34198439
Jul 24  17.563888 499.9570  15.47915486
Aug 24 -14.418589 500.6067  19.91190642
Sep 24 -32.362366 501.2564  15.10595819
Oct 24 -48.791405 500.8004   4.99100092
Nov 24 -27.707545 500.3444   8.86314518
Dec 24 -26.374826 498.8344  -2.95955468
Jan 25  45.540853 497.3244   1.83478598
Feb 25  54.276511 495.6048  -8.68131247
Mar 25  30.847662 493.8852  -3.23290398
Apr 25  -9.102209 491.9029 -13.10068873
May 25 -30.668211 489.9206 -24.85234260
Jun 25  41.196236 486.2315  15.17227411
Jul 25  17.563888 482.5424  17.19368545
Aug 25 -14.418589 477.4330  22.68558575
Sep 25 -32.362366 472.3236  25.83878626
Oct 25 -48.791405 466.4993  29.29213370
Nov 25 -27.707545 460.6750  -6.36741733
Dec 25 -26.374826 454.6594 -16.68456274
Jan 26  45.540853 448.6438 -26.68466762
Feb 26  54.276511 443.4559 -13.23237349
Mar 26  30.847662 438.2679 -17.91557243
Apr 26  -9.102209 435.3528  -8.85058628
May 26 -30.668211 432.4377 -21.86946925
Jun 26  41.196236 432.5247  10.97911212
Jul 26  17.563888 432.6116   4.82448813
Aug 26 -14.418589 434.4856   0.73301922
Sep 26 -32.362366 436.3595  12.50285053
Oct 26 -48.791405 438.6452 -13.55382830
Nov 26 -27.707545 440.9310  -7.62340560
Dec 26 -26.374826 444.6076 -12.43282218
Jan 27  45.540853 448.2843   6.97480175
Feb 27  54.276511 454.8737   4.84976496
Mar 27  30.847662 461.4631 -16.81076489
Apr 27  -9.102209 472.3562 -33.15399927
May 27 -30.668211 483.2493 -38.18110276
Jun 27  41.196236 501.3092  -4.50542700
Jul 27  17.563888 519.3691 -10.93295660
Aug 27 -14.418589 545.5855 -42.66695071
Sep 27 -32.362366 571.8020 -19.23964460
Oct 27 -48.791405 602.0898 -48.89837647
Nov 27 -27.707545 632.3776 -36.17000681
Dec 27 -26.374826 661.2399 -24.26505871
Jan 28  45.540853 690.1022  82.35692992
Feb 28  54.276511 714.2111  62.41237669
Mar 28  30.847662 738.3200  66.73233040
Apr 28  -9.102209 754.3745  36.72772710
May 28 -30.668211 770.4290  22.53925468
Jun 28  41.196236 776.2484  39.45541091
Jul 28  17.563888 782.0678  21.26836177
Aug 28 -14.418589 779.2235   4.79505889
Sep 28 -32.362366 776.3793   8.18305622
Oct 28 -48.791405 768.8627   4.32873604
Nov 28 -27.707545 761.3460 -10.53848261
Dec 28 -26.374826 753.4215  -7.54668941
Jan 29  45.540853 745.4970  26.36214431
Feb 29  54.276511 740.1696   8.85385097
Mar 29  30.847662 734.8423 -13.18993543
Apr 29  -9.102209 731.9694 -33.86721183
May 29 -30.668211 729.0966 -68.02835736
Jun 29  41.196236 727.9250  -3.62124636
Jul 29  17.563888 726.7535  13.38265928
Aug 29 -14.418589 727.0888  19.52977512
Sep 29 -32.362366 727.4242   7.53819118
Oct 29 -48.791405 726.8241   5.26729505
Nov 29 -27.707545 726.2240  10.98350045
Dec 29 -26.374826 722.8468   5.72799109
Jan 30  45.540853 719.4696  19.78952226
Feb 30  54.276511 714.1252  42.49825019
Mar 30  30.847662 708.7809  15.97148507
Apr 30  -9.102209 702.5441 -36.64191803
May 30 -30.668211 696.3074 -50.53919026
Jun 30  41.196236 688.4428  15.66095121
Jul 30  17.563888 680.5782  -4.04211269
Aug 30 -14.418589 672.2565  17.86209495
Sep 30 -32.362366 663.9348  12.12760279
Oct 30 -48.791405 655.9433  14.94809012
Nov 30 -27.707545 647.9519  14.35567896
Dec 30 -26.374826 640.2238 -25.84901246
Jan 31  45.540853 632.4958  11.66333664
Feb 31  54.276511 625.6815  -6.05805930
Mar 31  30.847662 618.8673  -1.81494830
Apr 31  -9.102209 615.7290 -37.82675295
May 31 -30.668211 612.5906 -36.22242674
Jun 31  41.196236 610.9613 -19.55754012
Jul 31  17.563888 609.3320  16.90414113
Aug 31 -14.418589 608.4648  -0.94617377
Sep 31 -32.362366 607.5976   4.46481153
Oct 31 -48.791405 607.0357 -12.24431082
Nov 31 -27.707545 606.4739 -15.86633166
Dec 31 -26.374826 606.0261  -7.15126876
> m$win
   s    t    l 
3721   19   13 
> m$deg
s t l 
0 0 1 
> m$jump
  s   t   l 
373   2   2 
> m$inner
[1] 2
> m$outer
[1] 0
> postscript(file="/var/wessaorg/rcomp/tmp/18cjf1324304075.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/2dugw1324304075.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/3iwy21324304075.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/4w9wc1324304075.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/5p1on1324304076.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/6992j1324304076.tab") 
> 
> try(system("convert tmp/18cjf1324304075.ps tmp/18cjf1324304075.png",intern=TRUE))
character(0)
> try(system("convert tmp/2dugw1324304075.ps tmp/2dugw1324304075.png",intern=TRUE))
character(0)
> try(system("convert tmp/3iwy21324304075.ps tmp/3iwy21324304075.png",intern=TRUE))
character(0)
> try(system("convert tmp/4w9wc1324304075.ps tmp/4w9wc1324304075.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  4.885   0.301   5.187