wfopen g:\bootdata.wf1 equation eq1.ls y c x eq1.makeresids resid01 genr u=resid01 !alpha=@coefs(1) !beta=@coefs(2) !n=1 !draws=999 vector(!draws) tvec=0 vector(!draws) betavec=0 for !n=1 to !draws u.resample u_b eq1.fit(u) y_hat genr yhat=y_hat+u_b equation eqb.ls yhat c x betavec(!n)=(@coefs(2)-!beta)/@sqr(1/@sumsq(x-@mean(x))) tvec(!n)=(@coefs(2)-!beta)/@stderrs(2) next wfsave g:\bootdata.wf1 wfcreate(wf="tstatbeta") u 999 fetch(d=g:\bootdata.wf1) tvec betavec series tser=0 series betaser=0 mtos(betavec,betaser) mtos(tvec,tser) betaser.hist tser.hist sort(d) tser scalar lower=tser(950) scalar upper=tser(50) table(2,2) tvalue tvalue.setwidth(1:2) 30 tvalue.setlines(a1:b2) +a tvalue.setfillcolor(a1:b2) yellow tvalue(1,1)="t-value lower (5%)" tvalue(2,1)=lower tvalue(1,2)="t-value upper (5%)" tvalue(2,2)=upper show tvalue