Home » date » 2008 » Jun » 01 »

exponential smoothing model tafelbier Maarten Verhaegen

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Sun, 01 Jun 2008 07:11:49 -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/Jun/01/t12123260607ihpumv9n7hrpvj.htm/, Retrieved Sun, 01 Jun 2008 13:14:20 +0000
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
0,63 0,64 0,65 0,65 0,65 0,65 0,65 0,66 0,65 0,66 0,66 0,66 0,66 0,68 0,69 0,7 0,71 0,71 0,7 0,7 0,7 0,7 0,71 0,7 0,7 0,7 0,69 0,7 0,69 0,69 0,69 0,7 0,7 0,71 0,71 0,71 0,72 0,73 0,74 0,74 0,74 0,74 0,75 0,75 0,76 0,76 0,76 0,76 0,76 0,77 0,77 0,78 0,78 0,78 0,78 0,78 0,78 0,78 0,8 0,8 0,8 0,81 0,81 0,81 0,8 0,81 0,81 0,81 0,8 0,82 0,83 0,83
 
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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.692780535122199
beta0.000961035668292444
gamma1


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
130.660.6429313445268870.0170686554731130
140.680.6729687817564330.0070312182435669
150.690.6861514079916560.00384859200834398
160.70.6971259599479190.00287404005208125
170.710.7074153638721790.00258463612782078
180.710.7075271654994930.00247283450050673
190.70.6980224062428790.00197759375712137
200.70.6990004956351570.000999504364843085
210.70.700963021739085-0.00096302173908469
220.70.702396296022734-0.00239629602273406
230.710.713723183028984-0.00372318302898378
240.70.704902961986993-0.00490296198699336
250.70.6913751476959470.00862485230405297
260.70.713212493024162-0.0132124930241621
270.690.711571455565236-0.0215714555652361
280.70.704671632031416-0.00467163203141574
290.690.709631505095134-0.0196315050951338
300.690.69435415404041-0.00435415404041006
310.690.6802789897077190.00972101029228079
320.70.6863361174770760.0136638825229239
330.70.6964496143082980.00355038569170241
340.710.7005432701471160.00945672985288382
350.710.719758398003111-0.00975839800311118
360.710.7063294190348160.0036705809651838
370.720.7027569331845470.0172430668154527
380.730.7239405933654880.00605940663451188
390.740.7330523906593130.0069476093406865
400.740.751888581219213-0.011888581219213
410.740.74723752220306-0.00723752220305962
420.740.745328199187258-0.00532819918725769
430.750.7342361906490550.0157638093509449
440.750.7455344225876390.0044655774123612
450.760.7458716668307570.0141283331692433
460.760.759217145650020.000782854349980378
470.760.766839818961247-0.00683981896124686
480.760.7592402859786730.000759714021327018
490.760.7574680663461260.00253193365387361
500.770.7652286498035220.00477135019647845
510.770.773881751186837-0.00388175118683709
520.780.7796496821565550.000350317843445302
530.780.785059912331307-0.00505991233130698
540.780.785340211931307-0.00534021193130729
550.780.780489099574688-0.000489099574688234
560.780.7768457284387160.00315427156128445
570.780.779114686603410.000885313396589704
580.780.7791111675915260.000888832408474038
590.80.7845146241867430.0154853758132567
600.80.794603404261130.00539659573887052
610.80.7964042503624480.00359574963755216
620.810.8058309719366830.00416902806331698
630.810.811445674255775-0.00144567425577458
640.810.820616388807891-0.0106163888078910
650.80.816825298821848-0.0168252988218480
660.810.8089104558460660.00108954415393403
670.810.809939619154376.03808456306476e-05
680.810.8076320487642730.00236795123572742
690.80.80856017343223-0.00856017343223037
700.820.8019321413478130.0180678586521871
710.830.8239785853165080.00602141468349215
720.830.8241996148526780.00580038514732162


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
730.8255636408256940.808631163890520.842496117760868
740.8328301999159460.8121831866689910.8534772131629
750.8337968977782780.8100552924292130.857538503127342
760.841272977808160.8146760117989230.867869943817398
770.8428385843161360.8137642967220040.871912871910267
780.8524806525253850.8209101051745260.884051199876244
790.8523406311952740.8187047697236060.885976492666942
800.8505176335205360.814996851489560.886038415551513
810.8461325866775210.8089242068990660.883340966455977
820.8538546191529120.8146013046278780.893107933677945
830.8598335790116650.8186798315902430.900987326433087
840.855589151959457-7.084163981916868.79534228583578
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/1zx9l1212325904.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/1zx9l1212325904.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/2fpqb1212325904.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/2fpqb1212325904.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/3ifmh1212325904.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Jun/01/t12123260607ihpumv9n7hrpvj/3ifmh1212325904.ps (open in new window)


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