'「自由度とt分布」
't分布はそもそも少ない標本でも仮説検定を
'正しく行うための道具です。
'今回は自由度が小さい時と大きな時で、どこに問題が
'生じるのか見てみましょう。

'ファイルの拡張子を.prgに変更して実行してください。
't分布の曲線が3つ(自由度5,30,100)ができます。
'自由度30のものと100のものでは、あまり目立った
'差はないことが分かります。

'新しいワークファイルを作成します。データの個数は100個
wfcreate(wf=degree of freedom) u 100

't値(横軸)を適当に用意します(-0.5<=t<=4.9,刻み0.1)。
series t=-5
smpl 2 100
t=t(-1)+0.1
smpl @all

'密度(縦軸)をそれぞれ求めます。
series dof5=@dtdist(t,5)
series dof30=@dtdist(t,30)
series dof100=@dtdist(t,100)

'グループオブジェクトgroup01を作成します。
group group01 t dof5 dof30 dof100
show group01.xyline

'最後にグラフのマイナス領域で面積が2.5%となる時の
't値をtableオブジェクトに入れておきます。
'興味のある方はご覧ください。

table(2,4) tvalue
tvalue(1,1)="DOF"
tvalue(1,2)="5"
tvalue(1,3)="30"
tvalue(1,4)="100"
tvalue(2,1)="t-value"
tvalue(2,2)=@qtdist(0.025,5) 
tvalue(2,3)=@qtdist(0.025,30) 
tvalue(2,4)=@qtdist(0.025,100) 

show tvalue
'自由度が小さいほど、2.5%に相当するt値が
'大きく変化することが分かります。