home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / pascal / library / dos / math / nr / cntab2.dem < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.5 KB  |  60 lines
  1. PROGRAM d13r14(input,output,dfile);
  2. (* driver for routine CNTAB2 *)
  3. (* contingency table in file TABLE1.DAT *)
  4. CONST
  5.    ni=9;
  6.    nmon=12;
  7. TYPE
  8.    narray = ARRAY [1..ni,1..nmon] OF integer;
  9.    chr5 = string[5];
  10.    chr15 = string[15];
  11. VAR
  12.    h,hx,hxgy,hy,hygx : real;
  13.    uxgy,uxy,uygx : real;
  14.    i,j : integer;
  15.    nmbr : narray;
  16.    fate : ARRAY [1..ni] OF chr15;
  17.    mon : ARRAY [1..nmon] OF chr5;
  18.    txt : string[64];
  19.    dfile : text;
  20.  
  21. (*$I MODFILE.PAS *)
  22. (*$I CNTAB2.PAS *)
  23.  
  24. BEGIN
  25.    glopen(dfile,'table1.dat');
  26.    readln(dfile);
  27.    readln(dfile,txt);
  28.    read(dfile,fate[1]);
  29.    FOR i := 1 to 12 DO read(dfile,mon[i]);
  30.    readln(dfile);
  31.    readln(dfile);
  32.    FOR i := 1 to ni DO BEGIN
  33.       read(dfile,fate[i]);
  34.       FOR j := 1 to 12 DO read(dfile,nmbr[i,j]);
  35.       readln(dfile)
  36.    END;
  37.    close(dfile);
  38.    writeln;
  39.    writeln(txt);
  40.    writeln;
  41.    write(' ':15);
  42.    FOR i := 1 to 12 DO write(mon[i]:5);
  43.    writeln;
  44.    FOR i := 1 to ni DO BEGIN
  45.       write(fate[i]);
  46.       FOR j := 1 to 12 DO write(nmbr[i,j]:5);
  47.       writeln
  48.    END;
  49.    cntab2(nmbr,ni,nmon,h,hx,hy,hygx,hxgy,uygx,uxgy,uxy);
  50.    writeln;
  51.    writeln('entropy of table           ',h:10:4);
  52.    writeln('entropy of x-distribution  ',hx:10:4);
  53.    writeln('entropy of y-distribution  ',hy:10:4);
  54.    writeln('entropy of y given x       ',hygx:10:4);
  55.    writeln('entropy of x given y       ',hxgy:10:4);
  56.    writeln('dependency of y on x       ',uygx:10:4);
  57.    writeln('dependency of x on y       ',uxgy:10:4);
  58.    writeln('symmetrical dependency     ',uxy:10:4)
  59. END.
  60.