Welch's t-test: Difference between revisions
Content added Content deleted
(+SAS) |
|||
Line 944: | Line 944: | ||
<pre>(0.021378001462867013 0.14884169660532798 0.035972271029796624 0.09077332428567102 0.01075139991904718)</pre> |
<pre>(0.021378001462867013 0.14884169660532798 0.035972271029796624 0.09077332428567102 0.01075139991904718)</pre> |
||
=={{header|SAS}}== |
|||
{{trans|Stata}} |
|||
<lang sas>proc iml; |
|||
x = {3.0,4.0,1.0,2.1}; |
|||
y = {490.2,340.0,433.9}; |
|||
n1 = nrow(x); |
|||
n2 = nrow(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 = 2*probt(-abs(t), df); |
|||
print(t || df || p); |
|||
quit;</lang> |
|||
'''Output''' |
|||
<pre>-9.559498 2.0008523 0.0107516</pre> |
|||
=={{header|Scilab}}== |
=={{header|Scilab}}== |
||
{{trans|Stata}} |
{{trans|Stata}} |