Home » date » 2008 » May » 25 »

Tom Knaepen - Goudwaarde te Brussel EUR/kg - Exponential Smoothing

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Sun, 25 May 2008 10:00:33 -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/May/25/t12117313075a1xx7vflq3pqlx.htm/, Retrieved Sun, 25 May 2008 18:01:47 +0200
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
10236 10893 10756 10940 10997 10827 10166 10186 10457 10368 10244 10511 10812 10738 10171 9721 9897 9828 9924 10371 10846 10413 10709 10662 10570 10297 10635 10872 10296 10383 10431 10574 10653 10805 10872 10625 10407 10463 10556 10646 10702 11353 11346 11451 11964 12574 13031 13812 14544 14931 14886 16005 17064 15168 16050 15839 15137 14954 15648 15305 15579 16348 15928 16171 15937 15713 15594 15683 16438 17032 17696 17745
 
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 time6 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.847259609453347
beta0.0285189422376599
gamma1


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
131081210914.545-102.545000000002
141073810758.5325100920-20.5325100919545
151017110152.218097716218.7819022838121
1697219702.5003584747818.4996415252208
1798979875.8238103002621.1761896997414
1898289802.5100057877325.4899942122684
1999249927.92536033625-3.92536033624674
201037110403.9900914676-32.9900914675727
211084610853.2156453249-7.21564532493721
221041310349.771161759963.2288382401257
231070910639.247572283469.7524277165885
241066210617.769914286544.2300857134705
251057011023.6991726782-453.699172678183
261029710576.4797557018-279.479755701808
27106359744.30297144189890.697028558112
281087210041.876900282830.123099717992
291029610931.4724699731-635.4724699731
301038310314.806630017168.193369982906
311043110485.2827160569-54.2827160568722
321057410926.3983518576-352.398351857648
331065311113.3771581826-460.37715818264
341080510230.2353995896574.764600410377
351087210959.9604811899-87.9604811899335
361062510802.9986097014-177.998609701406
371040710941.2567172338-534.256717233828
381046310447.116136778115.8838632218649
391055610045.7807801648510.219219835191
401064610004.4041825719641.59581742813
411070210498.5222248473203.477775152729
421135310708.5244120152644.475587984793
431134611370.859951866-24.8599518660067
441145111814.3868015072-363.386801507182
451196411998.3140926721-34.3140926720862
461257411667.3125361300906.687463870016
471303112618.1040195344412.89598046557
481381212924.9137617263887.086238273743
491454413990.0644512583553.935548741722
501493114607.1320573339323.867942666113
511488614654.8840954089231.115904591148
521600514502.99703430611502.00296569392
531706415785.87094270921278.12905729083
541516817126.3926391424-1958.39263914245
551605015570.9482258911479.051774108901
561583916491.6481334537-652.64813345372
571513716575.7049980647-1438.70499806466
581495415259.5607950411-305.560795041072
591564815139.5619948300508.438005170025
601530515633.7778245057-328.777824505683
611557915622.5407604476-43.5407604476204
621634815688.4637410486659.536258951406
631592816004.7712741020-76.771274102035
641617115777.0244373676393.975562632419
651593716051.0284935064-114.028493506383
661571315648.156709631664.843290368377
671559416158.5749061471-564.57490614714
681568315976.3390119373-293.339011937251
691643816207.5865352134230.413464786594
701703216481.8518848546550.148115145385
711769617235.0236337073460.976366292685
721774517583.8361052263161.16389477365


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
7318065.7981719416958.401088161719173.1952557182
7418311.575892362716842.707571924319780.444212801
7517976.260905439416203.923621161819748.5981897170
7617906.956157766615862.799559241619951.1127562915
7717781.543117319715485.634501603820077.4517330356
7817517.334503668314983.469482541620051.1995247950
7917889.839700335115127.986174389320651.6932262808
8018254.179473225415271.778224070521236.5807223803
8118847.852862487915650.58045297922045.1252719969
8219004.060528222315596.309389055722411.8116673890
8319292.526614841215677.723703562722907.3295261197
8419208.873151084015389.700515280823028.0457868871
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/1pb601211731224.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/1pb601211731224.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/20m5p1211731224.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/20m5p1211731224.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/3nl6h1211731224.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/25/t12117313075a1xx7vflq3pqlx/3nl6h1211731224.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