home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / PEARSN.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.2 KB  |  43 lines
  1. PROGRAM d13r15(input,output);
  2. (* driver for routine PEARSN *)
  3. CONST
  4.    n=10;
  5. TYPE
  6.    nparray = ARRAY [1..n] OF real;
  7. VAR
  8.    prob,r,z : real;
  9.    i : integer;
  10.    dose,spore : nparray;
  11.  
  12. (*$I MODFILE.PAS *)
  13. (*$I GAMMLN.PAS *)
  14.  
  15. (*$I BETACF.PAS *)
  16.  
  17. (*$I BETAI.PAS *)
  18.  
  19. (*$I PEARSN.PAS *)
  20.  
  21. BEGIN
  22.    dose[1] := 56.1; dose[2] := 64.1; dose[3] := 70.0;
  23.    dose[4] := 66.6; dose[5] := 82.0; dose[6] := 91.3;
  24.    dose[7] := 90.0; dose[8] := 99.7; dose[9] := 115.3;
  25.    dose[10] := 110.0;
  26.    spore[1] := 0.11; spore[2] := 0.40; spore[3] := 0.37;
  27.    spore[4] := 0.48; spore[5] := 0.75; spore[6] := 0.66;
  28.    spore[7] := 0.71; spore[8] := 1.20; spore[9] := 1.01;
  29.    spore[10] := 0.95;
  30.    writeln;
  31.    writeln('Effect of Gamma Rays on Man-in-the-Moon Marigolds');
  32.    writeln('Count Rate (cpm)':16,'Pollen Index':23);
  33.    FOR i := 1 to 10 DO BEGIN
  34.       writeln(dose[i]:10:2,spore[i]:25:2)
  35.    END;
  36.    pearsn(dose,spore,10,r,prob,z);
  37.    writeln;
  38.    writeln('PEARSN':30,'Expected':16);
  39.    writeln('Corr. Coeff.',' ':8,r:13,'  ',0.9069586:13);
  40.    writeln('Probability',' ':9,prob:13,'  ',0.2926505e-3:13);
  41.    writeln('Fisher''s z',' ':10,z:13,'  ',1.510110:13);
  42. END.
  43.