Home » date » 2008 » May » 19 » attachments

Aardappelprijs - Time Series Analysis - Kristof Vermeulen

R Software Module: rwasp_decompose.wasp (opens new window with default values)
Title produced by software: Classical Decomposition
Date of computation: Mon, 19 May 2008 14:16:09 -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/19/t121122825003bietcrfm6w6q6.htm/, Retrieved Mon, 19 May 2008 22:17:35 +0200
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
0,36 0,35 0,35 0,35 0,33 0,78 0,71 0,62 0,52 0,46 0,43 0,43 0,42 0,42 0,42 0,42 0,43 0,99 1,03 0,83 0,64 0,6 0,58 0,58 0,58 0,57 0,57 0,56 0,56 0,88 0,84 0,69 0,59 0,54 0,52 0,52 0,51 0,52 0,51 0,51 0,53 0,95 0,98 0,88 0,81 0,77 0,76 0,75 0,73 0,74 0,73 0,75 0,77 1,09 1,03 0,9 0,76 0,66 0,63 0,61 0,61 0,61 0,61 0,61 0,62 0,76 0,83 0,81 0,77 0,75 0,76 0,76
 
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 time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Classical Decomposition by Moving Averages
tObservationsFitTrendSeasonalRandom
10.36NANA-0.095138888888889NA
20.35NANA-0.098888888888889NA
30.35NANA-0.109305555555556NA
40.35NANA-0.111388888888889NA
50.33NANA-0.103055555555556NA
60.78NANA0.297986111111111NA
70.710.7632986111111110.4766666666666670.286631944444445-0.0532986111111113
80.620.6197569444444440.4820833333333330.1376736111111110.000243055555555649
90.520.4966319444444440.4879166666666670.008715277777777810.0233680555555555
100.460.4410069444444440.49375-0.05274305555555560.0189930555555555
110.430.4241319444444440.500833333333333-0.07670138888888890.00586805555555558
120.430.4299652777777780.51375-0.08378472222222223.47222222221166e-05
130.420.4406944444444440.535833333333333-0.095138888888889-0.0206944444444443
140.420.4590277777777780.557916666666667-0.098888888888889-0.0390277777777777
150.420.4623611111111110.571666666666667-0.109305555555556-0.042361111111111
160.420.4711111111111110.5825-0.111388888888889-0.0511111111111111
170.430.4915277777777780.594583333333333-0.103055555555556-0.0615277777777778
180.990.9050694444444440.6070833333333330.2979861111111110.0849305555555556
191.030.9066319444444440.620.2866319444444450.123368055555556
200.830.7705902777777780.6329166666666670.1376736111111110.0594097222222222
210.640.6541319444444440.6454166666666670.00871527777777781-0.0141319444444443
220.60.6047569444444440.6575-0.0527430555555556-0.0047569444444443
230.580.5920486111111110.66875-0.0767013888888889-0.0120486111111112
240.580.5857986111111110.669583333333333-0.0837847222222222-0.00579861111111113
250.580.5619444444444440.657083333333333-0.0951388888888890.0180555555555557
260.570.5444444444444440.643333333333333-0.0988888888888890.0255555555555556
270.570.5261111111111110.635416666666667-0.1093055555555560.0438888888888889
280.560.5194444444444440.630833333333333-0.1113888888888890.0405555555555556
290.560.5227777777777780.625833333333333-0.1030555555555560.0372222222222223
300.880.9188194444444440.6208333333333330.297986111111111-0.0388194444444445
310.840.9020486111111110.6154166666666670.286631944444445-0.0620486111111112
320.690.7480902777777780.6104166666666670.137673611111111-0.0580902777777778
330.590.6145486111111110.6058333333333330.00871527777777781-0.0245486111111111
340.540.5485069444444440.60125-0.0527430555555556-0.00850694444444433
350.520.5212152777777780.597916666666667-0.0767013888888889-0.00121527777777763
360.520.5157986111111110.599583333333333-0.08378472222222220.0042013888888891
370.510.5131944444444440.608333333333333-0.095138888888889-0.00319444444444428
380.520.5231944444444440.622083333333333-0.098888888888889-0.00319444444444439
390.510.5298611111111110.639166666666667-0.109305555555556-0.0198611111111111
400.510.5465277777777780.657916666666667-0.111388888888889-0.0365277777777777
410.530.5744444444444440.6775-0.103055555555556-0.0444444444444444
420.950.9950694444444440.6970833333333330.297986111111111-0.0450694444444444
430.981.002465277777780.7158333333333330.286631944444445-0.0224652777777777
440.880.8718402777777780.7341666666666670.1376736111111110.00815972222222228
450.810.7612152777777780.75250.008715277777777810.0487847222222223
460.770.7189236111111110.771666666666667-0.05274305555555560.0510763888888889
470.760.7149652777777780.791666666666667-0.07670138888888890.0450347222222223
480.750.7237152777777780.8075-0.08378472222222220.0262847222222222
490.730.7202777777777780.815416666666667-0.0951388888888890.00972222222222208
500.740.7194444444444440.818333333333333-0.0988888888888890.0205555555555555
510.730.7077777777777780.817083333333333-0.1093055555555560.0222222222222221
520.750.6990277777777780.810416666666667-0.1113888888888890.0509722222222223
530.770.6973611111111110.800416666666667-0.1030555555555560.072638888888889
541.091.087152777777780.7891666666666670.2979861111111110.00284722222222245
551.031.064965277777780.7783333333333330.286631944444445-0.0349652777777777
560.90.9055902777777780.7679166666666670.137673611111111-0.00559027777777776
570.760.7662152777777780.75750.00871527777777781-0.00621527777777786
580.660.6939236111111110.746666666666667-0.0527430555555556-0.0339236111111111
590.630.6578819444444440.734583333333333-0.0767013888888889-0.0278819444444444
600.610.6307986111111110.714583333333333-0.0837847222222222-0.0207986111111110
610.610.5973611111111110.6925-0.0951388888888890.0126388888888890
620.610.5815277777777780.680416666666667-0.0988888888888890.0284722222222221
630.610.5677777777777780.677083333333333-0.1093055555555560.0422222222222223
640.610.5698611111111110.68125-0.1113888888888890.0401388888888889
650.620.5873611111111110.690416666666666-0.1030555555555560.0326388888888891
660.761.000069444444440.7020833333333330.297986111111111-0.240069444444444
670.83NANA0.286631944444445NA
680.81NANA0.137673611111111NA
690.77NANA0.00871527777777781NA
700.75NANA-0.0527430555555556NA
710.76NANA-0.0767013888888889NA
720.76NANA-0.0837847222222222NA
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/1w49c1211228167.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/1w49c1211228167.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/2clml1211228167.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/2clml1211228167.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/35zv01211228167.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/35zv01211228167.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/4tp9f1211228167.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/May/19/t121122825003bietcrfm6w6q6/4tp9f1211228167.ps (open in new window)


 
Parameters (Session):
par1 = additive ; par2 = 12 ;
 
Parameters (R input):
par1 = additive ; par2 = 12 ;
 
R code (references can be found in the software module):
par2 <- as.numeric(par2)
x <- ts(x,freq=par2)
m <- decompose(x,type=par1)
m$figure
bitmap(file='test1.png')
plot(m)
dev.off()
mylagmax <- length(x)/2
bitmap(file='test2.png')
op <- par(mfrow = c(2,2))
acf(as.numeric(x),lag.max = mylagmax,main='Observed')
acf(as.numeric(m$trend),na.action=na.pass,lag.max = mylagmax,main='Trend')
acf(as.numeric(m$seasonal),na.action=na.pass,lag.max = mylagmax,main='Seasonal')
acf(as.numeric(m$random),na.action=na.pass,lag.max = mylagmax,main='Random')
par(op)
dev.off()
bitmap(file='test3.png')
op <- par(mfrow = c(2,2))
spectrum(as.numeric(x),main='Observed')
spectrum(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
spectrum(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
spectrum(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
bitmap(file='test4.png')
op <- par(mfrow = c(2,2))
cpgram(as.numeric(x),main='Observed')
cpgram(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
cpgram(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
cpgram(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Classical Decomposition by Moving Averages',6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t',header=TRUE)
a<-table.element(a,'Observations',header=TRUE)
a<-table.element(a,'Fit',header=TRUE)
a<-table.element(a,'Trend',header=TRUE)
a<-table.element(a,'Seasonal',header=TRUE)
a<-table.element(a,'Random',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(m$trend)) {
a<-table.row.start(a)
a<-table.element(a,i,header=TRUE)
a<-table.element(a,x[i])
if (par1 == 'additive') a<-table.element(a,m$trend[i]+m$seasonal[i]) else a<-table.element(a,m$trend[i]*m$seasonal[i])
a<-table.element(a,m$trend[i])
a<-table.element(a,m$seasonal[i])
a<-table.element(a,m$random[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.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