非常受教!本人按照文章的介绍用r软件编了协整残差的临界值,与大家分享一下!但是确实没有搞明白该实验设计的思想是什么,请高手指教!
- result<-c()
- j<-1
- while(j<=10000){
- u2t<-rnorm(1000)##生成解释变量x
- xxt<-c()
- xt<-c()
- xxt[1]<-u2t[1]
- i<-2
- while(i<=1000){
- xxt[i]<-xxt[i-1]+u2t[i]
- i<-i+1
- }
- xt<-xxt[951:1000]##解释变量
- ut<-rnorm(1000)##生成残差
- eet<-c()
- et<-c()
- eet[1]<-ut[1]
- i<-2
- while(i<=1000){
- eet[i]<-eet[i-1]+ut[i]
- i<-i+1
- }
- et<-eet[951:1000]##生成被解释变量
- rt<-lm(et~xt)$residuals
- result[j]<-summary(lm(diff(rt)~(rt[1:49])-1))$coefficients[1,3]
- j<-j+1
- }
- plot(density(result))
- quantile(result,c(0.01,0.05,0.1))
复制代码