Home » date » 2009 » Jan » 16 »

Opgave 10-2 - Consumptieprijzen pakketreizen - Chloë De Rijck

*Unverified author*
R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Fri, 16 Jan 2009 07:07:31 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo.htm/, Retrieved Fri, 16 Jan 2009 15:09:20 +0100
 
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/Jan/16/t1232114947ew5hztryhjohrfo.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 «
105,6 110,2 104,9 102,9 102,6 103,6 107,8 106,6 106 105,2 107,9 107,5 107,5 113,3 107,8 104,5 105,1 104,2 106,6 103,8 107,7 106,4 110 113,2 113,9 112 113,9 113,1 111,7 110,7 113,5 114 112,7 112,2 115,8 118,4 118,8 123,9 118 120,2 118,7 119,8 124,8 121,3 120,2 118,3 129,6 130,2 127,19 133,1 129,12 123,28 123,36 124,13 126,96 127,14 123,7 123,67 130,19 134,01 124,96 129,96 128,32 132,38 126,25 128,91 131,42 129,44 126,86 126,71 131,63 132,78 126,61 132,84 123,14 128,13 125,49 126,48 130,86 127,32 126,56 126,64 129,26 126,47 135,38 135,5 132,22 122,62 125,16 128,5 133,86 128,87 125,07 125,25 132,16 130,24
 
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


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.208653452026953
beta0.135928514801665
gamma0.144824151716476


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
13107.5106.5763087606840.923691239316184
14113.3112.8448735955540.455126404446062
15107.8107.6077457008720.192254299128066
16104.5104.3545546748650.145445325135029
17105.1104.9790552417800.120944758220134
18104.2103.9143739104810.285626089518516
19106.6109.008821488077-2.40882148807660
20103.8107.172744393442-3.37274439344178
21107.7105.5982171696392.10178283036056
22106.4105.0880796169201.31192038307984
23110107.9670098299662.03299017003442
24113.2107.9957199832425.20428001675758
25113.9109.4937615703374.40623842966316
26112116.732364194594-4.73236419459447
27113.9110.5326870464083.36731295359249
28113.1108.1766345546094.92336544539111
29111.7110.1707840629691.52921593703110
30110.7109.8342855096650.86571449033454
31113.5115.172891517895-1.67289151789539
32114113.8326908178590.167309182141238
33112.7114.177418506156-1.47741850615607
34112.2113.281628386918-1.08162838691824
35115.8116.127572854054-0.327572854054125
36118.4116.3440398185162.05596018148393
37118.8117.3212679713731.47873202862695
38123.9123.0462217353550.853778264645186
39118119.243343752074-1.24334375207384
40120.2116.2757865943173.92421340568275
41118.7117.8163231398090.883676860191187
42119.8117.3946343439262.40536565607441
43124.8122.9327667054501.86723329454969
44121.3122.811740874061-1.51174087406100
45120.2122.839629412424-2.63962941242441
46118.3121.935733616598-3.63573361659770
47129.6124.4517688362185.14823116378203
48130.2126.3558473873623.84415261263807
49127.19127.9626496001-0.772649600100024
50133.1133.404983124011-0.30498312401059
51129.12129.345879697500-0.225879697500488
52123.28127.437607809812-4.15760780981222
53123.36126.968918556119-3.60891855611925
54124.13125.682348349915-1.55234834991494
55126.96130.118887597700-3.15888759769963
56127.14128.205213479842-1.06521347984237
57123.7127.852993205802-4.15299320580236
58123.67126.132231869226-2.46223186922607
59130.19129.5461636116340.643836388366424
60134.01129.8795213116874.13047868831259
61124.96130.543665569088-5.58366556908798
62129.96134.42602440002-4.46602440002013
63128.32128.780018790534-0.460018790534349
64132.38125.6379000435356.7420999564651
65126.25127.081099647821-0.83109964782112
66128.91126.2633548680172.64664513198268
67131.42131.1645295212250.255470478774669
68129.44130.072678784374-0.632678784374292
69126.86129.338555499994-2.47855549999376
70126.71128.090158747334-1.38015874733365
71131.63132.045748662242-0.415748662241555
72132.78132.4874643656550.292535634345285
73126.61131.058510508721-4.44851050872084
74132.84135.159009132355-2.31900913235461
75123.14130.334189151014-7.19418915101356
76128.13126.3354658017791.79453419822060
77125.49125.4611778877890.0288221122108894
78126.48124.8285925238271.65140747617325
79130.86128.8270003686892.03299963131099
80127.32127.633600256416-0.313600256415839
81126.56126.3929051286360.167094871364441
82126.64125.5358512703261.10414872967354
83129.26129.904223453292-0.64422345329163
84126.47130.156856472641-3.68685647264132
85135.38127.0187882815878.36121171841295
86135.5134.0639945886301.43600541136965
87132.22129.5983066605672.62169333943251
88122.62129.090621782214-6.47062178221364
89125.16126.467761450972-1.30776145097161
90128.5125.8826856319612.61731436803927
91133.86130.2941952276503.56580477234968
92128.87129.362988061814-0.49298806181406
93125.07128.346168357144-3.27616835714363
94125.25126.986614866068-1.73661486606794
95132.16130.5898652531571.57013474684348
96130.24131.046618437654-0.806618437653555


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
97130.062784015298124.369193041398135.756374989198
98134.505063640771128.653735231408140.356392050133
99129.770247209067123.726985301821135.813509116313
100127.493757833897121.223347830058133.764167837736
101126.816485376104120.283394375182133.349576377026
102126.994981823671120.163958754473133.82600489287
103130.935716094095123.772268005701138.099164182490
104128.660829449433121.131568694351136.190090204515
105127.307386475152119.380258013804135.234514936500
106126.780260054706118.424666996682135.135853112730
107131.146489594850122.333338676713139.959640512988
108130.980361074005121.682049433254140.278672714757
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/17k3f1232114848.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/17k3f1232114848.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/2h0f91232114848.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/2h0f91232114848.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/3skj31232114848.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/16/t1232114947ew5hztryhjohrfo/3skj31232114848.ps (open in new window)


 
Parameters (Session):
par1 = 12 ; par2 = Triple ; par3 = additive ;
 
Parameters (R input):
par1 = 12 ; par2 = Triple ; par3 = additive ;
 
R code (references can be found in the software module):
par1 <- as.numeric(par1)
if (par2 == 'Single') K <- 1
if (par2 == 'Double') K <- 2
if (par2 == 'Triple') K <- par1
nx <- length(x)
nxmK <- nx - K
x <- ts(x, frequency = par1)
if (par2 == 'Single') fit <- HoltWinters(x, gamma=0, beta=0)
if (par2 == 'Double') fit <- HoltWinters(x, gamma=0)
if (par2 == 'Triple') fit <- HoltWinters(x, seasonal=par3)
fit
myresid <- x - fit$fitted[,'xhat']
bitmap(file='test1.png')
op <- par(mfrow=c(2,1))
plot(fit,ylab='Observed (black) / Fitted (red)',main='Interpolation Fit of Exponential Smoothing')
plot(myresid,ylab='Residuals',main='Interpolation Prediction Errors')
par(op)
dev.off()
bitmap(file='test2.png')
p <- predict(fit, par1, prediction.interval=TRUE)
np <- length(p[,1])
plot(fit,p,ylab='Observed (black) / Fitted (red)',main='Extrapolation Fit of Exponential Smoothing')
dev.off()
bitmap(file='test3.png')
op <- par(mfrow = c(2,2))
acf(as.numeric(myresid),lag.max = nx/2,main='Residual ACF')
spectrum(myresid,main='Residals Periodogram')
cpgram(myresid,main='Residal Cumulative Periodogram')
qqnorm(myresid,main='Residual Normal QQ Plot')
qqline(myresid)
par(op)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Estimated Parameters of Exponential Smoothing',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'Value',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,fit$alpha)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,fit$beta)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'gamma',header=TRUE)
a<-table.element(a,fit$gamma)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Interpolation Forecasts of Exponential Smoothing',4,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,'Residuals',header=TRUE)
a<-table.row.end(a)
for (i in 1:nxmK) {
a<-table.row.start(a)
a<-table.element(a,i+K,header=TRUE)
a<-table.element(a,x[i+K])
a<-table.element(a,fit$fitted[i,'xhat'])
a<-table.element(a,myresid[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Extrapolation Forecasts of Exponential Smoothing',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t',header=TRUE)
a<-table.element(a,'Forecast',header=TRUE)
a<-table.element(a,'95% Lower Bound',header=TRUE)
a<-table.element(a,'95% Upper Bound',header=TRUE)
a<-table.row.end(a)
for (i in 1:np) {
a<-table.row.start(a)
a<-table.element(a,nx+i,header=TRUE)
a<-table.element(a,p[i,'fit'])
a<-table.element(a,p[i,'lwr'])
a<-table.element(a,p[i,'upr'])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.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