代码贴出来？

library(BMR)
data(BMRVARData)
irf.points<-c(1979,1988,1997,2006)
yearlab<-seq(1955.00,2010.75,0.25)
USMacroData<-USMacroData[3:226,2:4]
bvartvptest <- BVARTVP(USMacroData,timelab=yearlab,
coefprior=NULL,tau=80,p=4,
irf.periods=20,irf.points=irf.points,
keep=30000,burnin=40000,
XiBeta=4,XiQ=0.005,gammaQ=NULL,
XiSigma=1,gammaS=NULL)

谢谢您的回复，麻烦您看下这组代码问题出在哪里，感激不尽~

BMR包附上去