Welch's t-test: Difference between revisions
+octave
(clarified that this page is about integration, is partly a translation of an R function) |
(+octave) |
||
Line 665:
0.010751
</pre>
=={{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>
=={{header|Perl}}==
| |||