Welch's t-test: Difference between revisions

=={{header|Octave}}==

{{trans|Stata}}

<lang octave>x = [3.0,4.0,1.0,2.1];

y = [490.2,340.0,433.9];

n1 = length(x);

n2 = length(y);

v1 = var(x);

v2 = var(y);

t = (mean(x)-mean(y))/(sqrt(v1/n1+v2/n2));

df = (v1/n1+v2/n2)^2/(v1^2/(n1^2*(n1-1))+v2^2/(n2^2*(n2-1)));

p = betainc(df/(t^2+df),df/2,1/2);

[t df p]

ans =

-9.559498 2.000852 0.010752</lang>