Home » date » 2010 » Nov » 28 »

W8 Smoothing Model

*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: Sun, 28 Nov 2010 14:41:54 +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/Nov/28/t12909558880s4jxk0srnzpif0.htm/, Retrieved Sun, 28 Nov 2010 15:51:28 +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/Nov/28/t12909558880s4jxk0srnzpif0.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 «
10057 10900 11771 11992 11993 14504 11727 11477 13578 11555 11846 11397 10066 10269 14279 13870 13695 14420 11424 9704 12464 14301 13464 9893 11572 12380 16692 16052 16459 14761 13654 13480 18068 16560 14530 10650 11651 13735 13360 17818 20613 16231 13862 12004 17734 15034 12609 12320 10833 11350 13648 14890 16325 18045 15616 11926 16855 15083 12520 12355
 
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.680282906594522
beta0
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
21090010057843
31177110630.47849025921140.52150974082
41199211406.3557779392585.644222060762
51199311804.7595315530188.240468446978
61450411932.81630456682571.18369543315
71172713681.9486222846-1954.94862228456
81147712352.0304912739-875.030491273861
91357811756.76220531121821.23779468875
101155512995.7191458819-1440.71914588191
111184612015.6225377350-169.622537734987
121139711900.2312247407-503.231224740692
131006611557.8916244850-1491.89162448497
141026910542.9832538563-273.983253856313
151427910356.59712956473922.40287043528
161387013024.9407550991845.059244900873
171369513599.820114464995.1798855351353
181442013664.5693636460755.43063635396
191142414178.4759126755-2754.47591267546
20970412304.653032656-2600.653032656
211246410535.47322855691928.52677144308
221430111847.41702607962453.58297392043
231346413516.547583149-52.5475831489912
24989313480.8003605499-3587.80036054988
251157211040.0811029941531.918897005868
261238011401.9364363218978.063563678164
271669212067.29636025504624.70363974499
281605215213.403194439838.596805561003
291645915783.8862667869675.113733213082
301476116243.154599499-1482.15459949899
311365415234.8701605294-1580.87016052938
321348014159.4312127759-679.431212775904
331806813697.22577251774370.77422748233
341656016670.5887680578-110.588768057776
351453016595.3571194867-2065.35711948673
361065015190.3299750866-4540.32997508661
371165112101.6211027365-450.621102736457
381373511795.07126919411939.92873080593
391336013114.7716247730245.228375227049
401781813281.59629665194536.40370334814
412061316367.63419345174245.36580654831
421623119255.6839838874-3024.68398388737
431386217198.0431717986-3336.04317179857
441200414928.5900263626-2924.59002636263
451773412939.04142263134794.95857736869
461503416200.9697806440-1166.96978064402
471260915407.1001863595-2798.10018635953
481232013503.6004586402-1183.60045864020
491083312698.4172983898-1865.41729838983
501135011429.4057966295-79.4057966294968
511364811375.38739049792272.61260950207
521489012921.40690205341968.59309794664
531632514260.60713662642064.39286337358
541804515664.97831407522380.02168592481
551561617284.0663843341-1668.06638433411
561192616149.3093360067-4223.30933600669
571685513276.26418546033578.73581453972
581508315710.8169873093-627.816987309274
591252015283.7238223731-2763.72382237311
601235513403.6097474646-1048.60974746461


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
6112690.25846057608101.3215432534717279.1953778986
6212690.25846057607140.139242810118240.3776783420
6312690.25846057606322.424870440119058.0920507120
6412690.25846057605598.3767686756219782.1401524765
6512690.25846057604941.6929654562420438.8239556958
6612690.25846057604336.471981834221044.0449393179
6712690.25846057603772.2301004137621608.2868207383
6812690.25846057603241.6230494223522138.8938717297
6912690.25846057602739.2689402742322641.2479808778
7012690.25846057602261.0842970382023119.4326241139
7112690.25846057601803.8837180984523576.6332030536
7212690.25846057601365.1255106322424015.3914105198
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/1dse31290955310.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/1dse31290955310.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/26jdo1290955310.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/26jdo1290955310.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/36jdo1290955310.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Nov/28/t12909558880s4jxk0srnzpif0/36jdo1290955310.ps (open in new window)


 
Parameters (Session):
par1 = 12 ; par2 = Single ; par3 = additive ;
 
Parameters (R input):
par1 = 12 ; par2 = Single ; 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=0, beta=0)
if (par2 == 'Double') fit <- HoltWinters(x, gamma=0)
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