Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
98.6 98 106.8 96.6 100.1 107.7 91.5 97.8 107.4 117.5 105.6 97.4 99.5 98 104.3 100.6 101.1 103.9 96.9 95.5 108.4 117 103.8 100.8 110.6 104 112.6 107.3 98.9 109.8 104.9 102.2 123.9 124.9 112.7 121.9 100.6 104.3 120.4 107.5 102.9 125.6 107.5 108.8 128.4 121.1 119.5 128.7 108.7 105.5 119.8 111.3 110.6 120.1 97.5 107.7 127.3 117.2 119.8 116.2
Data Y:
467 460 448 443 436 431 484 510 513 503 471 471 476 475 470 461 455 456 517 525 523 519 509 512 519 517 510 509 501 507 569 580 578 565 547 555 562 561 555 544 537 543 594 611 613 611 594 595 591 589 584 573 567 569 621 629 628 612 595 597
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)
0
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) 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')
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