sample 1 28
read Q Y P QLAG
1.917	2560	1.363	1.86
1.981	2652	1.361	1.917
1.937	2840	1.454	1.981
1.924	2882	1.62	1.937
1.947	2900	1.688	1.924
2.039	3167	1.664	1.947
1.985	3220	1.76	2.039
2.018	3350	1.812	1.985
2.016	3443	1.984	2.018
2.114	3546	1.996	2.016
2.114	3826	1.926	2.114
2.223	4014	1.875	2.114
2.357	4109	1.869	2.223
2.574	4304	1.913	2.357
2.403	4526	2.008	2.574
2.456	4784	2.148	2.403
2.517	4869	2.075	2.456
2.593	4906	2.347	2.517
2.723	4786	1.875	2.593
2.347	4638	2.145	2.723
2.527	4714	2.279	2.347
2.413	4808	2.392	2.527
2.346	4844	2.23	2.413
2.206	5006	2.106	2.346
2.107	5132	2.714	2.206
2.014	5409	2.539	2.107
2.165	5672	2.967	2.014
2.082	5723	3.968	2.165
genr lnq=log(Q)
genr lnY=log(Y)
genr lnP=log(P)
genr lnqlag=log(qlag)
stat lnq /mean=lnmean
gen1 gmean=exp(lnmean)
ols q y p qlag
gen1 ssel=$sse
gen1 ssel=ssel/(gmean)**2
ols lnq lnY lnP lnqlag
gen1 ssell=$sse
gen1 l=($N/2)*abs(log(ssel)-log(ssell))
print l
stop