Directory of image this file is from
This file as a plain text file
C TEILPROGRAMM BARTLE VON VAR1
C MULTIPLER T-TEST UND BARTLETT-TEST
C FORTRAN II / OPTION (I)
C FS : PFT.FT / PCHI.FT
C CH : DUNCAN.SV
COMMON NG,N1,FM,VS,NZ,FKW,FMW
DIMENSION NG(49),FM(49),VS(49),NZ(49),NC(49)
SPACE=' '
C MULTIPLER T-TEST
DO 2 I=1,N1
2 NC(I)=NG(I)
WRITE(1,100)LEER,
100 FORMAT('MULTIPLER T-TEST ? (JA=1,NEIN=0) : 'I0)
READ(1,103)ID
103 FORMAT(I1)
IF(ID)12,12,9
9 WRITE(1,104)LEER
104 FORMAT(I0)
FMW=SQRT(FMW)
IF(FKW-2047.)6,7,7
6 KW=IFIX(FKW)
GOTO 8
7 KW=2047
8 DO 11 I=1,N1
DO 10 J=1,N1
IF(NC(J))10,10,3
3 IF(NC(J)-NG(I))5,4,5
4 NC(J)=0
GOTO 10
5 IG=NG(I)
S JMS NTD
C NAME TEMUXX
IN1=1285
IN2=0853
IN3=IE
CALL IOPEN('P1',IN1)
READ(4,101)IOPT,KZ
101 FORMAT(A2)
READ(4,102)EFA,EFA,PARA
102 FORMAT(A6)
IG=NG(J)
S JMS NTD
C NAME TEMUXX
IN1=1285
IN2=0853
IN3=IE
CALL IOPEN('P1',IN1)
READ(4,101)KZ,KZ
READ(4,102)EFB,EFB,PARB
FNZI=FLOAT(NZ(I))
FNZJ=FLOAT(NZ(J))
T=SQRT(FNZI*FNZJ/(FNZI+FNZJ))
T=T*(FM(I)-FM(J))/FMW
TW=ABS(T)
TT=100.*PFT(1,1,KW,TW)
WRITE(1,105)NG(I),EFA,PARA,NG(J),EFB,PARB,TT
105 FORMAT(' #'I2' ('A6'/'A6') = #'I2' ('A6'/'A6
1') P ='F7.2' %')
10 CONTINUE
WRITE(1,107)SPACE
107 FORMAT(A1)
11 CONTINUE
12 WRITE(1,107)SPACE
FN=0.
DO 15 I=1,N1
15 FN=FN+FLOAT(NZ(I))
FT=0.
DO 16 I=1,N1
16 FT=FT+1./FLOAT(NZ(I)-1)
FT=FT-1./(FN-FLOAT(N1))
C=1./(3.*FLOAT(N1-1))
C=1.+C*FT
UM=2.3025851
C=UM/C
FT=0.
DO 17 I=1,N1
17 FT=FT+(FLOAT(NZ(I)-1))*ALOG(VS(I))/UM
FMW=FMW*FMW
FT=C*((FN-FLOAT(N1))*ALOG(FMW)/UM-FT)
WRITE(1,111)
111 FORMAT('BARTLETT-TEST'/13('=')/)
M=N1-1
C=PCHI(M,FT)*100.
WRITE(1,112)C
112 FORMAT(' SIGNIFIKANZ-NIVEAU FUER GLEICHHEIT DER STREU
1UNGEN : 'F6.2' %'/)
CALL CHAIN('DUNCAN')
SNTD, 0 /NUMMER DER DATEI
IZ=(IG/10)*10
IE=IG-IZ
IZ=IZ/10
S CLA CLL
S TAD \IZ
S CLL RTL;RTL;RTL
S TAD \IE
S TAD (6060
S DCA \IE
S JMP I NTD
END