FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_cross.wasp
Title produced by softwareCross Correlation Function
Date of computationFri, 05 Dec 2008 11:45:45 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Dec/05/t12285028302pjp55qgtfssgbs.htm/, Retrieved Thu, 16 May 2024 10:41:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=29388, Retrieved Thu, 16 May 2024 10:41:18 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact221

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Cross Correlation Function] [Q7 - non stationa...] [2008-12-01 00:17:04] [57850c80fd59ccfb28f882be994e814e]
-   P     [Cross Correlation Function] [] [2008-12-05 18:45:45] [0831954c833179c36e9320daee0825b5] [Current]
-   P       [Cross Correlation Function] [Verbetering Q7 - ...] [2008-12-05 18:48:38] [57850c80fd59ccfb28f882be994e814e]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
15107
15024
12083
15761
16943
15070
13660
14769
14725
15998
15371
14957
15470
15102
11704
16284
16727
14969
14861
14583
15306
17904
16379
15420
17871
15913
13867
17823
17872
17422
16705
15991
16584
19124
17839
17209
18587
16258
15142
19202
17747
19090
18040
17516
17752
21073
17170
19440
19795
17575
16165
19465
19932
19961
17343
18924
18574
21351
18595
19823
20844
19640
17735
19814
22239
20682
17819
21872
22117
21866
Dataseries Y:
Download CSV
Histogram 
12055
12113
9617
12646
13581
12162
10970
11880
11888
12927
12300
12093
12381
12197
9455
13168
13428
11981
11885
11692
12234
14341
13131
12421
14286
12865
11160
14316
14389
14014
13419
12770
13316
15333
14243
13824
14963
13203
12199
15509
14200
15170
14058
13786
14148
16542
13588
15582
15803
14131
12923
15612
16034
16037
14038
15331
15038
17402
14993
16044
16930
15921
14417
15961
17852
16484
14216
17430
17840
17629

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=29388&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=29388&T=0

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

As an alternative you can also use a QR Code:  
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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24






Cross Correlation Function
ParameterValue
Box-Cox transformation parameter (lambda) of X series1
Degree of non-seasonal differencing (d) of X series1
Degree of seasonal differencing (D) of X series1
Seasonal Period (s)4
Box-Cox transformation parameter (lambda) of Y series1
Degree of non-seasonal differencing (d) of Y series0
Degree of seasonal differencing (D) of Y series0
krho(Y[t],X[t+k])
-15-0.0419666399714604
-140.0671585380670575
-13-0.212009064718862
-120.192307514869653
-110.245796674669391
-10-0.209129476426147
-90.0177965892701358
-80.0303734465133481
-7-0.208942990927941
-60.177319193189249
-50.284855839673097
-4-0.256754725411921
-30.0154108373769437
-20.0646371382599219
-1-0.284804758693755
00.289449925469107
10.187759495085048
2-0.196568994994037
30.0338776216006262
4-0.038099124089524
5-0.195956607974375
60.250423984782653
70.178828966984696
8-0.224169173100057
90.0709209871660156
10-0.0928972352155877
11-0.173792212403547
120.247755276898692
130.0856987124009408
14-0.10460717620366
150.0500150307488364

\begin{tabular}{lllllllll}
\hline
Cross Correlation Function \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) of X series & 1 \tabularnewline
Degree of non-seasonal differencing (d) of X series & 1 \tabularnewline
Degree of seasonal differencing (D) of X series & 1 \tabularnewline
Seasonal Period (s) & 4 \tabularnewline
Box-Cox transformation parameter (lambda) of Y series & 1 \tabularnewline
Degree of non-seasonal differencing (d) of Y series & 0 \tabularnewline
Degree of seasonal differencing (D) of Y series & 0 \tabularnewline
k & rho(Y[t],X[t+k]) \tabularnewline
-15 & -0.0419666399714604 \tabularnewline
-14 & 0.0671585380670575 \tabularnewline
-13 & -0.212009064718862 \tabularnewline
-12 & 0.192307514869653 \tabularnewline
-11 & 0.245796674669391 \tabularnewline
-10 & -0.209129476426147 \tabularnewline
-9 & 0.0177965892701358 \tabularnewline
-8 & 0.0303734465133481 \tabularnewline
-7 & -0.208942990927941 \tabularnewline
-6 & 0.177319193189249 \tabularnewline
-5 & 0.284855839673097 \tabularnewline
-4 & -0.256754725411921 \tabularnewline
-3 & 0.0154108373769437 \tabularnewline
-2 & 0.0646371382599219 \tabularnewline
-1 & -0.284804758693755 \tabularnewline
0 & 0.289449925469107 \tabularnewline
1 & 0.187759495085048 \tabularnewline
2 & -0.196568994994037 \tabularnewline
3 & 0.0338776216006262 \tabularnewline
4 & -0.038099124089524 \tabularnewline
5 & -0.195956607974375 \tabularnewline
6 & 0.250423984782653 \tabularnewline
7 & 0.178828966984696 \tabularnewline
8 & -0.224169173100057 \tabularnewline
9 & 0.0709209871660156 \tabularnewline
10 & -0.0928972352155877 \tabularnewline
11 & -0.173792212403547 \tabularnewline
12 & 0.247755276898692 \tabularnewline
13 & 0.0856987124009408 \tabularnewline
14 & -0.10460717620366 \tabularnewline
15 & 0.0500150307488364 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=29388&T=1

[TABLE]
[ROW][C]Cross Correlation Function[/C][/ROW]
[ROW][C]Parameter[/C][C]Value[/C][/ROW]
[ROW][C]Box-Cox transformation parameter (lambda) of X series[/C][C]1[/C][/ROW]
[ROW][C]Degree of non-seasonal differencing (d) of X series[/C][C]1[/C][/ROW]
[ROW][C]Degree of seasonal differencing (D) of X series[/C][C]1[/C][/ROW]
[ROW][C]Seasonal Period (s)[/C][C]4[/C][/ROW]
[ROW][C]Box-Cox transformation parameter (lambda) of Y series[/C][C]1[/C][/ROW]
[ROW][C]Degree of non-seasonal differencing (d) of Y series[/C][C]0[/C][/ROW]
[ROW][C]Degree of seasonal differencing (D) of Y series[/C][C]0[/C][/ROW]
[ROW][C]k[/C][C]rho(Y[t],X[t+k])[/C][/ROW]
[ROW][C]-15[/C][C]-0.0419666399714604[/C][/ROW]
[ROW][C]-14[/C][C]0.0671585380670575[/C][/ROW]
[ROW][C]-13[/C][C]-0.212009064718862[/C][/ROW]
[ROW][C]-12[/C][C]0.192307514869653[/C][/ROW]
[ROW][C]-11[/C][C]0.245796674669391[/C][/ROW]
[ROW][C]-10[/C][C]-0.209129476426147[/C][/ROW]
[ROW][C]-9[/C][C]0.0177965892701358[/C][/ROW]
[ROW][C]-8[/C][C]0.0303734465133481[/C][/ROW]
[ROW][C]-7[/C][C]-0.208942990927941[/C][/ROW]
[ROW][C]-6[/C][C]0.177319193189249[/C][/ROW]
[ROW][C]-5[/C][C]0.284855839673097[/C][/ROW]
[ROW][C]-4[/C][C]-0.256754725411921[/C][/ROW]
[ROW][C]-3[/C][C]0.0154108373769437[/C][/ROW]
[ROW][C]-2[/C][C]0.0646371382599219[/C][/ROW]
[ROW][C]-1[/C][C]-0.284804758693755[/C][/ROW]
[ROW][C]0[/C][C]0.289449925469107[/C][/ROW]
[ROW][C]1[/C][C]0.187759495085048[/C][/ROW]
[ROW][C]2[/C][C]-0.196568994994037[/C][/ROW]
[ROW][C]3[/C][C]0.0338776216006262[/C][/ROW]
[ROW][C]4[/C][C]-0.038099124089524[/C][/ROW]
[ROW][C]5[/C][C]-0.195956607974375[/C][/ROW]
[ROW][C]6[/C][C]0.250423984782653[/C][/ROW]
[ROW][C]7[/C][C]0.178828966984696[/C][/ROW]
[ROW][C]8[/C][C]-0.224169173100057[/C][/ROW]
[ROW][C]9[/C][C]0.0709209871660156[/C][/ROW]
[ROW][C]10[/C][C]-0.0928972352155877[/C][/ROW]
[ROW][C]11[/C][C]-0.173792212403547[/C][/ROW]
[ROW][C]12[/C][C]0.247755276898692[/C][/ROW]
[ROW][C]13[/C][C]0.0856987124009408[/C][/ROW]
[ROW][C]14[/C][C]-0.10460717620366[/C][/ROW]
[ROW][C]15[/C][C]0.0500150307488364[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=29388&T=1

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

As an alternative you can also use a QR Code:  
QR Code


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

Cross Correlation Function
ParameterValue
Box-Cox transformation parameter (lambda) of X series1
Degree of non-seasonal differencing (d) of X series1
Degree of seasonal differencing (D) of X series1
Seasonal Period (s)4
Box-Cox transformation parameter (lambda) of Y series1
Degree of non-seasonal differencing (d) of Y series0
Degree of seasonal differencing (D) of Y series0
krho(Y[t],X[t+k])
-15-0.0419666399714604
-140.0671585380670575
-13-0.212009064718862
-120.192307514869653
-110.245796674669391
-10-0.209129476426147
-90.0177965892701358
-80.0303734465133481
-7-0.208942990927941
-60.177319193189249
-50.284855839673097
-4-0.256754725411921
-30.0154108373769437
-20.0646371382599219
-1-0.284804758693755
00.289449925469107
10.187759495085048
2-0.196568994994037
30.0338776216006262
4-0.038099124089524
5-0.195956607974375
60.250423984782653
70.178828966984696
8-0.224169173100057
90.0709209871660156
10-0.0928972352155877
11-0.173792212403547
120.247755276898692
130.0856987124009408
14-0.10460717620366
150.0500150307488364


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1 ; par2 = 1 ; par3 = 1 ; par4 = 4 ; par5 = 1 ; par6 = 0 ; par7 = 0 ;
Parameters (R input):
par1 = 1 ; par2 = 1 ; par3 = 1 ; par4 = 4 ; par5 = 1 ; par6 = 0 ; par7 = 0 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
par5 <- as.numeric(par5)
par6 <- as.numeric(par6)
par7 <- as.numeric(par7)
if (par1 == 0) {
x <- log(x)
} else {
x <- (x ^ par1 - 1) / par1
}
if (par5 == 0) {
y <- log(y)
} else {
y <- (y ^ par5 - 1) / par5
}
if (par2 > 0) x <- diff(x,lag=1,difference=par2)
if (par6 > 0) y <- diff(y,lag=1,difference=par6)
if (par3 > 0) x <- diff(x,lag=par4,difference=par3)
if (par7 > 0) y <- diff(y,lag=par4,difference=par7)
x
y
bitmap(file='test1.png')
(r <- ccf(x,y,main='Cross Correlation Function',ylab='CCF',xlab='Lag (k)'))
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Cross Correlation Function',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'Value',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Box-Cox transformation parameter (lambda) of X series',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of non-seasonal differencing (d) of X series',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of seasonal differencing (D) of X series',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Seasonal Period (s)',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Box-Cox transformation parameter (lambda) of Y series',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of non-seasonal differencing (d) of Y series',header=TRUE)
a<-table.element(a,par6)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of seasonal differencing (D) of Y series',header=TRUE)
a<-table.element(a,par7)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'k',header=TRUE)
a<-table.element(a,'rho(Y[t],X[t+k])',header=TRUE)
a<-table.row.end(a)
mylength <- length(r$acf)
myhalf <- floor((mylength-1)/2)
for (i in 1:mylength) {
a<-table.row.start(a)
a<-table.element(a,i-myhalf-1,header=TRUE)
a<-table.element(a,r$acf[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')