home *** CD-ROM | disk | FTP | other *** search
/ C!T ROM 2 / ctrom_ii_b.zip / ctrom_ii_b / PROGRAM / PASCAL / NRPAS13 / TTEST.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  779b  |  24 lines
  1. PROCEDURE ttest(data1: glnarray; n1: integer;
  2.       data2: glmarray; n2: integer; VAR t,prob: real);
  3. (* Programs using routine TTEST must define the types
  4. TYPE
  5.    glnarray = ARRAY [1..n1] OF real;
  6.    glmarray = ARRAY [1..n2] OF real;
  7. in the main routine, with n2 less than or equal to n1. *)
  8. VAR
  9.    i: integer;
  10.    var2,var1,svar,df,ave2,ave1: real;
  11. BEGIN
  12.    avevar(data1,n1,ave1,var1);
  13.    IF (n2 > n1) THEN BEGIN
  14.       writeln('pause in routine TTEST');
  15.       writeln('The first array should be the larger'); readln
  16.    END;
  17.    FOR i := 1 TO n2 DO data1[i] := data2[i];
  18.    avevar(data1,n2,ave2,var2);
  19.    df := n1+n2-2;
  20.    svar := ((n1-1)*var1+(n2-1)*var2)/df;
  21.    t := (ave1-ave2)/sqrt(svar*(1.0/n1+1.0/n2));
  22.    prob := betai(0.5*df,0.5,df/(df+sqr(t)))
  23. END;
  24.