Home » date » 2010 » Aug » 16 »

Exponential_Smoothing_Reeks_A_Jeroen_Kinne

*Unverified author*
R Software Module: /rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Mon, 16 Aug 2010 14:20:22 +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/Aug/16/t1281968478dc0cvsounltxsx0.htm/, Retrieved Mon, 16 Aug 2010 16:21:22 +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/2010/Aug/16/t1281968478dc0cvsounltxsx0.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:
Jeroen_Kinne
 
Dataseries X:
» Textbox « » Textfile « » CSV «
36 35 34 32 52 51 36 26 27 27 28 30 28 29 25 28 55 53 42 32 37 41 37 38 39 32 36 39 83 83 66 53 72 77 69 72 81 71 63 66 114 116 109 97 111 120 110 106 115 110 103 112 163 166 156 140 166 176 163 162 171 167 163 168 222 216 197 178 204 220 196 195 213 218 216 225 280 272 252 230 248 259 240 237 252 250 255 255 313 291 271 247 268 283 259 259 267 270 279 269 334 326 301 276 301 313 291 287 289 298 320 312 385 380 351 322 350 363 344 345
 
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'Gwilym Jenkins' @ 72.249.127.135


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.124791631525621
beta0.0903911823166393
gamma1


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
132828.7814054374197-0.781405437419746
142929.4996210015149-0.499621001514875
152525.0751391320068-0.0751391320067611
162827.48421870194200.515781298057966
175553.0885562974241.91144370257596
185350.80919476501922.19080523498079
194237.63996100340244.36003899659757
203228.32682261777093.67317738222907
213730.78919256731536.21080743268472
224132.49842234355618.5015776564439
233735.3690130880351.63098691196502
243838.4449775716736-0.444977571673597
253935.13593308337443.86406691662564
263237.2471740287663-5.24717402876632
273631.75292923961074.24707076038930
283936.33769003008892.66230996991113
298372.25836704559910.7416329544010
308371.141481359444111.8585186405559
316657.25833154375388.74166845624622
325344.15419723892968.84580276107038
337251.527489087505720.4725109124943
347758.644130420186718.3558695798133
356955.232603417589313.7673965824107
367259.23743729398512.7625627060150
378162.363640461465918.6363595385341
387154.748892659417216.2511073405828
396363.8486403153705-0.848640315370545
406669.4028479886932-3.40284798869324
41114145.936784916613-31.9367849166129
42116140.301520865075-24.3015208650754
43109107.7194219404631.28057805953746
449784.871548280509712.1284517194903
45111112.196601480667-1.19660148066666
46120115.3197864863984.68021351360176
47110100.5091321717979.49086782820338
48106103.0023076243482.99769237565154
49115111.6151439069723.38485609302802
5011094.108700003878815.8912999961212
5110384.83495336609918.1650466339010
5211291.393647089650820.6063529103492
53163166.495244554399-3.49524455439897
54166172.614676298980-6.61467629898044
55156161.279978360779-5.27997836077907
56140140.462656150105-0.462656150105460
57166160.7247897865335.27521021346718
58176173.4874663192852.51253368071482
59163157.3405032062095.65949679379102
60162151.56479153887810.4352084611216
61171165.0848614573695.91513854263061
62167155.22505358599811.7749464140016
63163142.69740530820820.3025946917918
64168153.30411743571914.6958825642812
65222225.738019612001-3.73801961200067
66216229.909230344462-13.9092303444623
67197214.720204210965-17.7202042109653
68178190.210088540833-12.2100885408333
69204222.047100034280-18.0471000342803
70220231.616305991051-11.6163059910511
71196211.21750758362-15.2175075836199
72195205.135795250163-10.1357952501627
73213212.8941663937890.105833606210808
74218204.67105396464213.3289460353577
75216196.66238952055219.3376104794482
76225201.49931227307323.5006877269274
77280269.15639009319510.8436099068049
78272263.8769733358968.12302666410437
79252243.1229428033878.87705719661264
80230221.8259581666758.17404183332474
81248257.536031178919-9.53603117891947
82259277.802840768603-18.8028407686027
83240247.255624673105-7.25562467310519
84237246.344461665583-9.34446166558308
85252267.536872966414-15.5368729664141
86250269.242449940353-19.2424499403534
87255260.480463680905-5.4804636809053
88255265.746466884705-10.7464668847052
89313325.740518194468-12.7405181944683
90291311.916455303031-20.9164553030314
91271283.436620072139-12.4366200721386
92247254.309660286664-7.30966028666396
93268272.464183474588-4.46418347458797
94283284.370560790463-1.37056079046289
95259262.656050995746-3.65605099574577
96259258.5915861673050.40841383269526
97267275.497912500129-8.49791250012868
98270273.299518289933-3.29951828993308
99279277.8327484972631.16725150273680
100269278.235774979787-9.23577497978715
101334340.411527037191-6.41152703719052
102326317.2207229347818.77927706521916
103301297.2623783374083.73762166259178
104276271.7660736866914.23392631330893
105301295.570494472625.42950552738017
106313312.6284742282030.371525771796712
107291286.3236242155274.67637578447273
108287286.6111606989010.388839301098983
109289296.411979002331-7.41197900233067
110298299.032979654439-1.03297965443949
111320308.50426102737911.4957389726211
112312299.98112342008812.0188765799124
113385375.3772975953489.6227024046521
114380366.61770374416113.3822962558390
115351339.83677411009511.1632258899053
116322312.6378395135149.36216048648572
117350341.8790566346348.1209433653662
118363356.9699972653056.03000273469485
119344332.37530147953611.6246985204639
120345329.70185570165015.2981442983503


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
121335.636009681172314.261896244569357.010123117774
122347.004585361403325.424263598509368.584907124297
123371.717489459965349.853301219648393.581677700281
124361.26754001249339.187271515152383.347808509828
125444.974623900517422.121586314335467.8276614867
126437.671737480202414.513484180969460.829990779434
127402.894824386870379.669632770995426.120016002746
128368.356574187279345.101706443548391.611441931009
129399.218236369251375.215519815881423.22095292262
130413.09177712158388.451803952472437.731750290688
131389.625832274864364.893737164886414.357927384841
132388.234357368062372.560490417341403.908224318783
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/18wca1281968419.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/18wca1281968419.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/28wca1281968419.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/28wca1281968419.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/3j5cu1281968419.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Aug/16/t1281968478dc0cvsounltxsx0/3j5cu1281968419.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=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