'山本拓 著 「計量経済学」（新世社）
'p62.一致性のプログラム例。n=50の場合。
'rndは開区間(0,1)で乱数を作成するコマンドです。
'テキストは閉区間[0,1]となっています。
'グラフができたら、X軸やY軸の範囲を調整して
'比較してみましょう。

wfcreate(wf=test) u 50
series x1
series x2
series y1
series y2

'パラメータを格納する行列
matrix(100,2) beta

series vi=rnd-0.5

for !i=1 to 100
	x1=2*vi
	x2=5*vi

	series ui=nrnd

	y1=3+0.75*x1+ui
	y2=3+0.75*x2+ui

	equation eq01.ls y1 c x1
	equation eq02.ls y2 c x2

	beta(!i,1)=eq01.@coefs(2)
	beta(!i,2)=eq02.@coefs(2)
	
next

beta.distplot