Home » date » 2009 » Jan » 25 »

Dennis Collin oef 10 deel 2

*Unverified author*
R Software Module: rwasp_exponentialsmoothing.wasp (opens new window with default values)
Title produced by software: Exponential Smoothing
Date of computation: Sun, 25 Jan 2009 05:50:17 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd.htm/, Retrieved Sun, 25 Jan 2009 13:50:51 +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/2009/Jan/25/t1232887851sgafbbxvzo8n9gd.htm/},
    year = {2009},
}
@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 = {2009},
    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 «
1,29 1,29 1,3 1,3 1,3 1,3 1,31 1,31 1,31 1,31 1,31 1,32 1,32 1,32 1,32 1,33 1,33 1,33 1,34 1,34 1,34 1,34 1,34 1,34 1,34 1,35 1,36 1,36 1,36 1,37 1,37 1,37 1,37 1,37 1,37 1,37 1,38 1,38 1,38 1,39 1,4 1,4 1,4 1,4 1,41 1,42 1,43 1,43 1,43 1,44 1,45 1,45 1,46 1,46 1,47 1,47 1,47 1,48 1,49 1,49 1,49 1,5 1,51 1,51 1,51 1,52 1,52 1,52 1,52 1,53 1,53 1,53 1,53 1,54 1,54 1,55 1,55 1,55 1,56 1,56 1,58 1,58 1,58 1,58 1,58 1,58 1,59 1,59 1,6 1,6 1,6 1,61 1,62 1,62 1,63 1,63
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001


Estimated Parameters of Exponential Smoothing
ParameterValue
alpha0.693880469628355
beta0.087046697378327
gamma0


Interpolation Forecasts of Exponential Smoothing
tObservedFittedResiduals
31.31.290.01
41.31.297542804728850.00245719527115185
51.31.30000021917251-2.19172507698318e-07
61.31.30075246848991-0.000752468489908953
71.311.300937297598430.0090627024015666
81.311.308480069349180.00151993065081557
91.311.31088086291148-0.000880862911483815
101.311.31156258858657-0.00156258858657199
111.311.31167689777435-0.00167689777435487
121.321.311610605418200.00838939458179566
131.321.319035836189620.000964163810384111
141.321.32136707984355-0.00136707984354878
151.321.32199814742907-0.00199814742907023
161.331.322070641430760.00792935856923993
171.331.32951057124020.000489428759800514
181.331.33181764055815-0.00181764055814893
191.341.332414094038660.00758590596133701
201.341.339993673539066.32646094178568e-06
211.341.34231431297440-0.00231431297439522
221.341.34288492151731-0.00288492151731257
231.341.34288534666746-0.0028853466674561
241.341.34271120186611-0.00271120186610796
251.341.34249413613976-0.00249413613976057
261.351.342277042250590.007722957749412
271.361.349615856940690.0103841430593148
281.361.359428418417980.000571581582015757
291.361.36246675865857-0.00246675865857293
301.371.363247861715120.00675213828487942
311.371.37083360648664-0.000833606486641436
321.371.37310540127911-0.00310540127911074
331.371.37361327578671-0.00361327578671133
341.371.37355050422318-0.00355050422318159
351.371.37331683815499-0.00331683815499084
361.371.37304497137257-0.00304497137257020
371.381.372777831359990.00722216864000713
381.381.38007107829046-7.10782904578267e-05
391.381.38229939048672-0.00229939048672101
401.391.382842637176880.0071573628231234
411.41.390380045032650.00961995496734902
421.41.40020724279332-0.000207242793318629
431.41.40320302249129-0.00320302249129223
441.41.40392662659663-0.00392662659662824
451.411.403910967686310.00608903231369462
461.421.411212756455890.0087872435441092
471.431.420917530839870.0090824691601279
481.431.43137573767979-0.0013757376797916
491.431.43449410448582-0.00449410448581755
501.441.435177253542430.00482274645757341
511.451.442616477408840.00738352259116426
521.451.45227853860998-0.00227853860998461
531.461.455098660507840.00490133949216065
541.461.46319680051637-0.0031968005163725
551.471.465482712571070.00451728742893054
561.471.47339412376743-0.00339412376742532
571.471.47561095616083-0.00561095616082552
581.481.475950670081620.00404932991838014
591.491.483238047384050.0067619526159457
601.491.49281608255782-0.00281608255782362
611.491.49557801479144-0.00557801479143771
621.51.496086584078140.00391341592185745
631.511.503417442100060.00658255789994211
641.511.51299795212982-0.00299795212982357
651.511.51574965704245-0.00574965704245223
661.521.516244728353940.00375527164606027
671.521.52356190246725-0.0035619024672473
681.521.52558669344984-0.00558669344983964
691.521.52586908521197-0.00586908521196827
701.531.525601038079660.0043989619203435
711.531.53272348562799-0.00272348562799318
721.531.53473930738579-0.00473930738578709
731.531.53507013561406-0.005070135614061
741.541.534865172388330.00513482761166628
751.541.54205137744268-0.00205137744268025
761.551.544127311954850.00587268804514696
771.551.55205631112614-0.00205631112613647
781.551.55435933143005-0.00435933143005185
791.561.554801027291350.00519897270864855
801.561.56218906168580-0.00218906168579602
811.581.564318463972040.0156815360279579
821.581.579795089819600.000204910180397588
831.581.58454516383136-0.0045451638313585
841.581.58572472634685-0.0057247263468525
851.581.58574003998015-0.00574003998014727
861.581.58539802934889-0.00539802934888645
871.591.584967292224810.00503270777518838
881.591.59207821544097-0.00207821544096998
891.61.594129483697010.00587051630299285
901.61.60205080087231-0.00205080087231280
911.61.60435180238658-0.00435180238657606
921.611.604793335011190.00520666498881384
931.621.612181784049440.00781821595056109
941.621.62185455756436-0.00185455756436070
951.631.624703567167210.00529643283278647
961.631.63283441390549-0.00283441390549433


Extrapolation Forecasts of Exponential Smoothing
tForecast95% Lower Bound95% Upper Bound
971.635152226281421.625303460207191.64500099235564
981.639436783109211.627100484989461.65177308122896
991.643721339937001.629005281799911.65843739807409
1001.648005896764791.630951857052661.66505993647692
1011.652290453592581.632906808644881.67167409854028
1021.656575010420371.634851144284381.67829887655635
1031.660859567248161.636773263967181.68494587052913
1041.665144124075951.638665748768191.69162249938371
1051.669428680903741.640523713347731.69833364845975
1061.673713237731531.642343887927901.70508258753515
1071.677997794559321.644124073743081.71187151537556
1081.682282351387111.645862803835691.71870189893853
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/1xfa11232887810.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/1xfa11232887810.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/26auz1232887810.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/26auz1232887810.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/3ysyu1232887810.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2009/Jan/25/t1232887851sgafbbxvzo8n9gd/3ysyu1232887810.ps (open in new window)


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