Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
134.79 107.16 107.83 108.85 109.52 110.19 111.20 111.54 111.88 112.55 123.34 112.55 114.24 116.26 116.60 118.62 119.63 120.64 121.65 122.33 122.66 123.00 123.34 124.68 125.02 125.02 125.36 125.70 125.70 126.03 126.37 126.37 126.71 126.71 127.04 127.04 127.38 127.72 128.05 129.40 131.09 131.42 131.76 132.10 132.43 132.77 132.77 133.11
Data Y:
120.52 102.84 137.41 118.97 125.01 118.57 130.61 116.30 99.15 110.26 107.59 107.01 113.77 93.33 147.32 124.48 106.79 134.39 111.41 132.43 98.26 109.81 115.28 108.97 99.19 105.46 138.97 124.52 117.37 123.86 116.39 124.70 97.46 103.24 112.39 107.19 100.53 95.73 143.54 101.99 120.66 121.46 102.97 121.32 85.02 106.21 110.39 87.10
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
Number of non-seasonal time lags in test
6
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation