首先，需要加载package,分别为 sp,Matrix,spdep.
  给定线性公式，如同lm(fomula)；然后需要给定空间权重矩阵W, W的格式为listW，格式很重要。
  例如：
  fomula <- y~x
  res <-  lagsarlm(fomula,data,nb2listw(COL.nb, style="W"))
  impacts(res)
  就可以得到 平均直接效应和间接效应了。
  遇到有问题的地方请观看spdep程序包的手册。

谢谢你

大神，我按照您给的代码，写完还有错误啊，我提供了listw，为什么还出现错误呢？
res <- spml(fm,data = mydata,index = c("province","year"),
            listw = mat2listw(gal.mat),
            model = "random",effect = "individual",
            lag = T,spatial.error = "b")
> impacts(res)
Error in impacts.splm(res) : either listw or tr should be provided

大神，面版数据的杜宾模型怎么算呢，代码怎么写？谢谢啦！

请问解决了吗