Opgave 9; oefening 2: Weuts Koen

*Unverified author*

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

Title produced by software: Classical Decomposition

Date of computation: Thu, 12 May 2011 13:18:29 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz.htm/, Retrieved Thu, 12 May 2011 15:15:09 +0200

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

KDGP2W92

Dataseries X:

» Textbox « » Textfile « » CSV «

190,2 180,7 193,6 192,8 195,5 197,2 196,9 178,9 172,4 156,4 143,7 153,6 168,8 185,8 199,9 205,4 197,5 199,6 200,5 193,7 179,6 169,1 169,8 195,5 194,8 204,5 203,8 204,8 204,9 240 248,3 258,4 254,9 288,3 333,6 346,3 357,5 490,7 468,2 471,2 517,1 609,2 682 614 554,2 406,8 348,6 298,8 313,7 282,1 232,9 239,3 241,9 265,7 276 271,5 254,6 269,9 293,5 306,1 365,4 347,9 352,1 377,9 377,4 372,2 362,5 341,9 354,8 369,2 406,7 454,7

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 2 seconds R Server 'Herman Ole Andreas Wold' @ www.yougetit.org Classical Decomposition by Moving Averages t Observations Fit Trend Seasonal Random 1 190.2 NA NA -12.8691666666667 NA 2 180.7 NA NA 6.55249999999999 NA 3 193.6 NA NA -7.14583333333334 NA 4 192.8 NA NA -2.09916666666666 NA 5 195.5 NA NA 1.97583333333333 NA 6 197.2 NA NA 26.855 NA 7 196.9 223.565833333333 178.433333333333 45.1325 -26.6658333333333 8 178.9 202.593333333333 177.754166666667 24.8391666666667 -23.6933333333333 9 172.4 180.194166666667 178.229166666667 1.96500000000002 -7.79416666666663 10 156.4 153.078333333333 179.016666666667 -25.9383333333334 3.32166666666672 11 143.7 150.368333333333 179.625 -29.2566666666667 -6.66833333333335 12 153.6 149.7975 179.808333333333 -30.0108333333333 3.80250000000001 13 168.8 167.189166666667 180.058333333333 -12.8691666666667 1.61083333333337 14 185.8 187.3775 180.825 6.55249999999999 -1.57749999999996 15 199.9 174.595833333333 181.741666666667 -7.14583333333334 25.3041666666667 16 205.4 180.471666666667 182.570833333333 -2.09916666666666 24.9283333333333 17 197.5 186.163333333333 184.1875 1.97583333333333 11.3366666666667 18 199.6 213.875833333333 187.020833333333 26.855 -14.2758333333334 19 200.5 234.9825 189.85 45.1325 -34.4825 20 193.7 216.551666666667 191.7125 24.8391666666667 -22.8516666666667 21 179.6 194.619166666667 192.654166666667 1.96500000000002 -15.0191666666667 22 169.1 166.853333333333 192.791666666667 -25.9383333333334 2.24666666666667 23 169.8 163.818333333333 193.075 -29.2566666666667 5.98166666666668 24 195.5 165.055833333333 195.066666666667 -30.0108333333333 30.4441666666666 25 194.8 185.8725 198.741666666667 -12.8691666666667 8.92750000000007 26 204.5 209.981666666667 203.429166666667 6.55249999999999 -5.48166666666665 27 203.8 202.116666666667 209.2625 -7.14583333333334 1.68333333333337 28 204.8 215.2675 217.366666666667 -2.09916666666666 -10.4675 29 204.9 231.134166666667 229.158333333333 1.97583333333333 -26.2341666666666 30 240 269.121666666667 242.266666666667 26.855 -29.1216666666666 31 248.3 300.461666666667 255.329166666667 45.1325 -52.1616666666667 32 258.4 298.8725 274.033333333333 24.8391666666667 -40.4725 33 254.9 298.94 296.975 1.96500000000002 -44.04 34 288.3 293.153333333333 319.091666666667 -25.9383333333334 -4.8533333333333 35 333.6 313.943333333333 343.2 -29.2566666666667 19.6566666666667 36 346.3 341.580833333333 371.591666666667 -30.0108333333333 4.71916666666669 37 357.5 392.176666666667 405.045833333333 -12.8691666666667 -34.6766666666667 38 490.7 444.485833333333 437.933333333333 6.55249999999999 46.2141666666666 39 468.2 458.075 465.220833333333 -7.14583333333334 10.125 40 471.2 480.53 482.629166666667 -2.09916666666666 -9.33000000000004 41 517.1 490.1675 488.191666666667 1.97583333333333 26.9325000000001 42 609.2 513.6925 486.8375 26.855 95.5075000000001 43 682 528.165833333333 483.033333333333 45.1325 153.834166666667 44 614 497.355833333333 472.516666666667 24.8391666666667 116.644166666667 45 554.2 455.985833333333 454.020833333333 1.96500000000002 98.2141666666668 46 406.8 408.615833333333 434.554166666667 -25.9383333333334 -1.81583333333327 47 348.6 384.168333333333 413.425 -29.2566666666667 -35.5683333333333 48 298.8 357.635 387.645833333333 -30.0108333333333 -58.835 49 313.7 343.5475 356.416666666667 -12.8691666666667 -29.8475000000001 50 282.1 331.781666666667 325.229166666667 6.55249999999999 -49.6816666666667 51 232.9 291.329166666667 298.475 -7.14583333333334 -58.4291666666667 52 239.3 278.188333333333 280.2875 -2.09916666666666 -38.8883333333333 53 241.9 274.263333333333 272.2875 1.97583333333333 -32.3633333333333 54 265.7 297.150833333333 270.295833333333 26.855 -31.4508333333333 55 276 317.886666666667 272.754166666667 45.1325 -41.8866666666667 56 271.5 302.489166666667 277.65 24.8391666666667 -30.9891666666666 57 254.6 287.323333333333 285.358333333333 1.96500000000002 -32.7233333333333 58 269.9 270.161666666667 296.1 -25.9383333333334 -0.261666666666713 59 293.5 278.264166666667 307.520833333333 -29.2566666666667 15.2358333333333 60 306.1 287.593333333333 317.604166666667 -30.0108333333333 18.5066666666667 61 365.4 312.776666666667 325.645833333333 -12.8691666666667 52.6233333333334 62 347.9 338.735833333333 332.183333333333 6.55249999999999 9.16416666666669 63 352.1 332.145833333333 339.291666666667 -7.14583333333334 19.9541666666667 64 377.9 345.505 347.604166666667 -2.09916666666666 32.395 65 377.4 358.434166666667 356.458333333333 1.97583333333333 18.9658333333333 66 372.2 394.221666666667 367.366666666667 26.855 -22.0216666666667 67 362.5 NA NA 45.1325 NA 68 341.9 NA NA 24.8391666666667 NA 69 354.8 NA NA 1.96500000000002 NA 70 369.2 NA NA -25.9383333333334 NA 71 406.7 NA NA -29.2566666666667 NA 72 454.7 NA NA -30.0108333333333 NA

Charts produced by software:

http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/1tu151305206307.png (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/1tu151305206307.ps (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/2e48p1305206307.png (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/2e48p1305206307.ps (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/3qiae1305206307.png (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/3qiae1305206307.ps (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/47py91305206307.png (open in new window) http://www.freestatistics.org/blog/date/2011/May/12/t1305206105xu4wxp08csvl9tz/47py91305206307.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

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