Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_autocorrelation.wasp
Title produced by software(Partial) Autocorrelation Function
Date of computationTue, 28 Aug 2012 10:14:16 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Aug/28/t1346163283v6v17ib0kzeuwa6.htm/, Retrieved Fri, 03 May 2024 10:05:02 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=169523, Retrieved Fri, 03 May 2024 10:05:02 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact130
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [(Partial) Autocorrelation Function] [ACF1] [2012-08-28 14:14:16] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
293403	111	74	91256	123	119
277108	70	69	86997	64	64
264020	76	76	55709	101	100
260646	109	60	75741	104	104
246100	81	89	92046	135	135
244051	67	111	84607	130	124
241329	54	57	73586	93	93
234730	106	116	162365	159	155
234509	125	122	70817	125	120
233482	68	90	59635	81	78
233406	96	85	109104	117	117
228548	106	65	120087	205	198
223914	104	89	72631	115	110
223696	88	82	104911	115	114
223004	87	84	85224	147	137
213765	84	56	58233	150	150
210554	81	73	117986	126	124
202204	44	79	67271	61	56
199512	75	59	55071	82	82
195304	93	47	114425	152	145
191467	76	75	79194	109	104
191381	87	71	101653	210	212
191276	112	90	81493	151	141
190410	84	107	64664	96	94
188967	86	75	63717	98	94
188780	98	85	72369	98	98
185139	121	83	86281	128	126
185039	94	73	63958	100	98
184217	69	45	73795	74	74
181853	87	93	96750	92	91
181379	92	123	83038	101	96
181344	75	114	65196	109	108
179562	76	89	62932	116	116
178863	86	78	57637	88	87
178140	56	91	70111	83	78
176789	115	66	123328	149	149
176460	97	55	38885	122	122
175877	95	81	54628	96	95
175568	106	80	74482	105	102
174107	49	71	76168	95	91
173587	70	70	71170	97	95
173260	41	78	37238	16	15
172684	87	112	101773	103	102
167845	105	77	103646	145	145
167131	71	69	37048	56	56
167105	56	32	85903	75	71
166790	49	59	43460	46	46
164767	51	87	90257	81	80
162810	49	76	70027	83	80
162336	111	84	111436	153	151
161678	75	59	65911	87	83
158980	84	75	105965	123	122
157250	84	106	61704	104	104
156833	79	73	48204	85	85
155383	83	75	60029	99	99
154991	63	87	52295	99	98
154730	78	82	82204	98	98
151503	93	83	56316	99	98
146455	65	68	95556	127	128
143937	98	66	78792	140	139
142339	75	67	125410	144	142
142146	108	88	76013	152	139
142141	73	87	91939	61	61
142069	66	88	57231	83	82
141933	90	75	51370	100	99
139350	70	79	99518	89	88
139144	57	76	56530	75	75
137793	70	78	56699	77	77
136911	95	86	74349	117	103
136548	89	62	83042	158	157
135171	80	61	71181	82	82
134043	54	69	55901	57	54
131876	27	83	38417	36	36
131122	56	50	65724	89	89
130539	60	47	48821	66	66
130533	64	76	85168	78	79
130232	102	83	55027	107	105
129100	38	60	73713	87	87
128655	75	70	79774	111	108
128066	42	48	42564	80	80
127619	49	50	36311	52	50
127324	79	87	56733	104	101
126683	71	123	63262	72	71
126681	39	90	94137	67	66
125971	61	45	38439	71	71
125366	69	22	34497	68	68
122433	51	91	58425	66	66
121135	50	51	42051	69	68
119291	83	38	64102	123	120
118958	52	68	54506	61	58
118807	56	81	55827	70	70
118372	72	35	66477	142	145
116900	42	36	28340	58	57
116775	30	83	73087	124	112
115199	84	54	51360	87	87
114928	44	72	53009	96	91
114397	70	65	55064	87	85
113337	58	37	63016	68	68
111664	55	59	38650	98	98
108715	64	35	40671	80	78
107342	77	53	82043	116	111
107335	48	61	49319	65	64
106539	36	68	77411	63	63
105615	57	70	202316	51	48
105410	62	72	89041	88	86
105324	42	71	26982	46	46
103012	30	37	29467	28	26
102531	46	63	40001	64	63
101324	81	104	70780	103	100
100885	39	29	49288	49	48
100672	38	69	50466	55	55
99946	106	80	99501	125	119
99768	24	62	15430	27	27
99246	27	63	37361	52	51
98599	48	55	36252	46	44
98030	30	41	31701	35	35
94763	94	75	56979	100	99
93340	41	63	43448	60	60
93125	30	29	50838	37	36
91185	57	66	21067	67	67
90961	42	78	63785	49	49
90938	40	51	37137	43	42
89318	75	78	44970	82	81
88817	70	60	46765	56	56
84944	54	72	54565	90	89
84572	43	82	72571	84	84
84256	97	58	59155	76	75
80953	49	27	56622	59	58
78800	20	66	33032	21	21
78776	30	18	26998	34	34
75812	28	57	35606	30	30
75426	3	19	47261	36	33
74398	41	30	31258	51	51
74112	28	54	174949	52	52
73567	37	31	23238	18	18
69471	22	63	22618	26	25
68948	31	47	35838	45	43
67746	18	35	62832	58	56
67507	101	112	78956	49	49
65029	21	61	32551	21	21
64320	16	56	62147	24	23
61857	23	30	25162	31	28
61499	28	75	36990	15	15
50999	2	66	63989	8	8
46660	12	13	6179	13	13
43287	13	64	43750	49	49
38214	16	21	8773	16	16
35523	0	53	52491	33	33
32750	1	22	22807	5	5
31414	18	9	14116	39	39
24188	8	7	5950	7	7
22938	12	0	1168	11	11
21054	4	0	855	4	4
17547	0	4	3926	3	3
14688	4	0	6023	5	5
7199	7	0	1644	6	6
969	0	0	0	0	0
455	0	0	0	0	0
203	0	0	0	0	0
98	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 0 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=169523&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=169523&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=169523&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net



Parameters (Session):
par1 = Default ; par2 = 1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ; par6 = White Noise ; par7 = 0.95 ;
Parameters (R input):
par1 = Default ; par2 = 1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ; par6 = White Noise ; par7 = 0.95 ; par8 = ;
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 (par8 != '') par8 <- as.numeric(par8)
ox <- x
if (par8 == '') {
if (par2 == 0) {
x <- log(x)
} else {
x <- (x ^ par2 - 1) / par2
}
} else {
x <- log(x,base=par8)
}
if (par3 > 0) x <- diff(x,lag=1,difference=par3)
if (par4 > 0) x <- diff(x,lag=par5,difference=par4)
bitmap(file='picts.png')
op <- par(mfrow=c(2,1))
plot(ox,type='l',main='Original Time Series',xlab='time',ylab='value')
if (par8=='') {
mytitle <- paste('Working Time Series (lambda=',par2,', d=',par3,', D=',par4,')',sep='')
mysub <- paste('(lambda=',par2,', d=',par3,', D=',par4,', CI=', par7, ', CI type=',par6,')',sep='')
} else {
mytitle <- paste('Working Time Series (base=',par8,', d=',par3,', D=',par4,')',sep='')
mysub <- paste('(base=',par8,', d=',par3,', D=',par4,', CI=', par7, ', CI type=',par6,')',sep='')
}
plot(x,type='l', main=mytitle,xlab='time',ylab='value')
par(op)
dev.off()
bitmap(file='pic1.png')
racf <- acf(x, par1, main='Autocorrelation', xlab='time lag', ylab='ACF', ci.type=par6, ci=par7, sub=mysub)
dev.off()
bitmap(file='pic2.png')
rpacf <- pacf(x,par1,main='Partial Autocorrelation',xlab='lags',ylab='PACF',sub=mysub)
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('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('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')