Home » date » 2008 » May » 27 »

smoothing champagne Katrien Broux

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Tue, 27 May 2008 10:27:48 -0600
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/May/27/t1211905700aar27iwki509kau.htm/, Retrieved Tue, 27 May 2008 18:28:24 +0200
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
2815 2672 2755 2721 2946 3036 2282 2212 2922 4301 5764 7312 2541 2475 3031 3266 3776 3230 3028 1759 3595 4474 6838 8357 3113 3006 4047 3523 3937 3986 3260 1573 3528 5211 7614 9254 5375 3088 3718 4514 4520 4539 3663 1643 4739 5428 8314 10651 3633 4292 4154 4121 4647 4753 3965 1723 5048 6922 9858 11331 4016 4276 4968 4677 3523 1821 5222 6872 10803 13916 2639 2899 3370 3740 2927 3986 4217 1738 5221 6424 9842 13076
 
Text written by user:
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'George Udny Yule' @ 72.249.76.132


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0
beta0
gamma0.89588635004732


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
1325412368.49607099695172.503929003050
1424752299.23063554187175.769364458128
1530312808.66370689647222.336293103529
1632662997.43890828896268.561091711039
1737763417.25819820299358.741801797015
1832302855.58006291701374.419937082992
1930282630.60624337706397.393756622937
2017591510.29684667595248.70315332405
2135953036.94605462600558.053945374004
2244743730.0172037031743.982796296901
2368385676.504419523171161.49558047683
2483576872.948968091051484.05103190895
2531132523.03998632031589.960013679685
2630062456.7000099164549.299990083599
2740473007.851757008041039.14824299196
2835233238.03912450669284.960875493313
2939373738.65008162431198.349918375688
3039863191.01777373524794.982226264764
3132602986.62588552958273.374114470421
3215731733.10660695269-160.106606952692
3335283536.89896687662-8.89896687661849
3452114396.54123557553814.458764424472
3576146717.07245571265896.927544287351
3692548202.490030351921051.50996964808
3753753051.577109649672323.42289035033
3830882948.81037311343139.189626886575
3937183938.81048358019-220.810483580192
4045143493.331683158681020.66831684132
4145203916.34906603009603.650933969909
4245393903.23149877607635.768501223933
4336633231.5380231399431.461976860098
4416431589.6692832313853.330716768616
4547393528.926503922331210.07349607767
4654285126.20372529982301.796274700182
4783147520.61759962115793.382400378851
48106519144.523459098311506.47654090169
4936335133.09996250202-1500.09996250202
5042923073.508459909291218.49154009071
5141543740.98938539335413.01061460665
5241214407.73449614259-286.734496142591
5346474457.15169796705189.848302032951
5447534472.80782081263280.192179187369
5539653618.07891877330346.921081226703
5617231637.4475444226385.5524555773732
5750484613.01483161235434.985168387647
5869225396.578888298841525.42111170116
5998588231.398062488341626.60193751166
601133110494.1552287586836.844771241364
6140163789.18088238996226.819117610035
6242764165.13839832469110.861601675308
6349684110.9999574441857.0000425559
6446774150.85297496075526.147025039253
6535234627.23420033803-1104.23420033803
6618214723.82816953661-2902.82816953661
6752223928.880779987961293.11922001204
6868721714.092821587425157.90717841258
69108035002.712106443885800.28789355612
70139166763.182840345927152.81715965408
7126399688.64853526556-7049.64853526556
72289911243.8730364222-8344.87303642225
7333703992.38503378657-622.385033786572
7437404264.45779400998-524.457794009983
7529274878.7745975599-1951.77459755990
7639864622.22091281142-636.220912811418
7742173637.96585299977579.034147000228
7817382123.22403591592-385.224035915915
7952215087.36863818058133.631361819417
8064246334.9914575383389.0085424616655
81984210199.1108566255-357.110856625526
821307613171.2940980642-95.294098064247


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
833372.96463989007-115.8826753137256861.81195509386
843767.81519021363278.9678750098387256.66250541742
853434.79877754344-54.04853766034926923.64609274723
863794.60321518051305.7558999767187283.4505303843
873130.20637723689-358.6409379669066619.05369244068
884052.23928140902563.391966205237541.08659661281
894156.71464150857667.8673263047777645.56195671236
901778.10708042871-1710.740234775085266.9543956325
915207.087151172831718.239835969048695.93446637663
926414.732995767352925.885680563569903.58031097114
939879.1801147216390.3327995172213368.0274299248
9413085.9214163684NANA
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/1gz4q1211905662.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/1gz4q1211905662.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/2ub901211905662.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/2ub901211905662.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/3c18d1211905662.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/27/t1211905700aar27iwki509kau/3c18d1211905662.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=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