Home » date » 2008 » May » 26 »

Exponential Smoothing eigen reeks - Jesse De Meyer

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Mon, 26 May 2008 12:59:43 -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/26/t1211828426low7uv6dl673mnt.htm/, Retrieved Mon, 26 May 2008 21:00:26 +0200
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
36409 33163 34122 35225 28249 30374 26311 22069 23651 28628 23187 14727 43080 32519 39657 33614 28671 34243 27336 22916 24537 26128 22602 15744 41086 39690 43129 37863 35953 29133 24693 22205 21725 27192 21790 13253 37702 30364 32609 30212 29965 28352 25814 22414 20506 28806 22228 13971 36845 35338 35022 34777 26887 23970 22780 17351 21382 24561 17409 11514 31514 27071 29462 26105 22397 23843 21705 18089 20764 25316 17704 15548 28029 29383 36438 32034 22679 24319 18004 17537 20366 22782 19169 13807 29743 25591 29096 26482 22405 27044 17970 18730 19684 19785 18479 10698 31956 29506 34506 27165 26736 23691 18157 17328 18205 20995 17382 9367
 
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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.51705592571119
beta0
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
33412233163959
43522533658.8566327571566.14336724297
52824934468.6403413033-6219.64034130328
63037431252.7384470401-878.73844704005
72631130798.3815258477-4487.38152584774
82206928478.1543169812-6409.15431698124
92365125164.2630985886-1513.26309858864
102862824381.82144630334246.1785536967
112318726577.3332291200-3390.33322911995
121472724824.3413428679-10097.3413428679
134308019603.451167609523476.5488323905
143251931742.1398566451776.860143354876
153965732143.81999721567513.1800027844
163361436028.5542385901-2414.5542385901
172867134780.094661576-6109.09466157602
183424331621.35106607752621.64893392246
192733632976.8901824966-5640.89018249657
202291630060.2344873506-7144.23448735064
212453726366.2657109957-1829.26571099574
222612825420.4330354251707.5669645749
232260225786.2847272960-3184.28472729603
241574424139.831439896-8395.83143989598
254108619798.717042625421287.2829573746
263969030805.43283802688884.5671619732
274312935399.25093650417729.7490634959
283786339395.9634940452-1532.96349404518
293595338603.3356355502-2650.33563555019
302913337232.9638900654-8099.96389006543
312469333044.8295626604-8351.82956266044
322220528726.4665967570-6521.46659675696
332172525354.5036485762-3629.50364857618
342719223477.84727968953714.15272031052
352179025398.2719527224-3608.27195272237
361325323532.5935579898-10279.5935579898
373770218217.468794928619484.5312050714
383036428292.06111421542071.93888578464
393260929363.36939282183245.63060717825
403021231041.5419309329-829.541930932872
412996530612.6223599181-647.622359918128
422835230277.7653810994-1925.76538109939
432581429282.0369792725-3468.03697927249
442241427488.8679085541-5074.86790855411
452050624864.8773842346-4358.87738423465
462880622611.09400326766194.90599673237
472222825814.2068581019-3586.20685810189
481397123959.9373512942-9988.9373512942
493684518795.098102249718049.9018977503
503533828127.90683698717210.09316301285
513502231855.92823185273166.07176814732
523477733492.96440080021284.03559919984
532688734156.8826161906-7269.88261619055
542397030397.9467302645-6427.94673026446
552278027074.3387832253-4294.33878322535
561735124853.9254683473-7502.9254683473
572138220974.4933947689407.506605231083
582456121185.19709977013375.8029002299
591740922930.675993367-5521.67599336699
601151420075.6607011394-8561.66070113936
613151415648.803301686615865.1966983134
622707123851.99726712323219.00273287683
632946225516.40170503763945.59829496235
642610527556.4966839239-1451.49668392390
652239726805.9917223509-4408.99172235091
662384324526.2964258978-683.296425897785
672170524172.9939598701-2467.99395987006
681808922896.9030582998-4807.90305829982
692076420410.9482917609353.051708239058
702531620593.49576958844722.5042304116
711770423035.2945661189-5331.29456611889
721554820278.7171189952-4730.71711899525
732802917832.671799755410196.3282002446
742938323104.7437161886278.25628381202
753643826350.953330866510087.0466691335
763203431566.5205840673467.479415932699
772267931808.2335862233-9129.2335862233
782431927087.9092632649-2768.90926326493
791800425656.2283209372-7652.22832093719
801753721699.5983227016-4162.59832270162
812036619547.3021935933818.697806406715
822278219970.61474576262811.38525423737
831916921424.2581509231-2255.25815092312
841380720258.1635599799-6451.16355997986
852974316922.551213560212820.4487864398
862559123551.44022886572039.55977113427
872909624606.00669437294489.99330562714
882648226927.5843394510-445.58433945095
892240526697.1923163337-4292.19231633373
902704424477.88884488132566.11115511867
911797025804.7118236690-7834.71182366903
921873021753.7276490014-3023.72764900143
931968420190.2913503485-506.291350348474
941978519928.5104075145-143.510407514474
951847919854.3075009079-1375.30750090789
961069819143.1966078884-8445.19660788841
973195614776.557657983717179.4423420163
982950623659.29012133695846.70987866306
993450626682.36611001387823.63388998617
1002716530727.6223734261-3562.62237342607
1012673628885.5473641749-2149.54736417485
1022369127774.1111619314-4083.11116193137
1031815725662.9143403173-7505.91434031725
1041732821781.9368527756-4453.93685277562
1051820519479.0024103045-1274.00241030453
1062099518820.27191468622174.72808531376
1071738219944.7279580083-2562.72795800827
108936718619.6542813344-9252.65428133435


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
10913835.5145566134199.38010308887927471.6490101379
11013835.5145566134-1515.5683167100029186.5974299368
11113835.5145566134-3057.3042707075130728.3333839343
11213835.5145566134-4469.6462042685632140.6753174954
11313835.5145566134-5780.5628553152533451.5919685421
11413835.5145566134-7009.1988158420234680.2279290688
11513835.5145566134-8169.3407741736435840.3698874005
11613835.5145566134-9271.3078272175636942.3369404444
11713835.5145566134-10323.062041087037994.0911543138
11813835.5145566134-11330.899684911039001.9287981378
11913835.5145566134-12299.901803663839970.9309168906
12013835.5145566134-13234.239244891540905.2683581183
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/26/t1211828426low7uv6dl673mnt/15axu1211828377.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/26/t1211828426low7uv6dl673mnt/15axu1211828377.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/26/t1211828426low7uv6dl673mnt/27d3x1211828377.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/26/t1211828426low7uv6dl673mnt/27d3x1211828377.ps (open in new window)


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