utest.sas

SAS for Windows统计分析系统教程,里面有源码

options nodate nonumber;
data utest;
  input n1 n2 x1 x2 s1 s2@@;
  d1=sqrt((s1*s1)/n1+(s2*s2)/n2);
  u1=abs(x1-x2)/d1;
  p1=1-probnorm(u1);
  d2=sqrt((s1*s1)/n2+(s2*s2)/n1);
  u2=abs(x1-x2)/d2;
  p2=1-probnorm(u2);
cards;
120 153 139.9 143.7   7.5   6.3
 50 600  4.49  4.36  0.53  0.43
700  60 134.0 119.8 59.88 48.10
;
proc print;
run;