Home » date » 2009 » Jun » 05 »

Nick Vermeulen Exonential smoothing model 2

*Unverified author*
R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Fri, 05 Jun 2009 09:48:16 -0600
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl.htm/, Retrieved Fri, 05 Jun 2009 17:48:51 +0200
 
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/Jun/05/t1244216931hfmc38psqc5xkkl.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:
maandelijkse verkoop auto's
 
Dataseries X:
» Textbox « » Textfile « » CSV «
14.620 16.005 16.683 15.487 15.684 15.962 12.000 13.769 14.031 16.078 15.827 13.149 15.969 16.628 16.670 16.487 16.883 16.201 12.168 14.010 16.556 17.404 16.435 13.123 16.744 17.410 16.484 17.103 17.301 17.301 12.843 13.748 16.904 17.342 15.476 15.424 15.988 19.244 18.715 17.780 17.160 17.349 11.171 13.438 16.713 18.369 17.067 14.055 15.500 18.475 19.423 18.686 19.646 19.733 12.605 16.616 19.156 21.348 20.049 18.020 20.262 21.789 20.603 21.928 21.025 19.346 11.786 19.082 20.127 20.217 20.385 16.653 13.065 20.275 21.776 20.260 22.523 23.033 14.133 20.110 19.682 22.197 17.212 11.784 15.467 17.002 15.952 18.767 20.605 19.809 14.233 19.311 20.827 23.388 20.181 14.344
 
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.422497061600848
beta0.0909939654833242
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
316.68317.39-0.707000000000001
415.48718.4491141865388-2.96211418653884
515.68418.4415717146605-2.75757171466052
615.96218.4144338655929-2.45243386559287
71218.4219328184999-6.42193281849994
813.76916.505440956064-2.73644095606400
914.03116.0408569708472-2.00985697084721
1016.07815.80598427019060.272015729809430
1115.82716.545653639023-0.718653639022998
1213.14916.8391396991379-3.69013969913790
1315.96915.73531537915930.233684620840727
1416.62816.29827923578910.329720764210869
1516.6716.9144940910722-0.244494091072188
1616.48717.2787053594341-0.791705359434147
1716.88317.3812846132337-0.498284613233697
1816.20117.5886768760961-1.38767687609613
1912.16817.3669547236250-5.19895472362505
2014.0115.3351067131389-1.32510671313893
2116.55614.88900479659051.66699520340947
2217.40415.77114425001861.63285574998144
2316.43516.7016345064469-0.266634506446856
2413.12316.8193450319068-3.69634503190685
2516.74415.34590812505991.39809187494006
2617.4116.07860504085981.33139495914022
2716.48416.8343077629217-0.350307762921705
2817.10316.86602855551430.236971444485693
2917.30117.15498337964860.146016620351379
3017.30117.4111236205358-0.110123620535756
3112.84317.5548116946036-4.71181169460362
3213.74815.5731559718304-1.82515597183041
3316.90414.74093626694642.16306373305364
3417.34215.67688596800721.66511403199276
3515.47616.4664681647305-0.990468164730478
3615.42416.0959964518442-0.671996451844244
3715.98815.83424341126020.153756588739805
3819.24415.92727972725033.31672027274969
3918.71517.48416916530341.23083083469664
4017.7818.2070953462707-0.427095346270672
4117.1618.2131330422129-1.05313304221290
4217.34917.9141842851612-0.565184285161223
4311.17117.7996641134734-6.62866411347339
4413.43814.8685046404258-1.43050464042579
4516.71314.07855697298482.63444302701516
4618.36915.10731807775973.26168192224026
4717.06716.52648010036720.540519899632798
4814.05516.8167392805015-2.76173928050145
4915.515.6056294691062-0.105629469106244
5018.47515.51265735685522.96234264314479
5119.42316.82978077107552.59321922892450
5218.68618.09064611872920.595353881270832
5319.64618.53030741878721.11569258121283
5419.73319.23270273826280.500297261737202
5512.60519.6943291131904-7.08932911319039
5616.61616.6768136353391-0.0608136353391231
5719.15616.62648733336102.52951266663903
5821.34817.76781259527093.58018740472914
5920.04919.49068401672110.55831598327892
6018.0219.9582879230864-1.93828792308639
6120.26219.29656715018760.965432849812373
6221.78919.89877563145551.89022436854452
6320.60320.9643748887044-0.361374888704443
6421.92821.06478713280820.863212867191827
6521.02521.7157700505258-0.690770050525792
6619.34621.6836433160837-2.33764331608373
6711.78620.8658475097668-9.07984750976678
6819.08216.85041838334032.23158161665973
6920.12717.69982727326582.42717272673416
7020.21718.7251847186361.49181528136400
7120.38519.41270875743720.972291242562768
7216.65319.9181148450569-3.26511484505695
7313.06518.50770300654-5.44270300654
7420.27515.96802357645184.30697642354823
7521.77617.71313540059064.06286459940941
7620.2619.51130623808790.748693761912065
7722.52319.93803292947022.58496707052975
7823.03321.23995793773521.79304206226475
7914.13322.2762298907614-8.14322989076142
8020.1118.8013922482941.30860775170598
8119.68219.37023724673550.311762753264453
8222.19719.5299037826972.66709621730300
8317.21220.7872274550166-3.57522745501657
8411.78419.2697392319153-7.4857392319153
8515.46715.8122844016269-0.345284401626868
8617.00215.35837640624961.64362359375035
8715.95215.80796478281270.144035217187284
8818.76715.62951886493853.13748113506152
8920.60516.83641451920893.76858548079106
9019.80918.45483197943221.35416802056776
9114.23319.1052257175672-4.8722257175672
9219.31116.93767522355852.37332477644154
9320.82717.92259024173052.90440975826947
9423.38819.24354621676054.14445378323948
9520.18121.2477493606671-1.06674936066707
9614.34421.0092236478139-6.66522364781392


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
9718.149116588661612.372043191243423.9261899860799
9818.105046935623211.743790531206524.4663033400399
9918.060977282584811.076081997009325.0458725681603
10018.016907629546410.372120569826225.6616946892666
10117.9728379765089.634597215837526.3110787371785
10217.92876832346968.8657728826719326.9917637642672
10317.88469867043128.0675551884582427.7018421524041
10417.84062901739287.2415642866619328.4396937481236
10517.79655936435436.3891870881552529.2039316405534
10617.75248971131595.5116209424537529.9933584801781
10717.70842005827754.6099084673183130.8069316492367
10817.66435040523913.684965221337131.6437355891411
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/11zb91244216891.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/11zb91244216891.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/2jm7j1244216891.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/2jm7j1244216891.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/35qap1244216891.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jun/05/t1244216931hfmc38psqc5xkkl/35qap1244216891.ps (open in new window)


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