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

(Source file text) 

---

C	NPAR2.FOR - NONPARAMETRIC STATISTICS USING:
C  WTEST	- KENDALL COEFFICIENT OF CONCORDANCE
C  QTEST	- COCHRAN Q-TEST
C  MOMEN	- FIRST FOUR MOMENTS
C  ...USES RANK, TIE
	DIMENSION A(8,3),A1(24),TEN(10),R(8,3),UBO(3),ANS(4)
	EQUIVALENCE (A,A1)
	DATA A/4.,2.,1.,7.,6.,3.,5.,8.,7.,2.,1.,6.,4.,5.,3.,8.,
	1 7.,4.,2.,6.,5.,3.,1.,8./
C
C  OUTPUT CHANNEL = IOUT
	IOUT=2
C
C	WTEST
	CALL WTEST(A,R,8,3,TEN,WKEN,CHIS,NDF,0)
	WRITE(IOUT,100) ((A(I,J),J=1,3),I=1,8),WKEN,CHIS,NDF
100	FORMAT(/' KENDALL WTEST: 8 VARIABLES, 3 CASES:'/
	1 8(3F10.0/)/' COEF. OF CONCORDANCE =',F8.4,/
	2' CHI-SQUARE =',F8.3,' SHOULD BE 0'/' DEG OF FREEDOM =',I4)
C
C	MOMEN
	DO 1 I=1,10
1	TEN(I)=0.
	DO 2 J=1,24
	K=IFIX(A1(J))
2	TEN(K)=TEN(K)+1.
	WRITE(IOUT,101) (K,TEN(K),K=1,10)
101	FORMAT(//' MOMENTS - FREQUENCY DIST. OF ABOVE ARRAY:'/
	1 ' NUMBER  FREQUENCY'/(I5,6X,F4.0))
	UBO(1)=0.
	UBO(2)=1.
	UBO(3)=10.
	CALL MOMEN(TEN,UBO,5,ANS)
	WRITE(IOUT,102) (K,ANS(K),K=1,4)
102	FORMAT(/' MOMENT     VALUE'/(I5,3X,F9.3))
C
C	QTEST
	XMEAN=ANS(1)
	DO 3 I=1,24
	IF(A1(I)-XMEAN) 4,4,5
4	A1(I)=0.
	GO TO 3
5	A1(I)=1.
3	CONTINUE
	A(6,1)=1.
	WRITE(IOUT,103) ((A(I,J),J=1,3),I=1,8)
103	FORMAT(//' Q-TEST - DICHOTOMOUS DATA (8 SETS/3 GROUPS)'//
	1 8(3F8.0/)/)
	CALL QTEST(A1,8,3,Q,NDF)
	WRITE(IOUT,104) Q,NDF
104	FORMAT(' COCHRAN Q-STATISTIC ='F8.3,' DEG OF FREEDOM ='I2/)
	STOP
	END

---