Home » date » 2009 » Aug » 19 »

Klassieke decompositie - Gem. prijs gebakken tong of forel - Niels Braspennincx

*Unverified author*

R Software Module: rwasp_decompose.wasp (opens new window with default values)

Title produced by software: Classical Decomposition

Date of computation: Tue, 18 Aug 2009 16:30:39 -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/Aug/19/t1250634701f2ga5wk7wsr1a6a.htm/, Retrieved Wed, 19 Aug 2009 00:31:41 +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/Aug/19/t1250634701f2ga5wk7wsr1a6a.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 «

15,22 15,27 15,31 15,33 15,42 15,49 15,65 15,67 15,69 15,83 15,92 15,99 15,94 15,96 16,03 16,09 16,04 16,23 16,2 16,2 16,26 16,28 16,27 16,29 16,3 16,37 16,39 16,42 16,43 16,37 16,37 16,39 16,48 16,51 16,5 16,54 16,52 16,56 16,61 16,75 16,75 16,79 16,82 16,84 17,14 17,25 17,28 17,3 17,34 17,44 17,48 17,55 17,59 17,66 17,67 17,64 17,68 17,72 17,78 17,83 17,88 18,11 18,16 18,27 18,29 18,35 18,35 18,38 18,41 18,41 18,42 18,43 18,48 18,54 18,65 18,66 18,69 18,72 18,72 18,73 18,84 18,83 18,91 18,91 18,94 18,97 19 19,08 19,18 19,24 19,23 19,25 19,3 19,33 19,35 19,35 19,31 19,47 19,7 19,76 19,9 19,97 20,1 20,26 20,44 20,43 20,57 20,6 20,69 20,93 20,98 21,11 21,14 21,16 21,32 21,32 21,48 21,58 21,74 21,75 21,81 21,89 22,21 22,37 22,47 22,51 22,55 22,61 22,58 22,85 22,93 22,98

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 seconds
R Server	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

Classical Decomposition by Moving Averages
t	Observations	Fit	Trend	Seasonal	Random
1	15.22	NA	NA	-0.0659131944444443	NA
2	15.27	NA	NA	-0.0205798611111119	NA
3	15.31	NA	NA	0.0187951388888892	NA
4	15.33	NA	NA	0.0458368055555562	NA
5	15.42	NA	NA	0.0293784722222218	NA
6	15.49	NA	NA	0.0230451388888882	NA
7	15.65	15.5967118055556	15.5958333333333	0.000878472222221809	0.0532881944444448
8	15.67	15.6254201388889	15.6545833333333	-0.0291631944444429	0.0445798611111083
9	15.69	15.7318368055556	15.7133333333333	0.0185034722222223	-0.0418368055555582
10	15.83	15.7804201388889	15.775	0.00542013888888879	0.0495798611111127
11	15.92	15.8362118055556	15.8325	0.00371180555555594	0.0837881944444447
12	15.99	15.8592534722222	15.8891666666667	-0.0299131944444451	0.130746527777781
13	15.94	15.8770034722222	15.9429166666667	-0.0659131944444443	0.0629965277777771
14	15.96	15.9673368055556	15.9879166666667	-0.0205798611111119	-0.00733680555554983
15	16.03	16.0525451388889	16.03375	0.0187951388888892	-0.0225451388888871
16	16.09	16.1220868055556	16.07625	0.0458368055555562	-0.0320868055555543
17	16.04	16.1389618055556	16.1095833333333	0.0293784722222218	-0.0989618055555539
18	16.23	16.1597118055556	16.1366666666667	0.0230451388888882	0.0702881944444478
19	16.2	16.1650451388889	16.1641666666667	0.000878472222221809	0.0349548611111103
20	16.2	16.1670868055556	16.19625	-0.0291631944444429	0.0329131944444470
21	16.26	16.2468368055556	16.2283333333333	0.0185034722222223	0.0131631944444486
22	16.28	16.2625034722222	16.2570833333333	0.00542013888888879	0.0174965277777801
23	16.27	16.2907951388889	16.2870833333333	0.00371180555555594	-0.0207951388888930
24	16.29	16.2792534722222	16.3091666666667	-0.0299131944444451	0.0107465277777798
25	16.3	16.2561701388889	16.3220833333333	-0.0659131944444443	0.0438298611111101
26	16.37	16.3165034722222	16.3370833333333	-0.0205798611111119	0.0534965277777779
27	16.39	16.3729618055556	16.3541666666667	0.0187951388888892	0.0170381944444422
28	16.42	16.4187534722222	16.3729166666667	0.0458368055555562	0.00124652777778067
29	16.43	16.4214618055556	16.3920833333333	0.0293784722222218	0.0085381944444407
30	16.37	16.4351284722222	16.4120833333333	0.0230451388888882	-0.0651284722222236
31	16.37	16.4325451388889	16.4316666666667	0.000878472222221809	-0.0625451388888898
32	16.39	16.4195868055556	16.44875	-0.0291631944444429	-0.0295868055555530
33	16.48	16.4843368055556	16.4658333333333	0.0185034722222223	-0.00433680555555327
34	16.51	16.4941701388889	16.48875	0.00542013888888879	0.0158298611111150
35	16.5	16.5195451388889	16.5158333333333	0.00371180555555594	-0.019545138888887
36	16.54	16.5167534722222	16.5466666666667	-0.0299131944444451	0.0232465277777791
37	16.52	16.5170034722222	16.5829166666667	-0.0659131944444443	0.00299652777777837
38	16.56	16.5998368055556	16.6204166666667	-0.0205798611111119	-0.0398368055555558
39	16.61	16.6854618055556	16.6666666666667	0.0187951388888892	-0.075461805555559
40	16.75	16.7708368055556	16.725	0.0458368055555562	-0.0208368055555574
41	16.75	16.8177118055556	16.7883333333333	0.0293784722222218	-0.0677118055555574
42	16.79	16.8755451388889	16.8525	0.0230451388888882	-0.0855451388888895
43	16.82	16.9192118055556	16.9183333333333	0.000878472222221809	-0.099211805555555
44	16.84	16.9600034722222	16.9891666666667	-0.0291631944444429	-0.120003472222223
45	17.14	17.0805868055556	17.0620833333333	0.0185034722222223	0.0594131944444456
46	17.25	17.1370868055556	17.1316666666667	0.00542013888888879	0.112913194444449
47	17.28	17.2037118055556	17.2	0.00371180555555594	0.0762881944444445
48	17.3	17.2413368055556	17.27125	-0.0299131944444451	0.0586631944444491
49	17.34	17.2770034722222	17.3429166666667	-0.0659131944444443	0.0629965277777806
50	17.44	17.3910868055556	17.4116666666667	-0.0205798611111119	0.0489131944444452
51	17.48	17.4862951388889	17.4675	0.0187951388888892	-0.00629513888888766
52	17.55	17.5554201388889	17.5095833333333	0.0458368055555562	-0.00542013888889059
53	17.59	17.5793784722222	17.55	0.0293784722222218	0.0106215277777757
54	17.66	17.6159618055556	17.5929166666667	0.0230451388888882	0.0440381944444432
55	17.67	17.6383784722222	17.6375	0.000878472222221809	0.0316215277777765
56	17.64	17.6587534722222	17.6879166666667	-0.0291631944444429	-0.0187534722222189
57	17.68	17.7626701388889	17.7441666666667	0.0185034722222223	-0.0826701388888829
58	17.72	17.8079201388889	17.8025	0.00542013888888879	-0.0879201388888866
59	17.78	17.8653784722222	17.8616666666667	0.00371180555555594	-0.0853784722222244
60	17.83	17.8896701388889	17.9195833333333	-0.0299131944444451	-0.0596701388888867
61	17.88	17.9107534722222	17.9766666666667	-0.0659131944444443	-0.0307534722222229
62	18.11	18.0152534722222	18.0358333333333	-0.0205798611111119	0.0947465277777795
63	18.16	18.1158784722222	18.0970833333333	0.0187951388888892	0.0441215277777758
64	18.27	18.2020868055556	18.15625	0.0458368055555562	0.0679131944444435
65	18.29	18.2410451388889	18.2116666666667	0.0293784722222218	0.0489548611111061
66	18.35	18.2863784722222	18.2633333333333	0.0230451388888882	0.0636215277777801
67	18.35	18.3142118055556	18.3133333333333	0.000878472222221809	0.0357881944444465
68	18.38	18.3270868055556	18.35625	-0.0291631944444429	0.0529131944444465
69	18.41	18.4130868055556	18.3945833333333	0.0185034722222223	-0.00308680555555441
70	18.41	18.4366701388889	18.43125	0.00542013888888879	-0.0266701388888855
71	18.42	18.4678784722222	18.4641666666667	0.00371180555555594	-0.047878472222223
72	18.43	18.4663368055556	18.49625	-0.0299131944444451	-0.0363368055555569
73	18.48	18.4611701388889	18.5270833333333	-0.0659131944444443	0.0188298611111115
74	18.54	18.5365034722222	18.5570833333333	-0.0205798611111119	0.0034965277777772
75	18.65	18.6083784722222	18.5895833333333	0.0187951388888892	0.0416215277777745
76	18.66	18.6708368055556	18.625	0.0458368055555562	-0.0108368055555523
77	18.69	18.6922951388889	18.6629166666667	0.0293784722222218	-0.00229513888888988
78	18.72	18.7263784722222	18.7033333333333	0.0230451388888882	-0.00637847222222732
79	18.72	18.7433784722222	18.7425	0.000878472222221809	-0.0233784722222268
80	18.73	18.7504201388889	18.7795833333333	-0.0291631944444429	-0.0204201388888912
81	18.84	18.8305868055556	18.8120833333333	0.0185034722222223	0.00941319444444488
82	18.83	18.8495868055556	18.8441666666667	0.00542013888888879	-0.0195868055555550
83	18.91	18.8857951388889	18.8820833333333	0.00371180555555594	0.0242048611111123
84	18.91	18.8942534722222	18.9241666666667	-0.0299131944444451	0.0157465277777824
85	18.94	18.9011701388889	18.9670833333333	-0.0659131944444443	0.0388298611111111
86	18.97	18.9894201388889	19.01	-0.0205798611111119	-0.0194201388888899
87	19	19.0696284722222	19.0508333333333	0.0187951388888892	-0.0696284722222238
88	19.08	19.1366701388889	19.0908333333333	0.0458368055555562	-0.0566701388888866
89	19.18	19.1593784722222	19.13	0.0293784722222218	0.0206215277777773
90	19.24	19.1897118055556	19.1666666666667	0.0230451388888882	0.0502881944444411
91	19.23	19.2012951388889	19.2004166666667	0.000878472222221809	0.0287048611111125
92	19.25	19.2075034722222	19.2366666666667	-0.0291631944444429	0.0424965277777822
93	19.3	19.3051701388889	19.2866666666667	0.0185034722222223	-0.0051701388888894
94	19.33	19.3495868055556	19.3441666666667	0.00542013888888879	-0.0195868055555550
95	19.35	19.4062118055556	19.4025	0.00371180555555594	-0.0562118055555558
96	19.35	19.4330034722222	19.4629166666667	-0.0299131944444451	-0.0830034722222202
97	19.31	19.4636701388889	19.5295833333333	-0.0659131944444443	-0.153670138888891
98	19.47	19.5873368055556	19.6079166666667	-0.0205798611111119	-0.117336805555556
99	19.7	19.7162951388889	19.6975	0.0187951388888892	-0.0162951388888892
100	19.76	19.8366701388889	19.7908333333333	0.0458368055555562	-0.0766701388888862
101	19.9	19.9168784722222	19.8875	0.0293784722222218	-0.0168784722222206
102	19.97	20.0134618055556	19.9904166666667	0.0230451388888882	-0.0434618055555553
103	20.1	20.1008784722222	20.1	0.000878472222221809	-0.000878472222218818
104	20.26	20.1891701388889	20.2183333333333	-0.0291631944444429	0.0708298611111111
105	20.44	20.3510034722222	20.3325	0.0185034722222223	0.0889965277777769
106	20.43	20.4475034722222	20.4420833333333	0.00542013888888879	-0.0175034722222200
107	20.57	20.5537118055556	20.55	0.00371180555555594	0.0162881944444422
108	20.6	20.6213368055556	20.65125	-0.0299131944444451	-0.0213368055555527
109	20.69	20.6857534722222	20.7516666666667	-0.0659131944444443	0.00424652777778078
110	20.93	20.8260868055556	20.8466666666667	-0.0205798611111119	0.103913194444445
111	20.98	20.9529618055556	20.9341666666667	0.0187951388888892	0.0270381944444438
112	21.11	21.0712534722222	21.0254166666667	0.0458368055555562	0.038746527777775
113	21.14	21.1514618055556	21.1220833333333	0.0293784722222218	-0.0114618055555553
114	21.16	21.2417951388889	21.21875	0.0230451388888882	-0.0817951388888929
115	21.32	21.3142118055556	21.3133333333333	0.000878472222221809	0.00578819444444179
116	21.32	21.3708368055556	21.4	-0.0291631944444429	-0.0508368055555586
117	21.48	21.5097534722222	21.49125	0.0185034722222223	-0.0297534722222217
118	21.58	21.6004201388889	21.595	0.00542013888888879	-0.0204201388888912
119	21.74	21.7066284722222	21.7029166666667	0.00371180555555594	0.0333715277777742
120	21.75	21.7846701388889	21.8145833333333	-0.0299131944444451	-0.0346701388888917
121	21.81	21.8561701388889	21.9220833333333	-0.0659131944444443	-0.0461701388888862
122	21.89	22.0065034722222	22.0270833333333	-0.0205798611111119	-0.116503472222217
123	22.21	22.1454618055556	22.1266666666667	0.0187951388888892	0.0645381944444452
124	22.37	22.2712534722222	22.2254166666667	0.0458368055555562	0.0987465277777773
125	22.47	22.3572951388889	22.3279166666667	0.0293784722222218	0.112704861111109
126	22.51	22.4517951388889	22.42875	0.0230451388888882	0.0582048611111112
127	22.55	NA	NA	0.000878472222221809	NA
128	22.61	NA	NA	-0.0291631944444429	NA
129	22.58	NA	NA	0.0185034722222223	NA
130	22.85	NA	NA	0.00542013888888879	NA
131	22.93	NA	NA	0.00371180555555594	NA
132	22.98	NA	NA	-0.0299131944444451	NA

Charts produced by software:

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/1tkic1250634634.png (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/1tkic1250634634.ps (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/28z3g1250634634.png (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/28z3g1250634634.ps (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/3hp8r1250634634.png (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/3hp8r1250634634.ps (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/4n07u1250634634.png (open in new window)

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Aug/19/t1250634701f2ga5wk7wsr1a6a/4n07u1250634634.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')

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.

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