R怎么做系统GMM

mus08.rar (36.76 KB) 本附件包括：

• mus08.csv

2017-2-13 20:23:13 上传

如题 我做了差分GMM，code如下：
library(plm)
mydata<-read.csv("C:/Users/ADMIN/Desktop/mus08.csv",header=TRUE)
pdata<-pdata.frame(mydata,index=c("id","t"))

lwage.diffgmm<-pgmm(lwage~lag(lwage,1:2)+lag(wks,0:1)+occ+south+smsa+ind+ms+union
                   |lag(lwage,2:4)+lag(ms,2:3)+lag(union,2:3)+lag(wks,2:3),
                   data=pdata,effect="individual",model="twosteps",
                   transformation="d")
summary(lwage.diffgmm,robust = TRUE,dianostics=TRUE)

跟正确答案比对过，答案是正确的，

现需要用同组数据做系统GMM，个人认为只需要把transformation="d"改成transformation="ld",但不知道为什么结果不对，STATA做出的正确结果我也贴上来，

求大神帮助我解决困难！

请问一下楼主，括号里的1：2，0：1，2：4，2：3等是什么意思？“|”前边是回归方程，后边是什么呢？小白一只~麻烦楼主能帮忙解答一下~~