10 REM
11 PRINT"[147]"CHR$(142):POKE53281,0:POKE53280,6
12 PRINTSPC(8)"LOADSTAR 2X2 STATISTICS":PRINT""SPC(9)"WRITTEN BY AL VEKOVIUS
13 [153]" CONVERTED TO THE C-64 BY ALAN GARDNER"
14 [153]""[166]12)"COPYRIGHT 1985"
15 [153]"    LOADSTAR IS NOT PUBLIC DOMAIN."
16 [153]"  PLEASE RESPECT THE AUTHOR'S RIGHTS."
17 [153]"        PRESS A KEY TO CONTINUE.":[151]198,0:[146]198,1:[161]R$
20 [137]500:ND[178]0
100 [143]   PLOT ROUTINE
110 [151]783,0:[151]781,VT:[151]782,HT
120 [158]65520
130 [142]
500 [143]   MAIN MENU
510 [153]"LOAD":[151]53281,0:[151]53280,5
520 [153]"2X2 LOADSTAR STATISTICS"
530 [153]"1> INPUT DATA"
540 [153]"2> VIEW YOUR DATA"
550 [153]"3> VIEW EXPECTED VALUES"
560 [153]"4> CHI-SQUARE TEST"
570 [153]"5> FISHER'S EXACT TEST"
580 [153]"6> MCNEMAR'S TEST"
581 [153]"7> SEE EXAMPLES"
582 [153]"8> EXIT"
590 :
600 [153]"MAKE A SELECTION ";:[151]198,0:[146]198,1:[161]R$
605 [139]R$[179]"1"[176]R$[177]"8"[167][153]:[153]"ONONONON";:[151]53280,[187](1)[172]15:[137]600
606 [139]R$[178]"7"[167][141]8541:[137]700:[143]        SAVE THOSE VALUES
607 [139]R$[178]"8"[167]20000
610 [145][197](R$)[141]7000,4500,8200,6000,4000,9500
620 [137]500
700 [143]   EXAMPLE MENU
701 [153]"LOAD":[153]"EXAMPLES MENU"
702 [153]"1.  EXAMPLE OF CHI-SQUARE TEST"
704 [153]"2.  EXAMPLE OF FISHER'S TEST"
706 [153]"3.  EXAMPLE OF MCNEMAR'S TEST"
708 [153]"4.  RETURN TO MAIN MENU"
710 :
720 [153]"MAKE A SELECTION":[151]198,0:[146]198,1:[161]R$
721 [139]R$[179]"1"[176]R$[177]"4"[167][153]"ONONONON";:[151]53280,[187](1)[172]15:[137]720
730 [139]R$[178]"1"[167]FLAG[178]1:[141]1920:[137]12000
734 [139]R$[178]"2"[167]FLAG[178]2:[141]1910:[137]12000
736 [139]R$[178]"3"[167]FLAG[178]3:[141]1930:[137]12000
760 [141]8651:[141]7400:[141]8000:[141]8500
761 [143]  GET VALUES, MARGINS, DISPLAY VALUES, SAVE VALUES
770 [137]500
900 :
1000 [143]  COMPUTE FACTORIAL
1005 [139]N[178]1[176]N[178]0[167]S[178]1:[137]1050
1010 S[178]1
1020 [129]K[178]2[164]N
1030 S[178]K[172]S
1040 [130]
1050 [142]
1090 :
1900 [143]   EXAMPLE
1910 B$[178]"FISHER'S EXACT EXAMPLE":X$(1)[178]"ALIVE":X$(2)[178]"DIED"
1911 Y$(1)[178]"GAM GLOB":Y$(2)[178]"PASTEUR":N(1,1)[178]11:N(1,2)[178]1:N(2,1)[178]2:N(2,2)[178]3
1912 [142]
1913 :
1920 B$[178]"CHI-SQUARE EXAMPLE":X$(1)[178]"MALE":X$(2)[178]"FEMALE"
1921 Y$(1)[178]"DEMOCRAT":Y$(2)[178]"REPUBLICAN":N(1,1)[178]15:N(1,2)[178]10:N(2,1)[178]7:N(2,2)[178]10
1922 [142]
1923 :
1930 B$[178]"MCNEMAR'S TEST":X$(1)[178]"TEST2 +":X$(2)[178]"TEST2 -"
1931 Y$(1)[178]"TEST1 +":Y$(2)[178]"TEST1 -":N(1,1)[178]32:N(1,2)[178]10:N(2,1)[178]7:N(2,2)[178]26
1932 [142]
1933 :
1940 :
1999 [142]
2000 [143]   DRAW 2 X 2 TABLE
2005 [153]"LOAD"B$
2006 D$[178]""
2010 [129]I[178]6[164]19:[153][202](D$,1,I)[166]10)"PEEK";[166]10)"PEEK"[166]10)"PEEK"
2020 [130]
2030 [129]X[178]1[164]3:[153]"":[129]P[178]1[164]4[172]X[170]2:[153]"";:     [130]
2032 [153][166]9)"---------------------------"
2034 [130]
2040 [153]""Y$(1)
2050 [153]""Y$(2)
2060 [153]""X$(1)
2070 [153]""X$(2)
2080 [153]
2100 [142]
3000 VT[178]9:HT[178]14:[141]100
3002 [153]D(1,1);:HT[178]24:[141]100:[153]D(1,2);:HT[178]33:[141]100:[153]R(1)
3010 VT[178]13:HT[178]14:[141]100
3012 [153]D(2,1);:HT[178]24:[141]100:[153]D(2,2);:HT[178]33:[141]100:[153]R(2)
3020 VT[178]17:HT[178]14:[141]100
3022 [153]C(1);:HT[178]24:[141]100:[153]C(2);:HT[178]33:[141]100:[153]T
3090 [142]
3500 VT[178]9:HT[178]14:[141]100:[153]D(1,1):HT[178]24:[141]100:[153]""D(1,2)"WAIT";:HT[178]33:[141]100
3501 [153]R(1)
3510 VT[178]13:HT[178]14:[141]100:[153]""D(2,1)"WAIT":HT[178]24:[141]100:[153]D(2,2);:HT[178]33:[141]100
3511 [153]R(2)
3520 VT[178]17:HT[178]14:[141]100:[153]C(1);:HT[178]24:[141]100:[153]C(2):HT[178]33:[141]100:[153]T
3600 :
3999 [142]
4000 [143]  FISHER'S EXACT TEST
4001 [139]T[178]0[167][141]15000:[137]4110
4002 PROB[178]0:TAG[178]0
4005 [139]T[177]33[167][153]"LOAD CANNOT COMPUTE THIS FUNCTION WITH THE"
4006 [139]T[177]33[167][153]""[166]12)"CURRENT DATA.":[129]DL[178]1[164]3000:[130]:[137]4100
4010 [141]5100
4020 [141]5000
4030 PROB[178]PROB[170]X[173]P
4035 B$[178]"FISHER'S EXACT TEST"
4039 TAG[178]0
4040 [139]N(LI,LJ)[178]0[167][141]2000:[141]3000:X[178]PROB:TAG[178]1
4042 [139]TAG[178]1[167]VT[178]20:HT[178]0:[141]100:[153]" P=";:[141]11000:
4043 [139]TAG[178]1[167][153]"       PRESS ANY KEY TO CONTINUE":[151]198,0:[146]198,1
4044 [139]TAG[178]1[167][161]R$:[153]:[137]4100
4050 [141]5200
4060 [137]4020
4100 [141]8600
4110 [142]
4200 :
4500 [143] VIEW DATA
4501 [139]T[178]0[167][141]15000:[137]4540
4510 B$[178]"":[141]2000
4520 [141]3000
4530 [153]"       PRESS A KEY TO CONTINUE":[151]198,0:[146]198,1:[161]R$
4540 [142]
4600 :
5000 [143]  PRINT FISHER
5005 P[178]1
5010 [129]I[178]1[164]2:[129]J[178]1[164]2:N[178]N(I,J):[141]1000:P[178]P[172]S:[130]J,I
5020 [129]I[178]1[164]2
5030 N[178]C(I):[141]1000:CF(I)[178]S
5040 N[178]R(I):[141]1000:RF(I)[178]S
5050 [130]I
5060 N[178]T:[141]1000:TF[178]S
5070 X[178](CF(1)[172]CF(2)[173]TF)[172]RF(1)[172]RF(2)
5080 EXACT[178]X[173]P
5090 [142]
5095 :
5100 [143]  SET UP FISHER
5105 M[178]9999
5110 [139]M[177]R(1)[167] M[178]R(1):X[178]1
5120 [139]M[177]R(2)[167] M[178]R(2):X[178]2
5130 [139]M[177]C(1)[167] M[178]C(1):X[178]3
5140 [139]M[177]C(2)[167] M[178]C(2):X[178]4
5150 [139]X[179]3[167]A1[178]N(X,1)[173]C(1):A2[178]N(X,2)[173]C(2)
5160 [139]X[177]2[167]A1[178]N(1,X[171]2)[173]R(1):A2[178]N(2,X[171]2)[173]R(2)
5170 [139]X[179]3[167]LI[178]X:LJ[178]1:[139]A2[179]A1[167]LJ[178]2
5180 [139]X[177]2[167]LJ[178]X[171]2:LI[178]1:[139]A2[179]A1[167]LI[178]2
5190 [142]
5195 :
5200 [143]  SET UP NEXT PASS FISHER
5210 [129]I[178]1[164]2:[129]J[178]1[164]2
5230 [139]I[178]LI[175]J[178]LJ[167]N(I,J)[178]N(I,J)[171]1:[137]5290
5240 [139]I[179][177]LI[175]J[179][177]LJ[167]N(I,J)[178]N(I,J)[171]1:[137]5290
5250 N(I,J)[178]N(I,J)[170]1
5290 [130]J:[130]I
5300 [142]
5400 :
6000 [143] CHI SQUARE AND YATES
6005 [139]T[178]0[167][141]15000:[137]6900
6006 B$[178]"CHI-SQUARE TEST"
6010 CN[178]N(1,1)[172]N(2,2)[171]N(2,1)[172]N(1,2)
6020 YT[178][182](CN)[171]T[173]2
6030 PROD[178]C(1)[172]C(2)[172]R(1)[172]R(2)
6040 CSQ[178]CN[174]2[172]T[173]PROD
6050 YTS[178]YT[174]2[172]T[173]PROD
6060 [141]8000:[141]2000:[141]3000
6100 [153]"CHI-SQUARE :"CSQ
6110 A[178]CSQ:[141]10000
6120 [153]A$
6125 [153]
6130 [153]"YATES ="YTS
6140 A[178]YTS:[141]10000
6150 [153]A$
6160 [151]198,0:[146]198,1:[161]R$
6900 [153]:[142]
6999 :
7000 [143]   INPUT DATA
7010 X$(1)[178]" COL 1":X$(2)[178]" COL 2":          Y$(1)[178]" ROW 1":Y$(2)[178]" ROW 2"
7020 [153]"LOAD"
7025 [141]2000
7030 HT[178]10:VT[178]20:[141]100:[153]"DO YOU WANT NEW LABELS"
7031 VT[178]21:HT[178]10:[141]100:[153]"FOR THIS TABLE (Y/N)";:
7032 [151]198,0:[146]198,1:[161]R$:[153]
7040 [139]R$[178]"Y"[167][141]7300:[137]7020
7045 [139]R$[178]"N"[167]7200
7050 [153]"ONONONON";:[137]7032
7100 :
7200 [143]  GET DATA VALUES
7201 ND[178]1:VT[178]1:HT[178]1:[141]100:[153]"PUT YOUR DATA IN THE APPROPRIATE CELL"
7220 VT[178]9:HT[178]14:[141]100:[133]N$:N(1,1)[178][197](N$)
7230 VT[178]9:HT[178]24:[141]100:[133]N$:N(1,2)[178][197](N$)
7240 VT[178]13:HT[178]14:[141]100:[133]N$:N(2,1)[178][197](N$)
7250 VT[178]13:HT[178]24:[141]100:[133]N$:N(2,2)[178][197](N$)
7260 VT[178]18:HT[178]10:[141]100:[153]"SATISFIED (Y/N)";:[151]198,0:[146]198,1:[161]R$
7270 [139]N(1,1)[177]16[176]N(1,2)[177]16[176]N(2,1)[177]16[176]N(2,2)[177]16[167][141]16000:[137]7000
7272 [139]N(1,1)[179]0[176]N(1,2)[179]0[176]N(2,1)[179]0[176]N(2,2)[179]0[167][141]17000:[137]7000
7275 [141]7400:[141]8000:[141]8100
7288 [141]8500
7290 [142]
7295 :
7300 [143]  INPUT LABELS
7310 [153]"LOAD"
7320 [133] "LABEL FOR ROW 1 ";Y$(1):[139][195](Y$(1))[177]8[167]Y$(1)[178][200](Y$(1),8)
7330 [133] "LABEL FOR ROW 2 ";Y$(2):[139][195](Y$(2))[177]8[167]Y$(2)[178][200](Y$(2),8)
7340 [133] "LABEL FOR COL 1 ";X$(1):[139][195](X$(1))[177]8[167]X$(1)[178][200](X$(1),8)
7350 [133] "LABEL FOR COL 2 ";X$(2):[139][195](X$(2))[177]8[167]X$(2)[178][200](X$(2),8)
7360 [142]
7390 :
7400 [143]  CALCULATE MARGINS
7405 R(1)[178]N(1,1)[170]N(1,2):R(2)[178]N(2,1)[170]N(2,2)
7410 C(1)[178]N(1,1)[170]N(2,1):C(2)[178]N(1,2)[170]N(2,2)
7420 T[178]C(1)[170]C(2)
7499 [142]
7500 :
8000 [143]  GET DISPLAY VALUES
8010 [129]I[178]1[164]2:[129]J[178]1[164]2:D(I,J)[178]N(I,J):[130]J,I
8020 [142]
8030 :
8100 [143]   EXPECTED VALUES
8105 EX[178]0
8110 B$[178]"EXPECTED VALUES"
8120 [129]I[178]1[164]2:[129]J[178]1[164]2
8130 E(I,J)[178]R(I)[172]C(J)[173]T:E(I,J)[178][181](E(I,J)[172]100[170].5)[173]100
8135 [139]E(I,J)[179]5[167]EX[178]1
8140 [130]J,I
8150 [142]
8160 :
8200 [143]  PRINT EXPECTED VALUES
8205 [139]T[178]0[167][141]15000:[137]8280
8230 [129]I[178]1[164]2:[129]J[178]1[164]2:D(I,J)[178]E(I,J):[130]J,I
8240 [141]2000
8245 [141]3000
8250 [139]EX[178]1[167]VT[178]20:HT[178]0:[141]100
8251 [139]EX[178]1[167][153]"YOU HAVE AN EXPECTED VALUES LESS THAN 5"
8260 [153]"       PRESS A KEY TO CONTINUE":[151]198,0:[146]198,1:[161]R$:[153]
8270 [141]8000
8280 [142]
8300 :
8500 [143]   HOLD VALUES
8510 [129]I[178]1[164]2:[129]J[178]1[164]2
8520 T(I,J)[178]N(I,J)
8530 [130]J,I
8540 [142]
8541 [143]   HOLD VALUES FOR EXAMPLES
8542 [129]I[178]1[164]2:[129]J[178]1[164]2
8543 L(I,J)[178]N(I,J)
8544 [130]J,I
8545 L1$[178]Y$(1):L2$[178]Y$(2):L3$[178]X$(1):L4$[178]X$(2)
8546 [142]
8550 :
8600 [143]   RETRIEVE VALUES
8610 [129]I[178]1[164]2:[129]J[178]1[164]2
8620 N(I,J)[178]T(I,J)
8630 [130]J,I
8640 [142]
8650 :
8651 [143]   RETRIEVE VALUES FOR EXAMPLES
8652 [129]I[178]1[164]2:[129]J[178]1[164]2
8653 N(I,J)[178]L(I,J)
8654 [130]J,I
8655 Y$(1)[178]L1$:Y$(2)[178]L2$:X$(1)[178]L3$:X$(2)[178]L4$
8656 [142]
8700 :
9000 [143]  CORR BINOMIAL PROB
9010 M1[178]N(1,2):M2[178]N(2,1):M3[178]M1[170]M2
9020 LT[178]M1
9030 [139]M2[179]M1[167]LT[178]M2
9040 BI[178]1
9045 [139]LT[178]0[167]9110
9050 N[178]M3:[141]1000:M4[178]S
9060 [129]J[178]1[164]LT
9070 N[178]J:[141]1000:JF[178]S
9080 N[178]M3[171]J:[141]1000:JS[178]S
9090 BI[178]BI[170]M4[173](JF[172]JS)
9100 [130]J
9110 BI[178]BI[173]2[174]M3
9120 [142]
9130 :
9500 [143]  MCNEMAR'S TEST
9505 [139]T[178]0[167][141]15000:[137]9550
9510 MC[178](N(1,2)[171]N(2,1))[174]2[173](N(1,2)[170]N(2,1))
9520 MC[178][181](MC[172]10000[170].5)[173]10000:A[178]MC:[141]10000
9525 B$[178]"MCNEMAR'S TEST"
9530 [141]2000:[141]3000:VT[178]20:HT[178]0:[141]100:[153]"CHI SQUARE ="MC:[153]A$
9535 [141]9000:[153]"ONE TAIL BINOMIAL PROB. =";:X[178]BI:[141]11000
9540 [151]198,0:[146]198,1:[161]R$:[153]
9550 [142]
9600 :
10000 [143]  CHI SQUARE TABLE DF=1
10010 A$[178]"NOT SIGNIFICANT AT .05 LEVEL"
10020 [139]A[177]3.84[167]A$[178]"SIGNIFICANT AT THE .05 LEVEL"
10030 [139]A[177]5.02[167]A$[178]"SIGNIFICANT AT THE .025 LEVEL"
10040 [139]A[177]6.63[167]A$[178]"SIGNIFICANT AT THE .010 LEVEL"
10050 [139]A[177]7.88[167]A$[178]"SIGNIFICANT AT THE .005 LEVEL"
10060 [139]A[177]10.8[167]A$[178]"SIGNIFICANT AT THE .0001 LEVEL"
10070 [142]
10080 :
11000 [143]   PRINT ROUTINE
11005 [153]".";
11010 [129]I[178]1[164]6
11020 X[178]10[172]X:XX[178][181](X):S$[178][196](XX):S$[178][202](S$,2):X[178]X[171]XX
11030 [153]S$;
11040 [130]
11050 [153]
11060 [142]
11070 :
12000 [143]  EXAMPLE OF FISHER
12003 [141]8500: [143]   SAVE VALUES
12005 [153]"LOAD":[139]FLAG[179][177]2[167]12020
12010 [153][163]4)"AN EXAMPLE OF FISHER'S EXACT TEST"
12012 [153]"  SEVENTEEN PERSONS HAD BEEN BITTEN IN  THE HEAD OR NECK BY THE SAME";
12013 [153]" RABID WOLF IN IRAN (BULL. WORLD HEALTH ORGAN.1955).
12014 PRINT"  THE STANDARD PASTEUR VACCINE TREATMENTWAS GIVEN.  IN ADDITION, 12";
12015 PRINT" PERSONS ALSORECEIVED ONE OR MORE DOSES OF ANTIRABIES";
12016 PRINT"GAMMA GLOBULIN.  THE RESULTS ARE SUMMAR-IZED AS FOLLOWS:"
12018 VT=23:HT=5:GOSUB100:PRINT"PRESS ANY KEY TO CONTINUE";:POKE198,0:WAIT198,1
12019 GETR$:GOTO12050
12020 IFFLAG<>1THEN12031
12021 PRINTTAB(5)"AN EXAMPLE OF CHI-SQUARE TEST":PRINT"  SUPPOSE THAT 42 PEOPLE ARE";
12022 PRINT" CALLED AT  RANDOM TO DETERMINE WHETHER THERE IS ANYRELATIONSHIP";
12023 PRINT" BETWEEN A PERSON'S SEX AND WHETHER HE WILL VOTE DEMOCRATIC OR  "
12024 PRINT"REPUBLICAN IN THE NEXT ELECTON.":PRINT"RESULTS ARE AS FOLLOWS:  OF THE 25"
12025 PRINT"MEN CONTACTED, 15 WILL VOTE DEMOCRATIC  AND 10 WILL VOTE REPUBLICAN."
12026 PRINT"OF THE 17 WOMEN CONTACTED, 7 WILL VOTE  DEMOCRATIC[160]AND 10 WILL VOTE";
12027 PRINT" REPUBLICAN.":PRINT"  THE QUESTION IS WHETHER THERE IS A    DIFFERENCE";
12028 PRINT" BETWEEN MALES AND FEMALES    BASED ON THE FINDINGS REPRESENTED."
12029 VT=24:HT=5:GOSUB100:PRINT"PRESS ANY KEY TO CONTINUE";:POKE198,0:WAIT198,1
12030 GETR$:GOTO12050
12031 PRINTTAB(6)"AN EXAMPLE OF MCNEMAR'S TEST":PRINT"  SUPPOSE THAT TWO DIAGNOSTIC";
12032 PRINT" TESTS ARE BOTH APPLIED TO 75 MEN.":PRINT"WHEN THE TESTS DO NOT AGREE,";
12033 PRINT" DO THEY DO SO IN THE SAME WAY?  THAT IS, DO THE"
12034 PRINT"CASES WHEN TEST 1 IS POSITIVE AND TEST 2"
12035 PRINT"[145]IS NEGATIVE OCCUR IN EQUAL NUMBERS TO   THE CASES WHEN TEST 1 IS";
12036 PRINT" NEGATIVE AND   TEST 2 IS POSITIVE?"
12038 VT=23:HT=5:GOSUB100:PRINT"PRESS ANY KEY TO CONTINUE":POKE198,0:WAIT198,1
12039 GETR$
12050 GOSUB8000:GOSUB7400:GOSUB2000
12052 IFFLAG=3THENGOSUB3500:GOTO12070
12055 GOSUB3000
12070 VT=23:HT=1:GOSUB100:PRINT"PRESS A KEY TO EXAMINE EXPECTED VALUES"
12072 POKE198,0:WAIT198,1:GETR$:GOSUB8105
12080 FORI=1TO2:FORJ=1TO2:D(I,J)=E(I,J):NEXTJ,I
12090 GOSUB2000:GOSUB3000
12100 IFFLAG=2THENPRINT:PRINT"SINCE AT LEAST ONE EXPECTED VALUE IS    LESS THAN 5, ";
12101 IFFLAG=2THENPRINT" CHI-SQUARE IS NOT VALID.  FISHER'S EXACT TEST IS THE";
12102 IFFLAG=2THENPRINT" TEST TO USE."
12110 VT=23:HT=5:GOSUB100:PRINT"PRESS ANY KEY TO CONTINUE":POKE198,0:WAIT198,1
12111 GETR$:
12120 IFFLAG=2THENGOSUB8000:GOSUB4000:GOTO12200
12130 IFFLAG=1THENGOSUB6000:GOTO12200
12140 GOSUB9500
12200 GOTO700
15000 REM   NO DATA ENTERED YET!
15005 PRINT"[147] SORRY! THERE IS NO DATA ENTERED YET.."
15010 PRINT"PLEASE PRESS A KEY"
15020 POKE198,0:WAIT198,1:GETR$
15030 RETURN
16000 REM   VALUES TO BIG
16010 PRINT"[147]       YOUR VALUES ARE TOO BIG..."
16012 PRINT"     PLEASE KEEP THEM LESS THAN 16."
16020 PRINT"  PRESS A KEY TO RE-ENTER THE VALUES."
16030 POKE198,0:WAIT198,1:GETR$
16040 RETURN
17000 REM   VALUES TO SMALL
17020 PRINT"[147]      NO NEGATIVE VALUES PLEASE..."
17030 PRINT"  PRESS A KEY TO RE-ENTER THE VALUES."
17040 POKE198,0:WAIT198,1:GETR$
17050 RETURN
20000 REM   RETURN TO LS
20010 PRINT"[147]LOAD"CHR$(34)"HELLO"CHR$(34)",8":PRINT"RUN"
20020 POKE631,13:POKE632,13:POKE198,2:
20030 END