﻿功能:t分布表
 
格式:studis(x,n)
 
说明:x可以为数据可以为矩阵,n必须为正数数据.当n不是整数时自动转化为整数求解.本程序求解返回的是概率p.
 
原理:
(1)当n<200时,studis(x,n)=1-{BetaIncomplete<矩阵运算\BetaIncomplete>}(n/2,1/2,n/(t^2+n))/2
(2)当n≥200时,studis(x,n)={Erf<矩阵运算\erf>}(z/{Sqrt<矩阵运算\Sqrt>}(2))/2+0.5
 
例子:
 
a =
[ 1.84913685258903    1.57791667318806    1.30849970006780    1.73471255401835
  1.31455755108714    0.35007148531734    0.49339727242169    0.47069499570443
  0.92218821818110    1.99381968006204    1.01760933968127    1.66493145919635
  0.80783143770314    1.50513572409057    1.96614754198405    1.85006954327695 ]
p=StuDis(a,6)//回车得到如下结果
p =
[ 0.94303683306030    0.91716964638301    0.88070541814087    0.93325871847805
  0.88166495693174    0.63087684591754    0.68036533496769    0.67275854276002
  0.80399817582184    0.95338936375401    0.82593600920542    0.92650965250028
  0.77499009286333    0.90850234713508    0.95156693507231    0.94311042292766 ]
 
//修改于2014/04/24