File: NONPAR.FT of Tape: Various/ETH/eth11-2

(Source file text) 

---

C	NONPAR.FOR - SAMPLE PROGRAM (NONPARAMETRIC STATISTICS) USING:
C  CHISQ	- CHI-SQUARE TEST FOR A CONTINGENCY TABLE
C  UTEST	- MANN-WHITNEY U-TEST
C  TWOAV	- TWO WAY ANALYSIS OF VARIANCE
C  SRANK	- SPEARMAN RANK CORRELATION
C  KRANK	- KENDALL RANK CORRELATION
C  ...ALSO...
C  RANK		- RANK OBSERVATIONS
C  TIE		- CALCULATION OF TIES IN RANKED OBSERVATIONS
C  GAUSS	- NORMAL RANDOM GENERATOR
C		   WHICH USES "RANDU" IN FORTRAN LIBRARY
C
C  THIS PROGRAM MERELY TESTS THESE ROUTINES AND IS NOT TO BE CONSIDERED
C	A TYPICAL APPLICATION
C
	DIMENSION A2(8,3),R2(8,3),A(65),B(65),WRK(25),A3(4,3)
	EQUIVALENCE (A2,A3)
C
C  OUTPUT CHANNEL = IOUT, INPUT CHANNEL = IN
	IN = 1
	IOUT = 2
C
C	CHI-SQUARE TEST
	NROW = 4
	MCOL = 3
	READ(IN,99) JUNK
99	FORMAT(A1)
C	JUST TO CLEAR AWAY LABEL CARD
	READ(IN,100) (WRK(I),WRK(I+4),WRK(I+8),I=1,4)
100	FORMAT(3F10.0)
	CALL CHISQ(WRK,NROW,MCOL,CHIS,NDF,IERR,A,B)
	WRITE(IOUT,101) NROW,MCOL,(WRK(I),WRK(I+4),WRK(I+8),I=1,4)
101	FORMAT(/' CHI-SQUARE STATISTICS FOR',I3,' BY',I3,
	1' CONTINGENCY TABLE:'//4(3F10.1/))
	WRITE(IOUT,102) CHIS,NDF
102	FORMAT(' CHI-SQUARE =',F8.3,' ACTUAL(0.003),DEG OF FREEDOM =',
	1I3,/)
C
C	MANN-WHITNEY UTEST,GAUSS
C	GENERATE 2 NORMAL SAMPLES:
C	1)  MEAN=65, SIGMA=10.5, NO. OF POINTS=30
C	2)  MEAN=70, SIGMA= 9.5, NO. OF POINTS=35
C	SEED FOR RANDOM NO. GENERATOR
	IX=13107
	DO 1 I=1,30
1	CALL GAUSS(IX,10.5,65.,A(I))
	DO 2 I=31,65
2	CALL GAUSS(IX,9.5,70.,A(I))
	WRITE(IOUT,103) (A(I),I=1,30)
103	FORMAT(//' UTEST STATISTICS FOR TWO SAMPLES:'/4X'SAMPLE #1 -
	1 30 NOS., MEAN=65, SIGMA=10.5'/(5F10.0))
	WRITE(IOUT,104) (A(I),I=31,65)
104	FORMAT(/4X'SAMPLE #2 - 35 NOS., MEAN=70, SIGMA=9.5'/(5F10.0))
	CALL UTEST(A,B,30,35,U,Z)
	WRITE(IOUT,105) U,Z
105	FORMAT(/' UTEST HOMOGENEITY:',F8.3,'  SIGNIFICANCE:',F8.3,//)
C	CLEAR CARD
	READ(IN,99) JUNK
C
C	TWOAV TEST
	MCOL=3
	NROW=4
	READ(IN,100) ((A3(I,J),J=1,3),I=1,4)
C	DATA UNRANKED
	CALL TWOAV(A3,R2,4,3,WRK,FRI,NDF,0)
	WRITE(IOUT,106)
106	FORMAT(//' FRIEDMAN TWO-WAY ANALYSIS FOR 4 GROUPS (3 CASES):'/
	13X,'GROUP    CASE 1    CASE 2    CASE 3'/)
	WRITE(IOUT,107) (I,(A3(I,J),J=1,3),I=1,4)
107	FORMAT(5X,I2,3X,3F10.0)
	WRITE(IOUT,108) FRI,NDF
108	FORMAT(/' FRIEDMAN STATISTIC ='F8.3,3X,' DEG OF FREEDOM ='I3//)
C
C	SRANK AND KRANK TESTS
	READ(IN,99) JUNK
	N=10
	READ(IN,109) (A(I),B(I),I=1,N)
109	FORMAT(2F10.0)
	WRITE(IOUT,110) (A(I),B(I),I=1,N)
110	FORMAT(//' SPEARMAN AND KENDALL RANK CORRELATIONS'/
	1'   GROUP 1   GROUP 2'/(2F10.0))
	CALL SRANK(A,B,WRK,10,RS,T,NDF,0)
	CALL KRANK(A,B,WRK,N,TAU,SD,Z,0)
	WRITE(IOUT,111) RS,T,NDF,TAU,Z,SD
111	FORMAT(/' SPEARMAN COEF:',F8.3,' ITS SIGNIFICANCE;',F8.3,
	1' DEG OF FREEDOM=',I3,/' KENDALL COEF:',F8.3,
	2' ITS SIGNIFICANCE:',F8.3,' STANDARD DEV. =',F8.3,//)
	STOP
	END

---