Home » date » 2008 » Aug » 13 »

Raf Mattheussen exponential smoothing roze garnalen

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Wed, 13 Aug 2008 07:54:47 -0600
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj.htm/, Retrieved Wed, 13 Aug 2008 13:56:51 +0000
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
16,8 16,91 16,91 17,16 17,02 17,23 17,22 17,29 17,3 17,22 17,19 17,23 17,36 17,39 17,29 17,28 17,4 17,51 17,54 17,64 17,65 17,5 17,37 17,56 17,49 17,61 17,79 17,83 17,56 17,95 18,09 18,38 18,38 18,44 18,84 19,01 19,06 19,06 18,97 18,98 19,41 19,55 19,64 19,71 19,48 19,48 19,41 19,25 19,14 19,21 19,3 19,53 19,14 19,16 19,24 19,38 19,27 19,27 19,07 19,15 19,24 19,36 19,57 19,59 19,36 19,46 19,65 19,46 19,51 19,64 19,64 19,69 19,28 19,67 19,65 19,6 19,53 19,64 19,67 19,81 19,73 19,87 19,97 20,12 19,94 20,31 20,13 20,22 20,38 20,44 20,34 20,14 19,97 19,82 19,98 20,12
 
Text written by user:
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.914696945069777
beta0
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
316.9116.910
417.1616.910.25
517.0217.1386742362674-0.118674236267445
617.2317.03012327489510.199876725104875
717.2217.21294990473910.00705009526089384
817.2917.21939860533670.0706013946633028
917.317.28397748535290.0160225146471156
1017.2217.2986332305529-0.0786332305529385
1117.1917.2267076547852-0.0367076547851966
1217.2317.19313127509250.0368687249074995
1317.3617.2268549851340.133145014865992
1417.3917.34864232348320.0413576765167996
1517.2917.3864720638483-0.096472063848303
1617.2817.2982293617617-0.0182293617616835
1717.417.28155502024770.118444979752297
1817.5117.38989628138600.120103718614022
1917.5417.49975478589370.0402452141062533
2017.6417.53656696029040.103433039709586
2117.6517.63117684573210.0188231542679453
2217.517.6483943274375-0.148394327437519
2317.3717.5126584894647-0.142658489464736
2417.5617.38216920496310.177830795036925
2517.4917.5448304899227-0.0548304899226792
2617.6117.49467720829370.115322791706276
2717.7917.60016261356440.189837386435627
2817.8317.77380629099710.0561937090029261
2917.5617.8252065049542-0.265206504954190
3017.9517.58262292506000.367377074940041
3118.0917.91866161319630.171338386803715
3218.3818.07538431217880.304615687821173
3318.3818.35401535124920.0259846487508177
3418.4418.37778343008030.0622165699197375
3518.8418.43469273651860.405307263481433
3619.0118.80542605223960.204573947760377
3719.0618.99254921729690.0674507827030943
3819.0619.0542462421780.00575375782200993
3918.9719.0595091868805-0.089509186880452
4018.9818.97763540708520.00236459291478042
4119.4118.97979829300070.430201706999298
4219.5519.37330248015680.176697519843238
4319.6419.53492716175880.105072838241220
4419.7119.63103696590780.0789630340921654
4519.4819.7032642119654-0.223264211965379
4619.4819.4990451193372-0.0190451193372354
4719.4119.4816246068610-0.0716246068609756
4819.2519.4161097977734-0.166109797773416
4919.1419.2641696732039-0.124169673203912
5019.2119.15059205245400.0594079475460205
5119.319.20493232058720.0950676794128107
5219.5319.29189043652100.238109563479039
5319.1419.5096885268271-0.369688526827137
5419.1619.171535560711-0.0115355607110104
5519.2419.1609840185690.0790159814310165
5619.3819.23325969539560.146740304604375
5719.2719.3674826037359-0.0974826037358554
5819.2719.2783155639012-0.00831556390122046
5919.0719.2707093430042-0.200709343004242
6019.1519.08712112011130.0628788798887001
6119.2419.14463623945490.0953637605450979
6219.3619.23186517989590.128134820104133
6319.5719.34906970840220.220930291597817
6419.5919.55115397120010.0388460287999202
6519.3619.5866863150715-0.226686315071461
6619.4619.37933703518650.0806629648135306
6719.6519.45311920268170.196880797318322
6819.4619.6332054665316-0.173205466531648
6919.5119.47477495542580.0352250445742364
7019.6419.50699519608780.133004803912236
7119.6419.62865428390590.0113457160941088
7219.6919.63903217575680.0509678242431981
7319.2819.6856522888889-0.405652288888909
7419.6719.31460337948170.355396620518338
7519.6519.63968358255790.0103164174420876
7619.619.6491199780763-0.0491199780762521
7719.5319.604190084188-0.074190084188011
7819.6419.53632864082680.103671359173230
7919.6719.63115651635380.0388434836462466
8019.8119.66668653218080.143313467819155
8119.7319.7977749233824-0.0677749233823803
8219.8719.73578140801220.134218591987821
8319.9719.8585507440750.111449255924992
8420.1219.96049303799990.159506962000105
8519.9420.1063935688588-0.166393568858751
8620.3119.95419387974440.355806120255604
8720.1320.2796486509793-0.149648650979326
8820.2220.14276548709470.077234512905278
8920.3820.21341166010310.166588339896869
9020.4420.36578950569100.0742104943089608
9120.3420.4336696181276-0.0936696181275636
9220.1420.3479903045804-0.207990304580427
9319.9720.1577422083766-0.18774220837658
9419.8219.9860149839139-0.166014983913868
9519.9819.83416158529200.145838414707956
9620.1219.96755953769920.152440462300770


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
9720.106996362870819.794622926709020.4193697990326
9820.106996362870819.683655958441320.5303367673002
9920.106996362870819.596254713477120.6177380122645
10020.106996362870819.521763958754020.6922287669876
10120.106996362870819.455738425313620.758254300428
10220.106996362870819.395816473086020.8181762526555
10320.106996362870819.340565172197320.8734275535442
10420.106996362870819.289037498110420.9249552276311
10520.106996362870819.240568847816620.973423877925
10620.106996362870819.194671545031221.0193211807103
10720.106996362870819.150975177096721.0630175486449
10820.106996362870819.109190549071421.1048021766701
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/1nz4y1218635680.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/1nz4y1218635680.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/24ggx1218635680.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/24ggx1218635680.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/3v6ng1218635680.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Aug/13/t1218635811hpmgacyij5zn9cj/3v6ng1218635680.ps (open in new window)


 
Parameters (Session):
par1 = 12 ; par2 = Double ; par3 = additive ;
 
Parameters (R input):
par1 = 12 ; par2 = Double ; 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