Home » date » 2010 » Dec » 29 »

*The author of this computation has been verified*
R Software Module: /rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Wed, 29 Dec 2010 22:33:08 +0000
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp.htm/, Retrieved Wed, 29 Dec 2010 23:31:59 +0100
 
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/2010/Dec/29/t1293661916zjykw53n265dnvp.htm/},
    year = {2010},
}
@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 = {2010},
    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:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
206010 198112 194519 185705 180173 176142 203401 221902 197378 185001 176356 180449 180144 173666 165688 161570 156145 153730 182698 200765 176512 166618 158644 159585 163095 159044 155511 153745 150569 150605 179612 194690 189917 184128 175335 179566 181140 177876 175041 169292 166070 166972 206348 215706 202108 195411 193111 195198 198770 194163 190420 189733 186029 191531 232571 243477 227247 217859 208679 213188 216234 213586 209465 204045 200237 203666 241476 260307 243324 244460 233575 237217 235243 230354 227184 221678 217142 219452 256446 265845 248624 241114 229245 231805 219277 219313 212610 214771 211142 211457 240048 240636 230580 208795 197922 194596 194581 185686 178106 172608 167302 168053 202300 202388 182516 173476 166444 171297 169701 164182 161914 159612 151001 158114 186530 187069 174330 169362 166827 178037 186413 189226 191563 188906 186005 195309 223532 226899 etc...
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ www.wessa.org


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.786721671208313
beta0.155737227444672
gamma1


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
13180144192256.77590812-12112.7759081197
14173666174646.428106398-980.428106397972
15165688164180.8073711891507.19262881079
16161570159602.1658100851967.83418991519
17156145154188.6069701651956.39302983464
18153730152119.1230186181610.876981382
19182698179952.3484015222745.65159847779
20200765200442.561219332322.43878066825
21176512176164.427462601347.572537398606
22166618164082.6938745692535.30612543123
23158644157564.6031125531079.39688744742
24159585162699.575040561-3114.57504056091
25163095157033.5205933096061.47940669145
26159044158060.686855114983.31314488608
27155511151876.2857412733634.7142587268
28153745151536.0705339142208.92946608641
29150569148805.6999799911763.30002000884
30150605148982.9099334371622.09006656322
31179612179540.64902107171.3509789286763
32194690199556.121835498-4866.12183549814
33189917172711.69245961417205.307540386
34184128177934.6383742186193.36162578195
35175335178007.834878517-2672.83487851685
36179566182860.56067138-3294.56067137973
37181140182552.108246814-1412.10824681402
38177876179243.047588603-1367.04758860334
39175041174113.551861832927.448138168431
40169292173346.182184128-4054.18218412847
41166070166832.87427879-762.874278789735
42166972165922.4910065951049.50899340483
43206348196558.7952843899789.20471561115
44215706225216.873886494-9510.87388649379
45202108200907.00696911200.99303089993
46195411190710.8561367064700.14386329381
47193111187055.8402206396055.15977936087
48195198199049.339913742-3851.33991374171
49198770199042.997871195-272.997871194704
50194163197117.930776858-2954.93077685832
51190420191512.249043426-1092.24904342552
52189733188129.6788898931603.32111010666
53186029187498.596717085-1469.59671708476
54191531187061.5533571314469.44664286863
55232571223314.1931842789256.80681572173
56243477248433.712722975-4956.71272297547
57227247231545.874832777-4298.8748327769
58217859218649.86048071-790.860480710166
59208679211171.888635608-2492.88863560779
60213188213488.229628964-300.229628963629
61216234216634.512356744-400.512356744439
62213586213617.211916038-31.2119160378934
63209465210647.255498067-1182.25549806724
64204045207696.056936062-3651.05693606209
65200237201559.353700862-1322.35370086171
66203666201806.3581361041859.64186389593
67241476236008.6293294945467.3706705058
68260307253632.9740701646674.02592983609
69243324245978.104966728-2654.10496672837
70244460235268.2824712649191.71752873648
71233575236647.931624633-3072.93162463326
72237217240271.635200018-3054.63520001838
73235243242188.15269305-6945.15269304995
74230354234257.517648547-3903.51764854678
75227184227677.911670693-493.911670693138
76221678224508.313359515-2830.31335951484
77217142219381.134929036-2239.13492903588
78219452219340.37877907111.621220930218
79256446252477.5642737673968.43572623306
80265845268537.035194257-2692.03519425693
81248624249733.664875926-1109.66487592596
82241114241164.041968422-50.0419684217777
83229245229923.595156388-678.595156388124
84231805232994.616274153-1189.61627415262
85219277233336.866000632-14059.866000632
86219313217374.1832982211938.81670177865
87212610213750.413699221-1140.41369922119
88214771207127.0340285987643.96597140183
89211142209202.7514305291939.24856947092
90211457212298.996501291-841.996501291229
91240048244740.097239513-4692.09723951275
92240636250736.072567318-10100.072567318
93230580223704.945028596875.05497141043
94208795219884.192553618-11089.1925536177
95197922196713.5364857461208.46351425373
96194596198279.950262193-3683.95026219264
97194581190729.0900505413851.9099494589
98185686191278.927945051-5592.92794505131
99178106179159.00117226-1053.00117226044
100172608172574.58126733.4187330002605
101167302164610.4379122582691.56208774194
102168053164961.7541554063091.2458445938
103202300197414.3760947654885.62390523477
104202388208703.725659621-6315.72565962127
105182516187645.695259936-5129.69525993612
106173476168453.7592840085022.24071599171
107166444160459.7414514655984.25854853509
108171297165203.6710899516093.32891004876
109169701168613.7145560091087.28544399096
110164182166297.126878904-2115.12687890409
111161914159630.5798763512283.4201236488
112159612158060.5395547431551.46044525696
113151001154201.424991103-3200.42499110274
114158114151624.5621344896489.43786551084
115186530189171.601050889-2641.60105088935
116187069193266.148679688-6197.1486796877
117174330173684.922388921645.077611079469
118169362163039.4116295686322.58837043165
119166827158271.002348838555.99765117007
120178037167373.95077015310663.049229847
121186413176183.81580852710229.1841914731
122189226184368.8393573164857.16064268383
123191563188972.4008920332590.59910796714
124188906192372.294027709-3466.2940277087
125186005187821.724755109-1816.72475510943
126195309192839.2157924012469.78420759932
127223532229223.087110721-5691.08711072078
128226899233733.22180034-6834.2218003402
129214126218605.044575765-4479.04457576541
130206903208006.298620989-1103.29862098879
131204442199829.4202577164612.57974228362
132220375207753.52926872812621.4707312718
133214320219725.686840237-5405.68684023729
134212588214263.16287187-1675.16287187007
135205816212242.323659017-6426.32365901687
136202196205149.962306198-2953.96230619837
137195722199310.402816273-3588.40281627342
138198563201587.356456262-3024.356456262
139229139228974.239624871164.760375128622
140229527235630.86822827-6103.8682282705
141211868219452.446170661-7584.4461706613
142203555204622.969306616-1067.96930661576
143195770195189.669246087580.330753913237


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
144198652.315847294188946.127403196208358.504291392
145192306.353232724179187.823603683205424.882861766
146188010.821329436171506.694169215204514.948489657
147182618.374637994162670.698972004202566.050303985
148178433.51109528154949.581220205201917.440970355
149172255.700582968145126.928102233199384.473063704
150175388.800333884144499.251040267206278.3496275
151204118.502515189169349.186597036238887.818433341
152207571.683921272168802.85614486246340.511697684
153194890.525319291152002.96560706237778.085031523
154187357.973538483140233.710876445234482.236200521
155179187.517956563127710.236010507230664.799902619
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/1eblz1293661987.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/1eblz1293661987.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/2eblz1293661987.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/2eblz1293661987.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/3eblz1293661987.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/29/t1293661916zjykw53n265dnvp/3eblz1293661987.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=F, beta=F)
if (par2 == 'Double') fit <- HoltWinters(x, gamma=F)
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