Home » date » 2009 » Jun » 02 »

Classical Deomposition eigen reeks

*Unverified author*
R Software Module: rwasp_decompose.wasp (opens new window with default values)
Title produced by software: Classical Decomposition
Date of computation: Tue, 02 Jun 2009 08:29:02 -0600
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv.htm/, Retrieved Tue, 02 Jun 2009 16:29:43 +0200
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv.htm/},
    year = {2009},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2009},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
310,95 312,97 315,07 315,43 315,73 315,77 315,77 315,77 313,99 314,57 314,63 314,65 314,65 314,93 315,27 316,26 316,98 317,01 317,07 317,07 317 317,08 317,04 317 317,05 321,59 325,59 326,23 326,28 326,35 326,35 326,35 326,39 326,74 326,9 326,9 326,91 336,93 348,5 349,43 349,26 349,26 349,28 349,61 349,66 349,68 349,91 349,91 350,89 355,52 356,36 357,04 360,28 360,63 360,79 360,97 361 361,01 361 361 361,58 363,19 363,61 364,14 365,51 365,51 365,5 365,5 364,59 364,63 364,54 363,67 365,22 369,05 370,45 370,46 370,46 370,58 370,58 370,22 370,21 370,29 370,29 370,2 370,2 372,55 374,51 375,58 375,75 375,75 375,75 375,69 375,76 377,5 377,51 377,74
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Classical Decomposition by Moving Averages
tObservationsFitTrendSeasonalRandom
1310.95NANA-3.08173115079366NA
2312.97NANA0.0988640873015803NA
3315.07NANA2.30737599206348NA
4315.43NANA2.25797123015873NA
5315.73NANA2.27767361111110NA
6315.77NANA1.60928075396824NA
7315.77315.786304563492314.76251.02380456349208-0.0163045634920991
8315.77315.336245039683314.9983333333330.3379117063492260.433754960317401
9313.99314.339221230159315.088333333333-0.74911210317461-0.349221230158719
10314.57313.836006944444315.13125-1.295243055555540.733993055555516
11314.63313.251661706349315.217916666667-1.966254960317461.37833829365076
12314.65312.501125992063315.321666666667-2.820540674603172.14887400793651
13314.65312.345768849206315.4275-3.081731150793662.30423115079367
14314.93315.634697420635315.5358333333330.0988640873015803-0.704697420634943
15315.27318.02279265873315.7154166666672.30737599206348-2.75279265873013
16316.26318.203387896825315.9454166666672.25797123015873-1.9433878968253
17316.98318.428090277778316.1504166666672.27767361111110-1.44809027777774
18317.01317.958030753968316.348751.60928075396824-0.948030753968226
19317.07317.570471230159316.5466666666671.02380456349208-0.500471230158723
20317.07317.262078373016316.9241666666670.337911706349226-0.192078373015875
21317316.882554563492317.631666666667-0.749112103174610.117445436507978
22317.08317.181840277778318.477083333333-1.29524305555554-0.101840277777853
23317.04317.313745039683319.28-1.96625496031746-0.273745039682467
24317317.236125992063320.056666666667-2.82054067460317-0.236125992063478
25317.05317.750768849206320.8325-3.08173115079366-0.700768849206383
26321.59321.704697420635321.6058333333330.0988640873015803-0.114697420634968
27325.59324.691125992064322.383752.307375992063480.898874007936456
28326.23325.435471230159323.17752.257971230158730.794528769841293
29326.28326.268506944444323.9908333333332.277673611111100.0114930555556043
30326.35326.423447420635324.8141666666671.60928075396824-0.0734474206348636
31326.35326.661304563492325.63751.02380456349208-0.311304563492058
32326.35327.025411706349326.68750.337911706349226-0.675411706349166
33326.39327.532137896825328.28125-0.74911210317461-1.14213789682537
34326.74328.907256944444330.2025-1.29524305555554-2.16725694444438
35326.9330.160411706349332.126666666667-1.96625496031746-3.26041170634926
36326.9331.218209325397334.03875-2.82054067460317-4.31820932539682
37326.91332.867018849206335.94875-3.08173115079366-5.95701884920635
38336.93337.972197420635337.8733333333330.0988640873015803-1.04219742063492
39348.5342.119459325397339.8120833333332.307375992063486.38054067460314
40349.43343.995471230159341.73752.257971230158735.43452876984128
41349.26345.929756944444343.6520833333332.277673611111103.33024305555557
42349.26347.178864087302345.5695833333331.609280753968242.08113591269847
43349.28348.551304563492347.52751.023804563492080.728695436507962
44349.61349.639161706349349.301250.337911706349226-0.029161706349214
45349.66349.654221230159350.403333333333-0.749112103174610.00577876984129944
46349.68349.752673611111351.047916666667-1.29524305555554-0.0726736111110995
47349.91349.857911706349351.824166666667-1.966254960317460.0520882936507974
48349.91349.93654265873352.757083333333-2.82054067460317-0.0265426587301363
49350.89350.628685515873353.710416666667-3.081731150793660.261314484126956
50355.52354.762197420635354.6633333333330.09886408730158030.757802579365034
51356.36357.91654265873355.6091666666672.30737599206348-1.55654265873017
52357.04358.811721230159356.553752.25797123015873-1.77172123015873
53360.28359.765590277778357.4879166666672.277673611111100.514409722222183
54360.63360.021364087302358.4120833333331.609280753968240.60863591269839
55360.79360.343387896825359.3195833333331.023804563492080.446612103174573
56360.97360.422495039683360.0845833333330.3379117063492260.54750496031744
57361359.957137896825360.70625-0.749112103174611.04286210317457
58361.01360.008923611111361.304166666667-1.295243055555541.00107638888886
59361359.851661706349361.817916666667-1.966254960317461.14833829365074
60361359.418625992063362.239166666667-2.820540674603171.58137400793657
61361.58359.557018849206362.63875-3.081731150793662.02298115079367
62363.19363.122614087301363.023750.09886408730158030.0673859126985121
63363.61365.669459325397363.3620833333332.30737599206348-2.0594593253968
64364.14365.920471230159363.66252.25797123015873-1.78047123015870
65365.51366.238506944444363.9608333333332.27767361111110-0.728506944444405
66365.51365.828864087302364.2195833333331.60928075396824-0.31886408730162
67365.5365.506304563492364.48251.02380456349208-0.00630456349210817
68365.5365.216245039683364.8783333333330.3379117063492260.283754960317424
69364.59364.658387896825365.4075-0.74911210317461-0.0683878968254135
70364.63364.660590277778365.955833333333-1.29524305555554-0.0305902777777192
71364.54364.459161706349366.425416666667-1.966254960317460.0808382936509133
72363.67364.022375992063366.842916666667-2.82054067460317-0.352375992063401
73365.22364.18410218254367.265833333333-3.081731150793661.03589781746035
74369.05367.773030753968367.6741666666670.09886408730158031.27696924603174
75370.45370.412375992064368.1052.307375992063480.037624007936472
76370.46370.832971230159368.5752.25797123015873-0.372971230158782
77370.46371.328090277778369.0504166666672.27767361111110-0.868090277777753
78370.58371.171364087302369.5620833333331.60928075396824-0.591364087301542
79370.58371.065471230159370.0416666666671.02380456349208-0.485471230158737
80370.22370.732911706349370.3950.337911706349226-0.512911706349257
81370.21369.960887896825370.71-0.749112103174610.249112103174582
82370.29369.797256944444371.0925-1.295243055555540.492743055555593
83370.29369.559995039683371.52625-1.966254960317460.730004960317444
84370.2369.14154265873371.962083333333-2.820540674603171.05845734126979
85370.2369.311185515873372.392916666667-3.081731150793660.888814484126897
86372.55372.935114087302372.836250.0988640873015803-0.385114087301645
87374.51375.60279265873373.2954166666672.30737599206348-1.09279265873022
88375.58376.085054563492373.8270833333332.25797123015873-0.505054563492081
89375.75376.706006944444374.4283333333332.27767361111110-0.956006944444425
90375.75376.652614087302375.0433333333331.60928075396824-0.902614087301572
91375.75NANA1.02380456349208NA
92375.69NANA0.337911706349226NA
93375.76NANA-0.74911210317461NA
94377.5NANA-1.29524305555554NA
95377.51NANA-1.96625496031746NA
96377.74NANA-2.82054067460317NA
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/1od6v1243952940.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/1od6v1243952940.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/2q5vt1243952940.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/2q5vt1243952940.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/3a5yy1243952940.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/3a5yy1243952940.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/47qjh1243952940.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/02/t1243952978sof3y4e4zk8l0jv/47qjh1243952940.ps (open in new window)


 
Parameters (Session):
par1 = additive ; par2 = 12 ;
 
Parameters (R input):
par1 = additive ; par2 = 12 ;
 
R code (references can be found in the software module):
par2 <- as.numeric(par2)
x <- ts(x,freq=par2)
m <- decompose(x,type=par1)
m$figure
bitmap(file='test1.png')
plot(m)
dev.off()
mylagmax <- length(x)/2
bitmap(file='test2.png')
op <- par(mfrow = c(2,2))
acf(as.numeric(x),lag.max = mylagmax,main='Observed')
acf(as.numeric(m$trend),na.action=na.pass,lag.max = mylagmax,main='Trend')
acf(as.numeric(m$seasonal),na.action=na.pass,lag.max = mylagmax,main='Seasonal')
acf(as.numeric(m$random),na.action=na.pass,lag.max = mylagmax,main='Random')
par(op)
dev.off()
bitmap(file='test3.png')
op <- par(mfrow = c(2,2))
spectrum(as.numeric(x),main='Observed')
spectrum(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
spectrum(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
spectrum(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
bitmap(file='test4.png')
op <- par(mfrow = c(2,2))
cpgram(as.numeric(x),main='Observed')
cpgram(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
cpgram(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
cpgram(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Classical Decomposition by Moving Averages',6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t',header=TRUE)
a<-table.element(a,'Observations',header=TRUE)
a<-table.element(a,'Fit',header=TRUE)
a<-table.element(a,'Trend',header=TRUE)
a<-table.element(a,'Seasonal',header=TRUE)
a<-table.element(a,'Random',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(m$trend)) {
a<-table.row.start(a)
a<-table.element(a,i,header=TRUE)
a<-table.element(a,x[i])
if (par1 == 'additive') a<-table.element(a,m$trend[i]+m$seasonal[i]) else a<-table.element(a,m$trend[i]*m$seasonal[i])
a<-table.element(a,m$trend[i])
a<-table.element(a,m$seasonal[i])
a<-table.element(a,m$random[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
 





Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.


FreeStatistics.org is powered by