Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap5_198111 - decus/20-0137/rmaov/rmaov.doc
There are 2 other files named rmaov.doc in the archive. Click here to see a list.
                         WESTERN MICHIGAN UNIVERSITY
                               COMPUTER CENTER
 
LIBRARY PROGRAM #1.9.4
CALLING NAME:  RMAOV
PREPARED BY:  BERENICE HOUCHARD
PROGRAMMED BY:  SAM ANEMA
STATISTICAL CONSULTANT:  DR. MICHAEL STOLINE
APPROVED BY:  JACK R. MEAGHER
DATE:  JANUARY, 1973
 
                 ANALYSIS OF VARIANCE WITH REPEATED MEASURES
 
1.0  GENERAL DESCRIPTION
 
THIS PROGRAM IS DESIGNED TO HANDLE THE FOLLOWING UNBALANCED 2-FACTOR REPEATED
MODEL:
 
                                                 B
 
                                            B(1)  ...       B(IQ)
        ------------------------------------------------------------------------
               A(1)    S(1)               X(1,1,1)   ...  X(1,1,IQ)
                       .                      .               .
                       .                      .               .
                       .                      .               .
                       S(N(1))            X(1,N(1),1) ... X(1, N(1), IQ)
        ------------------------------------------------------------------------
               A(2)    S(1)               X(2,1,1)   ...  X(2,1,IQ)
                       .                      .               .
                       .                      .               .
                       .                      .               .
                       S(N(2))            X(2,N(2),1) ... X(2,N(2),IQ)
   A
        ------------------------------------------------------------------------
                       .
                       .
                       .
        ------------------------------------------------------------------------
               A(IP)   S(1)               X(IP,1,1)   ... X(IP,1,IQ)
                       .                      .               .
                       .                      .               .
                       .                      .               .
                       S(N(IP))           S(IP,N(IP),1) ... X(IP,N(IP),IQ)
        ------------------------------------------------------------------------
 
*REFERENCE:  WINER, B.J.  "STATISTICAL PRINCIPLES IN EXPERIMENTAL DESIGN"
 (MCGRAW-HILL, 1962) PP. 306-307.
 
THERE ARE IP LEVELS ON FACTOR A WHICH IS BETWEEN SUBJECTS.
 
THERE ARE IQ LEVELS ON FACTOR B WHICH IS WITHIN SUBJECTS.
 
THERE ARE N(I) SUBJECTS FOR LEVEL (GROUP) I OF FACTOR A AND IQ REPEATED
OBSERVATIONS ON EACH SUBJECT.
 
1.1  PROGRAM OUTPUT
 
THE PROGRAM PRODUCES 2 TABLES:
 
   (A)  MEAN TABLE
   (B)  ANOVA TABLE
 
CASE 1  1-WAY REPEATED MEASURE ANOVA
 
--------------------------------------------------------------------------------
 
           SOURCE                 SS       DF       MS       F       PROB
 
--------------------------------------------------------------------------------
 
BETWEEN SUBJECTS - - -
 WITHIN SUBJECTS - - -
      B- - - - - - - -
 
 RESIDUAL  - - - - - -
   TOTAL - - - - - - -
 
--------------------------------------------------------------------------------
 
CASE 2  2-WAY REPEATED MEASURE ANOVA
 
--------------------------------------------------------------------------------
 
           SOURCE                 SS       DF       MS       F       PROB
 
--------------------------------------------------------------------------------
 
BETWEEN SUBJECTS - - - - - -
      A- - - - - - - - - - -
 
SUBJECTS WITHIN GROUPS - - -
    WITHIN SUBJECTS- - - - -
         B - - - - - - - - -
        AB - - - - - - - - -
 
B X SUBJECTS WITHIN GROUPS -
         TOTAL - - - - - - -
 
--------------------------------------------------------------------------------
 
1.2  LIMITATIONS
 
   (1)  MAXIMUM OF 30 LEVELS IN THE REPEATED FACTOR WHERE THE BREAKDOWN VARIABLE
        IS INCLUDED.
 
   (2)  MAXIMUM OF 30 GROUPS
 
   (3)  MAXIMUM OF 30 SUBJECTS PER GROUP
 
   (4)  2 LINES ARE ALLOTTED FOR OBJECT TIME FORMAT, EACH LINE BEING 80 COLUMNS
        LONG.
 
2.0  PROGRAM CAPABILITIES AND OPTIONS
 
THIS PROGRAM CAN BE PROCESSED BY BATCH OR BY TERMINAL.  THE PROGRAM IS IN
CONVERSATIONAL MODE TO ALLOW INTERACTION WITH THE USER ON TERMINAL JOBS.
 
2.1  PROGRAM OUTLINE
 
FOLLOWING IS AN OUTLINE OF THE PROGRAM OPTIONS.  THE QUESTIONS AND STATEMENTS
LISTED BELOW WILL BE TYPED OUT BY THE PROGRAM IN WHICH A RESPONSE IS IMPERATIVE
TO EACH LINE.  ON BATCH JOB, MAKE UP ONE CARD FOR EACH APPROPRIATE QUESTION OR
STATEMENT.
 
LINE 1  OUTPUT?
LINE 2  INPUT?
 
LINES 1 AND 2 DEFINE WHERE THE USER INTENDS TO WRITE HIS OUTPUT FILE (LINE 1)
AND FROM WHERE THE USER EXPECTS TO READ HIS INPUT DATA (LINE 2).  SEE NOTE (2)
BELOW FOR OTHER INPUT OPTIONS.
 
A PROPER RESPONSE TO EACH OF THESE QUESTIONS CONSISTS OF THREE BASIC PARTS:
A DEVICE, A FILENAME, AND A PROJECT-PROGRAMMER NUMBER.
 
THE GENERAL FORMAT FOR THESE THREE PARTS IS AS FOLLOWS:
 
                  DEV:FILE.EXT[PROJ,PROG]
 
1)  DEV:  ANY OF THE FOLLOWING DEVICES ARE APPROPRIATE WHERE INDICATED:
 
              DEVICE LIST             DEFINITION            STATEMENT USE
 
                  TTY:                TERMINAL              INPUT OR OUTPUT
                  DSK:                DISK                  INPUT OR OUTPUT
                  CDR:                CARD READER           INPUT ONLY
                  LPT:                LINE PRINTER          OUTPUT ONLY
                  DTA0:               DECTAPE 0             INPUT OR OUTPUT
                  DTA1:               DECTAPE 1             INPUT OR OUTPUT
                  DTA2:               DECTAPE 2             INPUT OR OUTPUT
                  DTA3:               DECTAPE 3             INPUT OR OUTPUT
                  DTA4:               DECTAPE 4             INPUT OR OUTPUT
                  DTA5:               DECTAPE 5             INPUT OR OUTPUT
                  DTA6:               DECTAPE 6             INPUT OR OUTPUT
                  DTA7:               DECTAPE 7             INPUT OR OUTPUT
                  MTA0:               MAGNETIC TAPE 0       INPUT OR OUTPUT
                  MTA1:               MAGNETIC TAPE 1       INPUT OR OUTPUT
 
INPUT MAY NOT BE DONE FROM THE LINE PRINTER NOR MAY OUTPUT GO TO THE CARD
READER
 
2)  FILE.EXT IS THE NAME AND EXTENSION OF THE FILE TO BE USED.  THIS PART OF
THE SPECIFICATION IS USED ONLY IF DISK OR DECTAPE IS USED.
 
3)  [PROJ,PROG]  IF A DISK IS USED AND THE USER WISHES TO READ A FILE IN
ANOTHER PERSON'S DIRECTORY, HE MAY DO SO BY SPECIFYING THE PROJECT-PROGRAMMER
NUMBER OF THE DIRECTORY FROM WHICH HE WISHES TO READ.  THE PROJECT NUMBER
AND THE PROGRAMMER NUMBER MUST BE SEPARATED BY A COMMA AND ENCLOSED IN BRACKETS.
OUTPUT MUST GO TO YOUR OWN AREA.
 
EXAMPLE:
          OUTPUT?  LPT:/2
          INPUT?   DSK:DATA.DAT[71171,71026]
 
IN THE EXAMPLE, TWO COPIES OF THE OUTPUT ARE TO BE PRINTED BY THE HIGH SPEED
LINE PRINTER.  THE INPUT DATA IS A DISK FILE OF NAME DATA.DAT IN USER DIRECTORY
[71171,71026].
 
DEFAULTS:
 
1)  IF NO DEVICE IS SPECIFIED BUT A FILENAME IS SPECIFIED THE DEFAULT DEVICE
    WILL BE DSK:
 
2)  IF NO FILENAME IS SPECIFIED AND A DISK OR DECTAPE IS USED THE DEFAULT ON
    INPUT WILL BE FROM INPUT.DAT; ON OUTPUT IT WILL BE OUTPT.DAT.
 
3)  IF THE PROGRAM IS RUN FROM THE TERMINAL AND NO SPECIFICATION IS GIVEN (JUST
    A CARRIAGE RETURN) BOTH INPUT AND OUTPUT DEVICES WILL BE THE TERMINAL.
 
4)  IF THE PROGRAM IS RUN THROUGH BATCH AND NO SPECIFICATION IS GIVEN, (A
    BLANK CARD) THE INPUT DEVICE WILL BE CDR: AND THE OUTPUT DEVICE WILL BE
    LPT:.
 
5)  IF NO PROJECT-PROGRAMMER NUMBER IS GIVEN, THE USER'S OWN NUMBER WILL BE
    ASSUMED.
 
       (1)  IF LPT: IS USED AS AN OUTPUT DEVICE MULTIPLE COPIES MAY BE OBTAINED
            BY SPECIFYING LPT:/N WHERE N REFERS TO THE NUMBER OF COPIES DESIRED.
 
       (2)  THE FOLLOWING TWO OPTIONS ARE NOT APPLICABLE FOR THE FIRST DATA SET,
            I.E., IT IS APPLICABLE ONLY WHEN THE PROGRAM BRANCHES BACK TO LINE
            2 UPON FIRST COMPLETION OF LINES 1-8.
 
            (A)  SAME OPTION
 
                 UPON RETURNING FROM LINE 8, IF THE SAME DATA FILE IS TO BE USED
                 AGAIN, SIMPLY ENTER "SAME<CR>", OTHERWISE, EITHER USE THE
                 FINISH OPTION OR ENTER ANOTHER FILENAME, ETC.
 
            (B)  FINISH OPTION
 
                 THE USER MUST ENTER "FINISH<CR>" TO BRANCH OUT OF THE PROGRAM.
                 FAILURE TO DO SO MIGHT RESULT IN LOSING THE ENTIRE OUTPUT FILE.
 
LINE 3  FORMAT:  (F-TYPE ONLY)
 
THERE ARE 3 OPTIONS AVAILABLE FOR THE FORMAT, NAMELY:
 
       (A)  STANDARD FORMAT OPTION
 
            UNLESS OTHERWISE SPECIFIED, THE PROGRAM ASSUMES THE STANDARD OPTION.
            IN THIS OPTION, THE DATA ARE ARRANGED IN GROUPS OF 10 PER LINE, TWO
            VALUES BEING SEPARATED BY A COMMA.
 
            TO USE THIS OPTION, SIMPLY TYPE IN "<CR>" ON TERMINAL JOBS OR USE A
            BLANK CARD FOR BATCH JOBS.
 
       (B)  OBJECT TIME FORMAT OPTION
 
            IF THE DATA IS SUCH THAT A USER'S OWN FORMAT IS REQUIRED, SIMPLY
            ENTER A LEFT PARENTHESIS FOLLOWED BY THE FIRST FORMAT SPECIFICATION,
            A COMMA AND THE SECOND SPECIFICATION, ETC.  WHEN YOU FINISH ENTER
            A RIGHT PARENTHESIS, AND THEN A CARRIAGE RETURN.  THERE CAN BE A
            MAXIMUM OF 2 LINES FOR THE FORMAT, EACH LINE BEING 80 COLUMNS LONG.
            NOTE THAT THE FORMAT SPECIFICATION LIST MUST USE THE FLOATING
            POINT (F-TYPE) NOTATION AND MUST CONTAIN SPECIFICATION FOR EACH OF
            THE VARIABLES.
 
            IF A BREAKDOWN VARIABLE IS USED BE SURE TO INCLUDE A SPECIFICATION
            FOR IT IN THE FORMAT.  THE SPECIFICATIONS FOR THE FORMAT ITSELF ARE
            THE SAME AS FOR THE FORTRAN IV FORMAT STATEMENT.
 
       (C)  SAME OPTION
 
            IF AN OBJECT TIME FORMAT WAS USED ON A DATA SET AND THE SUCCEEDING
            DATA SET UTILIZES THE SAME FORMAT, SIMPLY ENTER "SAME<CR>".
 
LINE 4  ENTER IDENTIFICATION
 
ENTER A LINE OF UP TO 50 CHARACTERS TO BE PRINTED ABOVE YOUR OUTPUT.  IF YOU
WANT NO HEADING ENTER A "<CR>" ON TERMINAL JOBS OR A BLANK CARD ON BATCH JOBS.
 
LINE 5  ENTER NO. OF LEVELS IN THE REPEATED FACTOR
 
ENTER A NUMBER BETWEEN 1 AND 30 TO INDICATE THE NUMBER OF LEVELS IN THE REPEATED
FACTOR.
 
LINE 6  WOULD YOU LIKE TO USE A BREAKDOWN VARIABLE?
 
THE LEVELS IN FACTOR A (THE NON-REPEATED FACTOR) MAY BE DEFINED IN TWO WAYS.
 
       (1)  THE LEVELS MAY BE DEFINED BY AN ADDITIONAL VARIABLE CALLED
            A BREAKDOWN VARIABLE.  THE LEVELS ARE DEFINED BY SPECIFYING THE
            RANGE OF VALUES THE BREAKDOWN VARIABLE MAY ASSUME FOR EACH LEVEL.
            FOR EXAMPLE:  IF THE VALUES ENTERED IN LINE 6.2 ARE 3, 7, 9 ALL
            THOSE SUBJECTS FOR WHICH THE BREAKDOWN VARIABLE IS 3 OR LESS WOULD
            BE IN GROUP 1, ALL THOSE WHOSE BREAKDOWN VARIABLE IS GREATER THAN
            3 AND LESS THAN OR EQUAL TO 7 WOULD BE IN GROUP 2 AND ALL THOSE
            WHOSE BREAKDOWN VARIABLE IS GREATER THAN 7 AND LESS THAN OR EQUAL
            TO 9 WOULD BE IN GROUP 3.
 
       (2)  THE LEVELS MAY ALSO BE DEFINED BY SPECIFYING THE SIZES OF THE
            GROUPS AND SUBMITTING THE DATA IN ORDER OF THE GROUP, I.E., SUBMIT
            ALL SUBJECTS IN GROUP 1 BEFORE SUBMITTING GROUP 2 SUBJECTS, ETC.
 
            A YES RESPONSE TO LINE 6 WILL YIELD THE FIRST OPTION.  A NO WILL
            YIELD THE SECOND OPTION.
 
            IF THE ANSWER IS "NO", PROCEED TO LINE 7, OTHERWISE PROCEED TO
            LINE 6.1.
 
LINE 6.1  ENTER BREAKDOWN VARIABLE NUMBER
 
INDICATE THE RELATIVE POSITION OF THE BREAKDOWN VARIABLE TO THE DATA VARIABLES.
IF IT APPEARS FIRST ENTER 1.  IF IT APPEARS SECOND ENTER 2, ETC.
 
LINE 6.2  ENTER UPPER BOUNDS ON THE BREAKDOWN VARIABLE
 
ENTER THE UPPER LIMITS ON THE BREAKDOWN VARIABLE FOR EACH GROUP.  FOR A MORE
COMPLETE DESCRIPTION SEE (1) UNDER THE DESCRIPTION FOR LINE 6.
 
PROCEED TO LINE 8.
 
LINE 7  ENTER THE NUMBER OF SUBJECTS IN EACH GROUP
 
THE NUMBER OF SUBJECTS IN EACH LEVEL OF FACTOR A MUST BE SUBMITTED HERE.  SEE
(2) IN THE DESCRIPTION FOR LINE 6.
 
LINE 8  (A)  ENTER DATA.
        (B)  YOUR DATA IS BEING READ.
 
THE PROGRAM WILL TYPE OUT MESSAGE A IF IT IS A TERMINAL JOB AND THE INPUT DATA
IS TO BE ENTERED DURING EXECUTION TIME (I.E. WHILE RUNNING THE PROGRAM).
OTHERWISE, MESSAGE B WILL BE TYPED OUT.
 
IMPORTANT:  YOU MUST ENTER THE DATA ACCORDING TO THE FORMAT YOU CHOSE.
 
AFTER ANALYZING THE DATA, THE PROGRAM WILL BRANCH TO LINE 2, THE USER CAN EITHER
PROCESS ANOTHER DATA SET, IN THIS CASE REPEAT LINES 2 TO 8; OR EXIT FROM THE
PROGRAM IN THIS CASE ENTER "FINISH<CR>".
 
3.0  METHOD OF USE
 
PROCEDURE AND SETUP FOR TERMINAL AND BATCH JOBS DIFFER SLIGHTLY.
 
3.1  TERMINAL JOB
 
FOLLOWING A SUCCESSFUL LOGIN, THE USER TYPES "R RMAOV<CR>".  INTERACTION
BETWEEN USER AND PROGRAM BEGINS AT THIS POINT.  THE PROGRAM WILL TYPE OUT
THE APPROPRIATE QUESTION OR STATEMENT AS OUTLINED IN SECTION 3.1.
 
EXAMPLE:
 
THE FOLLOWING EXAMPLE IS TAKEN FROM *WINER'S, "STATISTICAL PRINCIPLES IN
EXPERIMENTAL DESIGN"  PAGES 306-307.  EXCEPT FOR OUTPUT? AND INPUT?
ALL LINES TERMINATED BY THE <CR> ARE ENTERED BY THE USER.
TERMINAL IS USED FOR BOTH INPUT AND OUTPUT AND STANDARD
FORMAT OPTION IS SELECTED.
 
.R RMAOV<CR>
---WMU REPEATED MEASURE ANOVA
OUTPUT?<CR>
INPUT?<CR>
FORMAT:  (F-TYPE ONLY)
<CR>
ENTER IDENTIFICATION
WINER PAGE 306 (SAMPLE TERMINAL RUN)<CR>
ENTER NO OF LEVELS IN THE REPEATED FACTOR
4<CR>
WOULD YOU LIKE TO USE A BREAKDOWN VARIABLE?
NO<CR>
ENTER THE NUMBER OF SUBJECTS IN EACH GROUP
3,3<CR>
ENTER DATA
0,0,5,3<CR>
3,1,5,4<CR>
4,3,6,2<CR>
4,2,7,8<CR>
5,4,6,6<CR>
7,5,8,9<CR>
WINER PAGE 306 (SAMPLE TERMINAL RUN)
 



                                   MEANS
 FACTOR                         FACTOR B
   A  SIZE       1         2         3         4
   1     3      2.333     1.333     5.333     3.000
   2     3      5.333     3.667     7.000     7.667



                 2 - WAY    REPEATED  MEASURE  ANOVA

          SOURCE                     SS      DF     MS        F     PROB

      BETWEEN SUBJECTS  - - -      68.2083    5    13.6417
            A - - - - - - - -      51.0417    1    51.0417  11.893 0.026
   SUBJECTS WITHIN GROUPS - -      17.1667    4     4.2917
       WITHIN SUBJECTS  - - -      69.7500   18     3.8750
            B - - - - - - - -      47.4583    3    15.8194  12.798 0.000
           AB - - - - - - - -       7.4583    3     2.4861   2.011 0.166
 B X SUBJECTS WITHIN GROUPS -      14.8333   12     1.2361
          TOTAL - - - - - - -     137.9583   23
 
INPUT? FINISH<CR>

3.2  BATCH JOB
 
THE FOLLOWING IS A BATCH JOB SETUP:  (EACH LINE REPRESENTS ONE CARD, EACH CARD
STARTING IN COLUMN 1; DO NOT INCLUDE THE COMMENTS AT THE RIGHT).  IF INPUT DATA
ARE NOT IN CARDS, SKIP THE SECTION FROM $DATA TO $EOD INCLUSIVE.
 
--------------------------------------------------------------------------------
 
                                          COMMENTS
 
   $JOB [#,#]                             JOB CARD; INSERT USER'S PROJECT-
                                          PROGRAMMER NUMBER WITHIN THE BRACKET.
 
   $PASSWORD ######                       IN PLACE OF THE 6#'S, PUT IN THE
                                          PASSWORD.
 
   $DATA                                  SIGNIFY BEGINNING OF DATA DECK.
 
         (DATA CARDS)                     INSERT THE DATA CARD DECK TO BE
                                          ANALYZED.
 
   $EOD                                   SIGNIFY THE END OF DATA CARD DECK.
 
   .R RMAOV                               START THE EXECUTION.
 
        (RESPONSES TO LINES 1-8 IN
        SECTION 2.1 REPEATED OR NOT)
 
   (EOF)                                  AN END-OF-FILE CARD.
 
--------------------------------------------------------------------------------
 
EXAMPLE:  FOLLOWING IS A BATCH JOB SETUP USING THE SAME EXAMPLE IN SECTION 3.1.
THE INPUT DEVICE IS CDR:, WHILE OUTPUT DEVICE IS LPT:
 
                                          COMMENTS
 
   $JOB [420,420]                         JOB CARD
   $PASSWORD                              PASSWORD
   $DATA                                  START OF INPUT DATA
   0,0,5,3
   3,1,5,4
   4,3,6,2
   4,2,7,8                                INPUT DATA CARDS
   5,4,6,6
   7,5,8,9
   $EOD                                   END OF DATA DECK
   .R RMAOV                               START EXECUTION
   LPT:                                   OUTPUT DEVICE IS LINE PRINTER
   CDR:                                   INPUT DEVICE IS CARD READER
     (A BLANK CARD)                       STANDARD FORMAT
   SAMPLE BATCH RUN                       IDENTIFICATION
   4                                      NO. OF LEVELS IN THE REPEATED FACTOR
   3,3                                    NO. OF SUBJECTS IN EACH GROUP
   FINISH                                 TO BRANCH OUT OF THE PROGRAM
   (EOF)				  AN END OF FILE CARD