Home » date » 2007 » Nov » 27 » attachments

Workshop 4 Q2 Autocorrelation function 2

R Software Module: rwasp_autocorrelation.wasp (opens new window with default values)
Title produced by software: (Partial) Autocorrelation Function
Date of computation: Tue, 27 Nov 2007 03:51:59 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2007/Nov/27/t1196160546wmkbvl4rqx064q0.htm/, Retrieved Tue, 27 Nov 2007 11:49:09 +0100
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
15761.3 16943.0 15070.3 13659.6 14768.9 14725.1 15998.1 15370.6 14956.9 15469.7 15101.8 11703.7 16283.6 16726.5 14968.9 14861.0 14583.3 15305.8 17903.9 16379.4 15420.3 17870.5 15912.8 13866.5 17823.2 17872.0 17420.4 16704.4 15991.2 16583.6 19123.5 17838.7 17209.4 18586.5 16258.1 15141.6 19202.1 17746.5 19090.1 18040.3 17515.5 17751.8 21072.4 17170.0 19439.5 19795.4 17574.9 16165.4 19464.6 19932.1 19961.2 17343.4 18924.2 18574.1 21350.6 18840.1 20304.8 21132.4 19753.9 18009.9 20390.4
 
Text written by user:
 
Output produced by software:


Summary of compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Autocorrelation Function
Time lag kACF(k)T-STATP-value
017.81020
10.5598934.37292.4e-05
20.4742453.7040.000229
30.5646884.41042.1e-05
40.5251144.10136.2e-05
50.454513.54980.000375
60.5871494.58581.1e-05
70.3685892.87880.002749
80.4427343.45790.000499
90.352072.74980.003918
100.1912121.49340.070242
110.3065252.3940.009877
120.5100763.98389.2e-05
130.1885661.47280.07298
140.1583981.23710.110391
150.1794981.40190.083003
160.1290791.00810.158684
170.1035560.80880.210886
180.191121.49270.070335
190.0234870.18340.427531
200.0790510.61740.269631
21-0.006155-0.04810.519091
22-0.145772-1.13850.870321
23-0.022859-0.17850.570553
240.0896150.69990.243321
25-0.129389-1.01060.84189
26-0.124106-0.96930.831888
27-0.117826-0.92030.819467
28-0.158617-1.23880.889925
29-0.148552-1.16020.874761
30-0.133441-1.04220.849286
31-0.210425-1.64350.947285
32-0.132163-1.03220.846979
33-0.236687-1.84860.965317
34-0.308532-2.40970.990501
35-0.225493-1.76120.958389
36-0.170675-1.3330.906258
37-0.287837-2.24810.985902
38-0.288443-2.25280.98606
39-0.274535-2.14420.981994
40-0.277406-2.16660.98291
41-0.280015-2.1870.983706
42-0.272736-2.13010.981399
43-0.286772-2.23980.985619
44-0.239756-1.87260.967038
45-0.277657-2.16860.982988
46-0.298-2.32750.988362
47-0.242246-1.8920.96838


Partial Autocorrelation Function
Time lag kPACF(k)T-STATP-value
00.5598934.37292.4e-05
10.2341741.8290.036148
20.3548132.77120.003697
30.1628081.27160.104176
40.0575940.44980.327216
50.3008652.34980.011017
6-0.252115-1.96910.973254
70.2015091.57380.060349
8-0.293922-2.29560.987425
9-0.25271-1.97370.973527
100.1970251.53880.064511
110.3088232.4120.009445
12-0.178959-1.39770.916369
13-0.198128-1.54740.936534
14-0.071463-0.55810.710607
150.0037190.0290.488461
16-0.068302-0.53350.702173
170.095280.74420.229818
18-0.080587-0.62940.734282
19-0.078153-0.61040.728067
200.0328970.25690.399048
21-0.093652-0.73140.766347
220.0444330.3470.36488
23-0.058713-0.45860.67591
240.0392740.30670.380043
25-0.039107-0.30540.619464
26-0.019057-0.14880.558915
270.0855320.6680.253318
28-0.125766-0.98230.835075
29-0.098623-0.77030.777945
300.0024510.01910.492394
310.0397930.31080.37851
32-0.026506-0.2070.581657
33-0.002258-0.01760.507006
34-0.099949-0.78060.780979
35-0.069388-0.54190.705081
360.0200250.15640.438118
37-0.034084-0.26620.604512
380.0231810.1810.428465
39-0.012273-0.09590.538026
40-0.002488-0.01940.507722
410.0029390.0230.49088
42-0.027236-0.21270.583873
43-0.043558-0.34020.632563
440.0308690.24110.405144
450.0410210.32040.374885
460.0064330.05020.480047
47-0.043205-0.33740.631528
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/27/t1196160546wmkbvl4rqx064q0/1dvdk1196160717.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/27/t1196160546wmkbvl4rqx064q0/1dvdk1196160717.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/27/t1196160546wmkbvl4rqx064q0/2owsy1196160717.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/27/t1196160546wmkbvl4rqx064q0/2owsy1196160717.ps (open in new window)


 
Parameters:
par1 = 48 ; par2 = 1.1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ;
 
R code (references can be found in the software module):
if (par1 == 'Default') {
par1 = 10*log10(length(x))
} else {
par1 <- as.numeric(par1)
}
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
par5 <- as.numeric(par5)
if (par2 == 0) {
x <- log(x)
} else {
x <- (x ^ par2 - 1) / par2
}
if (par3 > 0) x <- diff(x,lag=1,difference=par3)
if (par4 > 0) x <- diff(x,lag=par5,difference=par4)
bitmap(file='pic1.png')
racf <- acf(x,par1,main='Autocorrelation',xlab='lags',ylab='ACF')
dev.off()
bitmap(file='pic2.png')
rpacf <- pacf(x,par1,main='Partial Autocorrelation',xlab='lags',ylab='PACF')
dev.off()
(myacf <- c(racf$acf))
(mypacf <- c(rpacf$acf))
lengthx <- length(x)
sqrtn <- sqrt(lengthx)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Autocorrelation Function',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Time lag k',header=TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/basics.htm','ACF(k)','click here for more information about the Autocorrelation Function'),header=TRUE)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,'P-value',header=TRUE)
a<-table.row.end(a)
for (i in 1:par1) {
a<-table.row.start(a)
a<-table.element(a,i-1,header=TRUE)
a<-table.element(a,round(myacf[i],6))
mytstat <- myacf[i]*sqrtn
a<-table.element(a,round(mytstat,4))
a<-table.element(a,round(1-pt(mytstat,lengthx),6))
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,'Partial Autocorrelation Function',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Time lag k',header=TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/basics.htm','PACF(k)','click here for more information about the Partial Autocorrelation Function'),header=TRUE)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,'P-value',header=TRUE)
a<-table.row.end(a)
for (i in 1:par1) {
a<-table.row.start(a)
a<-table.element(a,i-1,header=TRUE)
a<-table.element(a,round(mypacf[i],6))
mytstat <- mypacf[i]*sqrtn
a<-table.element(a,round(mytstat,4))
a<-table.element(a,round(1-pt(mytstat,lengthx),6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.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:

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