如题 我做了差分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做出的正确结果我也贴上来，

求大神帮助我解决困难！