Home » date » 2008 » May » 05 » attachments

double

R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Mon, 05 May 2008 07:14:57 -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/05/t1209993322z8wmwn1k9cixfth.htm/, Retrieved Mon, 05 May 2008 15:15:22 +0200
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
56421 53152 53536 52408 41454 38271 35306 26414 31917 38030 27534 18387 50556 43901 48572 43899 37532 40357 35489 29027 34485 42598 30306 26451 47460 50104 61465 53726 39477 43895 31481 29896 33842 39120 33702 25094 51442 45594 52518 48564 41745 49585 32747 33379 35645 37034 35681 20972 58552 54955 65540 51570 51145 46641 35704 33253 35193 41668 34865 21210 56126 49231 59723 48103 47472 50497 40059 34149 36860 46356 36577
 
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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.580952155663673
beta0
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
35353653152384
45240853375.0856277749-967.085627774853
54145452813.2551476077-11359.2551476077
63827146214.0713828713-7943.07138287133
73530641599.5269404018-6293.5269404018
82641437943.288897648-11529.2888976480
93191731245.3236592901671.676340709866
103803031635.53547733386394.46452266618
112753435350.4134260916-7816.41342609162
121838730809.4511966452-12422.4511966452
135055623592.601395327426963.3986046726
144390139257.04593873084643.95406126919
154857241954.96106142826617.03893857179
164389945799.1440969020-1900.14409690195
173753244695.2512877352-7163.25128773516
184035740533.7450105648-176.745010564839
193548940431.0646156744-4942.0646156744
202902737559.9615237692-8532.9615237692
213448532602.71913234031882.28086765969
224259833696.23425997178901.7657400283
233030638867.7342558542-8561.73425585416
242645133893.7762836962-7442.77628369618
254746029569.879357560417890.1206424396
265010439963.183509868910140.8164901311
276146545854.512710000315610.4872899997
285372654923.458952086-1197.45895208598
293947754227.7925925529-14750.7925925529
304389545658.2878381615-1763.28783816154
313148144633.9019675261-13152.9019675261
322989636992.6952162588-7096.69521625883
333384232869.8548322852972.145167714807
343912033434.62466308715685.37533691287
353370236737.5557208237-3035.55572082374
362509434974.043081174-9880.043081174
375144229234.21075511622207.789244884
384559442135.87378945593458.12621054411
395251844144.87966602858373.12033397146
404856449009.2619736806-445.261973680586
414174548750.5860702358-7005.58607023579
424958544680.67574104494904.32425895509
433274747529.8534913585-14782.8534913585
443337938941.7228886935-5562.72288869353
453564535710.0470351474-65.0470351473705
463703435672.2578198591361.74218014102
473568136463.3648748701-782.36487487005
482097236008.8483142988-15036.8483142988
495855227273.158871719231278.8411282808
505495545444.66905185559510.33094814451
516554050969.71631725514570.2836827450
525157059434.354031377-7864.35403137694
535114554865.5406039462-3720.54060394621
544664152704.0845198494-6063.08451984944
553570449181.7224980719-13477.7224980719
563325341351.8105593802-8098.81055938023
573519336646.7891065966-1453.78910659657
584166835802.20719123895865.79280876107
593486539209.9521681651-4344.95216816515
602121036685.7428398141-15475.7428398141
615612627695.076676527428430.9233234726
624923144212.08286880745018.91713119259
635972347127.833595271112595.1664047289
644810354445.022669041-6342.02266904102
654747250760.6109281938-3288.61092819376
665049748850.08532032051646.91467967952
674005949806.8639536744-9747.86395367444
683414944143.8213766711-9994.82137667106
693686038337.3083524207-1477.30835242065
704635637479.06288050198876.93711949808
713657742636.1386357652-6059.1386357652


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
7239116.068983852418610.572225801459621.5657419033
7339116.068983852415401.345340862662830.7926268422
7439116.068983852412577.402144082265654.7358236226
7539116.068983852410026.318936778168205.8190309267
7639116.06898385247681.5930383208670550.5449293839
7739116.06898385245500.0162985780472732.1216691267
7839116.06898385243451.6368657311274780.5011019736
7939116.06898385241514.6801999654276717.4577677393
8039116.0689838524-327.27213535251978559.4101030572
8139116.0689838524-2086.9634144420480319.1013821468
8239116.0689838524-3774.5197155795882006.6576832843
8339116.0689838524-5398.1458149372383630.283782642
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/1leoi1209993291.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/1leoi1209993291.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/26xro1209993291.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/26xro1209993291.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/3y6cy1209993291.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/05/t1209993322z8wmwn1k9cixfth/3y6cy1209993291.ps (open in new window)


 
Parameters (Session):
par1 = 12 ; par2 = Double ; par3 = additive ;
 
Parameters (R input):
par1 = 12 ; par2 = Double ; 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