Home » date » 2010 » Jun » 05 »

*Unverified author*
R Software Module: /rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Sat, 05 Jun 2010 10:28:11 +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/Jun/05/t1275733736vfusfwwr7aiyabm.htm/, Retrieved Sat, 05 Jun 2010 12:28:56 +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/Jun/05/t1275733736vfusfwwr7aiyabm.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:
KDGP2W62
 
Dataseries X:
» Textbox « » Textfile « » CSV «
7089 77102 35123 23109 11115 48105 7896 79102 2494 76101 9484 8183 9092 49109 27127 31116 77120 80115 67117 18110 3099 3105 5786 7387 5089 288 65108 44100 68102 8798 8183 288 3680 9374 2271 4671 3466 8068 2780 9883 6479 3373 6468 460 5553 6157 2345 3744 6144 9948 2049 9549 6049 745 9643 9938 5535 6932 1827 9726 7624 9529 7631 6132 7534 2137 8738 8836 3230 633 9126 123
 
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
beta0
gamma0.769366005522916


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
1390926099.769578111472992.23042188853
144910933007.012029395516101.9879706045
152712719441.17666573927685.82333426077
163111624225.11393938896890.8860606111
177712066085.344384459411034.6556155406
188011569059.170076988411055.8299230116
19671176678.9521176144360438.0478823856
201811069241.2740108332-51131.2740108332
2130992293.78721472168805.212785278321
22310569990.6959818707-66885.6959818707
2357867898.92425550947-2112.92425550947
2473876057.870654458711329.12934554129
2550898401.889945404-3312.88994540400
2628845395.3341953175-45107.3341953175
276510825354.387863574339753.6121364257
284410029526.727422354814573.2725776452
296810274575.0332977089-6473.03329770888
30879877565.1497825966-68767.1497825966
31818353177.9315984881-44994.9315984881
3228829902.6099678208-29614.6099678208
3336802913.29055892724766.709441072758
34937418531.1152376787-9157.11523767867
3522716273.31216107567-4002.31216107567
3646717080.4575898611-2409.4575898611
3734665853.06504137149-2387.06504137149
38806810691.2846656788-2623.28466567885
39278055939.4656380834-53159.4656380834
40988340738.9079328143-30855.9079328143
41647969594.9015258338-63115.9015258338
42337324658.0424431642-21285.0424431642
43646818560.3608057825-12092.3608057825
444607118.13579175937-6658.13579175937
4555533503.17073900212049.82926099790
46615711485.9420651528-5328.9420651528
4723453194.06924085309-849.069240853092
4837445226.70282847279-1482.70282847279
4961444016.538345568112127.46165443189
5099488673.0186210961274.98137890401
51204915040.3799043785-12991.3799043785
52954916999.4212997621-7450.42129976212
53604921035.6724839253-14986.6724839253
547458282.05436128123-7537.05436128123
5596439256.90947529575386.090524704254
5699381995.592453424317942.4075465757
5755355080.23968954004454.760310459956
5869327386.03519482315-454.035194823153
5918272540.82423060557-713.824230605574
6097264085.961675953155640.03832404685
6176245653.335020541551970.66497945845
6295299653.94595169947-124.945951699468
6376315045.253841116122585.74615888388
64613211267.3204249013-5135.32042490129
6575349505.4361388875-1971.4361388875
6621372483.30095393322-346.300953933217
6787389553.9544000577-815.954400057704
6888368106.21082176831729.789178231689
6932305430.11681306898-2200.11681306898
706337036.71595061524-6403.71595061524
7191261991.632133659097134.3678663409
721238425.21543232124-8302.21543232124


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
737169.4976640114-26616.704934043240955.700262066
749557.81678393419-24228.385814120443344.0193819888
757034.63903467283-26751.563563381840820.8416327274
767316.37946251474-26469.823135539941102.5820605693
777988.6801915681-25797.522406486541774.8827896227
782216.86877229684-31569.333825757836003.0713703514
798926.18682259646-24860.015775458142712.3894206511
808667.68580669828-25118.516791356342453.8884047529
813737.42172891429-30048.780869140337523.6243269689
822109.91458918701-31676.288008867635896.1171872416
837480.57224091685-26305.630357137841266.7748389714
842037.77310816554NANA
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/1icwt1275733687.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/1icwt1275733687.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/2blve1275733687.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/2blve1275733687.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/3blve1275733687.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/05/t1275733736vfusfwwr7aiyabm/3blve1275733687.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