Home » date » 2010 » Aug » 16 »

Aantal bezoekers

*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 12:20:14 +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/t12819611841umy7yd4oejoziy.htm/, Retrieved Mon, 16 Aug 2010 14:19:44 +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/t12819611841umy7yd4oejoziy.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 «
556 555 554 552 572 571 556 546 547 547 548 550 555 549 555 550 566 573 543 535 542 541 535 536 548 546 548 548 561 563 527 527 541 534 522 527 539 533 532 519 538 542 503 502 522 511 492 500 509 511 505 493 518 518 474 471 483 461 439 446 461 449 441 424 447 448 404 403 411 386 359 370 385 369 368 352 378 383 334 323 330 303 275 284 301 281 284 272 297 300 240 236 247 218 192 201 223 197 195 175 197 204 142 142 151 127 100 114 139 112 123 108 132 140 76 71 81 57 38 46
 
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.312700997165616
beta0.158746465680972
gamma1


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
13555556.096125968013-1.09612596801264
14549550.206838318706-1.20683831870599
15555555.909186031654-0.909186031654144
16550550.436065507093-0.436065507092508
17566566.440247225086-0.440247225086296
18573573.766014542637-0.766014542637436
19543550.225015114608-7.22501511460837
20535537.330053349687-2.33005334968720
21542536.6047280020035.3952719979975
22541537.4184025855073.58159741449333
23535539.109847151401-4.10984715140148
24536539.004461470953-3.00446147095295
25548541.5657069316396.43429306836106
26546537.8152689377338.18473106226702
27548546.7753467402641.22465325973565
28548542.6824894474845.31751055251561
29561560.9425504645390.0574495354612736
30563568.79555489801-5.79555489801044
31527539.90321780154-12.9032178015398
32527528.78717895982-1.78717895982038
33541533.5864776915737.41352230842665
34534534.030031809365-0.0300318093646865
35522529.406827228539-7.40682722853853
36527528.875073865166-1.8750738651662
37539538.0462919207910.953708079208923
38533533.498495593853-0.498495593853477
39532534.147861080172-2.14786108017165
40519530.908081947645-11.9080819476446
41538537.8418594691210.158140530879336
42542539.6953019868832.30469801311654
43503508.300312948131-5.30031294813119
44502506.150266856587-4.15026685658694
45522514.841278725737.15872127427008
46511509.2291197227821.77088027721834
47492499.454474038915-7.45447403891473
48500501.341235430462-1.34123543046184
49509510.936593504713-1.9365935047133
50511503.5556926946397.44430730536124
51505504.7224571997590.277542800241463
52493495.230095705913-2.23009570591250
53518512.3055081893275.69449181067341
54518517.2135282654740.786471734525662
55474481.725044855555-7.72504485555481
56471479.376457379085-8.37645737908474
57483493.150719195499-10.150719195499
58461477.835306271544-16.8353062715440
59439454.961426817135-15.9614268171351
60446454.966726298304-8.96672629830408
61461457.6649757482683.33502425173208
62449455.441866952845-6.44186695284452
63441444.42357547944-3.42357547943982
64424429.696144133860-5.69614413385955
65447443.9173378412723.0826621587276
66448440.4064896894457.59351031055508
67404403.5881191122040.411880887795746
68403400.0688968229092.93110317709068
69411410.9580522531380.0419477468622063
70386394.209388588765-8.20938858876485
71359375.020366169593-16.0203661695929
72370375.931993643433-5.9319936434328
73385383.4274373523471.57256264765283
74369373.170317363881-4.17031736388128
75368363.7471369205174.25286307948346
76352350.4840115893171.51598841068346
77378367.40936974605910.5906302539411
78383368.16718842694914.8328115730511
79334335.211434335291-1.21143433529147
80323332.266994940328-9.26699494032829
81330334.233761039406-4.23376103940632
82303312.879468444041-9.8794684440407
83275290.078212639396-15.0782126393955
84284293.336608130076-9.33660813007629
85301299.1344144423371.86558555766288
86281285.645488727339-4.64548872733883
87284279.6064483397054.39355166029526
88272265.7032184676726.29678153232817
89297282.20714796361814.7928520363818
90300284.61749479234315.3825052076568
91240250.715508243537-10.7155082435370
92236238.813634439328-2.8136344393275
93247241.6188081902925.38119180970767
94218223.721025347817-5.72102534781669
95192202.998784772789-10.9987847727890
96201206.224278226510-5.22427822650963
97223214.3439276798298.65607232017149
98197202.023412541456-5.02341254145594
99195199.753519815179-4.75351981517909
100175186.183052349202-11.1830523492015
101197192.6350373004884.36496269951246
102204188.57102078939415.4289792106057
103142153.54340498258-11.5434049825801
104142144.302293384114-2.30229338411394
105151145.1671743418045.83282565819627
106127127.096249244895-0.0962492448953896
107100110.522622860485-10.5226228604849
108114109.1339995098184.86600049018212
109139117.04143035189221.9585696481085
110112107.2690458402144.73095415978588
111123105.65045340779817.3495465922016
11210899.9957570241458.00424297585492
113132113.6730714323418.3269285676599
114140120.63178663317619.3682133668242
1157690.8084882621225-14.8084882621225
1167186.612711818103-15.6127118181030
1178184.7881741440841-3.78817414408411
1185768.9266573814876-11.9266573814876
1193850.9446403466986-12.9446403466986
1204649.5255992797587-3.52559927975867


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
12151.022895070708634.930370425497867.1154197159195
12235.683022751407818.850967149054752.5150783537609
12331.380663373630712.938361366365349.8229653808961
12420.87741969725471.5937341839944640.161105210515
12516.3552363945071-6.5336914196687539.244164208683
1267.4764579354967-18.297697872339133.2506137433325
127-1.39573003629813-22.510980873607419.7195208010111
128-7.5376989378109-30.784025796257715.7086279206359
129-16.8539920080947-45.634801567326811.9268175511373
130-19.5597083221400-46.35584889814347.23643225386336
131-20.4147863808383-45.43067590178074.60110314010403
132-35.6448778254921-66.1727852380443-5.11697041293999
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Aug/16/t12819611841umy7yd4oejoziy/1vtv61281961209.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Aug/16/t12819611841umy7yd4oejoziy/1vtv61281961209.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Aug/16/t12819611841umy7yd4oejoziy/2o2c91281961209.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Aug/16/t12819611841umy7yd4oejoziy/2o2c91281961209.ps (open in new window)


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