Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
100.00 99.97 101.03 101.00 101.30 101.43 101.49 102.14 102.58 102.34 102.07 101.83 102.14 102.58 102.63 102.74 103.32 103.27 102.48 102.14 102.07 101.69 101.15 100.90 101.77 101.93 102.27 102.49 102.80 102.82 102.83 102.89 102.87 102.67 102.96 103.22 103.53 104.63 104.63 104.17 103.93 104.01 104.16 105.22 105.85 106.21 105.77 105.63 106.49 107.51 110.43 111.42 111.58 111.34 111.08 111.66 112.36 112.31 111.52 110.87 111.13 112.71 113.25 113.09 112.55 112.87 113.59 115.14 116.38 116.50 116.25 116.73
Data Y:
100.00 99.95 102.07 102.02 102.63 102.88 103.01 104.33 105.21 104.73 104.17 103.69 104.33 105.21 105.31 105.54 106.70 106.60 105.01 104.33 104.17 103.42 102.33 101.82 103.57 103.90 104.58 105.03 105.67 105.69 105.72 105.84 105.79 105.39 105.97 106.50 107.13 109.36 109.36 108.42 107.94 108.10 108.40 110.55 111.81 112.55 111.66 111.38 113.10 115.18 121.07 123.07 123.40 122.92 122.39 123.55 124.97 124.87 123.27 121.96 122.49 125.68 126.76 126.44 125.35 126.01 127.45 130.58 133.09 133.34 132.84 133.80
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)
1
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)
1
0
1
2
Degree (D) of seasonal differencing (Y series)
0
0
1
2
Number of non-seasonal time lags in test
1
1
2
3
4
5
6
7
8
9
10
11
Chart options
Label y-axis:
Label x-axis:
R Code
library(lmtest) 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) par8 <- as.numeric(par8) ox <- x oy <- y 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 (gyx <- grangertest(y ~ x, order=par8)) (gxy <- grangertest(x ~ y, order=par8)) bitmap(file='test1.png') op <- par(mfrow=c(2,1)) (r <- ccf(ox,oy,main='Cross Correlation Function (raw data)',ylab='CCF',xlab='Lag (k)')) (r <- ccf(x,y,main='Cross Correlation Function (transformed and differenced)',ylab='CCF',xlab='Lag (k)')) par(op) dev.off() bitmap(file='test2.png') op <- par(mfrow=c(2,1)) acf(ox,lag.max=round(length(x)/2),main='ACF of x (raw)') acf(x,lag.max=round(length(x)/2),main='ACF of x (transformed and differenced)') par(op) dev.off() bitmap(file='test3.png') op <- par(mfrow=c(2,1)) acf(oy,lag.max=round(length(y)/2),main='ACF of y (raw)') acf(y,lag.max=round(length(y)/2),main='ACF of y (transformed and differenced)') par(op) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Granger Causality Test: Y = f(X)',5,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Model',header=TRUE) a<-table.element(a,'Res.DF',header=TRUE) a<-table.element(a,'Diff. DF',header=TRUE) a<-table.element(a,'F',header=TRUE) a<-table.element(a,'p-value',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Complete model',header=TRUE) a<-table.element(a,gyx$Res.Df[1]) a<-table.element(a,'') a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Reduced model',header=TRUE) a<-table.element(a,gyx$Res.Df[2]) a<-table.element(a,gyx$Df[2]) a<-table.element(a,gyx$F[2]) a<-table.element(a,gyx$Pr[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Granger Causality Test: X = f(Y)',5,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Model',header=TRUE) a<-table.element(a,'Res.DF',header=TRUE) a<-table.element(a,'Diff. DF',header=TRUE) a<-table.element(a,'F',header=TRUE) a<-table.element(a,'p-value',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Complete model',header=TRUE) a<-table.element(a,gxy$Res.Df[1]) a<-table.element(a,'') a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Reduced model',header=TRUE) a<-table.element(a,gxy$Res.Df[2]) a<-table.element(a,gxy$Df[2]) a<-table.element(a,gxy$F[2]) a<-table.element(a,gxy$Pr[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable2.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