基于马尔可夫的贝叶斯图向量自回归模型，可用于国际间通货膨胀、股票波动等仿真分析

library("bnlearn") library("qgraph") names=c("US","JP","UK","GE","FR","IT","CA") dag1=empty.graph(names) ma1=c(1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1) mat1=matrix(ma1,ncol=7,byrow=FALSE,dimnames =list(names,names)) qgraph(mat1,fade=TRUE,layout='circular',color="#61D04F") ma2=c(1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1) mat2=matrix(ma2,ncol=7,byrow=FALSE,dimnames =list(names,names)) qgraph(mat2,fade=TRUE,layout='circular',color="#61D04F") ma3=c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0) mat3=matrix(ma3,ncol=7,byrow=FALSE,dimnames =list(names,names)) qgraph(mat3,fade=TRUE,layout='circular',color="#61D04F") ma4=c(0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) mat4=matrix(ma4,ncol=7,byrow=FALSE,dimnames =list(names,names)) qgraph(mat4,fade=TRUE,layout='circular',color="#61D04F")