我用Eviews做PP检验的蒙特卡洛模拟，思路比较清楚，但是结果存在问题，请各位高手指教，或者给出一个正确的程序，多谢。下面是我的程序：

!N=10000
workfile dfvalue u 1 !N
series df
equation eq01
series e
scalar lambuda=0
 for !i=1 to !N
 smpl 1 100
 series u=@nrnd
 u(1)=0
series y
vector(24) gamma=0
y(1)=0
smpl 2 100
y=y(-1)+u
eq01.ls d(y) y(-1)^3
e=resid
!h=2
while !h<=25
 !ee=0
   for !t=!h+1 to 100
        !ee=!ee+e(!t-!h)*e(!t)
   next
 gamma(!h-1)=!ee/100
 lambuda=lambuda+(1-!h/25)*gamma(!h-1)
!h=!h+1
wend
lambuda=2*lambuda+eq01.@ssr*0.99
smpl @all
df(!i)=@sqrt(eq01.@ssr*0.99/lambuda)*eq01.@tstats(1)-50*(lambuda-eq01.@ssr*0.99)*eq01.@se/@sqrt(lambuda)/@sqrt(eq01.@ssr)
next