R语言实现多种组合格兰杰因果检验

如何用R语言实现格兰杰因果检验？包括之前的单位根等检验。主要是希望可以省略在EVIEWS中确定滞后项之类的步骤。