R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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. Natural language support but running in an English locale 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(10 + ,10 + ,10 + ,10 + ,10 + ,9.94 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,10.06 + ,9.94 + ,9.94 + ,9.94 + ,9.94 + ,9.94 + ,9.94 + ,10.06 + ,10.06 + ,9.94 + ,10.06 + ,10.06 + ,10.06 + ,10.18 + ,10.28 + ,10.28 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.28 + ,10.28 + ,10.18 + ,10.18 + ,10.18 + ,10.28 + ,10.28 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.18 + ,10.18 + ,10.28 + ,10.18 + ,10.18 + ,10.28 + ,10.18 + ,10.18 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.28 + ,10.18 + ,10.28 + ,10.28 + ,10.28 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.34 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.42 + ,10.34 + ,10.34 + ,10.42 + ,10.42 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.49 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.49 + ,10.55 + ,10.55 + ,10.55 + ,10.63 + ,10.63 + ,10.63 + ,10.63 + ,10.63 + ,10.57 + ,10.63 + ,10.63 + ,10.63 + ,10.63 + ,10.57 + ,10.57 + ,10.57 + ,10.63 + ,10.63 + ,10.63 + ,10.63 + ,10.57 + ,10.63 + ,10.63 + ,10.57 + ,10.63 + ,10.63 + ,10.63 + ,10.63 + ,10.57 + ,10.63 + ,10.6 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.6 + ,10.6 + ,10.6 + ,10.6 + ,10.7 + ,10.7 + ,10.6 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.7 + ,10.8 + ,10.8 + ,10.8 + ,10.8 + ,10.7 + ,10.7 + ,10.7 + ,10.85 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.75 + ,10.75 + ,10.85 + ,10.85 + ,10.85 + ,10.85 + ,10.75 + ,10.75 + ,10.75 + ,10.75 + ,10.75) > par8 = 'FALSE' > par7 = '1' > par6 = '' > par5 = '1' > 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 0.004849491 9.991738 0.0034126600 Feb 1 -0.003473069 9.997631 0.0058420522 Mar 1 -0.012443821 10.003524 0.0089196361 Apr 1 -0.018211463 10.009447 0.0087644293 May 1 -0.012944627 10.015370 -0.0024252557 Jun 1 0.004163783 10.021115 -0.0852783780 Jul 1 0.009203268 10.026859 0.0239374255 Aug 1 0.003275244 10.032408 0.0243168870 Sep 1 0.000105832 10.037956 0.0219377373 Oct 1 0.012216833 10.043792 0.0039906809 Nov 1 0.005362292 10.049629 0.0050091670 Dec 1 0.007896237 10.053793 -0.0016890081 Jan 2 0.004849491 10.057957 -0.0028064917 Feb 2 -0.003473069 10.056371 0.0071023582 Mar 2 -0.012443821 10.054784 0.0176593999 Apr 2 -0.018211463 10.047121 0.0310903736 May 2 -0.012944627 10.039458 0.0334868691 Jun 2 0.004163783 10.028507 0.0273292786 Jul 2 0.009203268 10.017556 0.0332406138 Aug 2 0.003275244 10.008289 0.0484354465 Sep 2 0.000105832 9.999022 -0.0591283320 Oct 2 0.012216833 9.993575 -0.0657918902 Nov 2 0.005362292 9.988128 -0.0534899059 Dec 2 0.007896237 9.987578 -0.0554741219 Jan 3 0.004849491 9.987028 -0.0518776463 Feb 3 -0.003473069 9.996089 -0.0526162416 Mar 3 -0.012443821 10.005150 0.0672933548 Apr 3 -0.018211463 10.026593 0.0516181879 May 3 -0.012944627 10.048036 -0.0950914572 Jun 3 0.004163783 10.073987 -0.0181508257 Jul 3 0.009203268 10.099938 -0.0491412686 Aug 3 0.003275244 10.124598 -0.0678734916 Sep 3 0.000105832 10.149258 0.0306356742 Oct 3 0.012216833 10.170955 0.0968286496 Nov 3 0.005362292 10.192651 0.0819871675 Dec 3 0.007896237 10.209476 -0.0373725143 Jan 4 0.004849491 10.226302 0.0488484955 Feb 4 -0.003473069 10.232923 0.0505503060 Mar 4 -0.012443821 10.239544 0.0529003084 Apr 4 -0.018211463 10.239112 -0.0409003310 May 4 -0.012944627 10.238680 0.0542645514 Jun 4 0.004163783 10.234879 0.0409568488 Jul 4 0.009203268 10.231079 -0.0602819282 Aug 4 0.003275244 10.224937 -0.0482125986 Sep 4 0.000105832 10.218796 -0.0389018802 Oct 4 0.012216833 10.213948 0.0538354733 Nov 4 0.005362292 10.209099 0.0655383693 Dec 4 0.007896237 10.209498 -0.0373945028 Jan 5 0.004849491 10.209897 -0.0347466833 Feb 5 -0.003473069 10.212392 -0.0289188327 Mar 5 -0.012443821 10.214887 -0.0224427903 Apr 5 -0.018211463 10.214647 -0.0164360323 May 5 -0.012944627 10.214408 -0.0214637525 Jun 5 0.004163783 10.214559 0.0612769042 Jul 5 0.009203268 10.214710 0.0560864865 Aug 5 0.003275244 10.216557 0.0601681679 Sep 5 0.000105832 10.218403 -0.0385087619 Oct 5 0.012216833 10.218013 -0.0502301482 Nov 5 0.005362292 10.217624 -0.0429859920 Dec 5 0.007896237 10.215928 0.0561761623 Jan 6 0.004849491 10.214232 -0.0390809917 Feb 6 -0.003473069 10.216337 -0.0328637102 Mar 6 -0.012443821 10.218442 0.0740017631 Apr 6 -0.018211463 10.224681 -0.0264697160 May 6 -0.012944627 10.230920 -0.0379756732 Jun 6 0.004163783 10.234713 -0.0588772075 Jul 6 0.009203268 10.238507 0.0322901839 Aug 6 0.003275244 10.238394 0.0383304533 Sep 6 0.000105832 10.238282 0.0416121114 Oct 6 0.012216833 10.235805 0.0319785026 Nov 6 0.005362292 10.233327 0.0413104363 Dec 6 0.007896237 10.228486 0.0436182161 Jan 7 0.004849491 10.223644 -0.0484933125 Feb 7 -0.003473069 10.216242 -0.0327690195 Mar 7 -0.012443821 10.208840 -0.0163965347 Apr 7 -0.018211463 10.205299 -0.0070875901 May 7 -0.012944627 10.201758 -0.0088131238 Jun 7 0.004163783 10.204900 -0.0290641478 Jul 7 0.009203268 10.208043 -0.0372462460 Aug 7 0.003275244 10.216150 -0.0394252844 Sep 7 0.000105832 10.224257 -0.0443629340 Oct 7 0.012216833 10.233196 0.0345874845 Nov 7 0.005362292 10.242134 0.0325034455 Dec 7 0.007896237 10.249511 0.0225925039 Jan 8 0.004849491 10.256888 0.0182622539 Feb 8 -0.003473069 10.261925 0.0215478352 Mar 8 -0.012443821 10.266962 0.0254816082 Apr 8 -0.018211463 10.267294 0.0309178029 May 8 -0.012944627 10.267625 0.0253195195 Jun 8 0.004163783 10.265269 -0.0894326595 Jul 8 0.009203268 10.262913 0.0078840872 Aug 8 0.003275244 10.259668 0.0170565724 Sep 8 0.000105832 10.256424 0.0234704465 Oct 8 0.012216833 10.253904 0.0138791285 Nov 8 0.005362292 10.251384 -0.0767466470 Dec 8 0.007896237 10.251465 0.0206386141 Jan 9 0.004849491 10.251546 0.0236045668 Feb 9 -0.003473069 10.252954 0.0305193142 Mar 9 -0.012443821 10.254362 -0.0619177466 Apr 9 -0.018211463 10.255763 -0.0575512867 May 9 -0.012944627 10.257164 0.0357806950 Jun 9 0.004163783 10.257802 0.0180342695 Jul 9 0.009203268 10.258440 0.0123567696 Aug 9 0.003275244 10.261234 0.0154911972 Sep 9 0.000105832 10.264027 0.0158670135 Oct 9 0.012216833 10.273331 -0.0055481748 Nov 9 0.005362292 10.282636 -0.0079978207 Dec 9 0.007896237 10.295042 -0.1229380218 Jan 10 0.004849491 10.307448 -0.0322975313 Feb 10 -0.003473069 10.321483 -0.0380103623 Mar 10 -0.012443821 10.335519 -0.0430750016 Apr 10 -0.018211463 10.350476 0.0877351453 May 10 -0.012944627 10.365434 0.0675108139 Jun 10 0.004163783 10.378899 0.0369369743 Jul 10 0.009203268 10.392365 0.0184320604 Aug 10 0.003275244 10.398101 0.0186235754 Sep 10 0.000105832 10.403838 0.0160564791 Oct 10 0.012216833 10.402811 0.0049726459 Nov 10 0.005362292 10.401783 0.0128543552 Dec 10 0.007896237 10.399882 0.0122218075 Jan 11 0.004849491 10.397981 0.0171699514 Feb 11 -0.003473069 10.396568 -0.0530953024 Mar 11 -0.012443821 10.395156 -0.0427123645 Apr 11 -0.018211463 10.394086 -0.0358744807 May 11 -0.012944627 10.393016 0.0399289248 Jun 11 0.004163783 10.394271 0.0215650521 Jul 11 0.009203268 10.395527 0.0152701050 Aug 11 0.003275244 10.399359 0.0173660908 Sep 11 0.000105832 10.403191 -0.0632965347 Oct 11 0.012216833 10.406807 0.0009766145 Nov 11 0.005362292 10.410422 0.0042153062 Dec 11 0.007896237 10.412357 -0.0002537340 Jan 12 0.004849491 10.414293 0.0008579174 Feb 12 -0.003473069 10.416215 0.0072582051 Mar 12 -0.012443821 10.418137 0.0143066846 Apr 12 -0.018211463 10.419361 0.0188509263 May 12 -0.012944627 10.420584 0.0123606899 Jun 12 0.004163783 10.420199 -0.0043627071 Jul 12 0.009203268 10.419814 -0.0090171785 Aug 12 0.003275244 10.417426 -0.0007013545 Sep 12 0.000105832 10.415038 0.0048558583 Oct 12 0.012216833 10.410109 -0.0023259637 Nov 12 0.005362292 10.405180 0.0094577569 Dec 12 0.007896237 10.400305 0.0117983817 Jan 13 0.004849491 10.395431 0.0197196982 Feb 13 -0.003473069 10.391302 0.0321712905 Mar 13 -0.012443821 10.387173 -0.0347289255 Apr 13 -0.018211463 10.382335 -0.0241239355 May 13 -0.012944627 10.377498 -0.0245534236 Jun 13 0.004163783 10.374710 -0.0388738284 Jul 13 0.009203268 10.371922 0.0388746925 Aug 13 0.003275244 10.372231 0.0444934746 Sep 13 0.000105832 10.372541 -0.0326463545 Oct 13 0.012216833 10.373183 -0.0454003286 Nov 13 0.005362292 10.373826 -0.0391887600 Dec 13 0.007896237 10.372660 0.0394441650 Jan 14 0.004849491 10.371493 0.0436577817 Feb 14 -0.003473069 10.371888 0.0515853042 Mar 14 -0.012443821 10.372283 -0.0198389816 Apr 14 -0.018211463 10.374265 -0.0160532914 May 14 -0.012944627 10.376247 -0.0233020795 Jun 14 0.004163783 10.377389 -0.0415526454 Jul 14 0.009203268 10.378531 -0.0477342856 Aug 14 0.003275244 10.381496 0.0352292070 Sep 14 0.000105832 10.384460 0.0354340883 Oct 14 0.012216833 10.390878 0.0169051312 Nov 14 0.005362292 10.397296 -0.0626582834 Dec 14 0.007896237 10.402984 0.0091195648 Jan 15 0.004849491 10.408672 0.0064781046 Feb 15 -0.003473069 10.411726 0.0117465720 Mar 15 -0.012443821 10.414781 0.0176632312 Apr 15 -0.018211463 10.416344 0.0218675744 May 15 -0.012944627 10.417907 0.0150374394 Jun 15 0.004163783 10.416851 -0.0010145108 Jul 15 0.009203268 10.415794 -0.0049975353 Aug 15 0.003275244 10.410988 0.0057367910 Sep 15 0.000105832 10.406182 0.0137125061 Oct 15 0.012216833 10.402629 0.0051545859 Nov 15 0.005362292 10.399075 0.0155622083 Dec 15 0.007896237 10.398218 0.0138855672 Jan 16 0.004849491 10.397361 -0.0622103823 Feb 16 -0.003473069 10.397706 -0.0542331201 Mar 16 -0.012443821 10.398051 -0.0456076661 Apr 16 -0.018211463 10.399358 0.0388529711 May 16 -0.012944627 10.400665 0.0322791302 Jun 16 0.004163783 10.403903 0.0119333596 Jul 16 0.009203268 10.407140 0.0036565146 Aug 16 0.003275244 10.411089 0.0056354921 Sep 16 0.000105832 10.415038 0.0048558584 Oct 16 0.012216833 10.416345 -0.0085615112 Nov 16 0.005362292 10.417651 -0.0030133382 Dec 16 0.007896237 10.415364 -0.0032606092 Jan 17 0.004849491 10.413078 0.0020728115 Feb 17 -0.003473069 10.411692 0.0117810741 Mar 17 -0.012443821 10.410306 0.0221375285 Apr 17 -0.018211463 10.414822 0.0233899227 May 17 -0.012944627 10.419337 0.0136078387 Jun 17 0.004163783 10.429137 -0.0933009577 Jul 17 0.009203268 10.438938 -0.1081408284 Aug 17 0.003275244 10.451232 -0.0345072778 Sep 17 0.000105832 10.463527 -0.0436323385 Oct 17 0.012216833 10.478098 0.0596854040 Nov 17 0.005362292 10.492669 0.0519686890 Dec 17 0.007896237 10.507431 0.0346728035 Jan 18 0.004849491 10.522193 0.0229576096 Feb 18 -0.003473069 10.530494 0.0229790327 Mar 18 -0.012443821 10.538795 0.0236486476 Apr 18 -0.018211463 10.539868 0.0283438931 May 18 -0.012944627 10.540940 0.0220046603 Jun 18 0.004163783 10.539823 0.0060133680 Jul 18 0.009203268 10.538706 0.0020910014 Aug 18 0.003275244 10.537705 -0.0509802915 Sep 18 0.000105832 10.536704 -0.0468101956 Oct 18 0.012216833 10.535953 0.0018299096 Nov 18 0.005362292 10.535202 0.0094355575 Dec 18 0.007896237 10.536266 0.0058380319 Jan 19 0.004849491 10.537329 0.0078211980 Feb 19 -0.003473069 10.539696 0.0137775098 Mar 19 -0.012443821 10.542062 0.0203820134 Apr 19 -0.018211463 10.542048 -0.0338366690 May 19 -0.012944627 10.542034 0.0209101703 Jun 19 0.004163783 10.539628 0.0062077237 Jul 19 0.009203268 10.537223 0.0035742028 Aug 19 0.003275244 10.535715 0.0110093192 Sep 19 0.000105832 10.534208 0.0156858244 Oct 19 0.012216833 10.533887 0.0038958368 Nov 19 0.005362292 10.533566 -0.0489286082 Dec 19 0.007896237 10.533617 -0.0515134901 Jan 20 0.004849491 10.533668 0.0114823196 Feb 20 -0.003473069 10.534614 0.0188593778 Mar 20 -0.012443821 10.535559 0.0268846278 Apr 20 -0.018211463 10.537662 -0.0294506308 May 20 -0.012944627 10.539765 0.0231796324 Jun 20 0.004163783 10.541311 0.0045247970 Jul 20 0.009203268 10.542858 -0.0020611127 Aug 20 0.003275244 10.541250 0.0054748813 Sep 20 0.000105832 10.539642 0.0102522640 Oct 20 0.012216833 10.537595 0.0001885464 Nov 20 0.005362292 10.535547 0.0090903713 Dec 20 0.007896237 10.534180 0.0079236106 Jan 21 0.004849491 10.532813 0.0123375415 Feb 21 -0.003473069 10.531557 -0.0380841848 Mar 21 -0.012443821 10.530302 -0.0278577194 Apr 21 -0.018211463 10.530416 -0.0222043095 May 21 -0.012944627 10.530530 0.0324146222 Jun 21 0.004163783 10.535202 0.0106343794 Jul 21 0.009203268 10.539874 0.0009230624 Aug 21 0.003275244 10.548991 -0.0622658413 Sep 21 0.000105832 10.558108 -0.0082133561 Oct 21 0.012216833 10.567652 -0.0298689960 Nov 21 0.005362292 10.577197 -0.0325590934 Dec 21 0.007896237 10.585580 0.0365233072 Jan 22 0.004849491 10.593964 0.0311863993 Feb 22 -0.003473069 10.601561 0.0319115757 Mar 22 -0.012443821 10.609159 0.0332849438 Apr 22 -0.018211463 10.612013 0.0361987974 May 22 -0.012944627 10.614866 -0.0319218271 Jun 22 0.004163783 10.613156 0.0126797398 Jul 22 0.009203268 10.611446 0.0093502325 Aug 22 0.003275244 10.609762 0.0169625660 Sep 22 0.000105832 10.608078 0.0218162883 Oct 22 0.012216833 10.608012 -0.0502286233 Nov 22 0.005362292 10.607946 -0.0433079922 Dec 22 0.007896237 10.608117 -0.0460130852 Jan 23 0.004849491 10.608288 0.0168625136 Feb 23 -0.003473069 10.608729 0.0247441501 Mar 23 -0.012443821 10.609170 0.0332739784 Apr 23 -0.018211463 10.611380 0.0368316289 May 23 -0.012944627 10.613590 -0.0306451989 Jun 23 0.004163783 10.614602 0.0112337917 Jul 23 0.009203268 10.615615 0.0051817080 Aug 23 0.003275244 10.615180 -0.0484556615 Sep 23 0.000105832 10.614746 0.0151483578 Oct 23 0.012216833 10.618185 -0.0004020555 Nov 23 0.005362292 10.621625 0.0030130736 Dec 23 0.007896237 10.629108 -0.0070037891 Jan 24 0.004849491 10.636590 -0.0714399602 Feb 24 -0.003473069 10.644867 -0.0113942284 Mar 24 -0.012443821 10.653144 -0.0407003047 Apr 24 -0.018211463 10.660913 0.0572989195 May 24 -0.012944627 10.668681 0.0442636656 Jun 24 0.004163783 10.675191 0.0206447206 Jul 24 0.009203268 10.681702 0.0090947013 Aug 24 0.003275244 10.682801 0.0139240266 Sep 24 0.000105832 10.683899 0.0159947406 Oct 24 0.012216833 10.679087 0.0086964256 Nov 24 0.005362292 10.674274 0.0203636531 Dec 24 0.007896237 10.669129 0.0229745455 Jan 25 0.004849491 10.663984 0.0311661294 Feb 25 -0.003473069 10.660661 -0.0571879158 Mar 25 -0.012443821 10.657338 -0.0448937693 Apr 25 -0.018211463 10.656262 -0.0380502285 May 25 -0.012944627 10.655186 -0.0422411660 Jun 25 0.004163783 10.657514 0.0383223634 Jul 25 0.009203268 10.659842 0.0309548185 Aug 25 0.003275244 10.665908 -0.0691836891 Sep 25 0.000105832 10.671975 0.0279191921 Oct 25 0.012216833 10.678787 0.0089966229 Nov 25 0.005362292 10.685598 0.0090395962 Dec 25 0.007896237 10.690004 0.0021002122 Jan 26 0.004849491 10.694409 0.0007415198 Feb 26 -0.003473069 10.697025 0.0064476928 Mar 26 -0.012443821 10.699642 0.0128020576 Apr 26 -0.018211463 10.702726 0.0154850369 May 26 -0.012944627 10.705811 0.0071335379 Jun 26 0.004163783 10.711674 -0.0158379488 Jul 26 0.009203268 10.717537 -0.0267405099 Aug 26 0.003275244 10.723604 -0.0268790174 Sep 26 0.000105832 10.729670 -0.0297761362 Oct 26 0.012216833 10.734079 -0.0462960285 Nov 26 0.005362292 10.738488 0.0561496217 Dec 26 0.007896237 10.743665 0.0484385245 Jan 27 0.004849491 10.748842 0.0463081190 Feb 27 -0.003473069 10.753059 0.0504142616 Mar 27 -0.012443821 10.757275 -0.0448314040 Apr 27 -0.018211463 10.757446 -0.0392349435 May 27 -0.012944627 10.757618 -0.0446729612 Jun 27 0.004163783 10.754932 0.0909044692 Jul 27 0.009203268 10.752246 -0.0114491747 Aug 27 0.003275244 10.751335 -0.0046099233 Sep 27 0.000105832 10.750423 -0.0005292831 Oct 27 0.012216833 10.752124 -0.0143405246 Nov 27 0.005362292 10.753824 -0.0091862237 Dec 27 0.007896237 10.757866 -0.0157620254 Jan 28 0.004849491 10.761908 -0.0167571356 Feb 28 -0.003473069 10.769535 -0.0160616224 Mar 28 -0.012443821 10.777162 -0.0147179175 Apr 28 -0.018211463 10.786574 -0.0183622481 May 28 -0.012944627 10.795986 -0.0330410569 Jun 28 0.004163783 10.805183 0.0406535733 Jul 28 0.009203268 10.814380 0.0264171292 Aug 28 0.003275244 10.821992 0.0247322770 Sep 28 0.000105832 10.829605 0.0202888135 Oct 28 0.012216833 10.832485 0.0052979931 Nov 28 0.005362292 10.835365 0.0092727152 Dec 28 0.007896237 10.834686 0.0074173671 Jan 29 0.004849491 10.834008 0.0111427105 Feb 29 -0.003473069 10.832146 0.0213271836 Mar 29 -0.012443821 10.830284 0.0321598485 Apr 29 -0.018211463 10.825710 -0.0574990139 May 29 -0.012944627 10.821137 -0.0581923546 Jun 29 0.004163783 10.812441 0.0333948074 Jul 29 0.009203268 10.803746 0.0370508951 Aug 29 0.003275244 10.795749 0.0509760353 Sep 29 0.000105832 10.787752 0.0621425642 Oct 29 0.012216833 10.779732 -0.0419484723 Nov 29 0.005362292 10.771712 -0.0270739664 Dec 29 0.007896237 10.763208 -0.0211043359 Jan 30 0.004849491 10.754705 -0.0095540139 Feb 30 -0.003473069 10.745844 0.0076290510 > m$win s t l 3501 19 13 > m$deg s t l 0 1 1 > m$jump s t l 351 2 2 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/www/html/freestat/rcomp/tmp/1lt071292236987.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/www/html/freestat/rcomp/tmp/2lt071292236987.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/www/html/freestat/rcomp/tmp/3ekzs1292236987.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/www/html/freestat/rcomp/tmp/4ekzs1292236987.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/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/www/html/freestat/rcomp/tmp/52lem1292236987.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/www/html/freestat/rcomp/tmp/6vud71292236987.tab") > > try(system("convert tmp/1lt071292236987.ps tmp/1lt071292236987.png",intern=TRUE)) character(0) > try(system("convert tmp/2lt071292236987.ps tmp/2lt071292236987.png",intern=TRUE)) character(0) > try(system("convert tmp/3ekzs1292236987.ps tmp/3ekzs1292236987.png",intern=TRUE)) character(0) > try(system("convert tmp/4ekzs1292236987.ps tmp/4ekzs1292236987.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.543 1.018 4.835