###test causality
library(vars)
data(Canada)
myVar<-VAR(Canada, p=2, season=12)
causality(myVar)
causality(myVar, cause="e")
causality(myVar, cause="prod")
causality(myVar, cause="rw")
causality(myVar, cause="U")
causality(myVar, cause=c("e", "prod"))
causality(myVar, cause=c("prod","e"))
causality(myVar, cause=c("e", "prod","rw"))
myVar2<-VAR(Canada, p=3, type="trend")
causality(myVar2, cause="e")
causality(myVar2, cause="prod")
causality(myVar2, cause="rw")
causality(myVar2, cause="U")
causality(myVar2, cause=c("e", "prod"))
causality(myVar, cause=c("prod","e"))
causality(myVar, cause=c("e", "prod","rw"))
myVar3<-VAR(Canada[,1:3], p=1, exogen=Canada[,4], type="none")
causality(myVar3)
causality(myVar3, cause="e")
causality(myVar3, cause="prod")
causality(myVar3, cause="rw")
causality(myVar3, cause=c("e", "prod"))
causality(myVar3, cause=c("prod","e"))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.