Home » date » 2010 » Jun » 02 »

Triple exponential smoothing - USA huisprijzen - Yannick Geerts

*Unverified author*
R Software Module: /rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Wed, 02 Jun 2010 09:07:40 +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/02/t1275469744h7yb2d6cxt1s6se.htm/, Retrieved Wed, 02 Jun 2010 11:09:08 +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/02/t1275469744h7yb2d6cxt1s6se.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 «
178600 181600 178500 176700 183500 175900 179800 186500 182700 182800 178600 183300 182900 191400 189300 192200 187900 193900 189100 193100 194800 200200 211500 202100 200300 199200 204900 207300 200000 197700 202200 200200 208300 215100 210700 208100 209000 211000 210200 205500 211400 211700 209300 207500 203300 207100 206900 228700 226900 265000 227100 228100 226500 225200 217800 221300 215300 231300 227100 237800 230200 233400 231100 237200 243700 239700 248400 241000 254500 242800 268300 253900 262100 264100 261000 269300 260400 263200 279200 272200 269200 289600 283200 284300 283000 289100 289600 289100 287400 279600 289300 295000 299600 293600 294400 290200 301000 307900 298800 310300 293900 305000 311300 317300 296200 306800 291800 301900 314600 321500 329400 311700 309700 306500 307100 301300 292200 310100 316800 284400 284600 301200 287600 314300 298200 299400 301900 265500 etc...
 
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'George Udny Yule' @ 72.249.76.132


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.458431564192018
beta0
gamma0.181610953212879


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
13182900176965.1555798805934.84442012027
14191400188445.3497983122954.65020168756
15189300187997.3065106011302.69348939907
16192200191306.444594224893.555405775842
17187900186404.0328902781495.96710972194
18193900191965.688018621934.31198137999
19189100192094.375243636-2994.37524363567
20193100198267.566327154-5167.56632715376
21194800192056.0887287082743.91127129228
22200200193310.9594464646889.04055353641
23211500192109.28448992919390.7155100714
24202100206287.940775720-4187.9407757204
25200300204369.896168068-4069.89616806776
26199200211879.599998405-12679.5999984053
27204900203847.1998726691052.80012733137
28207300207133.913531152166.086468847760
292e+05201451.743538031-1451.74353803066
30197700205983.873525062-8283.87352506205
31202200200833.0027570671366.99724293308
32200200209147.760854121-8947.76085412118
33208300201776.1027962206523.89720377969
34215100205122.7140658939977.28593410732
35210700206211.9877783464488.01222165412
36208100211286.821834192-3186.82183419238
37209000209806.216977001-806.216977001139
38211000218199.445553946-7199.44555394645
39210200213934.779083474-3734.77908347431
40205500215002.296545961-9502.29654596126
41211400204613.5359893146786.46401068577
42211700212357.148924390-657.148924390378
43209300211574.596376989-2274.59637698889
44207500217435.390967366-9935.39096736582
45203300211003.170770206-7703.17077020579
46207100208142.069715966-1042.06971596589
47206900203707.8275827743192.1724172261
48228700207410.80314461821289.1968553815
49226900217354.8996349159545.10036508532
50265000230270.86912400734729.1308759932
51227100245313.692234208-18213.6922342077
52228100239290.855097122-11190.8550971222
53226500229081.595426891-2581.59542689088
54225200232065.152869932-6865.152869932
55217800228142.047851911-10342.0478519109
56221300229845.229287021-8545.22928702127
57215300224076.663401462-8776.6634014617
58231300221395.8098438619904.19015613873
59227100221989.0072319995110.99276800075
60237800228521.5648345219278.43516547879
61230200231661.132096924-1461.13209692389
62233400242150.644114789-8750.64411478877
63231100232321.415878742-1221.41587874194
64237200234726.9286337172473.07136628326
65243700231563.77426680412136.2257331958
66239700240975.883431095-1275.88343109511
67248400239135.6663569389264.33364306192
68241000250512.630073953-9512.63007395333
69254500243987.00081699510512.9991830052
70242800252250.820802007-9450.82080200658
71268300242989.10389933925310.8961006614
72253900259642.715979594-5742.71597959386
73262100254519.4192343197580.5807656812
74264100269522.054232833-5422.05423283263
75261000261190.57480932-190.574809320009
76269300264719.2761642654580.72383573541
77260400262870.43460679-2470.43460678984
78263200264412.595223461-1212.59522346064
79279200263496.92099156115703.0790084389
80272200276416.950694722-4216.95069472247
81269200274117.190944808-4917.19094480819
82289600273363.62058062016236.3794193803
83283200278818.0080156294381.99198437115
84284300282904.1506268151395.84937318455
85283000282118.236810350881.763189650315
86289100293608.791797489-4508.79179748870
87289600285611.0082003843988.99179961625
88289100291826.345180788-2726.34518078784
89287400285437.6091760441962.39082395629
90279600289301.819961965-9701.81996196514
91289300286134.3639940553165.63600594504
92295000291496.968991983503.03100801975
93299600292559.8532775547040.14672244573
94293600299548.963460424-5948.96346042375
95294400293449.270544507950.729455492808
96290200295701.405161618-5501.40516161849
97301000291618.6658334939381.33416650724
98307900306918.809861528981.190138471895
99298800301927.448016464-3127.4480164643
100310300304335.9441961785964.055803822
101293900302061.557556553-8161.5575565531
102305000300153.7643761674846.23562383285
103311300305016.9568939826283.04310601793
104317300312040.3303814875259.66961851279
105296200314158.891597974-17958.8915979736
106306800308449.967021670-1649.96702166955
107291800304854.911626939-13054.9116269391
108301900300049.3414128191850.65858718136
109314600300754.87055926713845.1294407328
110321500317589.5162172383910.48378276231
111329400313268.15159515516131.8484048450
112311700325644.926006498-13944.9260064982
113309700312563.644983304-2863.64498330437
114306500314455.793152962-7955.79315296229
115307100313631.033810826-6531.03381082579
116301300314695.890035702-13395.8900357022
117292200306005.737167801-13805.7371678005
118310100303759.1777990106340.82220099046
119316800302694.05509109814105.9449089023
120284400311867.546794029-27467.5467940287
121284600300284.482903393-15684.4829033932
122301200302182.366317398-982.366317397682
123287600297122.015972044-9522.01597204432
124314300294666.78774632519633.2122536746
125298200298364.205412652-164.205412652169
126299400300916.429426254-1516.42942625372
127301900303117.26724659-1217.26724659029
128265500305861.502323839-40361.5023238392
129287100284977.3679734022122.63202659832
130274000291801.285710803-17801.2857108031
131290100280733.0131062389366.98689376237
132263100283739.576001709-20639.5760017087
133245200276325.894226481-31125.8942264812
134258600271633.909418633-13033.9094186330
135259800260972.012121791-1172.0121217907
136269800264761.8802517375038.11974826292
137274600260865.93351603513734.0664839647
138274800269504.6694433125295.33055668755
139271100274681.365786238-3581.36578623787
140257800272415.507552164-14615.5075521643
141290300267465.02016579722834.9798342034
142262200281694.041620673-19494.0416206730
143270000272550.080019718-2550.08001971780
144290600267344.13809180223255.861908198


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
145278851.652072713266969.162478297290734.141667129
146290965.791206565275324.750422770306606.83199036
147286985.703436384268722.894252073305248.512620694
148292315.461143935271390.176525797313240.745762074
149286339.199557518263646.250254348309032.148860689
150287878.024990257263133.060148915312622.989831598
151289798.848287068263125.288177991316472.408396145
152287878.767987302259706.456534833316051.07943977
153293610.987214270263379.898151688323842.076276851
154293073.160486438261445.854646713324700.466326164
155294564.959400048261415.384454246327714.534345851
156292751.610929177260609.419203813324893.80265454
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/1bmb51275469657.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/1bmb51275469657.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/24vaq1275469657.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/24vaq1275469657.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/34vaq1275469657.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Jun/02/t1275469744h7yb2d6cxt1s6se/34vaq1275469657.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