Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
15859.4 15258.9 15498.6 15106.5 15023.6 12083.0 15761.3 16942.6 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 17422.0 16704.5 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
Data Y:
12710.3 12120.8 12469.5 12054.6 12112.9 9617.2 12645.8 13581.3 12162.3 10969.7 11880.0 11887.6 12926.9 12300.0 12092.8 12380.8 12196.9 9455.0 13168.0 13427.9 11980.5 11884.8 11691.7 12233.8 14341.4 13130.7 12421.1 14285.8 12864.6 11160.2 14316.2 14388.7 14013.9 13419.0 12769.6 13315.5 15332.9 14243.0 13824.4 14962.9 13202.9 12199.0 15508.9 14199.8 15169.6 14058.0 13786.2 14147.9 16541.7 13587.5 15582.4 15802.8 14130.5 12923.2 15612.2 16033.7 16036.6 14037.8 15330.6 15038.3 17401.8
Sample Range:
(leave blank to include all observations)
From:
To:
Box-Cox transformation parameter (X series)
1
1
-2.0
-1.9
-1.8
-1.7
-1.6
-1.5
-1.4
-1.3
-1.2
-1.1
-1.0
-0.9
-0.8
-0.7
-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2.0
Degree (d) of non-seasonal differencing (X series)
0
0
1
2
Degree (D) of seasonal differencing (X series)
1
0
1
2
Seasonal Period
12
1
2
3
4
12
Box-Cox transformation parameter (Y series)
1
1
-2.0
-1.9
-1.8
-1.7
-1.6
-1.5
-1.4
-1.3
-1.2
-1.1
-1.0
-0.9
-0.8
-0.7
-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2.0
Degree (d) of non-seasonal differencing (Y series)
0
0
1
2
Degree (D) of seasonal differencing (Y series)
0
0
1
2
Treatment of missing data
(?)
na.fail
na.pass
Chart options
Label y-axis:
Label x-axis:
R Code
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) x <- diff(y,lag=par4,difference=par7) x y bitmap(file='test1.png') (r <- ccf(x,y,main='Cross Correlation Function',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')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation