Home » date » 2009 » Nov » 24 »

ws4

*The author of this computation has been verified*
R Software Module: /rwasp_autocorrelation.wasp (opens new window with default values)
Title produced by software: (Partial) Autocorrelation Function
Date of computation: Tue, 24 Nov 2009 03:14:44 -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/Nov/24/t1259058067wwy8oczc30xpc5v.htm/, Retrieved Tue, 24 Nov 2009 11:21:10 +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/Nov/24/t1259058067wwy8oczc30xpc5v.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 «
395.3 395.1 403.5 403.3 405.7 406.7 407.2 412.4 415.9 414 411.8 409.9 412.4 415.9 416.3 417.2 421.8 421.4 415.1 412.4 411.8 408.8 404.5 402.5 409.4 410.7 413.4 415.2 417.7 417.8 417.9 418.4 418.2 416.6 418.9 421 423.5 432.3 432.3 428.6 426.7 427.3 428.5 437 442 444.9 441.4 440.3 447.1 455.3 478.6 486.5 487.8 485.9 483.8 488.4 494 493.6 487.3 482.1 484.2 496.8 501.1 499.8 495.5 498.1 503.8 516.2 526.1 527.1 525.1 528.9 540.1 549 556 568.9 589.1 590.3 603.3 638.8 643 656.7 656.1 654.1 659.9 662.1 669.2 673.1 678.3 677.4 678.5 672.4 665.3 667.9 672.1 662.5 682.3 692.1 702.7 721.4 733.2 747.7 737.6 729.3 706.1 674.3 659 645.7
 
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


Autocorrelation Function
Time lag kACF(k)T-STATP-value
10.98502610.23670
20.96531810.03190
30.9425179.79490
40.9150659.50960
50.8852729.20
60.8542398.87750
70.8221538.54410
80.79178.22760
90.7637.92930
100.7362497.65130
110.7099587.37810
120.6835187.10330
130.6579486.83760
140.6307126.55460
150.6033396.27010
160.5757495.98340
170.547925.69410
180.5185785.38920
190.4876265.06761e-06
200.4557944.73683e-06
210.4234954.40111.3e-05
220.3906984.06034.7e-05
230.3572553.71270.000163
240.3232123.35890.000541
250.2897193.01080.001622
260.2558512.65890.004516
270.2226162.31350.011294
280.1913471.98850.024641
290.1609251.67240.048672
300.1338491.3910.083542
310.1081151.12360.131844
320.0823980.85630.196863
330.0587960.6110.271233
340.0364010.37830.352978
350.0148280.15410.438909
36-0.006041-0.06280.475028
37-0.025995-0.27010.393782
38-0.044801-0.46560.321227
39-0.063799-0.6630.254366
40-0.082817-0.86070.195667
41-0.100866-1.04820.148437
42-0.117687-1.2230.111989
43-0.133947-1.3920.083389
44-0.149249-1.5510.061908
45-0.164615-1.71070.045001
46-0.179951-1.87010.032088
47-0.195627-2.0330.022252
48-0.210658-2.18920.015366
49-0.225576-2.34430.010447
50-0.240824-2.50270.006911
51-0.254929-2.64930.004638
52-0.268814-2.79360.003084
53-0.282116-2.93180.002057
54-0.295114-3.06690.001367
55-0.307961-3.20040.000901
56-0.319959-3.32510.000604
57-0.330759-3.43730.000418
58-0.340233-3.53583e-04
59-0.347705-3.61350.00023
60-0.354931-3.68860.000177


Partial Autocorrelation Function
Time lag kPACF(k)T-STATP-value
10.98502610.23670
2-0.166766-1.73310.042968
3-0.091679-0.95280.171422
4-0.145774-1.51490.066356
5-0.048709-0.50620.306874
6-0.021972-0.22830.409907
7-0.020131-0.20920.417339
80.0636890.66190.254732
90.039530.41080.341015
100.0371050.38560.350273
11-0.038445-0.39950.345144
12-0.051669-0.5370.296199
13-0.006139-0.06380.474625
14-0.091126-0.9470.172874
15-0.003826-0.03980.484179
16-0.015689-0.1630.435394
170.0041320.04290.482914
18-0.058678-0.60980.271638
19-0.0634-0.65890.255691
20-0.032646-0.33930.367534
21-0.026435-0.27470.392028
22-0.022507-0.23390.407752
23-0.038946-0.40470.343235
24-0.0305-0.3170.375939
250.0056490.05870.476646
26-0.050721-0.52710.2996
27-0.004584-0.04760.481044
280.0271520.28220.389177
29-0.010246-0.10650.457701
300.0748080.77740.219304
31-0.023827-0.24760.402451
32-0.049788-0.51740.302962
330.0169660.17630.430187
34-0.013247-0.13770.445382
35-0.00241-0.0250.490034
36-0.012296-0.12780.44928
370.0223240.2320.408488
380.0130740.13590.446088
39-0.039457-0.410.341293
40-0.029876-0.31050.378399
41-0.003417-0.03550.485869
420.0284780.2960.383917
43-0.020677-0.21490.415133
440.0033350.03470.486207
45-0.033805-0.35130.36302
46-0.036777-0.38220.351534
47-0.056949-0.59180.277599
48-0.013187-0.1370.445626
49-0.013945-0.14490.442521
50-0.030816-0.32020.3747
510.017620.18310.427526
52-0.028753-0.29880.38283
53-0.013761-0.1430.443274
54-0.035975-0.37390.354618
55-0.04505-0.46820.3203
56-0.00225-0.02340.490696
570.0214360.22280.41207
580.0217030.22550.41099
590.0337040.35030.363411
60-0.036328-0.37750.35326
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Nov/24/t1259058067wwy8oczc30xpc5v/1m6qc1259057683.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/24/t1259058067wwy8oczc30xpc5v/1m6qc1259057683.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/24/t1259058067wwy8oczc30xpc5v/22tlp1259057683.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/24/t1259058067wwy8oczc30xpc5v/22tlp1259057683.ps (open in new window)


 
Parameters (Session):
par1 = 60 ; par2 = 1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ; par6 = MA ; par7 = 0.95 ;
 
Parameters (R input):
par1 = 60 ; par2 = 1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ; par6 = MA ; par7 = 0.95 ;
 
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 (par6 == 'White Noise') par6 <- 'white' else par6 <- 'ma'
par7 <- as.numeric(par7)
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='time lag', ylab='ACF', ci.type=par6, ci=par7, sub=paste('(lambda=',par2,', d=',par3,', D=',par4,', CI=', par7, ', CI type=',par6,')',sep=''))
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 2:(par1+1)) {
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(abs(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,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(abs(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:

  • 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