tq.m

来自「统计分析的软件包」· M 代码 · 共 13 行

function q = tq(p,v);
%TQ	t distribution quantiles.
%	Q = TQ(P,V) satisfies Pr(T < Q) = P where T follows a
%	t-distribution on V degrees of freedom.
%	V must be a scalar.

%	Gordon Smyth, University of Queensland, [e-mail]
%	2 August 1998

if v <= 0, error('Degrees of freedom must be positive.'); end;

q = sign(p-0.5)*sqrt( fq(2*abs(p-0.5),1,v) );