Home » date » 2009 » Jun » 01 »

Annelies Reul cons vertrouwen oef 6 BIS deel 2 inschr personenwagen

*Unverified author*

R Software Module: rwasp_autocorrelation.wasp (opens new window with default values)

Title produced by software: (Partial) Autocorrelation Function

Date of computation: Mon, 01 Jun 2009 15:13:26 -0600

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Jun/01/t1243890954f9xaeo6uidhbqoc.htm/, Retrieved Mon, 20 May 2013 00:47:59 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

System-generated keywords (parent):

(pk = 0)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

Autocorrelation Function
Time lag k	ACF(k)	T-STAT	P-value
1	-0.29901	-2.5195	0.007002
2	-0.125252	-1.0554	0.147412
3	0.14109	1.1888	0.11923
4	-0.134876	-1.1365	0.129788
5	0.039215	0.3304	0.371025
6	-0.092414	-0.7787	0.219375
7	-0.070657	-0.5954	0.276744
8	-0.027595	-0.2325	0.408401
9	0.081689	0.6883	0.246746
10	-0.124057	-1.0453	0.14971
11	-0.220857	-1.861	0.033444
12	0.694259	5.8499	0
13	-0.198323	-1.6711	0.049552
14	-0.098196	-0.8274	0.205388
15	0.131515	1.1082	0.135766
16	-0.173331	-1.4605	0.07428
17	0.064948	0.5473	0.292959
18	-0.061006	-0.514	0.304408
19	-0.10727	-0.9039	0.18456
20	0.077846	0.6559	0.256989
21	0.000546	0.0046	0.49817
22	-0.132975	-1.1205	0.133146
23	-0.10805	-0.9104	0.182834
24	0.504135	4.2479	3.2e-05
25	-0.158231	-1.3333	0.093352
26	-0.00396	-0.0334	0.486739
27	0.01758	0.1481	0.441328
28	-0.106521	-0.8976	0.186225
29	0.107413	0.9051	0.184242
30	-0.115098	-0.9698	0.167711
31	-0.068804	-0.5798	0.281959
32	0.057463	0.4842	0.314869
33	0.013123	0.1106	0.456133
34	-0.075043	-0.6323	0.264603
35	-0.080564	-0.6788	0.249721
36	0.301513	2.5406	0.006627
37	-0.097891	-0.8248	0.206112
38	0.022569	0.1902	0.424859
39	0.011443	0.0964	0.461731
40	-0.107815	-0.9085	0.183352
41	0.119872	1.0101	0.157949
42	-0.127614	-1.0753	0.142943
43	-0.003439	-0.029	0.488481
44	0.047851	0.4032	0.344005
45	-0.051805	-0.4365	0.331892
46	0.015052	0.1268	0.449717
47	-0.10213	-0.8606	0.196187
48	0.201693	1.6995	0.046802

Partial Autocorrelation Function
Time lag k	PACF(k)	T-STAT	P-value
1	-0.29901	-2.5195	0.007002
2	-0.235736	-1.9863	0.025428
3	0.028283	0.2383	0.406162
4	-0.120382	-1.0144	0.156929
5	-0.013668	-0.1152	0.454318
6	-0.152886	-1.2882	0.100922
7	-0.154438	-1.3013	0.098678
8	-0.196303	-1.6541	0.051263
9	-0.02271	-0.1914	0.424396
10	-0.196904	-1.6591	0.050749
11	-0.448937	-3.7828	0.00016
12	0.520609	4.3867	2e-05
13	0.138832	1.1698	0.122993
14	0.041534	0.35	0.363696
15	-0.041634	-0.3508	0.363384
16	-0.060008	-0.5056	0.307339
17	-0.121125	-1.0206	0.15545
18	-0.059157	-0.4985	0.30985
19	-0.010749	-0.0906	0.464044
20	0.0854	0.7196	0.237068
21	-0.060382	-0.5088	0.306239
22	-0.111788	-0.9419	0.174707
23	-0.023929	-0.2016	0.420392
24	0.078042	0.6576	0.256463
25	-0.044485	-0.3748	0.354451
26	0.126236	1.0637	0.145539
27	-0.143154	-1.2062	0.115865
28	0.031498	0.2654	0.395732
29	0.038083	0.3209	0.374618
30	-0.038024	-0.3204	0.374803
31	0.022141	0.1866	0.426267
32	-0.132996	-1.1206	0.133108
33	0.034598	0.2915	0.385747
34	0.111193	0.9369	0.175986
35	0.056881	0.4793	0.316603
36	-0.135844	-1.1446	0.128101
37	-0.060346	-0.5085	0.306344
38	-0.127665	-1.0757	0.142845
39	0.081697	0.6884	0.246726
40	-0.005407	-0.0456	0.481894
41	0.011957	0.1008	0.460016
42	-0.027509	-0.2318	0.408682
43	0.028352	0.2389	0.405937
44	-0.018327	-0.1544	0.438857
45	-0.025789	-0.2173	0.414299
46	-0.015715	-0.1324	0.447516
47	-0.065201	-0.5494	0.29223
48	-0.033106	-0.279	0.390543

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Jun/01/t1243890954f9xaeo6uidhbqoc/10pix1243890802.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Jun/01/t1243890954f9xaeo6uidhbqoc/10pix1243890802.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Jun/01/t1243890954f9xaeo6uidhbqoc/2e75e1243890802.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Jun/01/t1243890954f9xaeo6uidhbqoc/2e75e1243890802.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 48 ; par2 = 1 ; par3 = 1 ; par4 = 0 ; par5 = 12 ;

Parameters (R input):

par1 = 48 ; par2 = 1 ; par3 = 1 ; 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 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')