DataMuseum.dk

top - download

    Length: 3744 (0xea0)
    Types: TextFile
    Notes: Mikados_K
    Names: »VARIANS.K«

└─⟦eb89399bc⟧ Bits:30008990 SOM ISFORIG, MEN KUN K-FILER
    └─⟦this⟧ »VARIANS.K« 

Mikados K File

PROGRAM VARIANSANALYSE;
CONST  M=2;
       N=6;
       SS=4;
TYPE
SAMPLE=RECORD
       X:ARRAY (1..SS) OF REAL;
       MEAN:REAL;
       KVADSUM:REAL
END;
VAR    FIL:TEXT;
       F:STRING(20);
       I,J,K:INTEGER;
       R,S,T:REAL;
       A:ARRAY (1..M,1..N) OF SAMPLE;
       RMEAN:ARRAY (1..N) OF REAL;
       SMEAN:ARRAY (1..M) OF REAL;
       TMEAN:REAL;
PROCEDURE INDLÆS;
BEGIN
  TMEAN:=0;
  F:='ØRREDER:P2:0:K';
  RESET(FIL,F);
  FOR I:=1 TO M DO SMEAN(I):=0;
  FOR J:=1 TO N DO
  BEGIN
    RMEAN(J):=0;
    FOR I:=1 TO M DO
    BEGIN
      READLN(FIL);
      R:=0;
      FOR K:=1 TO SS DO
      BEGIN
        READ(FIL,A(I,J).X(K));
        R:=R+A(I,J).X(K)
      END;
      A(I,J).MEAN:=R/SS;
      RMEAN(J):=RMEAN(J)+A(I,J).MEAN;
      SMEAN(I):=SMEAN(I)+A(I,J).MEAN;
      TMEAN:=TMEAN+A(I,J).MEAN;
      R:=0;
      FOR K:=1 TO SS DO R:=R+SQR(A(I,J).X(K)-A(I,J).MEAN);
      A(I,J).KVADSUM:=R;
    END;
    RMEAN(J):=RMEAN(J)/M
  END;
  FOR I:=1 TO M DO SMEAN(I):=SMEAN(I)/N;
  TMEAN:=TMEAN/M/N;
END;
PROCEDURE PRINTOBS;
BEGIN
FOR I:=1 TO M DO WRITE(LIST,I:12,'!':17);
WRITELN(LIST);
WRITELN(LIST,'!':29,'!':29);
  FOR J:=1 TO N DO
  BEGIN
    K:=1;
    REPEAT
      FOR I:=1 TO M DO
        WRITE(LIST,A(I,J).X(K):12:3,A(I,J).X(K+1):12:3,'!':5);
      WRITELN(LIST);
      K:=K+2
    UNTIL K>SS;
    WRITELN(LIST,J,'Mean':11,'Kvadsum':12,'!':5,
'Mean':12,'Kvadsum':12,'!':5,'R-Mean':12);
    FOR I:=1 TO M DO WRITE(LIST,A(I,J).MEAN:12:3,A(I,J).KVADSUM:12:3,'!':5);
    WRITELN(LIST,RMEAN(J):12:3);
WRITELN(LIST,'_____________________________________________________________',
'__________');
    WRITELN(LIST);
  END;
  WRITELN(LIST,'S-Mean':12,' ':12,'!':5,'S-Mean':12,'!':17,'T-Mean':12);
  FOR I:=1 TO M DO WRITE(LIST,SMEAN(I):12:3,' ':12,'!':5);
  WRITELN(LIST,TMEAN:12:3);
WRITELN(LIST,'_____________________________________________________________',
'__________');
WRITELN(LIST);
WRITELN(LIST);
WRITELN(LIST);
END;
PROCEDURE ANALYSE;
BEGIN
R:=0;
FOR I:=1 TO M DO R:=R+SQR(SMEAN(I)-TMEAN);
R:=R*N*SS;
WRITELN(LIST,'Variation','!':9,' Kvadratsum','!':7,
' Frihedsgrader','!':4,' Varians');
WRITELN(LIST,'_____________________________________________________________',
'_________');
WRITELN(LIST);
WRITELN(LIST,'Mellem søjler     ',R:12:3,' ':6,M-1:12,' ':6,R/(M-1):12:3);
WRITELN(LIST);
R:=0;
FOR J:=1 TO N DO R:=R+SQR(RMEAN(J)-TMEAN);
R:=R*M*SS;
WRITELN(LIST,'Mellem rækker     ',R:12:3,' ':6,N-1:12,' ':6,R/(N-1):12:3);
WRITELN(LIST);
R:=0;
FOR J:=1 TO N DO
FOR I:=1 TO M DO R:=R+SQR(A(I,J).MEAN-RMEAN(J)-SMEAN(I)+TMEAN);
R:=R*SS;
WRITELN(LIST,'Vekselvirkning    ',R:12:3,' ':6,(N-1)*(M-1):12,' ':6,
                                               R/(N-1)/(M-1):12:3);
WRITELN(LIST);
R:=0;
FOR J:=1 TO N DO
FOR I:=1 TO M DO
                 R:=R+A(I,J).KVADSUM;
WRITELN(LIST,'Indenfor sæt      ',R:12:3,' ':6,N*M*(SS-1):12,' ':6,
                                               R/N/M/(SS-1):12:3);
WRITELN(LIST);
WRITELN(LIST,'_____________________________________________________________',
'__________');
WRITELN(LIST);
R:=0;
FOR J:=1 TO N DO
FOR I:=1 TO M DO
FOR K:=1 TO SS DO R:=R+SQR(A(I,J).X(K)-TMEAN);
WRITELN(LIST,'Total             ',R:12:3,' ':6,N*M*SS-1:12);
END;
BEGIN
  INDLÆS;
  PRINTOBS;
  ANALYSE;
END.