有没有大佬能帮忙看个R的编程，伍德里奇的计量经济学导论第五版的案例18.1，我做出来的结果和课本不一样，用stata做出来的和课本一样，感谢，有偿。
#R
a<-read.csv("HSEINV.csv")
b<-lm(linvpc~t,data=a)
c<-residuals(b)#除趋势之后，残差
library(dynlm)
a1<-read.csv("HSEINV.csv")
df<-data.frame(c,a1)#将除趋势之后的y放入数据框
tdf<-ts(df)#时间序列化数据
d<-dynlm(c~gprice+L(c),data=tdf)#做GDL回归
summary(d)