Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
2.86 2.55 2.28 2.26 2.57 3.08 2.76 2.51 2.87 3.14 3.12 3.16 2.48 2.57 2.88 2.63 2.38 1.69 1.96 2.19 1.87 1.6 1.63 1.22 1.21 1.49 1.64 1.66 1.77 1.82 1.78 1.28 1.29 1.37 1.12 1.51 2.24 2.94 3.09 3.46 3.64 4.39 4.15 5.21 5.8 5.91 5.39 5.46 4.72 3.14 2.63 2.32 1.93 0.62 0.6 -0.37 -1.1 -1.68 -0.78 -1.19
Data Y:
9 8.6 8.8 8.5 8.3 8.2 8 7.9 8 9.3 9.6 9 8.7 8.3 8.4 7.8 7.8 7.6 7.7 7.6 7.6 8.6 8.6 8.2 7.5 7.1 7 6.9 6.6 6.3 6.1 5.9 6 7.2 7.2 6.4 6.1 5.9 6.1 5.9 5.8 5.7 5.6 5.3 5.5 6.5 6.5 6.1 5.9 5.8 6.2 6.5 6.6 6.7 6.6 6.5 6.8 7.8 7.9 7.4
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation