R version 2.9.0 (2009-04-17) Copyright (C) 2009 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. 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(1.3866 + ,1.3582 + ,1.3332 + ,1.3595 + ,1.3617 + ,1.3684 + ,1.3394 + ,1.3262 + ,1.3173 + ,1.3085 + ,1.327 + ,1.3182 + ,1.293 + ,1.291 + ,1.2984 + ,1.2795 + ,1.299 + ,1.3174 + ,1.326 + ,1.3111 + ,1.2816 + ,1.276 + ,1.2849 + ,1.2818 + ,1.2829 + ,1.2796 + ,1.3008 + ,1.2967 + ,1.2938 + ,1.2833 + ,1.2823 + ,1.2765 + ,1.2634 + ,1.2596 + ,1.2705 + ,1.2591 + ,1.2798 + ,1.2763 + ,1.2795 + ,1.2782 + ,1.2644 + ,1.2596 + ,1.2615 + ,1.2555 + ,1.2555 + ,1.2658 + ,1.2565 + ,1.2783 + ,1.2786 + ,1.2782 + ,1.2905 + ,1.3042 + ,1.2942 + ,1.313 + ,1.3671 + ,1.3549 + ,1.3558 + ,1.3507 + ,1.3494 + ,1.3607 + ,1.3295 + ,1.3193 + ,1.3308 + ,1.3246 + ,1.3392 + ,1.3425 + ,1.3496 + ,1.3255 + ,1.3231 + ,1.3273 + ,1.3276 + ,1.3173 + ,1.3196 + ,1.3058 + ,1.2966 + ,1.2932 + ,1.2947 + ,1.305 + ,1.3232 + ,1.3125 + ,1.2992 + ,1.3266 + ,1.3275 + ,1.3223 + ,1.3403 + ,1.3322 + ,1.3363 + ,1.3425 + ,1.3574 + ,1.3683 + ,1.3623 + ,1.3563 + ,1.3518 + ,1.3494 + ,1.3612 + ,1.369 + ,1.3771 + ,1.3972 + ,1.401 + ,1.3908 + ,1.3901 + ,1.3856 + ,1.4098 + ,1.422 + ,1.4238 + ,1.4207 + ,1.4095 + ,1.4177 + ,1.3866 + ,1.3959 + ,1.4102 + ,1.3969 + ,1.4004 + ,1.385 + ,1.389 + ,1.384 + ,1.392 + ,1.3932 + ,1.3858 + ,1.3978 + ,1.4029 + ,1.394 + ,1.4096 + ,1.4058 + ,1.4134 + ,1.4096 + ,1.4049 + ,1.4009 + ,1.3897 + ,1.4019 + ,1.3901 + ,1.399 + ,1.3901 + ,1.3975 + ,1.3991 + ,1.4089 + ,1.413 + ,1.409 + ,1.4217 + ,1.4223 + ,1.4191 + ,1.4229 + ,1.4227 + ,1.4269 + ,1.4229 + ,1.4104 + ,1.4053 + ,1.4138 + ,1.4303 + ,1.4384 + ,1.441 + ,1.437 + ,1.4357 + ,1.4202 + ,1.4166 + ,1.417 + ,1.4293 + ,1.4294 + ,1.4072 + ,1.4101 + ,1.4112 + ,1.4243 + ,1.433 + ,1.4323 + ,1.4324 + ,1.427 + ,1.4268 + ,1.4364 + ,1.4272 + ,1.4314 + ,1.422 + ,1.4335 + ,1.4262 + ,1.433 + ,1.4473 + ,1.4522 + ,1.4545 + ,1.4594 + ,1.4561 + ,1.4611 + ,1.4671 + ,1.4712 + ,1.4705 + ,1.4658 + ,1.478 + ,1.4783 + ,1.4768 + ,1.467 + ,1.465 + ,1.4549 + ,1.4643 + ,1.4539 + ,1.4537 + ,1.4616 + ,1.4722 + ,1.4694 + ,1.4763 + ,1.475 + ,1.4765 + ,1.4864 + ,1.4881 + ,1.4864 + ,1.4869 + ,1.4918 + ,1.4971 + ,1.4921 + ,1.5 + ,1.502 + ,1.5019 + ,1.4874 + ,1.4785 + ,1.4788 + ,1.48 + ,1.4772 + ,1.4658 + ,1.4761 + ,1.4867 + ,1.4862 + ,1.4984 + ,1.4966 + ,1.5037 + ,1.4922 + ,1.4868 + ,1.4965 + ,1.4875 + ,1.4957 + ,1.4863 + ,1.4815 + ,1.4968 + ,1.4969 + ,1.5083 + ,1.5071 + ,1.4918 + ,1.5023 + ,1.5074 + ,1.509 + ,1.512 + ,1.5068 + ,1.4787 + ,1.4774 + ,1.4768 + ,1.473 + ,1.4757 + ,1.4647 + ,1.4541 + ,1.456 + ,1.4343 + ,1.4337 + ,1.4368 + ,1.4279 + ,1.4276 + ,1.4398 + ,1.4405 + ,1.4433 + ,1.4338 + ,1.4406 + ,1.4389 + ,1.4442 + ,1.435 + ,1.4304 + ,1.4273 + ,1.4528 + ,1.4481 + ,1.4563 + ,1.4486 + ,1.4374 + ,1.4369 + ,1.4279 + ,1.4132 + ,1.4064 + ,1.4135 + ,1.4151 + ,1.4085 + ,1.4072 + ,1.3999 + ,1.3966 + ,1.3913 + ,1.3937 + ,1.3984 + ,1.3847 + ,1.3691 + ,1.3675 + ,1.376 + ,1.374 + ,1.3718 + ,1.3572 + ,1.3607 + ,1.3649 + ,1.3726 + ,1.3567 + ,1.3519 + ,1.3626 + ,1.3577 + ,1.3547 + ,1.3489 + ,1.357 + ,1.3525 + ,1.3548 + ,1.3641 + ,1.3668 + ,1.3582 + ,1.3662 + ,1.3557 + ,1.361 + ,1.3657 + ,1.3765 + ,1.3705 + ,1.3723 + ,1.3756 + ,1.366 + ,1.3548 + ,1.3471 + ,1.3519 + ,1.3338 + ,1.3356 + ,1.3353 + ,1.3471 + ,1.3482 + ,1.3479 + ,1.3468 + ,1.3396 + ,1.334 + ,1.3296 + ,1.3384 + ,1.3585 + ,1.3583 + ,1.3615 + ,1.3544 + ,1.3535 + ,1.3432 + ,1.3486 + ,1.3373 + ,1.3339 + ,1.3311 + ,1.3321 + ,1.329 + ,1.3245 + ,1.3256 + ,1.3315 + ,1.3238 + ,1.3089 + ,1.2924 + ,1.2727 + ,1.2746 + ,1.2969 + ,1.2698 + ,1.2686 + ,1.2587 + ,1.2492 + ,1.2349 + ,1.2428 + ,1.227 + ,1.2334 + ,1.2497 + ,1.236 + ,1.2223 + ,1.2309 + ,1.2255 + ,1.2384 + ,1.2307 + ,1.2155 + ,1.2218 + ,1.2268 + ,1.206 + ,1.1959 + ,1.1942 + ,1.201 + ,1.2045 + ,1.2127 + ,1.2249 + ,1.2258 + ,1.2277 + ,1.2363 + ,1.2372 + ,1.2391 + ,1.2258 + ,1.2271 + ,1.2262 + ,1.2294 + ,1.2339 + ,1.2198 + ,1.2271 + ,1.2328 + ,1.2548 + ,1.2531 + ,1.2579 + ,1.2567 + ,1.266 + ,1.2637 + ,1.2572 + ,1.2569 + ,1.2703 + ,1.2828 + ,1.3 + ,1.2957 + ,1.2844 + ,1.2817 + ,1.285 + ,1.2897 + ,1.2931 + ,1.3033 + ,1.2992 + ,1.3069 + ,1.3028 + ,1.3073 + ,1.3221 + ,1.3206 + ,1.3184 + ,1.3176 + ,1.3253 + ,1.3133 + ,1.3016 + ,1.279 + ,1.2799 + ,1.282 + ,1.286 + ,1.288 + ,1.2836 + ,1.2711 + ,1.2704 + ,1.2611 + ,1.2613 + ,1.2693 + ,1.2713 + ,1.27 + ,1.268 + ,1.28 + ,1.2818 + ,1.2834 + ,1.2874 + ,1.2744 + ,1.2697 + ,1.2715 + ,1.2725 + ,1.2801 + ,1.285 + ,1.2989 + ,1.3078 + ,1.306 + ,1.3074 + ,1.312 + ,1.3364 + ,1.3323 + ,1.3412 + ,1.3477 + ,1.346 + ,1.3611 + ,1.3648 + ,1.3726 + ,1.3705 + ,1.378 + ,1.3856 + ,1.397 + ,1.3874 + ,1.3936 + ,1.3833 + ,1.3958 + ,1.4101 + ,1.4089 + ,1.3896 + ,1.3859 + ,1.3861 + ,1.4016 + ,1.3934 + ,1.4031 + ,1.3912 + ,1.3803 + ,1.3857 + ,1.3857 + ,1.3926 + ,1.4018 + ,1.4014 + ,1.4244 + ,1.4084 + ,1.3917 + ,1.3945 + ,1.377 + ,1.37 + ,1.3711 + ,1.3626 + ,1.3612 + ,1.3481 + ,1.3647 + ,1.3674 + ,1.3647 + ,1.3496 + ,1.3339 + ,1.3321 + ,1.3225 + ,1.3146 + ,1.2998) > par8 = 'FALSE' > par7 = '1' > par6 = '' > par5 = '1' > par4 = '' > par3 = '0' > par2 = 'periodic' > par1 = '5' > main = 'Seasonal Decomposition by Loess' > par1 <- 5 > 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 Time Series: Start = c(1, 1) End = c(99, 1) Frequency = 5 seasonal trend remainder 1.0 -0.0011337034 1.368933 1.880096e-02 1.2 -0.0006619511 1.364769 -5.906719e-03 1.4 0.0002066172 1.360908 -2.791431e-02 1.6 0.0004998966 1.356627 2.373178e-03 1.8 0.0010891406 1.353636 6.974528e-03 2.0 -0.0011337034 1.350091 1.944284e-02 2.2 -0.0006619511 1.342269 -2.206565e-03 2.4 0.0002066172 1.332138 -6.144289e-03 2.6 0.0004998966 1.324046 -7.245838e-03 2.8 0.0010891406 1.318703 -1.129202e-02 3.0 -0.0011337034 1.313554 1.457999e-02 3.2 -0.0006619511 1.308664 1.019764e-02 3.4 0.0002066172 1.303552 -1.075845e-02 3.6 0.0004998966 1.296999 -6.499030e-03 3.8 0.0010891406 1.293891 3.420217e-03 4.0 -0.0011337034 1.297073 -1.643902e-02 4.2 -0.0006619511 1.302799 -3.136839e-03 4.4 0.0002066172 1.306705 1.048838e-02 4.6 0.0004998966 1.306777 1.872327e-02 4.8 0.0010891406 1.302314 7.696380e-03 5.0 -0.0011337034 1.294942 -1.220787e-02 5.2 -0.0006619511 1.287308 -1.064575e-02 5.4 0.0002066172 1.282532 2.161093e-03 5.6 0.0004998966 1.282252 -9.516702e-04 5.8 0.0010891406 1.285139 -3.328281e-03 6.0 -0.0011337034 1.288271 -7.537677e-03 6.2 -0.0006619511 1.290508 1.095362e-02 6.4 0.0002066172 1.291359 5.134452e-03 6.6 0.0004998966 1.290231 3.068614e-03 6.8 0.0010891406 1.286090 -3.878907e-03 7.0 -0.0011337034 1.279921 3.512490e-03 7.2 -0.0006619511 1.274046 3.115815e-03 7.4 0.0002066172 1.269721 -6.527192e-03 7.6 0.0004998966 1.266745 -7.645300e-03 7.8 0.0010891406 1.266700 2.710729e-03 8.0 -0.0011337034 1.269321 -9.087566e-03 8.2 -0.0006619511 1.272560 7.902377e-03 8.4 0.0002066172 1.274726 1.367109e-03 8.6 0.0004998966 1.274714 4.286489e-03 8.8 0.0010891406 1.272014 5.096838e-03 9.0 -0.0011337034 1.268166 -2.632759e-03 9.2 -0.0006619511 1.263772 -3.509892e-03 9.4 0.0002066172 1.260298 9.949418e-04 9.6 0.0004998966 1.259189 -4.189363e-03 9.8 0.0010891406 1.259734 -5.323138e-03 10.0 -0.0011337034 1.262463 4.470536e-03 10.2 -0.0006619511 1.266853 -9.690795e-03 10.4 0.0002066172 1.271608 6.485358e-03 10.6 0.0004998966 1.277475 6.251668e-04 10.8 0.0010891406 1.284330 -7.219036e-03 11.0 -0.0011337034 1.290057 1.576376e-03 11.2 -0.0006619511 1.298428 6.433847e-03 11.4 0.0002066172 1.311993 -1.799970e-02 11.6 0.0004998966 1.325858 -1.335746e-02 11.8 0.0010891406 1.337954 2.805678e-02 12.0 -0.0011337034 1.348197 7.836996e-03 12.2 -0.0006619511 1.353731 2.731239e-03 12.4 0.0002066172 1.353332 -2.838571e-03 12.6 0.0004998966 1.348976 -7.592539e-05 12.8 0.0010891406 1.343050 1.656052e-02 13.0 -0.0011337034 1.337592 -6.958490e-03 13.2 -0.0006619511 1.332641 -1.267856e-02 13.4 0.0002066172 1.329915 6.786084e-04 13.6 0.0004998966 1.331910 -7.809492e-03 13.8 0.0010891406 1.335742 2.368487e-03 14.0 -0.0011337034 1.336965 6.668564e-03 14.2 -0.0006619511 1.335909 1.435319e-02 14.4 0.0002066172 1.333340 -8.047043e-03 14.6 0.0004998966 1.329467 -6.866951e-03 14.8 0.0010891406 1.325380 8.304900e-04 15.0 -0.0011337034 1.322640 6.094118e-03 15.2 -0.0006619511 1.319093 -1.131161e-03 15.4 0.0002066172 1.313335 6.058542e-03 15.6 0.0004998966 1.306994 -1.693765e-03 15.8 0.0010891406 1.302038 -6.526991e-03 16.0 -0.0011337034 1.300113 -5.779177e-03 16.2 -0.0006619511 1.302149 -6.786997e-03 16.4 0.0002066172 1.305409 -6.160692e-04 16.6 0.0004998966 1.308512 1.418784e-02 16.8 0.0010891406 1.312782 -1.371589e-03 17.0 -0.0011337034 1.316369 -1.603490e-02 17.2 -0.0006619511 1.318952 8.309503e-03 17.4 0.0002066172 1.323571 3.722188e-03 17.6 0.0004998966 1.328709 -6.908611e-03 17.8 0.0010891406 1.332000 7.211227e-03 18.0 -0.0011337034 1.335649 -2.315260e-03 18.2 -0.0006619511 1.341149 -4.187178e-03 18.4 0.0002066172 1.347263 -4.969135e-03 18.6 0.0004998966 1.353098 3.802200e-03 18.8 0.0010891406 1.357315 9.896220e-03 19.0 -0.0011337034 1.358735 4.698416e-03 19.2 -0.0006619511 1.357634 -6.719749e-04 19.4 0.0002066172 1.356511 -4.918025e-03 19.6 0.0004998966 1.357807 -8.906599e-03 19.8 0.0010891406 1.362593 -2.482502e-03 20.0 -0.0011337034 1.370947 -8.135595e-04 20.2 -0.0006619511 1.380181 -2.418568e-03 20.4 0.0002066172 1.387043 9.950631e-03 20.6 0.0004998966 1.391071 9.429333e-03 20.8 0.0010891406 1.393130 -3.419318e-03 21.0 -0.0011337034 1.395393 -4.159166e-03 21.2 -0.0006619511 1.399794 -1.353224e-02 21.4 0.0002066172 1.406165 3.428145e-03 21.6 0.0004998966 1.412518 8.982467e-03 21.8 0.0010891406 1.416833 5.877653e-03 22.0 -0.0011337034 1.416936 4.898019e-03 22.2 -0.0006619511 1.412242 -2.079861e-03 22.4 0.0002066172 1.406831 1.066197e-02 22.6 0.0004998966 1.403449 -1.734933e-02 22.8 0.0010891406 1.400960 -6.149412e-03 23.0 -0.0011337034 1.398977 1.235628e-02 23.2 -0.0006619511 1.397772 -2.100202e-04 23.4 0.0002066172 1.395404 4.789431e-03 23.6 0.0004998966 1.391872 -7.371460e-03 23.8 0.0010891406 1.389714 -1.803137e-03 24.0 -0.0011337034 1.388934 -3.800779e-03 24.2 -0.0006619511 1.389258 3.404003e-03 24.4 0.0002066172 1.390991 2.002189e-03 24.6 0.0004998966 1.393388 -8.088292e-03 24.8 0.0010891406 1.395396 1.315075e-03 25.0 -0.0011337034 1.398256 5.777352e-03 25.2 -0.0006619511 1.401803 -7.141023e-03 25.4 0.0002066172 1.404693 4.700840e-03 25.6 0.0004998966 1.406791 -1.491341e-03 25.8 0.0010891406 1.407925 4.385487e-03 26.0 -0.0011337034 1.406803 3.930700e-03 26.2 -0.0006619511 1.404026 1.536330e-03 26.4 0.0002066172 1.400888 -1.948364e-04 26.6 0.0004998966 1.397949 -8.749182e-03 26.8 0.0010891406 1.396009 4.801476e-03 27.0 -0.0011337034 1.395087 -3.853326e-03 27.2 -0.0006619511 1.395193 4.469394e-03 27.4 0.0002066172 1.396029 -6.135886e-03 27.6 0.0004998966 1.398469 -1.468569e-03 27.8 0.0010891406 1.401885 -3.874264e-03 28.0 -0.0011337034 1.405868 4.165648e-03 28.2 -0.0006619511 1.410360 3.301813e-03 28.4 0.0002066172 1.414351 -5.557800e-03 28.6 0.0004998966 1.417052 4.148160e-03 28.8 0.0010891406 1.419291 1.919905e-03 29.0 -0.0011337034 1.421366 -1.132689e-03 29.2 -0.0006619511 1.422556 1.006118e-03 29.4 0.0002066172 1.422582 -8.829184e-05 29.6 0.0004998966 1.420982 5.417618e-03 29.8 0.0010891406 1.418178 3.632972e-03 30.0 -0.0011337034 1.416168 -4.634146e-03 30.2 -0.0006619511 1.416576 -1.061452e-02 30.4 0.0002066172 1.420022 -6.428819e-03 30.6 0.0004998966 1.425913 3.886975e-03 30.8 0.0010891406 1.431997 5.314225e-03 31.0 -0.0011337034 1.435352 6.781656e-03 31.2 -0.0006619511 1.434237 3.424818e-03 31.4 0.0002066172 1.430112 5.381388e-03 31.6 0.0004998966 1.425822 -6.122100e-03 31.8 0.0010891406 1.423385 -7.874304e-03 32.0 -0.0011337034 1.422101 -3.967232e-03 32.2 -0.0006619511 1.420616 9.346432e-03 32.4 0.0002066172 1.418932 1.026188e-02 32.6 0.0004998966 1.417221 -1.052115e-02 32.8 0.0010891406 1.416308 -7.297314e-03 33.0 -0.0011337034 1.417907 -5.573463e-03 33.2 -0.0006619511 1.422160 2.801532e-03 33.4 0.0002066172 1.426619 6.174535e-03 33.6 0.0004998966 1.429393 2.407445e-03 33.8 0.0010891406 1.430354 9.565191e-04 34.0 -0.0011337034 1.430479 -2.345425e-03 34.2 -0.0006619511 1.430143 -2.681499e-03 34.4 0.0002066172 1.429678 6.515056e-03 34.6 0.0004998966 1.429363 -2.662781e-03 34.8 0.0010891406 1.429199 1.112226e-03 35.0 -0.0011337034 1.428745 -5.611468e-03 35.2 -0.0006619511 1.429605 4.557152e-03 35.4 0.0002066172 1.432844 -6.850137e-03 35.6 0.0004998966 1.437789 -5.289340e-03 35.8 0.0010891406 1.443265 2.945856e-03 36.0 -0.0011337034 1.449034 4.300075e-03 36.2 -0.0006619511 1.453642 1.519941e-03 36.4 0.0002066172 1.456731 2.462768e-03 36.6 0.0004998966 1.459594 -3.994163e-03 36.8 0.0010891406 1.462617 -2.606332e-03 37.0 -0.0011337034 1.465260 2.973660e-03 37.2 -0.0006619511 1.467663 4.199148e-03 37.4 0.0002066172 1.470200 9.365685e-05 37.6 0.0004998966 1.472339 -7.039052e-03 37.8 0.0010891406 1.473531 3.380033e-03 38.0 -0.0011337034 1.473634 5.799593e-03 38.2 -0.0006619511 1.472225 5.236961e-03 38.4 0.0002066172 1.468833 -2.039659e-03 38.6 0.0004998966 1.464989 -4.889068e-04 38.8 0.0010891406 1.461362 -7.551047e-03 39.0 -0.0011337034 1.458818 6.615385e-03 39.2 -0.0006619511 1.458521 -3.958924e-03 39.4 0.0002066172 1.460276 -6.782548e-03 39.6 0.0004998966 1.462856 -1.755525e-03 39.8 0.0010891406 1.466694 4.416915e-03 40.0 -0.0011337034 1.470745 -2.113991e-04 40.2 -0.0006619511 1.473894 3.067745e-03 40.4 0.0002066172 1.476841 -2.048003e-03 40.6 0.0004998966 1.479954 -3.953621e-03 40.8 0.0010891406 1.482644 2.666843e-03 41.0 -0.0011337034 1.485183 4.050842e-03 41.2 -0.0006619511 1.487705 -6.426052e-04 41.4 0.0002066172 1.489647 -2.953411e-03 41.6 0.0004998966 1.491293 7.545174e-06 41.8 0.0010891406 1.493663 2.347930e-03 42.0 -0.0011337034 1.496360 -3.126372e-03 42.2 -0.0006619511 1.497811 2.850572e-03 42.4 0.0002066172 1.496839 4.954659e-03 42.6 0.0004998966 1.493866 7.533674e-03 42.8 0.0010891406 1.489617 -3.306300e-03 43.0 -0.0011337034 1.484944 -5.310514e-03 43.2 -0.0006619511 1.480456 -9.941082e-04 43.4 0.0002066172 1.477050 2.743418e-03 43.6 0.0004998966 1.475912 7.878204e-04 43.8 0.0010891406 1.476714 -1.200305e-02 44.0 -0.0011337034 1.478898 -1.664784e-03 44.2 -0.0006619511 1.483140 4.222390e-03 44.4 0.0002066172 1.488734 -2.740277e-03 44.6 0.0004998966 1.493365 4.535562e-03 44.8 0.0010891406 1.495451 5.968426e-05 45.0 -0.0011337034 1.495657 9.176474e-03 45.2 -0.0006619511 1.494814 -1.951631e-03 45.4 0.0002066172 1.493229 -6.635954e-03 45.6 0.0004998966 1.491748 4.251781e-03 45.8 0.0010891406 1.490692 -4.280912e-03 46.0 -0.0011337034 1.489840 6.993764e-03 46.2 -0.0006619511 1.489843 -2.881334e-03 46.4 0.0002066172 1.491281 -9.987875e-03 46.6 0.0004998966 1.494159 2.140784e-03 46.8 0.0010891406 1.497800 -1.988956e-03 47.0 -0.0011337034 1.500316 9.117681e-03 47.2 -0.0006619511 1.501650 6.112043e-03 47.4 0.0002066172 1.502678 -1.108475e-02 47.6 0.0004998966 1.503466 -1.665410e-03 47.8 0.0010891406 1.505032 1.278811e-03 48.0 -0.0011337034 1.506046 4.087680e-03 48.2 -0.0006619511 1.503024 9.638395e-03 48.4 0.0002066172 1.496925 9.668116e-03 48.6 0.0004998966 1.489741 -1.154137e-02 48.8 0.0010891406 1.482510 -6.199477e-03 49.0 -0.0011337034 1.477028 9.052292e-04 49.2 -0.0006619511 1.473261 4.005242e-04 49.4 0.0002066172 1.469163 6.330460e-03 49.6 0.0004998966 1.463921 2.787833e-04 49.8 0.0010891406 1.456851 -3.840030e-03 50.0 -0.0011337034 1.449222 7.911617e-03 50.2 -0.0006619511 1.442905 -7.943497e-03 50.4 0.0002066172 1.437440 -3.946909e-03 50.6 0.0004998966 1.433524 2.776031e-03 50.8 0.0010891406 1.433002 -6.190667e-03 51.0 -0.0011337034 1.434302 -5.568314e-03 51.2 -0.0006619511 1.435881 4.581028e-03 51.4 0.0002066172 1.437575 2.718775e-03 51.6 0.0004998966 1.438994 3.806467e-03 51.8 0.0010891406 1.439525 -6.814264e-03 52.0 -0.0011337034 1.439549 2.185113e-03 52.2 -0.0006619511 1.438907 6.554484e-04 52.4 0.0002066172 1.437559 6.434760e-03 52.6 0.0004998966 1.436333 -1.832991e-03 52.8 0.0010891406 1.437015 -7.703921e-03 53.0 -0.0011337034 1.439246 -1.081244e-02 53.2 -0.0006619511 1.442972 1.049021e-02 53.4 0.0002066172 1.446633 1.259911e-03 53.6 0.0004998966 1.447731 8.069498e-03 53.8 0.0010891406 1.445408 2.102770e-03 54.0 -0.0011337034 1.440424 -1.889893e-03 54.2 -0.0006619511 1.432990 4.571518e-03 54.4 0.0002066172 1.425252 2.441694e-03 54.6 0.0004998966 1.419344 -6.643790e-03 54.8 0.0010891406 1.414967 -9.655775e-03 55.0 -0.0011337034 1.411965 2.669147e-03 55.2 -0.0006619511 1.410337 5.424957e-03 55.4 0.0002066172 1.408484 -1.904250e-04 55.6 0.0004998966 1.405134 1.565882e-03 55.8 0.0010891406 1.401079 -2.268011e-03 56.0 -0.0011337034 1.397898 -1.641569e-04 56.2 -0.0006619511 1.395580 -3.618178e-03 56.4 0.0002066172 1.392521 9.719739e-04 56.6 0.0004998966 1.387947 9.953376e-03 56.8 0.0010891406 1.383016 5.947203e-04 57.0 -0.0011337034 1.378470 -8.236530e-03 57.2 -0.0006619511 1.374511 -6.348684e-03 57.4 0.0002066172 1.371810 3.983444e-03 57.6 0.0004998966 1.369791 3.708691e-03 57.8 0.0010891406 1.367763 2.947529e-03 58.0 -0.0011337034 1.365919 -7.585015e-03 58.2 -0.0006619511 1.364595 -3.232648e-03 58.4 0.0002066172 1.363065 1.628843e-03 58.6 0.0004998966 1.361930 1.016999e-02 58.8 0.0010891406 1.361212 -5.601059e-03 59.0 -0.0011337034 1.359605 -6.571737e-03 59.2 -0.0006619511 1.357277 5.984627e-03 59.4 0.0002066172 1.355907 1.586778e-03 59.6 0.0004998966 1.355350 -1.149516e-03 59.8 0.0010891406 1.354351 -6.540424e-03 60.0 -0.0011337034 1.354106 4.027452e-03 60.2 -0.0006619511 1.355782 -2.619982e-03 60.4 0.0002066172 1.358187 -3.593154e-03 60.6 0.0004998966 1.359963 3.637162e-03 60.8 0.0010891406 1.361574 4.137239e-03 61.0 -0.0011337034 1.362109 -2.774834e-03 61.2 -0.0006619511 1.361656 5.205677e-03 61.4 0.0002066172 1.362087 -6.593156e-03 61.6 0.0004998966 1.364108 -3.608168e-03 61.8 0.0010891406 1.366480 -1.869299e-03 62.0 -0.0011337034 1.369314 8.319456e-03 62.2 -0.0006619511 1.371603 -4.408870e-04 62.4 0.0002066172 1.371297 7.960534e-04 62.6 0.0004998966 1.368009 7.090662e-03 62.8 0.0010891406 1.363429 1.481792e-03 63.0 -0.0011337034 1.357987 -2.053650e-03 63.2 -0.0006619511 1.351229 -3.467148e-03 63.4 0.0002066172 1.345142 6.551438e-03 63.6 0.0004998966 1.341426 -8.125751e-03 63.8 0.0010891406 1.340172 -5.660981e-03 64.0 -0.0011337034 1.340681 -4.247678e-03 64.2 -0.0006619511 1.342821 4.940668e-03 64.4 0.0002066172 1.344958 3.035773e-03 64.6 0.0004998966 1.345275 2.125437e-03 64.8 0.0010891406 1.343172 2.538974e-03 65.0 -0.0011337034 1.340017 7.164454e-04 65.2 -0.0006619511 1.338372 -3.710475e-03 65.4 0.0002066172 1.339852 -1.045830e-02 65.6 0.0004998966 1.343940 -6.039767e-03 65.8 0.0010891406 1.349365 8.045948e-03 66.0 -0.0011337034 1.354187 5.246282e-03 66.2 -0.0006619511 1.356101 6.060891e-03 66.4 0.0002066172 1.354547 -3.536604e-04 66.6 0.0004998966 1.351498 1.502301e-03 66.8 0.0010891406 1.347553 -5.442518e-03 67.0 -0.0011337034 1.343283 6.450689e-03 67.2 -0.0006619511 1.339379 -1.417166e-03 67.4 0.0002066172 1.336069 -2.375567e-03 67.6 0.0004998966 1.332889 -2.289179e-03 67.8 0.0010891406 1.330344 6.672459e-04 68.0 -0.0011337034 1.328858 1.276132e-03 68.2 -0.0006619511 1.328047 -2.885124e-03 68.4 0.0002066172 1.326391 -9.977017e-04 68.6 0.0004998966 1.322647 8.353409e-03 68.8 0.0010891406 1.315917 6.793370e-03 69.0 -0.0011337034 1.305896 4.137904e-03 69.2 -0.0006619511 1.295563 -2.501175e-03 69.4 0.0002066172 1.287885 -1.539142e-02 69.6 0.0004998966 1.281818 -7.718266e-03 69.8 0.0010891406 1.277342 1.846924e-02 70.0 -0.0011337034 1.273508 -2.573806e-03 70.2 -0.0006619511 1.266979 2.283421e-03 70.4 0.0002066172 1.257873 6.201218e-04 70.6 0.0004998966 1.249985 -1.284397e-03 70.8 0.0010891406 1.243091 -9.280250e-03 71.0 -0.0011337034 1.238464 5.469942e-03 71.2 -0.0006619511 1.237332 -9.670043e-03 71.4 0.0002066172 1.236769 -3.575986e-03 71.6 0.0004998966 1.235170 1.402978e-02 71.8 0.0010891406 1.234054 8.568794e-04 72.0 -0.0011337034 1.232473 -9.039509e-03 72.2 -0.0006619511 1.230653 9.088762e-04 72.4 0.0002066172 1.229574 -4.280503e-03 72.6 0.0004998966 1.228361 9.538821e-03 72.8 0.0010891406 1.226944 2.667071e-03 73.0 -0.0011337034 1.224916 -8.281958e-03 73.2 -0.0006619511 1.220185 2.277140e-03 73.4 0.0002066172 1.214189 1.240436e-02 73.6 0.0004998966 1.208928 -3.428140e-03 73.8 0.0010891406 1.204252 -9.441378e-03 74.0 -0.0011337034 1.201300 -5.966698e-03 74.2 -0.0006619511 1.202518 -8.556737e-04 74.4 0.0002066172 1.207454 -3.160804e-03 74.6 0.0004998966 1.213540 -1.339574e-03 74.8 0.0010891406 1.219489 4.322175e-03 75.0 -0.0011337034 1.225242 1.691423e-03 75.2 -0.0006619511 1.229935 -1.573344e-03 75.4 0.0002066172 1.232686 3.407602e-03 75.6 0.0004998966 1.233448 3.252328e-03 75.8 0.0010891406 1.232799 5.212186e-03 76.0 -0.0011337034 1.231124 -4.190470e-03 76.2 -0.0006619511 1.229457 -1.694889e-03 76.4 0.0002066172 1.228367 -2.373991e-03 76.6 0.0004998966 1.227661 1.238998e-03 76.8 0.0010891406 1.227626 5.184828e-03 77.0 -0.0011337034 1.229224 -8.289826e-03 77.2 -0.0006619511 1.232976 -5.213987e-03 77.4 0.0002066172 1.238300 -5.706653e-03 77.6 0.0004998966 1.244945 9.355111e-03 77.8 0.0010891406 1.251458 5.531684e-04 78.0 -0.0011337034 1.256608 2.425538e-03 78.2 -0.0006619511 1.259332 -1.970193e-03 78.4 0.0002066172 1.260249 5.544552e-03 78.6 0.0004998966 1.260834 2.366411e-03 78.8 0.0010891406 1.262669 -6.558191e-03 79.0 -0.0011337034 1.266649 -8.615328e-03 79.2 -0.0006619511 1.273336 -2.374207e-03 79.4 0.0002066172 1.280884 1.708979e-03 79.6 0.0004998966 1.286425 1.307483e-02 79.8 0.0010891406 1.288845 5.765600e-03 80.0 -0.0011337034 1.288692 -3.157882e-03 80.2 -0.0006619511 1.287392 -5.029650e-03 80.4 0.0002066172 1.287579 -2.785333e-03 80.6 0.0004998966 1.290377 -1.176882e-03 80.8 0.0010891406 1.294291 -2.280021e-03 81.0 -0.0011337034 1.298152 6.281642e-03 81.2 -0.0006619511 1.301247 -1.384840e-03 81.4 0.0002066172 1.304064 2.629438e-03 81.6 0.0004998966 1.307563 -5.262733e-03 81.8 0.0010891406 1.311390 -5.179514e-03 82.0 -0.0011337034 1.314623 8.610211e-03 82.2 -0.0006619511 1.317591 3.671013e-03 82.4 0.0002066172 1.319542 -1.348767e-03 82.6 0.0004998966 1.318642 -1.541935e-03 82.8 0.0010891406 1.314675 9.535929e-03 83.0 -0.0011337034 1.307870 6.563865e-03 83.2 -0.0006619511 1.299533 2.728866e-03 83.4 0.0002066172 1.291449 -1.265568e-02 83.6 0.0004998966 1.285856 -6.456168e-03 83.8 0.0010891406 1.283689 -2.778159e-03 84.0 -0.0011337034 1.283471 3.662248e-03 84.2 -0.0006619511 1.282338 6.323816e-03 84.4 0.0002066172 1.279364 4.029828e-03 84.6 0.0004998966 1.274699 -4.098823e-03 84.8 0.0010891406 1.269891 -5.799064e-04 85.0 -0.0011337034 1.267075 -4.841494e-03 85.2 -0.0006619511 1.266476 -4.513904e-03 85.4 0.0002066172 1.266969 2.123931e-03 85.6 0.0004998966 1.268595 2.204778e-03 85.8 0.0010891406 1.271310 -2.399027e-03 86.0 -0.0011337034 1.274075 -4.941584e-03 86.2 -0.0006619511 1.276930 3.732148e-03 86.4 0.0002066172 1.279770 1.823518e-03 86.6 0.0004998966 1.280848 2.052240e-03 86.8 0.0010891406 1.279541 6.769813e-03 87.0 -0.0011337034 1.277245 -1.711064e-03 87.2 -0.0006619511 1.275044 -4.682214e-03 87.4 0.0002066172 1.274221 -2.927434e-03 87.6 0.0004998966 1.276401 -4.400632e-03 87.8 0.0010891406 1.281720 -2.709564e-03 88.0 -0.0011337034 1.288679 -2.545725e-03 88.2 -0.0006619511 1.295471 4.091257e-03 88.4 0.0002066172 1.301261 6.332614e-03 88.6 0.0004998966 1.306851 -1.351383e-03 88.8 0.0010891406 1.312973 -6.662159e-03 89.0 -0.0011337034 1.319210 -6.076746e-03 89.2 -0.0006619511 1.326276 1.078557e-02 89.4 0.0002066172 1.333798 -1.704192e-03 89.6 0.0004998966 1.340278 4.220279e-04 89.8 0.0010891406 1.345945 6.657766e-04 90.0 -0.0011337034 1.352045 -4.911324e-03 90.2 -0.0006619511 1.358079 3.683307e-03 90.4 0.0002066172 1.363585 1.008821e-03 90.6 0.0004998966 1.369073 3.026967e-03 90.8 0.0010891406 1.374528 -5.117378e-03 91.0 -0.0011337034 1.379853 -7.196567e-04 91.2 -0.0006619511 1.384289 1.973293e-03 91.4 0.0002066172 1.387760 9.033438e-03 91.6 0.0004998966 1.389717 -2.817368e-03 91.8 0.0010891406 1.391348 1.162943e-03 92.0 -0.0011337034 1.394124 -9.690338e-03 92.2 -0.0006619511 1.397175 -7.132468e-04 92.4 0.0002066172 1.398371 1.152189e-02 92.6 0.0004998966 1.397897 1.050291e-02 92.8 0.0010891406 1.396008 -7.496979e-03 93.0 -0.0011337034 1.393644 -6.610720e-03 93.2 -0.0006619511 1.392508 -5.745941e-03 93.4 0.0002066172 1.393732 7.661541e-03 93.6 0.0004998966 1.394812 -1.911919e-03 93.8 0.0010891406 1.393687 8.324075e-03 94.0 -0.0011337034 1.391158 1.175497e-03 94.2 -0.0006619511 1.388829 -7.866795e-03 94.4 0.0002066172 1.387775 -2.281649e-03 94.6 0.0004998966 1.389511 -4.310898e-03 94.8 0.0010891406 1.394141 -2.630586e-03 95.0 -0.0011337034 1.400393 2.541069e-03 95.2 -0.0006619511 1.404935 -2.873019e-03 95.4 0.0002066172 1.405802 1.839126e-02 95.6 0.0004998966 1.403610 4.290260e-03 95.8 0.0010891406 1.397796 -7.185052e-03 96.0 -0.0011337034 1.389085 6.548465e-03 96.2 -0.0006619511 1.381180 -3.518085e-03 96.4 0.0002066172 1.374674 -4.880958e-03 96.6 0.0004998966 1.368522 2.077863e-03 96.8 0.0010891406 1.363703 -2.192417e-03 97.0 -0.0011337034 1.361379 9.542330e-04 97.2 -0.0006619511 1.360789 -1.202670e-02 97.4 0.0002066172 1.360613 3.880640e-03 97.6 0.0004998966 1.359230 7.670563e-03 97.8 0.0010891406 1.355624 7.987170e-03 98.0 -0.0011337034 1.349045 1.689084e-03 98.2 -0.0006619511 1.340196 -5.633587e-03 98.4 0.0002066172 1.330720 1.173241e-03 98.6 0.0004998966 1.321169 8.313340e-04 98.8 0.0010891406 1.311724 1.786513e-03 99.0 -0.0011337034 1.302343 -1.409096e-03 > m$win s t l 4911 9 5 > m$deg s t l 0 1 1 > m$jump s t l 492 1 1 > m$inner [1] 2 > m$outer [1] 0 > postscript(file="/var/www/html/rcomp/tmp/148d61292335871.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/rcomp/tmp/2ficr1292335871.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/rcomp/tmp/3ficr1292335871.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/rcomp/tmp/4qrtc1292335871.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/rcomp/tmp/5eaq51292335871.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/rcomp/tmp/6717q1292335871.tab") > > try(system("convert tmp/148d61292335871.ps tmp/148d61292335871.png",intern=TRUE)) character(0) > try(system("convert tmp/2ficr1292335871.ps tmp/2ficr1292335871.png",intern=TRUE)) character(0) > try(system("convert tmp/3ficr1292335871.ps tmp/3ficr1292335871.png",intern=TRUE)) character(0) > try(system("convert tmp/4qrtc1292335871.ps tmp/4qrtc1292335871.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.510 1.031 15.727