如何用R语言做面板格兰杰因果检验？

Panel.set<-plm.data(Pdata_2, index = c("city", "year"))
#grangertest
library(lmtest)
grangertest(A~B)