Source Code 1992 March

home *** CD-ROM | disk | FTP | other *** search
/ Source Code 1992 March / Source_Code_CD-ROM_Walnut_Creek_March_1992.iso / unix_c / utils / tot_info.tar < prev next >
Tape Archive | 1989-07-14 | 56.0 KB
open in:
MacOS 8.1 extracted |
Win98 extracted |
DOS extracted
view JSON data |
view as text

This file was processed as: Tape Archive (archive/tar).

You can browse this item here: tot_info.tar
Confidence	Program	Detection	Match Type	Support
`100%`	dexvert	Tape Archive `(archive/tar)`	magic	Supported
`100%`	file	tar archive	default
`99%`	file	tar archive (V7), file tot-info/, mode 755 , uid 243 , gid 0 , size 0 , seconds 4457210476	default
`98%`	file	data	default
`100%`	siegfried	x-fmt/265 Tape Archive Format	default
`100%`	lsar	Tar	default
`100%`	gt2	TAR Archiv gefunden (Auflistung ist deaktiviert)	default
`100%`	disktype	Pre-POSIX tar archive	default
hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 74 6f 74 2d 69 6e 66 6f | 2f 00 00 00 00 00 00 00 |tot-info|/.......|
|00000010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000050| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000060| 00 00 00 00 20 20 20 37 | 35 35 20 00 20 20 20 32 |....   7|55 .   2|
|00000070| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00000080| 20 20 20 20 20 20 30 20 | 20 34 34 35 37 32 31 30 |      0 | 4457210|
|00000090| 34 37 36 20 20 20 35 35 | 34 31 00 20 00 00 00 00 |476   55|41. ....|
|000000a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000100| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000110| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000120| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000130| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000140| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000150| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000160| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000170| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000180| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000190| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000001f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000200| 74 6f 74 2d 69 6e 66 6f | 2f 53 68 61 72 46 69 6c |tot-info|/SharFil|
|00000210| 65 2e 48 65 61 64 65 72 | 73 00 00 00 00 00 00 00 |e.Header|s.......|
|00000220| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000230| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000240| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000250| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000260| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00000270| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00000280| 20 20 20 20 33 32 30 20 | 20 34 34 35 37 32 31 30 |    320 | 4457210|
|00000290| 34 37 37 20 20 31 30 35 | 37 34 00 20 00 00 00 00 |477  105|74. ....|
|000002a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000002b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000002c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000002d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000002e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000002f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000300| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000310| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000320| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000330| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000340| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000350| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000360| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000370| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000380| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000390| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000003f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000400| 3d 3d 3d 20 74 6f 74 2d | 69 6e 66 6f 2e 73 68 61 |=== tot-|info.sha|
|00000410| 72 20 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |r ======|========|
|00000420| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 0a 46 72 6f 6d |========|===.From|
|00000430| 20 72 69 61 63 73 21 65 | 6f 73 21 61 6d 65 73 21 | riacs!e|os!ames!|
|00000440| 78 61 6e 74 68 21 67 69 | 6e 6f 73 6b 6f 21 75 75 |xanth!gi|nosko!uu|
|00000450| 6e 65 74 21 61 6c 6c 62 | 65 72 79 20 4d 6f 6e 20 |net!allb|ery Mon |
|00000460| 4a 75 6c 20 31 30 20 31 | 31 3a 32 34 3a 34 37 20 |Jul 10 1|1:24:47 |
|00000470| 50 44 54 20 31 39 38 39 | 0a 0a 50 6f 73 74 69 6e |PDT 1989|..Postin|
|00000480| 67 2d 6e 75 6d 62 65 72 | 3a 20 56 6f 6c 75 6d 65 |g-number|: Volume|
|00000490| 20 37 2c 20 49 73 73 75 | 65 20 36 33 0a 53 75 62 | 7, Issu|e 63.Sub|
|000004a0| 6d 69 74 74 65 64 2d 62 | 79 3a 20 67 77 79 6e 40 |mitted-b|y: gwyn@|
|000004b0| 42 52 4c 2e 4d 49 4c 0a | 41 72 63 68 69 76 65 2d |BRL.MIL.|Archive-|
|000004c0| 6e 61 6d 65 3a 20 74 6f | 74 5f 69 6e 66 6f 0a 0a |name: to|t_info..|
|000004d0| 61 2e 63 00 00 00 00 bc | 00 00 78 f3 00 14 00 07 |a.c.....|..x.....|
|000004e0| 67 61 6d 6d 61 2e 68 00 | 00 00 00 cc 00 00 78 f4 |gamma.h.|......x.|
|000004f0| 00 14 00 06 69 6e 66 6f | 2e 63 00 ee 00 00 00 dc |....info|.c......|
|00000500| 00 00 78 f5 00 14 00 05 | 73 74 64 2e 68 00 10 c7 |..x.....|std.h...|
|00000510| 00 00 00 f0 00 00 78 f6 | 00 18 00 0a 74 6f 74 5f |......x.|....tot_|
|00000520| 69 6e 66 6f 2e 31 00 07 | 00 00 01 04 00 00 78 f7 |info.1..|......x.|
|00000530| 00 18 00 0a 74 6f 74 5f | 69 6e 66 6f 2e 63 00 00 |....tot_|info.c..|
|00000540| 00 00 01 1c 00 00 78 f8 | 00 1c 00 0c 74 6f 74 5f |......x.|....tot_|
|00000550| 69 6e 66 6f 2e 65 78 70 | 00 14 00 07 00 00 02 00 |info.exp|........|
|00000560| 00 00 78 f9 00 18 00 0b | 74 6f 74 5f 69 6e 66 6f |..x.....|tot_info|
|00000570| 2e 69 6e 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.in.....|........|
|00000580| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000590| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000005f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000600| 74 6f 74 2d 69 6e 66 6f | 2f 4d 61 6b 65 66 69 6c |tot-info|/Makefil|
|00000610| 65 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |e.......|........|
|00000620| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000630| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000640| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000650| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000660| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00000670| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00000680| 20 20 20 33 35 31 31 20 | 20 34 34 35 37 32 31 30 |   3511 | 4457210|
|00000690| 34 36 31 20 20 20 37 32 | 36 30 00 20 00 00 00 00 |461   72|60. ....|
|000006a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000700| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000710| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000720| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000730| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000740| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000750| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000760| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000770| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000780| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000790| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000007f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000800| 23 09 4d 61 6b 65 66 69 | 6c 65 20 2d 2d 20 2d 2d |#.Makefi|le -- --|
|00000810| 20 6d 61 6b 65 66 69 6c | 65 20 66 6f 72 20 22 74 | makefil|e for "t|
|00000820| 6f 74 5f 69 6e 66 6f 22 | 20 75 74 69 6c 69 74 79 |ot_info"| utility|
|00000830| 0a 0a 23 09 6c 61 73 74 | 20 65 64 69 74 3a 09 38 |..#.last| edit:.8|
|00000840| 39 2f 30 32 2f 30 36 09 | 44 20 41 20 47 77 79 6e |9/02/06.|D A Gwyn|
|00000850| 0a 0a 23 09 53 43 43 53 | 20 49 44 3a 09 40 28 23 |..#.SCCS| ID:.@(#|
|00000860| 29 74 6f 74 5f 69 6e 66 | 6f 2e 6d 6b 09 31 2e 31 |)tot_inf|o.mk.1.1|
|00000870| 20 28 65 64 69 74 65 64 | 20 66 6f 72 20 70 75 62 | (edited| for pub|
|00000880| 6c 69 63 61 74 69 6f 6e | 29 0a 0a 50 52 4f 44 55 |lication|)..PRODU|
|00000890| 43 54 20 3d 20 74 6f 74 | 5f 69 6e 66 6f 0a 0a 4d |CT = tot|_info..M|
|000008a0| 41 4b 45 46 49 4c 09 3d | 20 4d 61 6b 65 66 69 6c |AKEFIL.=| Makefil|
|000008b0| 65 0a 48 45 41 44 45 52 | 53 20 3d 20 63 68 69 73 |e.HEADER|S = chis|
|000008c0| 71 2e 68 20 67 61 6d 6d | 61 2e 68 0a 43 46 49 4c |q.h gamm|a.h.CFIL|
|000008d0| 45 53 09 3d 20 24 28 50 | 52 4f 44 55 43 54 29 2e |ES.= $(P|RODUCT).|
|000008e0| 63 20 69 6e 66 6f 2e 63 | 20 67 61 6d 6d 61 2e 63 |c info.c| gamma.c|
|000008f0| 20 63 68 69 73 71 2e 63 | 20 23 20 63 68 69 73 71 | chisq.c| # chisq|
|00000900| 2e 63 20 69 73 20 61 20 | 22 66 72 65 65 62 69 65 |.c is a |"freebie|
|00000910| 22 0a 4f 42 4a 45 43 54 | 53 09 3d 20 24 28 50 52 |".OBJECT|S.= $(PR|
|00000920| 4f 44 55 43 54 29 2e 6f | 20 69 6e 66 6f 2e 6f 20 |ODUCT).o| info.o |
|00000930| 67 61 6d 6d 61 2e 6f 20 | 23 20 63 68 69 73 71 2e |gamma.o |# chisq.|
|00000940| 6f 0a 4c 49 42 54 45 53 | 54 09 3d 20 67 5f 74 65 |o.LIBTES|T.= g_te|
|00000950| 73 74 0a 54 53 52 43 09 | 3d 20 24 28 4c 49 42 54 |st.TSRC.|= $(LIBT|
|00000960| 45 53 54 29 2e 63 0a 54 | 4f 42 4a 09 3d 20 24 28 |EST).c.T|OBJ.= $(|
|00000970| 4c 49 42 54 45 53 54 29 | 2e 6f 0a 54 46 49 4c 45 |LIBTEST)|.o.TFILE|
|00000980| 53 09 3d 20 24 28 50 52 | 4f 44 55 43 54 29 2e 69 |S.= $(PR|ODUCT).i|
|00000990| 6e 20 24 28 50 52 4f 44 | 55 43 54 29 2e 65 78 70 |n $(PROD|UCT).exp|
|000009a0| 0a 4c 49 42 45 53 09 3d | 20 2d 6c 6d 0a 4c 4c 49 |.LIBES.=| -lm.LLI|
|000009b0| 42 45 53 09 3d 20 24 28 | 4c 49 42 45 53 29 0a 4d |BES.= $(|LIBES).M|
|000009c0| 41 4e 53 52 43 53 09 3d | 20 24 28 50 52 4f 44 55 |ANSRCS.=| $(PRODU|
|000009d0| 43 54 29 2e 31 20 63 68 | 69 73 71 2e 33 20 67 61 |CT).1 ch|isq.3 ga|
|000009e0| 6d 6d 61 2e 33 0a 0a 23 | 20 22 73 74 61 6e 64 61 |mma.3..#| "standa|
|000009f0| 72 64 22 20 55 4e 49 58 | 20 75 74 69 6c 69 74 69 |rd" UNIX| utiliti|
|00000a00| 65 73 20 74 68 61 74 20 | 6d 61 79 20 6e 65 65 64 |es that |may need|
|00000a10| 20 69 6e 76 6f 63 61 74 | 69 6f 6e 73 20 74 77 65 | invocat|ions twe|
|00000a20| 61 6b 65 64 3a 0a 49 4e | 53 09 3d 20 63 70 0a 4c |aked:.IN|S.= cp.L|
|00000a30| 49 4e 54 09 3d 20 6c 69 | 6e 74 0a 46 4c 4f 57 09 |INT.= li|nt.FLOW.|
|00000a40| 3d 20 63 66 6c 6f 77 0a | 58 52 45 46 09 3d 20 63 |= cflow.|XREF.= c|
|00000a50| 78 72 65 66 20 2d 63 20 | 2d 73 20 2d 77 31 33 32 |xref -c |-s -w132|
|00000a60| 0a 0a 23 20 70 72 69 6e | 74 65 72 20 63 6f 6e 66 |..# prin|ter conf|
|00000a70| 69 67 75 72 61 74 69 6f | 6e 3a 0a 4c 50 52 09 3d |iguratio|n:.LPR.=|
|00000a80| 20 6c 70 0a 50 52 09 3d | 20 70 72 20 2d 66 0a 0a | lp.PR.=| pr -f..|
|00000a90| 23 20 44 57 42 20 63 6f | 6e 66 69 67 75 72 61 74 |# DWB co|nfigurat|
|00000aa0| 69 6f 6e 3a 0a 54 4f 50 | 54 09 3d 20 2d 54 69 33 |ion:.TOP|T.= -Ti3|
|00000ab0| 30 30 0a 54 4d 41 43 09 | 3d 20 2d 6d 61 6e 0a 45 |00.TMAC.|= -man.E|
|00000ac0| 51 4e 09 3d 20 65 71 6e | 20 24 28 54 4f 50 54 29 |QN.= eqn| $(TOPT)|
|00000ad0| 0a 54 52 4f 46 46 09 3d | 20 74 72 6f 66 66 20 24 |.TROFF.=| troff $|
|00000ae0| 28 54 4f 50 54 29 20 24 | 28 54 4d 41 43 29 0a 54 |(TOPT) $|(TMAC).T|
|00000af0| 50 4f 53 54 09 3d 20 7c | 20 64 69 6d 70 0a 0a 0a |POST.= || dimp...|
|00000b00| 24 28 50 52 4f 44 55 43 | 54 29 3a 09 24 28 4f 42 |$(PRODUC|T):.$(OB|
|00000b10| 4a 45 43 54 53 29 0a 09 | 24 28 43 43 29 20 24 28 |JECTS)..|$(CC) $(|
|00000b20| 4c 44 46 4c 41 47 53 29 | 20 2d 6f 20 24 40 20 24 |LDFLAGS)| -o $@ $|
|00000b30| 28 4f 42 4a 45 43 54 53 | 29 20 24 28 4c 49 42 45 |(OBJECTS|) $(LIBE|
|00000b40| 53 29 0a 0a 63 68 69 73 | 71 2e 6f 20 69 6e 66 6f |S)..chis|q.o info|
|00000b50| 2e 6f 3a 09 63 68 69 73 | 71 2e 68 0a 0a 67 61 6d |.o:.chis|q.h..gam|
|00000b60| 6d 61 2e 6f 3a 09 67 61 | 6d 6d 61 2e 68 0a 0a 70 |ma.o:.ga|mma.h..p|
|00000b70| 72 69 6e 74 3a 09 09 24 | 28 4d 41 4b 45 46 49 4c |rint:..$|(MAKEFIL|
|00000b80| 29 20 24 28 48 45 41 44 | 45 52 53 29 20 24 28 43 |) $(HEAD|ERS) $(C|
|00000b90| 46 49 4c 45 53 29 20 24 | 28 54 46 49 4c 45 53 29 |FILES) $|(TFILES)|
|00000ba0| 20 24 28 4d 41 4e 53 52 | 43 53 29 0a 09 24 28 50 | $(MANSR|CS)..$(P|
|00000bb0| 52 29 20 24 28 4d 41 4b | 45 46 49 4c 29 20 24 28 |R) $(MAK|EFIL) $(|
|00000bc0| 48 45 41 44 45 52 53 29 | 20 24 28 43 46 49 4c 45 |HEADERS)| $(CFILE|
|00000bd0| 53 29 20 24 28 54 46 49 | 4c 45 53 29 20 24 28 4d |S) $(TFI|LES) $(M|
|00000be0| 41 4e 53 52 43 53 29 20 | 7c 20 24 28 4c 50 52 29 |ANSRCS) || $(LPR)|
|00000bf0| 0a 0a 74 79 70 65 73 65 | 74 3a 09 24 28 4d 41 4e |..typese|t:.$(MAN|
|00000c00| 53 52 43 53 29 0a 09 24 | 28 45 51 4e 29 20 24 28 |SRCS)..$|(EQN) $(|
|00000c10| 50 52 4f 44 55 43 54 29 | 2e 31 20 7c 20 24 28 54 |PRODUCT)|.1 | $(T|
|00000c20| 52 4f 46 46 29 20 24 28 | 54 50 4f 53 54 29 0a 09 |ROFF) $(|TPOST)..|
|00000c30| 24 28 45 51 4e 29 20 63 | 68 69 73 71 2e 33 20 7c |$(EQN) c|hisq.3 ||
|00000c40| 20 24 28 54 52 4f 46 46 | 29 20 24 28 54 50 4f 53 | $(TROFF|) $(TPOS|
|00000c50| 54 29 0a 09 24 28 45 51 | 4e 29 20 67 61 6d 6d 61 |T)..$(EQ|N) gamma|
|00000c60| 2e 33 20 7c 20 24 28 54 | 52 4f 46 46 29 20 24 28 |.3 | $(T|ROFF) $(|
|00000c70| 54 50 4f 53 54 29 0a 0a | 6c 69 6e 74 3a 09 09 24 |TPOST)..|lint:..$|
|00000c80| 28 48 45 41 44 45 52 53 | 29 20 24 28 43 46 49 4c |(HEADERS|) $(CFIL|
|00000c90| 45 53 29 0a 09 24 28 4c | 49 4e 54 29 20 24 28 43 |ES)..$(L|INT) $(C|
|00000ca0| 46 49 4c 45 53 29 20 24 | 28 4c 4c 49 42 45 53 29 |FILES) $|(LLIBES)|
|00000cb0| 20 3e 20 24 40 0a 0a 66 | 6c 6f 77 3a 09 09 24 28 | > $@..f|low:..$(|
|00000cc0| 48 45 41 44 45 52 53 29 | 20 24 28 43 46 49 4c 45 |HEADERS)| $(CFILE|
|00000cd0| 53 29 0a 09 24 28 46 4c | 4f 57 29 20 24 28 43 46 |S)..$(FL|OW) $(CF|
|00000ce0| 49 4c 45 53 29 20 3e 20 | 24 40 0a 0a 78 72 65 66 |ILES) > |$@..xref|
|00000cf0| 3a 09 09 24 28 48 45 41 | 44 45 52 53 29 20 24 28 |:..$(HEA|DERS) $(|
|00000d00| 43 46 49 4c 45 53 29 0a | 09 24 28 58 52 45 46 29 |CFILES).|.$(XREF)|
|00000d10| 20 24 28 43 46 49 4c 45 | 53 29 20 3e 20 24 40 0a | $(CFILE|S) > $@.|
|00000d20| 0a 74 65 73 74 3a 09 09 | 67 61 6d 6d 61 5f 74 65 |.test:..|gamma_te|
|00000d30| 73 74 20 74 6f 74 5f 74 | 65 73 74 0a 0a 74 6f 74 |st tot_t|est..tot|
|00000d40| 5f 74 65 73 74 3a 09 24 | 28 50 52 4f 44 55 43 54 |_test:.$|(PRODUCT|
|00000d50| 29 20 24 28 54 46 49 4c | 45 53 29 0a 09 2e 2f 24 |) $(TFIL|ES).../$|
|00000d60| 28 50 52 4f 44 55 43 54 | 29 20 3c 20 24 28 50 52 |(PRODUCT|) < $(PR|
|00000d70| 4f 44 55 43 54 29 2e 69 | 6e 20 3e 20 24 28 50 52 |ODUCT).i|n > $(PR|
|00000d80| 4f 44 55 43 54 29 2e 6f | 75 74 0a 09 40 69 66 20 |ODUCT).o|ut..@if |
|00000d90| 63 6d 70 20 2d 73 20 24 | 28 50 52 4f 44 55 43 54 |cmp -s $|(PRODUCT|
|00000da0| 29 2e 65 78 70 20 24 28 | 50 52 4f 44 55 43 54 29 |).exp $(|PRODUCT)|
|00000db0| 2e 6f 75 74 20 3b 20 5c | 0a 09 74 68 65 6e 09 65 |.out ; \|..then.e|
|00000dc0| 63 68 6f 20 27 54 65 73 | 74 65 64 20 6f 6b 61 79 |cho 'Tes|ted okay|
|00000dd0| 27 20 3b 20 5c 0a 09 65 | 6c 73 65 09 65 63 68 6f |' ; \..e|lse.echo|
|00000de0| 20 27 2a 2a 2a 20 54 65 | 73 74 20 66 61 69 6c 65 | '*** Te|st faile|
|00000df0| 64 3b 20 64 69 66 66 65 | 72 65 6e 63 65 73 3a 27 |d; diffe|rences:'|
|00000e00| 20 3b 20 5c 0a 09 09 64 | 69 66 66 20 24 28 50 52 | ; \...d|iff $(PR|
|00000e10| 4f 44 55 43 54 29 2e 65 | 78 70 20 24 28 50 52 4f |ODUCT).e|xp $(PRO|
|00000e20| 44 55 43 54 29 2e 6f 75 | 74 20 3b 20 5c 0a 09 09 |DUCT).ou|t ; \...|
|00000e30| 65 78 69 74 20 31 20 3b | 20 5c 0a 09 66 69 0a 0a |exit 1 ;| \..fi..|
|00000e40| 67 61 6d 6d 61 5f 74 65 | 73 74 3a 09 24 28 4c 49 |gamma_te|st:.$(LI|
|00000e50| 42 54 45 53 54 29 0a 09 | 2e 2f 24 28 4c 49 42 54 |BTEST)..|./$(LIBT|
|00000e60| 45 53 54 29 0a 0a 24 28 | 4c 49 42 54 45 53 54 29 |EST)..$(|LIBTEST)|
|00000e70| 3a 09 24 28 54 4f 42 4a | 29 20 67 61 6d 6d 61 2e |:.$(TOBJ|) gamma.|
|00000e80| 6f 0a 09 24 28 43 43 29 | 20 24 28 4c 44 46 4c 41 |o..$(CC)| $(LDFLA|
|00000e90| 47 53 29 20 2d 6f 20 24 | 40 20 24 28 54 4f 42 4a |GS) -o $|@ $(TOBJ|
|00000ea0| 29 20 67 61 6d 6d 61 2e | 6f 20 24 28 4c 49 42 45 |) gamma.|o $(LIBE|
|00000eb0| 53 29 0a 0a 24 28 4c 49 | 42 54 45 53 54 29 2e 6f |S)..$(LI|BTEST).o|
|00000ec0| 3a 09 67 61 6d 6d 61 2e | 68 0a 0a 63 6c 65 61 6e |:.gamma.|h..clean|
|00000ed0| 3a 0a 09 2d 72 6d 20 2d | 66 20 24 28 4f 42 4a 45 |:..-rm -|f $(OBJE|
|00000ee0| 43 54 53 29 20 24 28 54 | 4f 42 4a 29 20 24 28 4c |CTS) $(T|OBJ) $(L|
|00000ef0| 49 42 54 45 53 54 29 20 | 24 28 50 52 4f 44 55 43 |IBTEST) |$(PRODUC|
|00000f00| 54 29 2e 6f 75 74 20 5c | 0a 09 6c 69 6e 74 20 66 |T).out \|..lint f|
|00000f10| 6c 6f 77 20 78 72 65 66 | 20 63 6f 72 65 20 6d 6f |low xref| core mo|
|00000f20| 6e 2e 6f 75 74 0a 0a 63 | 6c 6f 62 62 65 72 3a 09 |n.out..c|lobber:.|
|00000f30| 63 6c 65 61 6e 0a 09 2d | 72 6d 20 2d 66 20 24 28 |clean..-|rm -f $(|
|00000f40| 50 52 4f 44 55 43 54 29 | 0a 00 00 00 00 00 00 00 |PRODUCT)|........|
|00000f50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000f60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000f70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000f80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000f90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000fa0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000fb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000fc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000fd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000fe0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000ff0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001000| 74 6f 74 2d 69 6e 66 6f | 2f 63 68 69 73 71 2e 33 |tot-info|/chisq.3|
|00001010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001040| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001050| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001060| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00001070| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00001080| 20 20 31 34 33 34 33 20 | 20 34 34 35 37 32 31 30 |  14343 | 4457210|
|00001090| 34 36 32 20 20 20 37 30 | 34 31 00 20 00 00 00 00 |462   70|41. ....|
|000010a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000010b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000010c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000010d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000010e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000010f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001100| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001110| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001120| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001130| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001140| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001150| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001160| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001170| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001180| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001190| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001200| 27 5c 22 20 65 0a 2e 54 | 48 20 43 48 49 53 51 20 |'\" e..T|H CHISQ |
|00001210| 33 56 20 56 4d 42 0a 27 | 5c 22 09 6c 61 73 74 20 |3V VMB.'|\".last |
|00001220| 65 64 69 74 3a 09 38 38 | 2f 30 39 2f 31 39 09 44 |edit:.88|/09/19.D|
|00001230| 20 41 20 47 77 79 6e 0a | 27 5c 22 09 53 43 43 53 | A Gwyn.|'\".SCCS|
|00001240| 20 49 44 3a 09 40 28 23 | 29 63 68 69 73 71 2e 33 | ID:.@(#|)chisq.3|
|00001250| 09 31 2e 32 20 28 65 64 | 69 74 65 64 20 66 6f 72 |.1.2 (ed|ited for|
|00001260| 20 70 75 62 6c 69 63 61 | 74 69 6f 6e 29 0a 27 5c | publica|tion).'\|
|00001270| 22 0a 27 5c 22 09 54 68 | 69 73 20 73 6f 75 72 63 |".'\".Th|is sourc|
|00001280| 65 20 6d 75 73 74 20 62 | 65 20 74 79 70 65 73 65 |e must b|e typese|
|00001290| 74 20 76 69 61 20 22 65 | 71 6e 20 7c 20 74 72 6f |t via "e|qn | tro|
|000012a0| 66 66 20 2d 6d 61 6e 22 | 0a 27 5c 22 0a 2e 45 51 |ff -man"|.'\"..EQ|
|000012b0| 0a 64 65 6c 69 6d 20 40 | 40 0a 2e 45 4e 0a 2e 53 |.delim @|@..EN..S|
|000012c0| 48 20 4e 41 4d 45 0a 63 | 68 69 73 71 20 5c 2d 20 |H NAME.c|hisq \- |
|000012d0| 69 6e 64 65 70 65 6e 64 | 65 6e 63 65 20 74 65 73 |independ|ence tes|
|000012e0| 74 73 20 66 6f 72 20 32 | 2d 77 61 79 20 63 6f 6e |ts for 2|-way con|
|000012f0| 74 69 6e 67 65 6e 63 79 | 20 74 61 62 6c 65 73 0a |tingency| tables.|
|00001300| 2e 53 48 20 53 59 4e 4f | 50 53 49 53 0a 2e 64 73 |.SH SYNO|PSIS..ds|
|00001310| 20 63 57 20 28 43 57 5c | 22 20 63 68 61 6e 67 65 | cW (CW\|" change|
|00001320| 20 74 6f 20 42 20 28 77 | 69 74 68 6f 75 74 20 74 | to B (w|ithout t|
|00001330| 68 65 20 70 61 72 65 6e | 29 20 69 66 20 79 6f 75 |he paren|) if you|
|00001340| 20 64 6f 6e 27 74 20 68 | 61 76 65 20 61 20 43 57 | don't h|ave a CW|
|00001350| 20 66 6f 6e 74 0a 5c 66 | 5c 2a 28 63 57 0a 23 69 | font.\f|\*(cW.#i|
|00001360| 6e 63 6c 75 64 65 20 22 | 63 68 69 73 71 2e 68 22 |nclude "|chisq.h"|
|00001370| 0a 2e 62 72 0a 2f 2a 20 | 5c 66 50 54 68 65 20 66 |..br./* |\fPThe f|
|00001380| 6f 6c 6c 6f 77 69 6e 67 | 20 66 75 6e 63 74 69 6f |ollowing| functio|
|00001390| 6e 73 20 61 72 65 20 64 | 65 63 6c 61 72 65 64 20 |ns are d|eclared |
|000013a0| 69 6e 20 74 68 69 73 20 | 68 65 61 64 65 72 20 66 |in this |header f|
|000013b0| 69 6c 65 2e 5c 66 5c 2a | 28 63 57 20 2a 2f 0a 2e |ile.\f\*|(cW */..|
|000013c0| 73 70 0a 64 6f 75 62 6c | 65 20 43 68 69 53 71 54 |sp.doubl|e ChiSqT|
|000013d0| 62 6c 28 69 6e 74 20 72 | 2c 20 69 6e 74 20 63 2c |bl(int r|, int c,|
|000013e0| 20 63 6f 6e 73 74 20 6c | 6f 6e 67 20 2a 66 2c 20 | const l|ong *f, |
|000013f0| 69 6e 74 20 2a 64 66 29 | 3b 0a 2e 62 72 0a 64 6f |int *df)|;..br.do|
|00001400| 75 62 6c 65 20 49 6e 66 | 6f 54 62 6c 28 69 6e 74 |uble Inf|oTbl(int|
|00001410| 20 72 2c 20 69 6e 74 20 | 63 2c 20 63 6f 6e 73 74 | r, int |c, const|
|00001420| 20 6c 6f 6e 67 20 2a 66 | 2c 20 69 6e 74 20 2a 64 | long *f|, int *d|
|00001430| 66 29 3b 5c 66 50 0a 2e | 53 48 20 44 45 53 43 52 |f);\fP..|SH DESCR|
|00001440| 49 50 54 49 4f 4e 0a 5c | 66 5c 2a 28 63 57 43 68 |IPTION.\|f\*(cWCh|
|00001450| 69 53 71 54 62 6c 5c 66 | 50 0a 72 65 74 75 72 6e |iSqTbl\f|P.return|
|00001460| 73 20 50 65 61 72 73 6f | 6e 27 73 20 40 63 68 69 |s Pearso|n's @chi|
|00001470| 20 73 75 70 20 32 40 20 | 73 74 61 74 69 73 74 69 | sup 2@ |statisti|
|00001480| 63 0a 66 6f 72 20 74 68 | 65 20 5c 66 5c 2a 28 63 |c.for th|e \f\*(c|
|00001490| 57 72 5c 66 50 2d 62 79 | 2d 5c 66 5c 2a 28 63 57 |Wr\fP-by|-\f\*(cW|
|000014a0| 63 5c 66 50 20 63 6f 6e | 74 69 6e 67 65 6e 63 79 |c\fP con|tingency|
|000014b0| 20 74 61 62 6c 65 0a 6f | 66 20 6f 62 73 65 72 76 | table.o|f observ|
|000014c0| 65 64 20 66 72 65 71 75 | 65 6e 63 69 65 73 20 6f |ed frequ|encies o|
|000014d0| 66 20 6f 63 63 75 72 72 | 65 6e 63 65 20 66 6f 72 |f occurr|ence for|
|000014e0| 20 65 61 63 68 20 63 6f | 6d 62 69 6e 61 74 69 6f | each co|mbinatio|
|000014f0| 6e 20 6f 66 0a 72 6f 77 | 20 61 6e 64 20 63 6f 6c |n of.row| and col|
|00001500| 75 6d 6e 20 63 61 74 65 | 67 6f 72 69 65 73 20 73 |umn cate|gories s|
|00001510| 74 6f 72 65 64 20 69 6e | 20 72 6f 77 2d 6d 61 6a |tored in| row-maj|
|00001520| 6f 72 20 6f 72 64 65 72 | 0a 69 6e 20 74 68 65 20 |or order|.in the |
|00001530| 76 65 63 74 6f 72 20 73 | 74 61 72 74 69 6e 67 20 |vector s|tarting |
|00001540| 61 74 20 5c 66 5c 2a 28 | 63 57 66 5c 66 50 2e 0a |at \f\*(|cWf\fP..|
|00001550| 54 68 65 20 63 6f 6d 70 | 75 74 65 64 20 6e 75 6d |The comp|uted num|
|00001560| 62 65 72 20 6f 66 20 64 | 65 67 72 65 65 73 20 6f |ber of d|egrees o|
|00001570| 66 20 66 72 65 65 64 6f | 6d 20 69 73 20 73 74 6f |f freedo|m is sto|
|00001580| 72 65 64 20 69 6e 74 6f | 20 74 68 65 20 6c 6f 63 |red into| the loc|
|00001590| 61 74 69 6f 6e 0a 70 6f | 69 6e 74 65 64 20 74 6f |ation.po|inted to|
|000015a0| 20 62 79 20 5c 66 5c 2a | 28 63 57 64 66 5c 66 50 | by \f\*|(cWdf\fP|
|000015b0| 2e 0a 2e 62 72 0a 2e 45 | 51 0a 63 68 69 20 73 75 |...br..E|Q.chi su|
|000015c0| 70 20 32 20 7e 3d 3d 7e | 20 72 6f 6d 61 6e 20 73 |p 2 ~==~| roman s|
|000015d0| 75 6d 20 66 72 6f 6d 20 | 69 3d 30 20 74 6f 20 72 |um from |i=0 to r|
|000015e0| 2d 31 20 72 6f 6d 61 6e | 20 73 75 6d 20 66 72 6f |-1 roman| sum fro|
|000015f0| 6d 20 6a 3d 30 20 74 6f | 20 63 2d 31 0a 7b 20 6c |m j=0 to| c-1.{ l|
|00001600| 65 66 74 20 28 20 7b 20 | 66 20 73 75 62 20 69 6a |eft ( { |f sub ij|
|00001610| 20 7e 2d 7e 20 7b 20 66 | 20 73 75 62 20 7b 20 69 | ~-~ { f| sub { i|
|00001620| 2e 20 7d 20 66 20 73 75 | 62 20 7b 20 2e 6a 20 7d |. } f su|b { .j }|
|00001630| 7d 20 6f 76 65 72 20 4e | 20 7d 20 72 69 67 68 74 |} over N| } right|
|00001640| 20 29 20 73 75 70 20 32 | 20 7d 0a 6f 76 65 72 20 | ) sup 2| }.over |
|00001650| 7b 20 7b 20 66 20 73 75 | 62 20 7b 20 69 2e 20 7d |{ { f su|b { i. }|
|00001660| 20 66 20 73 75 62 20 7b | 20 2e 6a 20 7d 7d 20 6f | f sub {| .j }} o|
|00001670| 76 65 72 20 4e 20 7d 0a | 2e 45 4e 0a 77 68 65 72 |ver N }.|.EN.wher|
|00001680| 65 20 74 68 65 20 72 6f | 77 20 61 6e 64 20 63 6f |e the ro|w and co|
|00001690| 6c 75 6d 6e 20 73 75 6d | 73 20 61 72 65 20 64 65 |lumn sum|s are de|
|000016a0| 6e 6f 74 65 64 0a 40 66 | 20 73 75 62 20 7b 20 69 |noted.@f| sub { i|
|000016b0| 2e 20 7d 7e 3d 3d 7e 72 | 6f 6d 61 6e 20 73 75 6d |. }~==~r|oman sum|
|000016c0| 20 66 72 6f 6d 20 6a 3d | 30 20 74 6f 20 63 2d 31 | from j=|0 to c-1|
|000016d0| 20 66 20 73 75 62 20 69 | 6a 40 0a 61 6e 64 0a 40 | f sub i|j@.and.@|
|000016e0| 66 20 73 75 62 20 7b 20 | 2e 6a 20 7d 7e 3d 3d 7e |f sub { |.j }~==~|
|000016f0| 72 6f 6d 61 6e 20 73 75 | 6d 20 66 72 6f 6d 20 69 |roman su|m from i|
|00001700| 3d 30 20 74 6f 20 72 2d | 31 20 66 20 73 75 62 20 |=0 to r-|1 f sub |
|00001710| 69 6a 40 0a 61 6e 64 20 | 74 68 65 20 74 6f 74 61 |ij@.and |the tota|
|00001720| 6c 20 6e 75 6d 62 65 72 | 20 6f 66 20 6f 63 63 75 |l number| of occu|
|00001730| 72 72 65 6e 63 65 73 20 | 69 73 0a 40 4e 7e 3d 3d |rrences |is.@N~==|
|00001740| 7e 72 6f 6d 61 6e 20 73 | 75 6d 20 66 72 6f 6d 20 |~roman s|um from |
|00001750| 69 3d 30 20 74 6f 20 72 | 2d 31 20 72 6f 6d 61 6e |i=0 to r|-1 roman|
|00001760| 20 73 75 6d 20 66 72 6f | 6d 20 6a 3d 30 20 74 6f | sum fro|m j=0 to|
|00001770| 20 63 2d 31 20 66 20 73 | 75 62 20 69 6a 40 2e 0a | c-1 f s|ub ij@..|
|00001780| 54 65 72 6d 73 20 69 6e | 76 6f 6c 76 69 6e 67 20 |Terms in|volving |
|00001790| 7a 65 72 6f 20 72 6f 77 | 20 6f 72 20 63 6f 6c 75 |zero row| or colu|
|000017a0| 6d 6e 20 73 75 6d 20 61 | 72 65 20 6f 6d 69 74 74 |mn sum a|re omitt|
|000017b0| 65 64 20 61 6e 64 20 74 | 68 65 20 72 65 74 75 72 |ed and t|he retur|
|000017c0| 6e 65 64 0a 6e 75 6d 62 | 65 72 20 6f 66 20 64 65 |ned.numb|er of de|
|000017d0| 67 72 65 65 73 20 6f 66 | 20 66 72 65 65 64 6f 6d |grees of| freedom|
|000017e0| 20 61 72 65 0a 63 6f 72 | 72 65 73 70 6f 6e 64 69 | are.cor|respondi|
|000017f0| 6e 67 6c 79 20 72 65 64 | 75 63 65 64 20 66 72 6f |ngly red|uced fro|
|00001800| 6d 20 74 68 65 20 6e 6f | 6d 69 6e 61 6c 20 76 61 |m the no|minal va|
|00001810| 6c 75 65 0a 40 28 72 2d | 31 29 28 63 2d 31 29 40 |lue.@(r-|1)(c-1)@|
|00001820| 2e 0a 2e 50 0a 5c 66 5c | 2a 28 63 57 49 6e 66 6f |...P.\f\|*(cWInfo|
|00001830| 54 62 6c 5c 66 50 0a 72 | 65 74 75 72 6e 73 20 32 |Tbl\fP.r|eturns 2|
|00001840| 20 74 69 6d 65 73 20 4b | 75 6c 6c 62 61 63 6b 27 | times K|ullback'|
|00001850| 73 0a 40 22 5c 76 27 2d | 30 2e 32 27 5e 5c 76 27 |s.@"\v'-|0.2'^\v'|
|00001860| 30 2e 32 27 5c 68 27 2d | 5c 77 3b 5e 3b 75 27 49 |0.2'\h'-|\w;^;u'I|
|00001870| 5c 5e 22 20 28 20 48 20 | 73 75 62 20 31 20 3a 20 |\^" ( H |sub 1 : |
|00001880| 48 20 73 75 62 20 32 20 | 29 40 20 73 74 61 74 69 |H sub 2 |)@ stati|
|00001890| 73 74 69 63 0a 66 6f 72 | 20 61 20 73 69 6d 69 6c |stic.for| a simil|
|000018a0| 61 72 20 63 6f 6e 74 69 | 6e 67 65 6e 63 79 20 74 |ar conti|ngency t|
|000018b0| 61 62 6c 65 2e 0a 54 68 | 65 20 6e 75 6d 62 65 72 |able..Th|e number|
|000018c0| 20 6f 66 20 64 65 67 72 | 65 65 73 20 6f 66 20 66 | of degr|ees of f|
|000018d0| 72 65 65 64 6f 6d 20 74 | 6f 20 62 65 20 75 73 65 |reedom t|o be use|
|000018e0| 64 20 77 68 65 6e 20 72 | 65 6c 61 74 69 6e 67 20 |d when r|elating |
|000018f0| 74 68 69 73 20 73 74 61 | 74 69 73 74 69 63 0a 74 |this sta|tistic.t|
|00001900| 6f 20 74 68 65 20 40 63 | 68 69 20 73 75 70 20 32 |o the @c|hi sup 2|
|00001910| 40 20 64 69 73 74 72 69 | 62 75 74 69 6f 6e 20 69 |@ distri|bution i|
|00001920| 73 20 73 74 6f 72 65 64 | 20 69 6e 74 6f 20 74 68 |s stored| into th|
|00001930| 65 20 6c 6f 63 61 74 69 | 6f 6e 0a 70 6f 69 6e 74 |e locati|on.point|
|00001940| 65 64 20 74 6f 20 62 79 | 20 5c 66 5c 2a 28 63 57 |ed to by| \f\*(cW|
|00001950| 64 66 5c 66 50 2e 0a 28 | 53 65 65 20 5c 73 2d 31 |df\fP..(|See \s-1|
|00001960| 44 49 53 43 55 53 53 49 | 4f 4e 5c 73 30 20 62 65 |DISCUSSI|ON\s0 be|
|00001970| 6c 6f 77 2e 29 0a 2e 62 | 72 0a 2e 45 51 0a 32 20 |low.)..b|r..EQ.2 |
|00001980| 22 5c 76 27 2d 30 2e 32 | 27 5e 5c 76 27 30 2e 32 |"\v'-0.2|'^\v'0.2|
|00001990| 27 5c 68 27 2d 5c 77 3b | 5e 3b 75 27 49 5c 5e 22 |'\h'-\w;|^;u'I\^"|
|000019a0| 20 28 20 48 20 73 75 62 | 20 31 20 3a 20 48 20 73 | ( H sub| 1 : H s|
|000019b0| 75 62 20 32 20 29 20 7e | 3d 3d 7e 0a 32 20 72 6f |ub 2 ) ~|==~.2 ro|
|000019c0| 6d 61 6e 20 73 75 6d 20 | 66 72 6f 6d 20 69 3d 30 |man sum |from i=0|
|000019d0| 20 74 6f 20 72 2d 31 20 | 72 6f 6d 61 6e 20 73 75 | to r-1 |roman su|
|000019e0| 6d 20 66 72 6f 6d 20 6a | 3d 30 20 74 6f 20 63 2d |m from j|=0 to c-|
|000019f0| 31 0a 7b 20 66 20 73 75 | 62 20 69 6a 20 6c 6f 67 |1.{ f su|b ij log|
|00001a00| 20 7b 20 7b 20 4e 20 66 | 20 73 75 62 20 69 6a 20 | { { N f| sub ij |
|00001a10| 7d 0a 6f 76 65 72 20 7b | 20 66 20 73 75 62 20 7b |}.over {| f sub {|
|00001a20| 20 69 2e 20 7d 20 66 20 | 73 75 62 20 7b 20 2e 6a | i. } f |sub { .j|
|00001a30| 20 7d 20 7d 20 7d 20 7d | 0a 2e 45 4e 0a 77 68 65 | } } } }|..EN.whe|
|00001a40| 72 65 20 74 68 65 20 72 | 6f 77 20 61 6e 64 20 63 |re the r|ow and c|
|00001a50| 6f 6c 75 6d 6e 20 73 75 | 6d 73 20 40 66 20 73 75 |olumn su|ms @f su|
|00001a60| 62 20 7b 20 69 2e 20 7d | 40 20 61 6e 64 20 40 66 |b { i. }|@ and @f|
|00001a70| 20 73 75 62 20 7b 20 2e | 6a 20 7d 40 0a 61 6e 64 | sub { .|j }@.and|
|00001a80| 20 74 68 65 20 74 6f 74 | 61 6c 20 6e 75 6d 62 65 | the tot|al numbe|
|00001a90| 72 20 6f 66 20 6f 62 73 | 65 72 76 61 74 69 6f 6e |r of obs|ervation|
|00001aa0| 73 20 40 4e 40 0a 61 72 | 65 20 74 68 65 20 73 61 |s @N@.ar|e the sa|
|00001ab0| 6d 65 20 61 73 20 66 6f | 72 20 50 65 61 72 73 6f |me as fo|r Pearso|
|00001ac0| 6e 27 73 20 40 63 68 69 | 20 73 75 70 20 32 40 2e |n's @chi| sup 2@.|
|00001ad0| 0a 2e 53 48 20 44 49 53 | 43 55 53 53 49 4f 4e 0a |..SH DIS|CUSSION.|
|00001ae0| 40 32 20 22 5c 76 27 2d | 30 2e 32 27 5e 5c 76 27 |@2 "\v'-|0.2'^\v'|
|00001af0| 30 2e 32 27 5c 68 27 2d | 5c 77 3b 5e 3b 75 27 49 |0.2'\h'-|\w;^;u'I|
|00001b00| 5c 5e 22 20 28 20 48 20 | 73 75 62 20 31 20 3a 20 |\^" ( H |sub 1 : |
|00001b10| 48 20 73 75 62 20 32 20 | 29 40 0a 72 65 70 72 65 |H sub 2 |)@.repre|
|00001b20| 73 65 6e 74 73 20 74 77 | 69 63 65 20 74 68 65 20 |sents tw|ice the |
|00001b30| 65 73 74 69 6d 61 74 65 | 64 20 69 6e 66 6f 72 6d |estimate|d inform|
|00001b40| 61 74 69 6f 6e 0a 69 6e | 20 66 61 76 6f 72 20 6f |ation.in| favor o|
|00001b50| 66 20 68 79 70 6f 74 68 | 65 73 69 73 20 40 48 20 |f hypoth|esis @H |
|00001b60| 73 75 62 20 31 40 20 6f | 76 65 72 20 68 79 70 6f |sub 1@ o|ver hypo|
|00001b70| 74 68 65 73 69 73 20 40 | 48 20 73 75 62 20 32 40 |thesis @|H sub 2@|
|00001b80| 0a 63 6f 6e 74 61 69 6e | 65 64 20 69 6e 20 74 68 |.contain|ed in th|
|00001b90| 65 20 6f 62 73 65 72 76 | 65 64 20 66 72 65 71 75 |e observ|ed frequ|
|00001ba0| 65 6e 63 69 65 73 2c 0a | 77 68 65 72 65 20 74 68 |encies,.|where th|
|00001bb0| 65 20 5c 66 49 6e 75 6c | 6c 20 68 79 70 6f 74 68 |e \fInul|l hypoth|
|00001bc0| 65 73 69 73 5c 66 50 20 | 40 48 20 73 75 62 20 32 |esis\fP |@H sub 2|
|00001bd0| 40 20 69 73 20 74 68 61 | 74 20 74 68 65 20 72 6f |@ is tha|t the ro|
|00001be0| 77 20 61 6e 64 20 63 6f | 6c 75 6d 6e 0a 63 61 74 |w and co|lumn.cat|
|00001bf0| 65 67 6f 72 69 7a 61 74 | 69 6f 6e 73 20 61 72 65 |egorizat|ions are|
|00001c00| 20 73 74 61 74 69 73 74 | 69 63 61 6c 6c 79 20 69 | statist|ically i|
|00001c10| 6e 64 65 70 65 6e 64 65 | 6e 74 20 61 6e 64 0a 74 |ndepende|nt and.t|
|00001c20| 68 65 20 5c 66 49 61 6c | 74 65 72 6e 61 74 69 76 |he \fIal|ternativ|
|00001c30| 65 20 68 79 70 6f 74 68 | 65 73 69 73 5c 66 50 20 |e hypoth|esis\fP |
|00001c40| 40 48 20 73 75 62 20 31 | 40 20 69 73 20 74 68 61 |@H sub 1|@ is tha|
|00001c50| 74 20 74 68 65 79 20 61 | 72 65 6e 27 74 2e 0a 54 |t they a|ren't..T|
|00001c60| 68 69 73 20 73 74 61 74 | 69 73 74 69 63 20 69 73 |his stat|istic is|
|00001c70| 20 61 73 79 6d 70 74 6f | 74 69 63 61 6c 6c 79 20 | asympto|tically |
|00001c80| 64 69 73 74 72 69 62 75 | 74 65 64 20 61 73 20 40 |distribu|ted as @|
|00001c90| 63 68 69 20 73 75 70 20 | 32 40 0a 77 69 74 68 20 |chi sup |2@.with |
|00001ca0| 40 28 72 2d 31 29 28 63 | 2d 31 29 40 20 64 65 67 |@(r-1)(c|-1)@ deg|
|00001cb0| 72 65 65 73 20 6f 66 20 | 66 72 65 65 64 6f 6d 3b |rees of |freedom;|
|00001cc0| 0a 74 68 65 72 65 66 6f | 72 65 20 62 79 20 75 73 |.therefo|re by us|
|00001cd0| 65 20 6f 66 20 74 68 65 | 20 5c 66 5c 2a 28 63 57 |e of the| \f\*(cW|
|00001ce0| 51 43 68 69 53 71 5c 66 | 50 20 66 75 6e 63 74 69 |QChiSq\f|P functi|
|00001cf0| 6f 6e 20 28 73 65 65 20 | 5c 66 49 67 61 6d 6d 61 |on (see |\fIgamma|
|00001d00| 5c 66 50 28 33 56 29 29 | 0a 69 74 20 63 61 6e 20 |\fP(3V))|.it can |
|00001d10| 62 65 20 63 6f 6e 76 65 | 72 74 65 64 20 74 6f 20 |be conve|rted to |
|00001d20| 74 68 65 20 70 72 6f 62 | 61 62 69 6c 69 74 79 20 |the prob|ability |
|00001d30| 74 68 61 74 0a 74 68 65 | 20 73 74 61 74 69 73 74 |that.the| statist|
|00001d40| 69 63 20 77 6f 75 6c 64 | 20 62 65 20 61 73 20 6c |ic would| be as l|
|00001d50| 61 72 67 65 20 61 73 20 | 77 61 73 20 6f 62 73 65 |arge as |was obse|
|00001d60| 72 76 65 64 0a 69 66 20 | 74 68 65 20 63 61 74 65 |rved.if |the cate|
|00001d70| 67 6f 72 69 7a 61 74 69 | 6f 6e 73 20 72 65 61 6c |gorizati|ons real|
|00001d80| 6c 79 20 61 72 65 20 69 | 6e 64 65 70 65 6e 64 65 |ly are i|ndepende|
|00001d90| 6e 74 2e 0a 54 68 69 73 | 20 69 73 20 6f 66 20 63 |nt..This| is of c|
|00001da0| 6f 75 72 73 65 20 74 68 | 65 20 74 72 61 64 69 74 |ourse th|e tradit|
|00001db0| 69 6f 6e 61 6c 20 75 73 | 65 20 6f 66 20 50 65 61 |ional us|e of Pea|
|00001dc0| 72 73 6f 6e 27 73 20 40 | 63 68 69 20 73 75 70 20 |rson's @|chi sup |
|00001dd0| 32 40 20 73 74 61 74 69 | 73 74 69 63 2c 0a 77 68 |2@ stati|stic,.wh|
|00001de0| 69 63 68 20 40 32 20 22 | 5c 76 27 2d 30 2e 32 27 |ich @2 "|\v'-0.2'|
|00001df0| 5e 5c 76 27 30 2e 32 27 | 5c 68 27 2d 5c 77 3b 5e |^\v'0.2'|\h'-\w;^|
|00001e00| 3b 75 27 49 5c 5e 22 40 | 20 61 70 70 72 6f 61 63 |;u'I\^"@| approac|
|00001e10| 68 65 73 20 66 6f 72 20 | 6c 61 72 67 65 20 73 61 |hes for |large sa|
|00001e20| 6d 70 6c 65 73 0a 77 68 | 65 6e 20 40 48 20 73 75 |mples.wh|en @H su|
|00001e30| 62 20 32 40 20 69 73 20 | 74 72 75 65 2e 0a 54 68 |b 2@ is |true..Th|
|00001e40| 65 20 6d 61 69 6e 20 64 | 69 66 66 65 72 65 6e 63 |e main d|ifferenc|
|00001e50| 65 20 69 73 20 74 68 61 | 74 20 50 65 61 72 73 6f |e is tha|t Pearso|
|00001e60| 6e 27 73 20 40 63 68 69 | 20 73 75 70 20 32 40 20 |n's @chi| sup 2@ |
|00001e70| 69 73 20 6e 6f 74 20 75 | 73 65 66 75 6c 0a 66 6f |is not u|seful.fo|
|00001e80| 72 20 73 6d 61 6c 6c 20 | 73 61 6d 70 6c 65 20 73 |r small |sample s|
|00001e90| 69 7a 65 73 2c 0a 77 68 | 65 72 65 61 73 20 40 32 |izes,.wh|ereas @2|
|00001ea0| 20 22 5c 76 27 2d 30 2e | 32 27 5e 5c 76 27 30 2e | "\v'-0.|2'^\v'0.|
|00001eb0| 32 27 5c 68 27 2d 5c 77 | 3b 5e 3b 75 27 49 5c 5e |2'\h'-\w|;^;u'I\^|
|00001ec0| 22 40 20 66 75 6c 6c 79 | 20 74 61 6b 65 73 20 69 |"@ fully| takes i|
|00001ed0| 6e 74 6f 20 61 63 63 6f | 75 6e 74 0a 61 6c 6c 20 |nto acco|unt.all |
|00001ee0| 61 76 61 69 6c 61 62 6c | 65 20 69 6e 66 6f 72 6d |availabl|e inform|
|00001ef0| 61 74 69 6f 6e 20 66 6f | 72 20 65 76 65 6e 20 74 |ation fo|r even t|
|00001f00| 68 65 20 73 6d 61 6c 6c | 65 73 74 20 73 61 6d 70 |he small|est samp|
|00001f10| 6c 65 73 2e 0a 2e 50 0a | 40 32 20 22 5c 76 27 2d |les...P.|@2 "\v'-|
|00001f20| 30 2e 32 27 5e 5c 76 27 | 30 2e 32 27 5c 68 27 2d |0.2'^\v'|0.2'\h'-|
|00001f30| 5c 77 3b 5e 3b 75 27 49 | 5c 5e 22 20 28 20 48 20 |\w;^;u'I|\^" ( H |
|00001f40| 73 75 62 20 31 20 3a 20 | 48 20 73 75 62 20 32 20 |sub 1 : |H sub 2 |
|00001f50| 29 40 0a 28 61 6c 6f 6e | 67 20 77 69 74 68 20 69 |)@.(alon|g with i|
|00001f60| 74 73 20 63 6f 72 72 65 | 73 70 6f 6e 64 69 6e 67 |ts corre|sponding|
|00001f70| 20 64 65 67 72 65 65 73 | 20 6f 66 20 66 72 65 65 | degrees| of free|
|00001f80| 64 6f 6d 0a 66 6f 72 20 | 75 73 65 20 77 69 74 68 |dom.for |use with|
|00001f90| 20 5c 66 5c 2a 28 63 57 | 51 43 68 69 53 71 5c 66 | \f\*(cW|QChiSq\f|
|00001fa0| 50 29 20 69 73 20 5c 66 | 49 61 64 64 69 74 69 76 |P) is \f|Iadditiv|
|00001fb0| 65 5c 66 50 20 66 6f 72 | 20 69 6e 64 65 70 65 6e |e\fP for| indepen|
|00001fc0| 64 65 6e 74 20 73 61 6d | 70 6c 65 73 2c 0a 73 6f |dent sam|ples,.so|
|00001fd0| 20 74 68 61 74 20 74 68 | 65 20 69 6e 66 6f 72 6d | that th|e inform|
|00001fe0| 61 74 69 6f 6e 20 63 61 | 6e 20 62 65 20 61 63 63 |ation ca|n be acc|
|00001ff0| 75 6d 75 6c 61 74 65 64 | 20 6f 76 65 72 20 74 68 |umulated| over th|
|00002000| 65 20 63 6f 75 72 73 65 | 20 6f 66 20 6d 61 6e 79 |e course| of many|
|00002010| 0a 69 6e 64 65 70 65 6e | 64 65 6e 74 20 65 78 70 |.indepen|dent exp|
|00002020| 65 72 69 6d 65 6e 74 73 | 20 69 6e 20 6f 72 64 65 |eriments| in orde|
|00002030| 72 20 74 6f 20 69 6d 70 | 72 6f 76 65 20 6f 6e 65 |r to imp|rove one|
|00002040| 27 73 20 61 62 69 6c 69 | 74 79 0a 74 6f 20 64 65 |'s abili|ty.to de|
|00002050| 74 65 63 74 20 61 20 76 | 69 6f 6c 61 74 69 6f 6e |tect a v|iolation|
|00002060| 20 6f 66 20 74 68 65 20 | 6e 75 6c 6c 20 68 79 70 | of the |null hyp|
|00002070| 6f 74 68 65 73 69 73 2e | 0a 2e 53 48 20 43 41 56 |othesis.|..SH CAV|
|00002080| 45 41 54 53 0a 50 65 61 | 72 73 6f 6e 27 73 20 40 |EATS.Pea|rson's @|
|00002090| 63 68 69 20 73 75 70 20 | 32 40 20 74 65 73 74 20 |chi sup |2@ test |
|000020a0| 69 73 20 75 6e 72 65 6c | 69 61 62 6c 65 20 66 6f |is unrel|iable fo|
|000020b0| 72 20 6c 6f 77 20 66 72 | 65 71 75 65 6e 63 69 65 |r low fr|equencie|
|000020c0| 73 3b 0a 63 6f 6e 73 65 | 71 75 65 6e 74 6c 79 20 |s;.conse|quently |
|000020d0| 69 74 20 69 73 20 75 73 | 75 61 6c 6c 79 20 72 65 |it is us|ually re|
|000020e0| 63 6f 6d 6d 65 6e 64 65 | 64 20 74 68 61 74 0a 63 |commende|d that.c|
|000020f0| 61 74 65 67 6f 72 69 65 | 73 20 62 65 20 63 68 6f |ategorie|s be cho|
|00002100| 73 65 6e 20 74 6f 20 74 | 61 6c 6c 79 20 61 74 20 |sen to t|ally at |
|00002110| 6c 65 61 73 74 0a 35 20 | 6f 63 63 75 72 72 65 6e |least.5 |occurren|
|00002120| 63 65 73 20 69 6e 20 65 | 61 63 68 20 62 69 6e 2e |ces in e|ach bin.|
|00002130| 0a 48 6f 77 65 76 65 72 | 2c 20 65 78 63 65 73 73 |.However|, excess|
|00002140| 69 76 65 20 63 6f 6d 62 | 69 6e 61 74 69 6f 6e 20 |ive comb|ination |
|00002150| 6f 66 20 62 69 6e 73 20 | 63 61 75 73 65 73 20 6c |of bins |causes l|
|00002160| 6f 73 73 20 6f 66 20 69 | 6e 66 6f 72 6d 61 74 69 |oss of i|nformati|
|00002170| 6f 6e 20 61 6e 64 0a 63 | 6f 6e 73 65 71 75 65 6e |on and.c|onsequen|
|00002180| 74 6c 79 20 6c 6f 73 73 | 20 6f 66 20 64 69 73 63 |tly loss| of disc|
|00002190| 72 69 6d 69 6e 61 74 69 | 6e 67 20 70 6f 77 65 72 |riminati|ng power|
|000021a0| 2e 0a 42 65 63 61 75 73 | 65 20 40 32 20 22 5c 76 |..Becaus|e @2 "\v|
|000021b0| 27 2d 30 2e 32 27 5e 5c | 76 27 30 2e 32 27 5c 68 |'-0.2'^\|v'0.2'\h|
|000021c0| 27 2d 5c 77 3b 5e 3b 75 | 27 49 5c 5e 22 40 20 69 |'-\w;^;u|'I\^"@ i|
|000021d0| 73 20 63 6f 72 72 65 63 | 74 20 66 6f 72 20 61 6e |s correc|t for an|
|000021e0| 79 20 73 61 6d 70 6c 65 | 20 73 69 7a 65 2c 0a 69 |y sample| size,.i|
|000021f0| 74 20 64 6f 65 73 20 6e | 6f 74 20 72 65 71 75 69 |t does n|ot requi|
|00002200| 72 65 20 63 6f 6d 62 69 | 6e 61 74 69 6f 6e 20 6f |re combi|nation o|
|00002210| 66 20 62 69 6e 73 0a 61 | 6e 64 20 69 73 20 74 68 |f bins.a|nd is th|
|00002220| 65 72 65 66 6f 72 65 20 | 69 6d 6d 75 6e 65 20 66 |erefore |immune f|
|00002230| 72 6f 6d 20 74 68 69 73 | 20 70 72 6f 62 6c 65 6d |rom this| problem|
|00002240| 2e 0a 2e 62 70 0a 2e 53 | 48 20 45 58 41 4d 50 4c |...bp..S|H EXAMPL|
|00002250| 45 0a 5c 66 5c 2a 28 63 | 57 0a 2e 74 61 20 38 6e |E.\f\*(c|W..ta 8n|
|00002260| 20 31 36 6e 20 32 34 6e | 20 33 32 6e 20 34 30 6e | 16n 24n| 32n 40n|
|00002270| 20 34 38 6e 20 35 36 6e | 20 36 34 6e 0a 2e 6e 66 | 48n 56n| 64n..nf|
|00002280| 0a 5c 26 0a 2f 2a 0a 09 | 45 78 61 6d 70 6c 65 20 |.\&./*..|Example |
|00002290| 70 72 6f 67 72 61 6d 20 | 74 6f 20 72 65 61 64 20 |program |to read |
|000022a0| 74 61 62 6c 65 20 64 69 | 6d 65 6e 73 69 6f 6e 73 |table di|mensions|
|000022b0| 2c 20 74 68 65 6e 20 66 | 72 65 71 75 65 6e 63 69 |, then f|requenci|
|000022c0| 65 73 2c 0a 09 61 6e 64 | 20 74 6f 20 70 72 69 6e |es,..and| to prin|
|000022d0| 74 20 62 6f 74 68 20 73 | 74 61 74 69 73 74 69 63 |t both s|tatistic|
|000022e0| 73 20 61 6c 6f 6e 67 20 | 77 69 74 68 20 73 69 67 |s along |with sig|
|000022f0| 6e 69 66 69 63 61 6e 63 | 65 20 6c 65 76 65 6c 73 |nificanc|e levels|
|00002300| 2e 0a 2a 2f 0a 23 69 6e | 63 6c 75 64 65 09 3c 73 |..*/.#in|clude.<s|
|00002310| 74 64 69 6f 2e 68 3e 0a | 23 69 6e 63 6c 75 64 65 |tdio.h>.|#include|
|00002320| 09 3c 73 74 64 6c 69 62 | 2e 68 3e 09 09 09 2f 2a |.<stdlib|.h>.../*|
|00002330| 20 66 6f 72 20 45 58 49 | 54 5f 53 55 43 43 45 53 | for EXI|T_SUCCES|
|00002340| 53 20 2a 2f 0a 23 69 6e | 63 6c 75 64 65 09 22 63 |S */.#in|clude."c|
|00002350| 68 69 73 71 2e 68 22 0a | 23 69 6e 63 6c 75 64 65 |hisq.h".|#include|
|00002360| 09 22 67 61 6d 6d 61 2e | 68 22 09 09 09 2f 2a 20 |."gamma.|h".../* |
|00002370| 66 6f 72 20 51 43 68 69 | 53 71 28 5c 7c 29 20 2a |for QChi|Sq(\|) *|
|00002380| 2f 0a 0a 23 64 65 66 69 | 6e 65 20 4d 41 58 54 42 |/..#defi|ne MAXTB|
|00002390| 4c 09 31 30 30 30 0a 73 | 74 61 74 69 63 20 6c 6f |L.1000.s|tatic lo|
|000023a0| 6e 67 09 66 5b 4d 41 58 | 54 42 4c 5d 3b 09 09 2f |ng.f[MAX|TBL];../|
|000023b0| 2a 20 66 72 65 71 75 65 | 6e 63 79 20 74 61 6c 6c |* freque|ncy tall|
|000023c0| 69 65 73 20 2a 2f 0a 73 | 74 61 74 69 63 20 69 6e |ies */.s|tatic in|
|000023d0| 74 09 72 3b 09 09 09 2f | 2a 20 23 20 6f 66 20 72 |t.r;.../|* # of r|
|000023e0| 6f 77 73 20 2a 2f 0a 73 | 74 61 74 69 63 20 69 6e |ows */.s|tatic in|
|000023f0| 74 09 63 3b 09 09 09 2f | 2a 20 23 20 6f 66 20 63 |t.c;.../|* # of c|
|00002400| 6f 6c 75 6d 6e 73 20 2a | 2f 0a 23 64 65 66 69 6e |olumns *|/.#defin|
|00002410| 65 20 78 28 69 2c 6a 29 | 09 66 5b 28 69 29 2a 63 |e x(i,j)|.f[(i)*c|
|00002420| 2b 28 6a 29 5d 09 2f 2a | 20 63 6f 6e 76 65 6e 69 |+(j)]./*| conveni|
|00002430| 65 6e 74 20 77 61 79 20 | 74 6f 20 61 63 63 65 73 |ent way |to acces|
|00002440| 73 20 66 72 65 71 73 20 | 2a 2f 0a 0a 73 74 61 74 |s freqs |*/..stat|
|00002450| 69 63 20 76 6f 69 64 20 | 43 61 6c 63 28 63 68 61 |ic void |Calc(cha|
|00002460| 72 20 2a 6e 61 6d 65 2c | 20 64 6f 75 62 6c 65 20 |r *name,| double |
|00002470| 28 2a 66 75 6e 63 29 28 | 69 6e 74 2c 20 69 6e 74 |(*func)(|int, int|
|00002480| 2c 20 63 6f 6e 73 74 20 | 6c 6f 6e 67 20 2a 2c 20 |, const |long *, |
|00002490| 69 6e 74 20 2a 29 29 0a | 7b 0a 09 69 6e 74 09 64 |int *)).|{..int.d|
|000024a0| 66 3b 09 09 09 2f 2a 20 | 64 65 67 72 65 65 73 20 |f;.../* |degrees |
|000024b0| 6f 66 20 66 72 65 65 64 | 6f 6d 20 66 6f 72 20 63 |of freed|om for c|
|000024c0| 68 69 2d 73 71 75 61 72 | 65 20 2a 2f 0a 09 64 6f |hi-squar|e */..do|
|000024d0| 75 62 6c 65 09 73 74 61 | 74 20 3d 20 28 2a 66 75 |uble.sta|t = (*fu|
|000024e0| 6e 63 29 28 72 2c 20 63 | 2c 20 66 2c 20 26 64 66 |nc)(r, c|, f, &df|
|000024f0| 29 3b 09 09 2f 2a 20 63 | 6f 6d 70 75 74 65 64 20 |);../* c|omputed |
|00002500| 73 74 61 74 69 73 74 69 | 63 20 2a 2f 0a 0a 09 2f |statisti|c */.../|
|00002510| 2a 20 70 72 69 6e 74 20 | 72 65 73 75 6c 74 73 20 |* print |results |
|00002520| 2a 2f 0a 09 69 66 20 28 | 20 73 74 61 74 20 3e 3d |*/..if (| stat >=|
|00002530| 20 30 2e 30 20 29 0a 09 | 09 28 76 6f 69 64 29 70 | 0.0 )..|.(void)p|
|00002540| 72 69 6e 74 66 28 22 25 | 73 20 3d 20 25 35 2e 32 |rintf("%|s = %5.2|
|00002550| 66 5c 65 74 64 66 20 3d | 20 25 32 64 5c 65 74 71 |f\etdf =| %2d\etq|
|00002560| 20 3d 20 25 37 2e 34 66 | 5c 65 6e 22 2c 0a 09 09 | = %7.4f|\en",...|
|00002570| 09 20 20 20 20 20 6e 61 | 6d 65 2c 20 73 74 61 74 |.     na|me, stat|
|00002580| 2c 20 64 66 2c 20 51 43 | 68 69 53 71 28 73 74 61 |, df, QC|hiSq(sta|
|00002590| 74 2c 20 64 66 29 0a 09 | 09 09 20 20 20 20 29 3b |t, df)..|..    );|
|000025a0| 0a 09 65 6c 73 65 0a 09 | 09 28 76 6f 69 64 29 70 |..else..|.(void)p|
|000025b0| 72 69 6e 74 66 28 73 74 | 61 74 20 3c 20 2d 33 2e |rintf(st|at < -3.|
|000025c0| 35 20 3f 20 22 6f 75 74 | 20 6f 66 20 6d 65 6d 6f |5 ? "out| of memo|
|000025d0| 72 79 5c 65 6e 22 0a 09 | 09 09 20 20 20 3a 20 73 |ry\en"..|..   : s|
|000025e0| 74 61 74 20 3c 20 2d 32 | 2e 35 20 3f 20 22 74 61 |tat < -2|.5 ? "ta|
|000025f0| 62 6c 65 20 74 6f 6f 20 | 73 6d 61 6c 6c 5c 65 6e |ble too |small\en|
|00002600| 22 0a 09 09 09 20 20 20 | 3a 20 73 74 61 74 20 3c |"....   |: stat <|
|00002610| 20 2d 31 2e 35 20 3f 20 | 22 6e 65 67 61 74 69 76 | -1.5 ? |"negativ|
|00002620| 65 20 66 72 65 71 75 65 | 6e 63 79 5c 65 6e 22 0a |e freque|ncy\en".|
|00002630| 09 09 09 20 20 20 3a 20 | 22 74 6f 6f 20 6d 61 6e |...   : |"too man|
|00002640| 79 20 7a 65 72 6f 73 5c | 65 6e 22 0a 09 09 09 20 |y zeros\|en".... |
|00002650| 20 20 20 29 3b 0a 7d 0a | 0a 69 6e 74 20 6d 61 69 |   );.}.|.int mai|
|00002660| 6e 28 69 6e 74 20 61 72 | 67 63 2c 20 63 68 61 72 |n(int ar|gc, char|
|00002670| 20 2a 61 72 67 76 5b 5d | 29 0a 7b 0a 09 72 65 67 | *argv[]|).{..reg|
|00002680| 69 73 74 65 72 20 69 6e | 74 09 69 3b 09 09 2f 2a |ister in|t.i;../*|
|00002690| 20 72 6f 77 20 69 6e 64 | 65 78 20 2a 2f 0a 09 72 | row ind|ex */..r|
|000026a0| 65 67 69 73 74 65 72 20 | 69 6e 74 09 6a 3b 09 09 |egister |int.j;..|
|000026b0| 2f 2a 20 63 6f 6c 75 6d | 6e 20 69 6e 64 65 78 20 |/* colum|n index |
|000026c0| 2a 2f 0a 0a 09 77 68 69 | 6c 65 20 28 20 73 63 61 |*/...whi|le ( sca|
|000026d0| 6e 66 28 22 25 64 20 25 | 64 5c 65 6e 22 2c 20 26 |nf("%d %|d\en", &|
|000026e0| 72 2c 20 26 63 29 20 3d | 3d 20 32 20 29 09 2f 2a |r, &c) =|= 2 )./*|
|000026f0| 20 73 74 61 72 74 20 6e | 65 77 20 74 61 62 6c 65 | start n|ew table|
|00002700| 20 2a 2f 0a 09 09 7b 0a | 09 09 2f 2a 20 69 6e 70 | */...{.|../* inp|
|00002710| 75 74 20 74 61 6c 6c 69 | 65 73 20 2a 2f 0a 09 09 |ut talli|es */...|
|00002720| 66 6f 72 20 28 20 69 20 | 3d 20 30 3b 20 69 20 3c |for ( i |= 0; i <|
|00002730| 20 72 3b 20 2b 2b 69 20 | 29 0a 09 09 09 66 6f 72 | r; ++i |)....for|
|00002740| 20 28 20 6a 20 3d 20 30 | 3b 20 6a 20 3c 20 63 3b | ( j = 0|; j < c;|
|00002750| 20 2b 2b 6a 20 29 0a 09 | 09 09 09 28 76 6f 69 64 | ++j )..|...(void|
|00002760| 29 73 63 61 6e 66 28 22 | 20 25 6c 64 22 2c 20 26 |)scanf("| %ld", &|
|00002770| 78 28 69 2c 6a 29 29 3b | 0a 0a 09 09 2f 2a 20 63 |x(i,j));|..../* c|
|00002780| 6f 6d 70 75 74 65 20 73 | 74 61 74 69 73 74 69 63 |ompute s|tatistic|
|00002790| 73 20 61 6e 64 20 70 72 | 69 6e 74 20 72 65 73 75 |s and pr|int resu|
|000027a0| 6c 74 73 20 2a 2f 0a 09 | 09 43 61 6c 63 28 22 63 |lts */..|.Calc("c|
|000027b0| 68 69 73 71 22 2c 20 43 | 68 69 53 71 54 62 6c 29 |hisq", C|hiSqTbl)|
|000027c0| 3b 0a 09 09 43 61 6c 63 | 28 22 32 69 6e 66 6f 22 |;...Calc|("2info"|
|000027d0| 2c 20 49 6e 66 6f 54 62 | 6c 29 3b 0a 09 09 7d 0a |, InfoTb|l);...}.|
|000027e0| 09 72 65 74 75 72 6e 20 | 45 58 49 54 5f 53 55 43 |.return |EXIT_SUC|
|000027f0| 43 45 53 53 3b 0a 7d 5c | 66 50 0a 2e 62 70 0a 2e |CESS;.}\|fP..bp..|
|00002800| 53 48 20 46 49 4c 45 53 | 0a 2e 74 61 20 33 69 0a |SH FILES|..ta 3i.|
|00002810| 63 68 69 73 71 2e 68 09 | 68 65 61 64 65 72 20 66 |chisq.h.|header f|
|00002820| 69 6c 65 0a 2e 62 72 0a | 63 68 69 73 71 2e 6f 2c |ile..br.|chisq.o,|
|00002830| 20 69 6e 66 6f 2e 6f 09 | 6f 62 6a 65 63 74 20 6d | info.o.|object m|
|00002840| 6f 64 75 6c 65 73 0a 2e | 53 48 20 52 45 46 45 52 |odules..|SH REFER|
|00002850| 45 4e 43 45 0a 53 6f 6c | 6f 6d 6f 6e 20 4b 75 6c |ENCE.Sol|omon Kul|
|00002860| 6c 62 61 63 6b 2c 20 5c | 66 49 49 6e 66 6f 72 6d |lback, \|fIInform|
|00002870| 61 74 69 6f 6e 20 54 68 | 65 6f 72 79 20 61 6e 64 |ation Th|eory and|
|00002880| 20 53 74 61 74 69 73 74 | 69 63 73 5c 66 50 20 28 | Statist|ics\fP (|
|00002890| 44 6f 76 65 72 2c 20 31 | 39 36 38 29 2e 0a 2e 53 |Dover, 1|968)...S|
|000028a0| 48 20 22 53 45 45 20 41 | 4c 53 4f 22 0a 67 61 6d |H "SEE A|LSO".gam|
|000028b0| 6d 61 28 33 56 29 2e 0a | 2e 53 48 20 44 49 41 47 |ma(3V)..|.SH DIAG|
|000028c0| 4e 4f 53 54 49 43 53 0a | 42 6f 74 68 20 74 68 65 |NOSTICS.|Both the|
|000028d0| 73 65 20 66 75 6e 63 74 | 69 6f 6e 73 20 61 20 72 |se funct|ions a r|
|000028e0| 65 74 75 72 6e 20 6e 65 | 67 61 74 69 76 65 20 76 |eturn ne|gative v|
|000028f0| 61 6c 75 65 20 66 6f 72 | 20 74 68 65 20 73 74 61 |alue for| the sta|
|00002900| 74 69 73 74 69 63 0a 77 | 68 65 6e 20 69 74 20 63 |tistic.w|hen it c|
|00002910| 61 6e 6e 6f 74 20 62 65 | 20 6d 65 61 6e 69 6e 67 |annot be| meaning|
|00002920| 66 75 6c 6c 79 20 63 6f | 6d 70 75 74 65 64 2c 20 |fully co|mputed, |
|00002930| 61 73 20 66 6f 6c 6c 6f | 77 73 3a 0a 2e 52 53 0a |as follo|ws:..RS.|
|00002940| 2e 50 44 20 30 0a 2e 54 | 50 0a 5c 2d 31 2e 30 0a |.PD 0..T|P.\-1.0.|
|00002950| 74 6f 6f 20 6d 61 6e 79 | 20 30 20 65 6e 74 72 69 |too many| 0 entri|
|00002960| 65 73 20 69 6e 20 74 68 | 65 20 74 61 62 6c 65 0a |es in th|e table.|
|00002970| 28 66 6f 72 20 5c 66 5c | 2a 28 63 57 43 68 69 53 |(for \f\|*(cWChiS|
|00002980| 71 54 62 6c 5c 66 50 2c | 20 74 68 69 73 20 6d 65 |qTbl\fP,| this me|
|00002990| 61 6e 73 20 69 6e 73 75 | 66 66 69 63 69 65 6e 74 |ans insu|fficient|
|000029a0| 20 64 65 67 72 65 65 73 | 20 6f 66 20 66 72 65 65 | degrees| of free|
|000029b0| 64 6f 6d 3b 0a 66 6f 72 | 20 5c 66 5c 2a 28 63 57 |dom;.for| \f\*(cW|
|000029c0| 49 6e 66 6f 54 62 6c 5c | 66 50 2c 20 74 68 69 73 |InfoTbl\|fP, this|
|000029d0| 20 6d 65 61 6e 73 20 5c | 66 49 65 76 65 72 79 5c | means \|fIevery\|
|000029e0| 66 50 20 65 6e 74 72 79 | 20 69 6e 20 74 68 65 20 |fP entry| in the |
|000029f0| 74 61 62 6c 65 20 77 61 | 73 20 30 29 0a 2e 54 50 |table wa|s 0)..TP|
|00002a00| 0a 5c 2d 32 2e 30 0a 73 | 6f 6d 65 20 66 72 65 71 |.\-2.0.s|ome freq|
|00002a10| 75 65 6e 63 79 20 77 61 | 73 20 73 70 65 63 69 66 |uency wa|s specif|
|00002a20| 69 65 64 20 61 73 20 6c | 65 73 73 20 74 68 61 6e |ied as l|ess than|
|00002a30| 20 30 0a 2e 54 50 0a 5c | 2d 33 2e 30 0a 73 70 65 | 0..TP.\|-3.0.spe|
|00002a40| 63 69 66 69 65 64 20 6e | 75 6d 62 65 72 20 6f 66 |cified n|umber of|
|00002a50| 20 72 6f 77 73 20 6f 72 | 20 63 6f 6c 75 6d 6e 73 | rows or| columns|
|00002a60| 20 77 61 73 20 6c 65 73 | 73 20 74 68 61 6e 20 32 | was les|s than 2|
|00002a70| 0a 2e 54 50 0a 5c 2d 34 | 2e 30 0a 75 6e 61 62 6c |..TP.\-4|.0.unabl|
|00002a80| 65 20 74 6f 20 64 79 6e | 61 6d 69 63 61 6c 6c 79 |e to dyn|amically|
|00002a90| 20 61 6c 6c 6f 63 61 74 | 65 20 65 6e 6f 75 67 68 | allocat|e enough|
|00002aa0| 20 77 6f 72 6b 69 6e 67 | 20 73 74 6f 72 61 67 65 | working| storage|
|00002ab0| 0a 2e 50 44 0a 2e 52 45 | 0a 2e 53 48 20 41 55 54 |..PD..RE|..SH AUT|
|00002ac0| 48 4f 52 0a 44 6f 75 67 | 6c 61 73 20 41 2e 5c 26 |HOR.Doug|las A.\&|
|00002ad0| 20 47 77 79 6e 2c 20 42 | 52 4c 2f 56 4c 44 2d 56 | Gwyn, B|RL/VLD-V|
|00002ae0| 4d 42 0a 69 53 71 5c 66 | 50 20 66 75 6e 63 74 69 |MB.iSq\f|P functi|
|00002af0| 6f 6e 20 28 73 65 65 20 | 5c 66 49 67 61 6d 6d 61 |on (see |\fIgamma|
|00002b00| 5c 66 50 28 33 56 29 29 | 0a 69 74 20 63 61 6e 20 |\fP(3V))|.it can |
|00002b10| 62 65 20 63 6f 6e 76 65 | 72 74 65 64 20 74 6f 20 |be conve|rted to |
|00002b20| 74 68 65 20 70 72 6f 62 | 61 62 69 6c 69 74 79 20 |the prob|ability |
|00002b30| 74 68 61 74 0a 74 68 65 | 20 73 74 61 74 69 73 74 |that.the| statist|
|00002b40| 69 63 20 77 6f 75 6c 64 | 20 62 65 20 61 73 20 6c |ic would| be as l|
|00002b50| 61 72 67 65 20 61 73 20 | 77 61 73 20 6f 62 73 65 |arge as |was obse|
|00002b60| 72 76 65 64 0a 69 66 20 | 74 68 65 20 63 61 74 65 |rved.if |the cate|
|00002b70| 67 6f 72 69 7a 61 74 69 | 6f 6e 73 20 72 65 61 6c |gorizati|ons real|
|00002b80| 6c 79 20 61 72 65 20 69 | 6e 64 65 70 65 6e 64 65 |ly are i|ndepende|
|00002b90| 6e 74 2e 0a 54 68 69 73 | 20 69 73 20 6f 66 20 63 |nt..This| is of c|
|00002ba0| 6f 75 72 73 65 20 74 68 | 65 20 74 72 61 64 69 74 |ourse th|e tradit|
|00002bb0| 69 6f 6e 61 6c 20 75 73 | 65 20 6f 66 20 50 65 61 |ional us|e of Pea|
|00002bc0| 72 73 6f 6e 27 73 20 40 | 63 68 69 20 73 75 70 20 |rson's @|chi sup |
|00002bd0| 32 40 20 73 74 61 74 69 | 73 74 69 63 2c 0a 77 68 |2@ stati|stic,.wh|
|00002be0| 69 63 68 20 40 32 20 22 | 5c 76 27 2d 30 2e 32 27 |ich @2 "|\v'-0.2'|
|00002bf0| 5e 5c 76 27 30 2e 32 27 | 5c 68 27 2d 5c 77 3b 5e |^\v'0.2'|\h'-\w;^|
|00002c00| 74 6f 74 2d 69 6e 66 6f | 2f 63 68 69 73 71 2e 63 |tot-info|/chisq.c|
|00002c10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c60| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00002c70| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00002c80| 20 20 20 34 37 31 34 20 | 20 34 34 35 37 32 31 30 |   4714 | 4457210|
|00002c90| 34 36 32 20 20 20 37 31 | 30 32 00 20 00 00 00 00 |462   71|02. ....|
|00002ca0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ce0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cf0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d00| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002da0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002db0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002dc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002dd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002de0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002df0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e00| 2f 2a 0a 09 43 68 69 53 | 71 54 62 6c 20 2d 2d 20 |/*..ChiS|qTbl -- |
|00002e10| 50 65 61 72 73 6f 6e 27 | 73 20 63 68 69 2d 73 71 |Pearson'|s chi-sq|
|00002e20| 75 61 72 65 20 74 65 73 | 74 20 66 6f 72 20 61 20 |uare tes|t for a |
|00002e30| 32 2d 77 61 79 20 63 6f | 6e 74 69 6e 67 65 6e 63 |2-way co|ntingenc|
|00002e40| 79 20 74 61 62 6c 65 0a | 0a 09 6c 61 73 74 20 65 |y table.|..last e|
|00002e50| 64 69 74 3a 09 38 38 2f | 30 39 2f 31 39 09 44 20 |dit:.88/|09/19.D |
|00002e60| 41 20 47 77 79 6e 0a 0a | 09 53 43 43 53 20 49 44 |A Gwyn..|.SCCS ID|
|00002e70| 3a 09 40 28 23 29 63 68 | 69 73 71 2e 63 09 31 2e |:.@(#)ch|isq.c.1.|
|00002e80| 31 20 28 65 64 69 74 65 | 64 20 66 6f 72 20 70 75 |1 (edite|d for pu|
|00002e90| 62 6c 69 63 61 74 69 6f | 6e 29 0a 0a 09 53 70 65 |blicatio|n)...Spe|
|00002ea0| 63 69 61 6c 20 72 65 74 | 75 72 6e 20 76 61 6c 75 |cial ret|urn valu|
|00002eb0| 65 73 3a 0a 09 09 2d 31 | 2e 30 09 69 6e 73 75 66 |es:...-1|.0.insuf|
|00002ec0| 66 69 63 69 65 6e 74 20 | 64 65 67 72 65 65 73 20 |ficient |degrees |
|00002ed0| 6f 66 20 66 72 65 65 64 | 6f 6d 20 28 74 6f 6f 20 |of freed|om (too |
|00002ee0| 6d 61 6e 79 20 30 20 65 | 6e 74 72 69 65 73 29 0a |many 0 e|ntries).|
|00002ef0| 09 09 2d 32 2e 30 09 69 | 6e 76 61 6c 69 64 20 74 |..-2.0.i|nvalid t|
|00002f00| 61 62 6c 65 20 65 6e 74 | 72 79 20 28 66 72 65 71 |able ent|ry (freq|
|00002f10| 75 65 6e 63 79 20 6c 65 | 73 73 20 74 68 61 6e 20 |uency le|ss than |
|00002f20| 30 29 0a 09 09 2d 33 2e | 30 09 69 6e 76 61 6c 69 |0)...-3.|0.invali|
|00002f30| 64 20 74 61 62 6c 65 20 | 64 69 6d 65 6e 73 69 6f |d table |dimensio|
|00002f40| 6e 73 20 28 72 20 6f 72 | 20 63 20 6c 65 73 73 20 |ns (r or| c less |
|00002f50| 74 68 61 6e 20 32 29 0a | 09 09 2d 34 2e 30 09 75 |than 2).|..-4.0.u|
|00002f60| 6e 61 62 6c 65 20 74 6f | 20 61 6c 6c 6f 63 61 74 |nable to| allocat|
|00002f70| 65 20 65 6e 6f 75 67 68 | 20 77 6f 72 6b 69 6e 67 |e enough| working|
|00002f80| 20 73 74 6f 72 61 67 65 | 0a 2a 2f 0a 0a 23 69 66 | storage|.*/..#if|
|00002f90| 64 65 66 20 5f 5f 53 54 | 44 43 5f 5f 0a 23 69 6e |def __ST|DC__.#in|
|00002fa0| 63 6c 75 64 65 09 3c 73 | 74 64 6c 69 62 2e 68 3e |clude.<s|tdlib.h>|
|00002fb0| 09 09 2f 2a 20 6d 61 6c | 6c 6f 63 2c 20 66 72 65 |../* mal|loc, fre|
|00002fc0| 65 20 2a 2f 0a 0a 23 69 | 6e 63 6c 75 64 65 09 22 |e */..#i|nclude."|
|00002fd0| 73 74 64 2e 68 22 0a 23 | 65 6c 73 65 0a 23 69 6e |std.h".#|else.#in|
|00002fe0| 63 6c 75 64 65 09 22 73 | 74 64 2e 68 22 0a 0a 65 |clude."s|td.h"..e|
|00002ff0| 78 74 65 72 6e 20 70 6f | 69 6e 74 65 72 09 6d 61 |xtern po|inter.ma|
|00003000| 6c 6c 6f 63 28 29 3b 0a | 65 78 74 65 72 6e 20 76 |lloc();.|extern v|
|00003010| 6f 69 64 09 66 72 65 65 | 28 29 3b 0a 23 65 6e 64 |oid.free|();.#end|
|00003020| 69 66 0a 0a 64 6f 75 62 | 6c 65 0a 43 68 69 53 71 |if..doub|le.ChiSq|
|00003030| 54 62 6c 28 20 72 2c 20 | 63 2c 20 66 2c 20 70 64 |Tbl( r, |c, f, pd|
|00003040| 66 20 29 0a 09 69 6e 74 | 09 09 72 3b 09 09 2f 2a |f )..int|..r;../*|
|00003050| 20 23 20 72 6f 77 73 20 | 69 6e 20 74 61 62 6c 65 | # rows |in table|
|00003060| 20 2a 2f 0a 09 69 6e 74 | 09 09 63 3b 09 09 2f 2a | */..int|..c;../*|
|00003070| 20 23 20 63 6f 6c 75 6d | 6e 73 20 69 6e 20 74 61 | # colum|ns in ta|
|00003080| 62 6c 65 20 2a 2f 0a 09 | 63 6f 6e 73 74 20 6c 6f |ble */..|const lo|
|00003090| 6e 67 09 2a 66 3b 09 09 | 2f 2a 20 2d 3e 20 72 2a |ng.*f;..|/* -> r*|
|000030a0| 63 20 66 72 65 71 75 65 | 6e 63 79 20 74 61 6c 6c |c freque|ncy tall|
|000030b0| 69 65 73 20 2a 2f 0a 09 | 69 6e 74 09 09 2a 70 64 |ies */..|int..*pd|
|000030c0| 66 3b 09 09 2f 2a 20 2d | 3e 20 72 65 74 75 72 6e |f;../* -|> return|
|000030d0| 20 23 20 64 65 67 72 65 | 65 73 20 6f 66 20 66 72 | # degre|es of fr|
|000030e0| 65 65 64 6f 6d 20 2a 2f | 0a 09 7b 0a 23 64 65 66 |eedom */|..{.#def|
|000030f0| 69 6e 65 09 78 28 69 2c | 6a 29 09 66 5b 28 69 29 |ine.x(i,|j).f[(i)|
|00003100| 2a 63 2b 28 6a 29 5d 09 | 09 2f 2a 20 63 6f 6e 76 |*c+(j)].|./* conv|
|00003110| 65 6e 69 65 6e 74 20 77 | 61 79 20 74 6f 20 61 63 |enient w|ay to ac|
|00003120| 63 65 73 73 20 66 72 65 | 71 73 20 2a 2f 0a 09 72 |cess fre|qs */..r|
|00003130| 65 67 69 73 74 65 72 20 | 69 6e 74 09 69 3b 09 09 |egister |int.i;..|
|00003140| 2f 2a 20 72 6f 77 20 69 | 6e 64 65 78 20 2a 2f 0a |/* row i|ndex */.|
|00003150| 09 72 65 67 69 73 74 65 | 72 20 69 6e 74 09 6a 3b |.registe|r int.j;|
|00003160| 09 09 2f 2a 20 63 6f 6c | 75 6d 6e 20 69 6e 64 65 |../* col|umn inde|
|00003170| 78 20 2a 2f 0a 09 64 6f | 75 62 6c 65 09 09 4e 3b |x */..do|uble..N;|
|00003180| 09 09 2f 2a 20 28 64 6f | 75 62 6c 65 29 6e 20 2a |../* (do|uble)n *|
|00003190| 2f 0a 09 64 6f 75 62 6c | 65 09 09 63 68 69 73 71 |/..doubl|e..chisq|
|000031a0| 3b 09 09 2f 2a 20 61 63 | 63 75 6d 75 6c 61 74 65 |;../* ac|cumulate|
|000031b0| 73 20 63 68 69 2d 73 71 | 75 61 72 65 20 2a 2f 0a |s chi-sq|uare */.|
|000031c0| 09 64 6f 75 62 6c 65 09 | 09 2a 78 69 3b 09 09 2f |.double.|.*xi;../|
|000031d0| 2a 20 72 6f 77 20 73 75 | 6d 73 20 2a 2f 0a 09 64 |* row su|ms */..d|
|000031e0| 6f 75 62 6c 65 09 09 2a | 78 6a 3b 09 09 2f 2a 20 |ouble..*|xj;../* |
|000031f0| 63 6f 6c 20 73 75 6d 73 | 20 2a 2f 0a 09 69 6e 74 |col sums| */..int|
|00003200| 09 09 72 64 66 20 3d 20 | 72 20 2d 20 31 3b 09 2f |..rdf = |r - 1;./|
|00003210| 2a 20 72 6f 77 20 64 65 | 67 72 65 65 73 20 6f 66 |* row de|grees of|
|00003220| 20 66 72 65 65 64 6f 6d | 20 2a 2f 0a 09 69 6e 74 | freedom| */..int|
|00003230| 09 09 63 64 66 20 3d 20 | 63 20 2d 20 31 3b 09 2f |..cdf = |c - 1;./|
|00003240| 2a 20 63 6f 6c 75 6d 6e | 20 64 65 67 72 65 65 73 |* column| degrees|
|00003250| 20 6f 66 20 66 72 65 65 | 64 6f 6d 20 2a 2f 0a 0a | of free|dom */..|
|00003260| 09 69 66 20 28 20 72 64 | 66 20 3c 3d 20 30 20 7c |.if ( rd|f <= 0 ||
|00003270| 7c 20 63 64 66 20 3c 3d | 20 30 20 29 0a 09 09 7b || cdf <=| 0 )...{|
|00003280| 0a 09 09 63 68 69 73 71 | 20 3d 20 2d 33 2e 30 3b |...chisq| = -3.0;|
|00003290| 0a 09 09 67 6f 74 6f 20 | 72 65 74 33 3b 0a 09 09 |...goto |ret3;...|
|000032a0| 7d 0a 0a 09 69 66 20 28 | 20 28 78 69 20 3d 20 28 |}...if (| (xi = (|
|000032b0| 64 6f 75 62 6c 65 20 2a | 29 6d 61 6c 6c 6f 63 28 |double *|)malloc(|
|000032c0| 20 72 20 2a 20 73 69 7a | 65 6f 66 28 64 6f 75 62 | r * siz|eof(doub|
|000032d0| 6c 65 29 20 29 29 20 3d | 3d 20 4e 55 4c 4c 20 29 |le) )) =|= NULL )|
|000032e0| 0a 09 09 7b 0a 09 09 63 | 68 69 73 71 20 3d 20 2d |...{...c|hisq = -|
|000032f0| 34 2e 30 3b 0a 09 09 67 | 6f 74 6f 20 72 65 74 33 |4.0;...g|oto ret3|
|00003300| 3b 0a 09 09 7d 0a 0a 09 | 69 66 20 28 20 28 78 6a |;...}...|if ( (xj|
|00003310| 20 3d 20 28 64 6f 75 62 | 6c 65 20 2a 29 6d 61 6c | = (doub|le *)mal|
|00003320| 6c 6f 63 28 20 63 20 2a | 20 73 69 7a 65 6f 66 28 |loc( c *| sizeof(|
|00003330| 64 6f 75 62 6c 65 29 20 | 29 29 20 3d 3d 20 4e 55 |double) |)) == NU|
|00003340| 4c 4c 20 29 0a 09 09 7b | 0a 09 09 63 68 69 73 71 |LL )...{|...chisq|
|00003350| 20 3d 20 2d 34 2e 30 3b | 0a 09 09 67 6f 74 6f 20 | = -4.0;|...goto |
|00003360| 72 65 74 32 3b 0a 09 09 | 7d 0a 0a 09 2f 2a 20 63 |ret2;...|}.../* c|
|00003370| 6f 6d 70 75 74 65 20 72 | 6f 77 20 73 75 6d 73 20 |ompute r|ow sums |
|00003380| 61 6e 64 20 74 6f 74 61 | 6c 20 2a 2f 0a 0a 09 4e |and tota|l */...N|
|00003390| 20 3d 20 30 2e 30 3b 0a | 0a 09 66 6f 72 20 28 20 | = 0.0;.|..for ( |
|000033a0| 69 20 3d 20 30 3b 20 69 | 20 3c 20 72 3b 20 2b 2b |i = 0; i| < r; ++|
|000033b0| 69 20 29 0a 09 09 7b 0a | 09 09 64 6f 75 62 6c 65 |i )...{.|..double|
|000033c0| 09 73 75 6d 20 3d 20 30 | 2e 30 3b 09 2f 2a 20 61 |.sum = 0|.0;./* a|
|000033d0| 63 63 75 6d 75 6c 61 74 | 6f 72 20 2a 2f 0a 0a 09 |ccumulat|or */...|
|000033e0| 09 66 6f 72 20 28 20 6a | 20 3d 20 30 3b 20 6a 20 |.for ( j| = 0; j |
|000033f0| 3c 20 63 3b 20 2b 2b 6a | 20 29 0a 09 09 09 7b 0a |< c; ++j| )....{.|
|00003400| 09 09 09 72 65 67 69 73 | 74 65 72 20 6c 6f 6e 67 |...regis|ter long|
|00003410| 09 6b 20 3d 20 78 28 69 | 2c 6a 29 3b 0a 0a 09 09 |.k = x(i|,j);....|
|00003420| 09 69 66 20 28 20 6b 20 | 3c 20 30 4c 20 29 0a 09 |.if ( k |< 0L )..|
|00003430| 09 09 09 7b 0a 09 09 09 | 09 63 68 69 73 71 20 3d |...{....|.chisq =|
|00003440| 20 2d 32 2e 30 3b 0a 09 | 09 09 09 67 6f 74 6f 20 | -2.0;..|...goto |
|00003450| 72 65 74 31 3b 0a 09 09 | 09 09 7d 0a 0a 09 09 09 |ret1;...|..}.....|
|00003460| 73 75 6d 20 2b 3d 20 28 | 64 6f 75 62 6c 65 29 6b |sum += (|double)k|
|00003470| 3b 0a 09 09 09 7d 0a 0a | 09 09 69 66 20 28 20 28 |;....}..|..if ( (|
|00003480| 78 69 5b 69 5d 20 3d 20 | 73 75 6d 29 20 3c 3d 20 |xi[i] = |sum) <= |
|00003490| 30 2e 30 20 29 0a 09 09 | 09 2d 2d 72 64 66 3b 0a |0.0 )...|.--rdf;.|
|000034a0| 09 09 65 6c 73 65 0a 09 | 09 09 4e 20 2b 3d 20 73 |..else..|..N += s|
|000034b0| 75 6d 3b 0a 09 09 7d 0a | 0a 09 2f 2a 20 63 6f 6d |um;...}.|../* com|
|000034c0| 70 75 74 65 20 63 6f 6c | 75 6d 6e 20 73 75 6d 73 |pute col|umn sums|
|000034d0| 20 2a 2f 0a 0a 09 66 6f | 72 20 28 20 6a 20 3d 20 | */...fo|r ( j = |
|000034e0| 30 3b 20 6a 20 3c 20 63 | 3b 20 2b 2b 6a 20 29 0a |0; j < c|; ++j ).|
|000034f0| 09 09 7b 0a 09 09 64 6f | 75 62 6c 65 09 73 75 6d |..{...do|uble.sum|
|00003500| 20 3d 20 30 2e 30 3b 09 | 2f 2a 20 61 63 63 75 6d | = 0.0;.|/* accum|
|00003510| 75 6c 61 74 6f 72 20 2a | 2f 0a 0a 09 09 66 6f 72 |ulator *|/....for|
|00003520| 20 28 20 69 20 3d 20 30 | 3b 20 69 20 3c 20 72 3b | ( i = 0|; i < r;|
|00003530| 20 2b 2b 69 20 29 0a 09 | 09 09 73 75 6d 20 2b 3d | ++i )..|..sum +=|
|00003540| 20 28 64 6f 75 62 6c 65 | 29 78 28 69 2c 6a 29 3b | (double|)x(i,j);|
|00003550| 0a 0a 09 09 69 66 20 28 | 20 28 78 6a 5b 6a 5d 20 |....if (| (xj[j] |
|00003560| 3d 20 73 75 6d 29 20 3c | 3d 20 30 2e 30 20 29 0a |= sum) <|= 0.0 ).|
|00003570| 09 09 09 2d 2d 63 64 66 | 3b 0a 09 09 7d 0a 0a 09 |...--cdf|;...}...|
|00003580| 69 66 20 28 20 72 64 66 | 20 3c 3d 20 30 20 7c 7c |if ( rdf| <= 0 |||
|00003590| 20 63 64 66 20 3c 3d 20 | 30 20 29 0a 09 09 7b 0a | cdf <= |0 )...{.|
|000035a0| 09 09 63 68 69 73 71 20 | 3d 20 2d 31 2e 30 3b 0a |..chisq |= -1.0;.|
|000035b0| 09 09 67 6f 74 6f 20 72 | 65 74 31 3b 0a 09 09 7d |..goto r|et1;...}|
|000035c0| 0a 0a 09 2a 70 64 66 20 | 3d 20 72 64 66 20 2a 20 |...*pdf |= rdf * |
|000035d0| 63 64 66 3b 09 09 2f 2a | 20 74 6f 74 61 6c 20 64 |cdf;../*| total d|
|000035e0| 65 67 72 65 65 73 20 6f | 66 20 66 72 65 65 64 6f |egrees o|f freedo|
|000035f0| 6d 20 2a 2f 0a 0a 09 2f | 2a 20 63 6f 6d 70 75 74 |m */.../|* comput|
|00003600| 65 20 63 68 69 2d 73 71 | 75 61 72 65 20 2a 2f 0a |e chi-sq|uare */.|
|00003610| 0a 09 63 68 69 73 71 20 | 3d 20 30 2e 30 3b 0a 0a |..chisq |= 0.0;..|
|00003620| 09 66 6f 72 20 28 20 69 | 20 3d 20 30 3b 20 69 20 |.for ( i| = 0; i |
|00003630| 3c 20 72 3b 20 2b 2b 69 | 20 29 0a 09 09 7b 0a 09 |< r; ++i| )...{..|
|00003640| 09 64 6f 75 62 6c 65 09 | 70 69 20 3d 20 78 69 5b |.double.|pi = xi[|
|00003650| 69 5d 3b 09 2f 2a 20 72 | 6f 77 20 73 75 6d 20 2a |i];./* r|ow sum *|
|00003660| 2f 0a 0a 09 09 69 66 20 | 28 20 70 69 20 3e 20 30 |/....if |( pi > 0|
|00003670| 2e 30 20 29 0a 09 09 09 | 66 6f 72 20 28 20 6a 20 |.0 )....|for ( j |
|00003680| 3d 20 30 3b 20 6a 20 3c | 20 63 3b 20 2b 2b 6a 20 |= 0; j <| c; ++j |
|00003690| 29 0a 09 09 09 09 7b 0a | 09 09 09 09 64 6f 75 62 |).....{.|....doub|
|000036a0| 6c 65 09 70 6a 20 3d 20 | 78 6a 5b 6a 5d 3b 09 2f |le.pj = |xj[j];./|
|000036b0| 2a 20 63 6f 6c 75 6d 6e | 20 73 75 6d 20 2a 2f 0a |* column| sum */.|
|000036c0| 0a 09 09 09 09 69 66 20 | 28 20 70 6a 20 3e 20 30 |.....if |( pj > 0|
|000036d0| 2e 30 20 29 0a 09 09 09 | 09 09 7b 0a 09 09 09 09 |.0 )....|..{.....|
|000036e0| 09 64 6f 75 62 6c 65 09 | 65 78 70 65 63 74 65 64 |.double.|expected|
|000036f0| 20 3d 20 70 69 20 2a 20 | 70 6a 20 2f 20 4e 3b 0a | = pi * |pj / N;.|
|00003700| 09 09 09 09 09 64 6f 75 | 62 6c 65 09 64 65 6c 74 |.....dou|ble.delt|
|00003710| 61 20 3d 0a 09 09 09 09 | 09 09 28 64 6f 75 62 6c |a =.....|..(doubl|
|00003720| 65 29 78 28 69 2c 6a 29 | 20 2d 20 65 78 70 65 63 |e)x(i,j)| - expec|
|00003730| 74 65 64 3b 0a 0a 09 09 | 09 09 09 63 68 69 73 71 |ted;....|...chisq|
|00003740| 20 2b 3d 20 64 65 6c 74 | 61 20 2a 20 64 65 6c 74 | += delt|a * delt|
|00003750| 61 20 2f 20 65 78 70 65 | 63 74 65 64 3b 0a 09 09 |a / expe|cted;...|
|00003760| 09 09 09 7d 0a 09 09 09 | 09 7d 0a 09 09 7d 0a 0a |...}....|.}...}..|
|00003770| 20 20 20 20 72 65 74 31 | 3a 0a 09 66 72 65 65 28 |    ret1|:..free(|
|00003780| 20 28 70 6f 69 6e 74 65 | 72 29 78 6a 20 29 3b 0a | (pointe|r)xj );.|
|00003790| 20 20 20 20 72 65 74 32 | 3a 0a 09 66 72 65 65 28 |    ret2|:..free(|
|000037a0| 20 28 70 6f 69 6e 74 65 | 72 29 78 69 20 29 3b 0a | (pointe|r)xi );.|
|000037b0| 20 20 20 20 72 65 74 33 | 3a 0a 09 72 65 74 75 72 |    ret3|:..retur|
|000037c0| 6e 20 63 68 69 73 71 3b | 0a 09 7d 0a 74 68 69 73 |n chisq;|..}.this|
|000037d0| 20 6d 65 61 6e 73 20 5c | 66 49 65 76 65 72 79 5c | means \|fIevery\|
|000037e0| 66 50 20 65 6e 74 72 79 | 20 69 6e 20 74 68 65 20 |fP entry| in the |
|000037f0| 74 61 62 6c 65 20 77 61 | 73 20 30 29 0a 2e 54 50 |table wa|s 0)..TP|
|00003800| 74 6f 74 2d 69 6e 66 6f | 2f 63 68 69 73 71 2e 68 |tot-info|/chisq.h|
|00003810| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003820| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003830| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003840| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003850| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003860| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00003870| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00003880| 20 20 20 20 35 37 33 20 | 20 34 34 35 37 32 31 30 |    573 | 4457210|
|00003890| 34 36 33 20 20 20 37 30 | 36 37 00 20 00 00 00 00 |463   70|67. ....|
|000038a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000038b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000038c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000038d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000038e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000038f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003900| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003910| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003920| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003930| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003940| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003950| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003960| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003970| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003980| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003990| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000039f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003a00| 2f 2a 0a 09 63 68 69 73 | 71 2e 68 20 2d 2d 20 64 |/*..chis|q.h -- d|
|00003a10| 65 66 69 6e 69 74 69 6f | 6e 73 20 66 6f 72 20 63 |efinitio|ns for c|
|00003a20| 6f 6e 74 69 6e 67 65 6e | 63 79 2d 74 61 62 6c 65 |ontingen|cy-table|
|00003a30| 20 72 6f 75 74 69 6e 65 | 73 0a 0a 09 6c 61 73 74 | routine|s...last|
|00003a40| 20 65 64 69 74 3a 09 38 | 38 2f 30 39 2f 31 39 09 | edit:.8|8/09/19.|
|00003a50| 44 20 41 20 47 77 79 6e | 0a 0a 09 53 43 43 53 20 |D A Gwyn|...SCCS |
|00003a60| 49 44 3a 09 40 28 23 29 | 63 68 69 73 71 2e 68 09 |ID:.@(#)|chisq.h.|
|00003a70| 31 2e 31 20 28 65 64 69 | 74 65 64 20 66 6f 72 20 |1.1 (edi|ted for |
|00003a80| 70 75 62 6c 69 63 61 74 | 69 6f 6e 29 0a 2a 2f 0a |publicat|ion).*/.|
|00003a90| 0a 2f 2a 20 6c 69 62 72 | 61 72 79 20 72 6f 75 74 |./* libr|ary rout|
|00003aa0| 69 6e 65 73 3a 20 2a 2f | 0a 0a 23 69 66 64 65 66 |ines: */|..#ifdef|
|00003ab0| 20 5f 5f 53 54 44 43 5f | 5f 0a 65 78 74 65 72 6e | __STDC_|_.extern|
|00003ac0| 20 64 6f 75 62 6c 65 09 | 43 68 69 53 71 54 62 6c | double.|ChiSqTbl|
|00003ad0| 28 20 69 6e 74 20 72 2c | 20 69 6e 74 20 63 2c 20 |( int r,| int c, |
|00003ae0| 63 6f 6e 73 74 20 6c 6f | 6e 67 20 2a 66 2c 20 69 |const lo|ng *f, i|
|00003af0| 6e 74 20 2a 70 64 66 20 | 29 3b 0a 65 78 74 65 72 |nt *pdf |);.exter|
|00003b00| 6e 20 64 6f 75 62 6c 65 | 09 49 6e 66 6f 54 62 6c |n double|.InfoTbl|
|00003b10| 28 20 69 6e 74 20 72 2c | 20 69 6e 74 20 63 2c 20 |( int r,| int c, |
|00003b20| 63 6f 6e 73 74 20 6c 6f | 6e 67 20 2a 66 2c 20 69 |const lo|ng *f, i|
|00003b30| 6e 74 20 2a 70 64 66 20 | 29 3b 0a 23 65 6c 73 65 |nt *pdf |);.#else|
|00003b40| 0a 65 78 74 65 72 6e 20 | 64 6f 75 62 6c 65 09 43 |.extern |double.C|
|00003b50| 68 69 53 71 54 62 6c 28 | 29 3b 0a 65 78 74 65 72 |hiSqTbl(|);.exter|
|00003b60| 6e 20 64 6f 75 62 6c 65 | 09 49 6e 66 6f 54 62 6c |n double|.InfoTbl|
|00003b70| 28 29 3b 0a 23 65 6e 64 | 69 66 0a 72 6b 69 6e 67 |();.#end|if.rking|
|00003b80| 20 73 74 6f 72 61 67 65 | 0a 2a 2f 0a 0a 23 69 66 | storage|.*/..#if|
|00003b90| 64 65 66 20 5f 5f 53 54 | 44 43 5f 5f 0a 23 69 6e |def __ST|DC__.#in|
|00003ba0| 63 6c 75 64 65 09 3c 73 | 74 64 6c 69 62 2e 68 3e |clude.<s|tdlib.h>|
|00003bb0| 09 09 2f 2a 20 6d 61 6c | 6c 6f 63 2c 20 66 72 65 |../* mal|loc, fre|
|00003bc0| 65 20 2a 2f 0a 0a 23 69 | 6e 63 6c 75 64 65 09 22 |e */..#i|nclude."|
|00003bd0| 73 74 64 2e 68 22 0a 23 | 65 6c 73 65 0a 23 69 6e |std.h".#|else.#in|
|00003be0| 63 6c 75 64 65 09 22 73 | 74 64 2e 68 22 0a 0a 65 |clude."s|td.h"..e|
|00003bf0| 78 74 65 72 6e 20 70 6f | 69 6e 74 65 72 09 6d 61 |xtern po|inter.ma|
|00003c00| 74 6f 74 2d 69 6e 66 6f | 2f 67 5f 74 65 73 74 2e |tot-info|/g_test.|
|00003c10| 63 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |c.......|........|
|00003c20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003c30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003c40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003c50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003c60| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00003c70| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00003c80| 20 20 31 30 37 34 33 20 | 20 34 34 35 37 32 31 30 |  10743 | 4457210|
|00003c90| 34 36 36 20 20 20 37 33 | 30 33 00 20 00 00 00 00 |466   73|03. ....|
|00003ca0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003cb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003cc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003cd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003ce0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003cf0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d00| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003d90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003da0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003db0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003dc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003dd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003de0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003df0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003e00| 2f 2a 0a 09 67 5f 74 65 | 73 74 20 2d 2d 20 74 65 |/*..g_te|st -- te|
|00003e10| 73 74 73 20 66 6f 72 20 | 67 61 6d 6d 61 20 61 6e |sts for |gamma an|
|00003e20| 64 20 72 65 6c 61 74 65 | 64 20 66 75 6e 63 74 69 |d relate|d functi|
|00003e30| 6f 6e 73 0a 0a 09 6c 61 | 73 74 20 65 64 69 74 3a |ons...la|st edit:|
|00003e40| 09 38 38 2f 30 39 2f 30 | 39 09 44 20 41 20 47 77 |.88/09/0|9.D A Gw|
|00003e50| 79 6e 0a 09 53 43 43 53 | 20 49 44 3a 09 40 28 23 |yn..SCCS| ID:.@(#|
|00003e60| 29 67 5f 74 65 73 74 2e | 63 09 31 2e 31 20 28 65 |)g_test.|c.1.1 (e|
|00003e70| 64 69 74 65 64 20 66 6f | 72 20 70 75 62 6c 69 63 |dited fo|r public|
|00003e80| 61 74 69 6f 6e 29 0a 2a | 2f 0a 0a 23 69 6e 63 6c |ation).*|/..#incl|
|00003e90| 75 64 65 09 3c 73 74 64 | 69 6f 2e 68 3e 0a 23 69 |ude.<std|io.h>.#i|
|00003ea0| 66 64 65 66 20 5f 5f 53 | 54 44 43 5f 5f 0a 23 69 |fdef __S|TDC__.#i|
|00003eb0| 6e 63 6c 75 64 65 09 3c | 73 74 64 6c 69 62 2e 68 |nclude.<|stdlib.h|
|00003ec0| 3e 09 09 2f 2a 20 66 6f | 72 20 4e 55 4c 4c 20 65 |>../* fo|r NULL e|
|00003ed0| 74 63 2e 20 2a 2f 0a 23 | 65 6e 64 69 66 0a 0a 23 |tc. */.#|endif..#|
|00003ee0| 69 6e 63 6c 75 64 65 09 | 22 73 74 64 2e 68 22 0a |include.|"std.h".|
|00003ef0| 0a 23 69 6e 63 6c 75 64 | 65 09 22 67 61 6d 6d 61 |.#includ|e."gamma|
|00003f00| 2e 68 22 0a 0a 23 64 65 | 66 69 6e 65 09 50 72 69 |.h"..#de|fine.Pri|
|00003f10| 6e 74 66 09 28 76 6f 69 | 64 29 70 72 69 6e 74 66 |ntf.(voi|d)printf|
|00003f20| 0a 0a 23 64 65 66 69 6e | 65 09 54 4f 4c 09 31 2e |..#defin|e.TOL.1.|
|00003f30| 30 65 2d 35 09 09 09 2f | 2a 20 74 6f 6c 65 72 61 |0e-5.../|* tolera|
|00003f40| 6e 63 65 20 66 6f 72 20 | 63 68 65 63 6b 73 20 2a |nce for |checks *|
|00003f50| 2f 0a 0a 73 74 61 74 69 | 63 20 69 6e 74 09 65 72 |/..stati|c int.er|
|00003f60| 72 73 20 3d 20 30 3b 09 | 09 2f 2a 20 74 61 6c 6c |rs = 0;.|./* tall|
|00003f70| 79 20 65 72 72 6f 72 73 | 20 2a 2f 0a 0a 73 74 61 |y errors| */..sta|
|00003f80| 74 69 63 20 64 6f 75 62 | 6c 65 0a 52 65 6c 44 69 |tic doub|le.RelDi|
|00003f90| 66 28 20 61 2c 20 62 20 | 29 09 09 09 2f 2a 20 72 |f( a, b |).../* r|
|00003fa0| 65 74 75 72 6e 73 20 72 | 65 6c 61 74 69 76 65 20 |eturns r|elative |
|00003fb0| 64 69 66 66 65 72 65 6e | 63 65 3a 09 2a 2f 0a 09 |differen|ce:.*/..|
|00003fc0| 64 6f 75 62 6c 65 09 61 | 2c 20 62 3b 09 09 2f 2a |double.a|, b;../*|
|00003fd0| 20 30 2e 30 20 69 66 20 | 65 78 61 63 74 6c 79 20 | 0.0 if |exactly |
|00003fe0| 74 68 65 20 73 61 6d 65 | 2c 0a 09 09 09 09 20 20 |the same|,.....  |
|00003ff0| 20 6f 74 68 65 72 77 69 | 73 65 20 72 61 74 69 6f | otherwi|se ratio|
|00004000| 20 6f 66 20 64 69 66 66 | 65 72 65 6e 63 65 0a 09 | of diff|erence..|
|00004010| 09 09 09 20 20 20 74 6f | 20 74 68 65 20 6c 61 72 |...   to| the lar|
|00004020| 67 65 72 20 6f 66 20 74 | 68 65 20 74 77 6f 09 2a |ger of t|he two.*|
|00004030| 2f 0a 09 7b 0a 09 64 6f | 75 62 6c 65 09 63 20 3d |/..{..do|uble.c =|
|00004040| 20 41 62 73 28 20 61 20 | 29 3b 0a 09 64 6f 75 62 | Abs( a |);..doub|
|00004050| 6c 65 09 64 20 3d 20 41 | 62 73 28 20 62 20 29 3b |le.d = A|bs( b );|
|00004060| 0a 0a 09 64 20 3d 20 4d | 61 78 28 20 63 2c 20 64 |...d = M|ax( c, d|
|00004070| 20 29 3b 0a 0a 09 72 65 | 74 75 72 6e 20 64 20 3d | );...re|turn d =|
|00004080| 3d 20 30 2e 30 20 3f 20 | 30 2e 30 20 3a 20 41 62 |= 0.0 ? |0.0 : Ab|
|00004090| 73 28 20 61 20 2d 20 62 | 20 29 20 2f 20 64 3b 0a |s( a - b| ) / d;.|
|000040a0| 09 7d 0a 0a 73 74 61 74 | 69 63 20 76 6f 69 64 0a |.}..stat|ic void.|
|000040b0| 52 43 68 65 63 6b 28 20 | 64 2c 20 72 20 29 09 09 |RCheck( |d, r )..|
|000040c0| 09 09 2f 2a 20 63 68 65 | 63 6b 20 72 65 61 6c 20 |../* che|ck real |
|000040d0| 6e 75 6d 62 65 72 20 2a | 2f 0a 09 64 6f 75 62 6c |number *|/..doubl|
|000040e0| 65 09 64 3b 09 09 09 2f | 2a 20 72 65 61 6c 20 74 |e.d;.../|* real t|
|000040f0| 6f 20 62 65 20 63 68 65 | 63 6b 65 64 20 2a 2f 0a |o be che|cked */.|
|00004100| 09 64 6f 75 62 6c 65 09 | 72 3b 09 09 09 2f 2a 20 |.double.|r;.../* |
|00004110| 65 78 70 65 63 74 65 64 | 20 76 61 6c 75 65 20 2a |expected| value *|
|00004120| 2f 0a 09 7b 0a 09 69 66 | 20 28 20 52 65 6c 44 69 |/..{..if| ( RelDi|
|00004130| 66 28 20 64 2c 20 72 20 | 29 20 3e 20 54 4f 4c 20 |f( d, r |) > TOL |
|00004140| 29 0a 09 09 7b 0a 09 09 | 65 72 72 73 20 3d 20 31 |)...{...|errs = 1|
|00004150| 3b 0a 09 09 50 72 69 6e | 74 66 28 20 22 76 61 6c |;...Prin|tf( "val|
|00004160| 75 65 20 73 2e 62 2e 20 | 25 67 2c 20 77 61 73 20 |ue s.b. |%g, was |
|00004170| 25 67 5c 6e 22 2c 20 72 | 2c 20 64 20 29 3b 0a 09 |%g\n", r|, d );..|
|00004180| 09 7d 0a 09 7d 0a 0a 76 | 6f 69 64 0a 47 54 65 73 |.}..}..v|oid.GTes|
|00004190| 74 28 29 0a 09 7b 0a 09 | 73 74 61 74 69 63 20 73 |t()..{..|static s|
|000041a0| 74 72 75 63 74 0a 09 09 | 7b 0a 09 09 64 6f 75 62 |truct...|{...doub|
|000041b0| 6c 65 09 69 6e 3b 09 09 | 09 2f 2a 20 69 6e 70 75 |le.in;..|./* inpu|
|000041c0| 74 20 76 61 6c 75 65 73 | 20 2a 2f 0a 09 09 64 6f |t values| */...do|
|000041d0| 75 62 6c 65 09 65 78 70 | 3b 09 09 09 2f 2a 20 65 |uble.exp|;.../* e|
|000041e0| 78 70 65 63 74 65 64 20 | 6f 75 74 70 75 74 20 76 |xpected |output v|
|000041f0| 61 6c 75 65 73 20 2a 2f | 0a 09 09 7d 09 74 62 6c |alues */|...}.tbl|
|00004200| 5b 5d 20 3d 20 09 09 2f | 2a 20 74 61 62 6c 65 20 |[] = ../|* table |
|00004210| 6f 66 20 74 65 73 74 20 | 76 61 6c 75 65 73 20 2a |of test |values *|
|00004220| 2f 0a 09 09 7b 0a 09 09 | 7b 09 31 2e 30 2c 09 09 |/...{...|{.1.0,..|
|00004230| 31 2e 30 09 09 09 7d 2c | 0a 09 09 7b 09 32 2e 30 |1.0...},|...{.2.0|
|00004240| 2c 09 09 31 2e 30 09 09 | 09 7d 2c 0a 09 09 7b 09 |,..1.0..|.},...{.|
|00004250| 33 2e 30 2c 09 09 32 2e | 30 09 09 09 7d 2c 0a 09 |3.0,..2.|0...},..|
|00004260| 09 7b 09 34 2e 30 2c 09 | 09 36 2e 30 09 09 09 7d |.{.4.0,.|.6.0...}|
|00004270| 2c 0a 09 09 7b 09 35 2e | 30 2c 09 09 32 34 2e 30 |,...{.5.|0,..24.0|
|00004280| 09 09 09 7d 2c 0a 09 09 | 7b 09 36 2e 30 2c 09 09 |...},...|{.6.0,..|
|00004290| 31 32 30 2e 30 09 09 09 | 7d 2c 0a 09 09 7b 09 30 |120.0...|},...{.0|
|000042a0| 2e 35 2c 09 09 31 2e 37 | 37 32 34 35 33 38 35 30 |.5,..1.7|72453850|
|000042b0| 39 09 09 7d 2c 0a 09 09 | 7b 09 31 2e 35 2c 09 09 |9..},...|{.1.5,..|
|000042c0| 30 2e 38 38 36 32 32 36 | 39 32 35 35 09 09 7d 2c |0.886226|9255..},|
|000042d0| 0a 09 09 7b 09 30 2e 32 | 35 2c 09 09 33 2e 36 32 |...{.0.2|5,..3.62|
|000042e0| 35 36 30 39 39 30 38 32 | 09 09 7d 2c 0a 09 09 7b |56099082|..},...{|
|000042f0| 09 30 2e 33 33 33 33 33 | 33 33 33 33 33 33 33 2c |.0.33333|3333333,|
|00004300| 09 32 2e 36 37 38 39 33 | 38 35 33 34 37 09 09 7d |.2.67893|85347..}|
|00004310| 2c 0a 09 09 7b 09 30 2e | 36 36 36 36 36 36 36 36 |,...{.0.|66666666|
|00004320| 36 36 36 37 2c 09 31 2e | 33 35 34 31 31 37 39 33 |6667,.1.|35411793|
|00004330| 39 34 09 09 7d 2c 0a 09 | 09 7b 09 30 2e 37 35 2c |94..},..|.{.0.75,|
|00004340| 09 09 31 2e 32 32 35 34 | 31 36 37 30 32 34 09 09 |..1.2254|167024..|
|00004350| 7d 2c 0a 09 09 7b 09 31 | 30 2e 30 2c 09 09 33 36 |},...{.1|0.0,..36|
|00004360| 32 38 38 30 2e 30 09 09 | 7d 2c 0a 09 09 7b 09 32 |2880.0..|},...{.2|
|00004370| 30 2e 30 2c 09 09 31 2e | 32 31 36 34 35 31 30 30 |0.0,..1.|21645100|
|00004380| 34 31 65 2b 31 37 09 7d | 2c 0a 09 09 7d 3b 0a 09 |41e+17.}|,...};..|
|00004390| 72 65 67 69 73 74 65 72 | 20 69 6e 74 09 69 3b 09 |register| int.i;.|
|000043a0| 09 2f 2a 20 69 6e 64 65 | 78 65 73 20 74 62 6c 5b |./* inde|xes tbl[|
|000043b0| 5d 20 74 65 73 74 20 76 | 61 6c 75 65 73 20 2a 2f |] test v|alues */|
|000043c0| 0a 0a 09 66 6f 72 20 28 | 20 69 20 3d 20 30 3b 20 |...for (| i = 0; |
|000043d0| 69 20 3c 20 73 69 7a 65 | 6f 66 20 74 62 6c 20 2f |i < size|of tbl /|
|000043e0| 20 73 69 7a 65 6f 66 20 | 74 62 6c 5b 30 5d 3b 20 | sizeof |tbl[0]; |
|000043f0| 2b 2b 69 20 29 0a 09 09 | 52 43 68 65 63 6b 28 20 |++i )...|RCheck( |
|00004400| 47 61 6d 6d 61 28 20 74 | 62 6c 5b 69 5d 2e 69 6e |Gamma( t|bl[i].in|
|00004410| 20 29 2c 20 74 62 6c 5b | 69 5d 2e 65 78 70 20 29 | ), tbl[|i].exp )|
|00004420| 3b 0a 09 7d 0a 0a 76 6f | 69 64 0a 46 54 65 73 74 |;..}..vo|id.FTest|
|00004430| 28 29 0a 09 7b 0a 09 73 | 74 61 74 69 63 20 73 74 |()..{..s|tatic st|
|00004440| 72 75 63 74 0a 09 09 7b | 0a 09 09 69 6e 74 09 69 |ruct...{|...int.i|
|00004450| 6e 3b 09 09 09 2f 2a 20 | 69 6e 70 75 74 20 76 61 |n;.../* |input va|
|00004460| 6c 75 65 73 20 2a 2f 0a | 09 09 64 6f 75 62 6c 65 |lues */.|..double|
|00004470| 09 65 78 70 3b 09 09 09 | 2f 2a 20 65 78 70 65 63 |.exp;...|/* expec|
|00004480| 74 65 64 20 6f 75 74 70 | 75 74 20 76 61 6c 75 65 |ted outp|ut value|
|00004490| 73 20 2a 2f 0a 09 09 7d | 09 74 62 6c 5b 5d 20 3d |s */...}|.tbl[] =|
|000044a0| 20 09 09 2f 2a 20 74 61 | 62 6c 65 20 6f 66 20 74 | ../* ta|ble of t|
|000044b0| 65 73 74 20 76 61 6c 75 | 65 73 20 2a 2f 0a 09 09 |est valu|es */...|
|000044c0| 7b 0a 09 09 7b 09 30 2c | 09 09 31 2e 30 09 09 09 |{...{.0,|..1.0...|
|000044d0| 7d 2c 0a 09 09 7b 09 31 | 2c 09 09 31 2e 30 09 09 |},...{.1|,..1.0..|
|000044e0| 09 7d 2c 0a 09 09 7b 09 | 32 2c 09 09 32 2e 30 09 |.},...{.|2,..2.0.|
|000044f0| 09 09 7d 2c 0a 09 09 7b | 09 33 2c 09 09 36 2e 30 |..},...{|.3,..6.0|
|00004500| 09 09 09 7d 2c 0a 09 09 | 7b 09 34 2c 09 09 32 34 |...},...|{.4,..24|
|00004510| 2e 30 09 09 09 7d 2c 0a | 09 09 7b 09 35 2c 09 09 |.0...},.|..{.5,..|
|00004520| 31 32 30 2e 30 09 09 09 | 7d 2c 0a 09 09 7b 09 36 |120.0...|},...{.6|
|00004530| 2c 09 09 37 32 30 2e 30 | 09 09 09 7d 2c 0a 09 09 |,..720.0|...},...|
|00004540| 7b 09 31 30 2c 09 09 33 | 36 32 38 38 30 30 2e 30 |{.10,..3|628800.0|
|00004550| 09 09 7d 2c 0a 09 09 7b | 09 32 30 2c 09 09 32 2e |..},...{|.20,..2.|
|00004560| 34 33 32 39 30 32 30 30 | 38 32 65 2b 31 38 09 7d |43290200|82e+18.}|
|00004570| 2c 0a 09 09 7d 3b 0a 09 | 72 65 67 69 73 74 65 72 |,...};..|register|
|00004580| 20 69 6e 74 09 69 3b 09 | 09 2f 2a 20 69 6e 64 65 | int.i;.|./* inde|
|00004590| 78 65 73 20 74 62 6c 5b | 5d 20 74 65 73 74 20 76 |xes tbl[|] test v|
|000045a0| 61 6c 75 65 73 20 2a 2f | 0a 0a 09 66 6f 72 20 28 |alues */|...for (|
|000045b0| 20 69 20 3d 20 30 3b 20 | 69 20 3c 20 73 69 7a 65 | i = 0; |i < size|
|000045c0| 6f 66 20 74 62 6c 20 2f | 20 73 69 7a 65 6f 66 20 |of tbl /| sizeof |
|000045d0| 74 62 6c 5b 30 5d 3b 20 | 2b 2b 69 20 29 0a 09 09 |tbl[0]; |++i )...|
|000045e0| 52 43 68 65 63 6b 28 20 | 46 61 63 74 6f 72 69 61 |RCheck( |Factoria|
|000045f0| 6c 28 20 74 62 6c 5b 69 | 5d 2e 69 6e 20 29 2c 20 |l( tbl[i|].in ), |
|00004600| 74 62 6c 5b 69 5d 2e 65 | 78 70 20 29 3b 0a 09 7d |tbl[i].e|xp );..}|
|00004610| 0a 0a 76 6f 69 64 0a 42 | 43 54 65 73 74 28 29 0a |..void.B|CTest().|
|00004620| 09 7b 0a 09 73 74 61 74 | 69 63 20 73 74 72 75 63 |.{..stat|ic struc|
|00004630| 74 0a 09 09 7b 0a 09 09 | 69 6e 74 09 6e 3b 09 09 |t...{...|int.n;..|
|00004640| 09 2f 2a 20 74 6f 70 20 | 70 61 72 74 73 20 6f 66 |./* top |parts of|
|00004650| 20 69 6e 70 75 74 73 20 | 2a 2f 0a 09 09 69 6e 74 | inputs |*/...int|
|00004660| 09 6b 3b 09 09 09 2f 2a | 20 62 6f 74 74 6f 6d 20 |.k;.../*| bottom |
|00004670| 70 61 72 74 73 20 6f 66 | 20 69 6e 70 75 74 73 20 |parts of| inputs |
|00004680| 2a 2f 0a 09 09 64 6f 75 | 62 6c 65 09 65 78 70 3b |*/...dou|ble.exp;|
|00004690| 09 09 09 2f 2a 20 65 78 | 70 65 63 74 65 64 20 6f |.../* ex|pected o|
|000046a0| 75 74 70 75 74 20 76 61 | 6c 75 65 73 20 2a 2f 0a |utput va|lues */.|
|000046b0| 09 09 7d 09 74 62 6c 5b | 5d 20 3d 20 09 09 2f 2a |..}.tbl[|] = ../*|
|000046c0| 20 74 61 62 6c 65 20 6f | 66 20 74 65 73 74 20 76 | table o|f test v|
|000046d0| 61 6c 75 65 73 20 2a 2f | 0a 09 09 7b 0a 09 09 7b |alues */|...{...{|
|000046e0| 09 31 2c 09 09 30 2c 09 | 09 31 2e 30 09 09 7d 2c |.1,..0,.|.1.0..},|
|000046f0| 0a 09 09 7b 09 31 2c 09 | 09 31 2c 09 09 31 2e 30 |...{.1,.|.1,..1.0|
|00004700| 09 09 7d 2c 0a 09 09 7b | 09 32 2c 09 09 30 2c 09 |..},...{|.2,..0,.|
|00004710| 09 31 2e 30 09 09 7d 2c | 0a 09 09 7b 09 32 2c 09 |.1.0..},|...{.2,.|
|00004720| 09 31 2c 09 09 32 2e 30 | 09 09 7d 2c 0a 09 09 7b |.1,..2.0|..},...{|
|00004730| 09 32 2c 09 09 32 2c 09 | 09 31 2e 30 09 09 7d 2c |.2,..2,.|.1.0..},|
|00004740| 0a 09 09 7b 09 33 2c 09 | 09 30 2c 09 09 31 2e 30 |...{.3,.|.0,..1.0|
|00004750| 09 09 7d 2c 0a 09 09 7b | 09 33 2c 09 09 31 2c 09 |..},...{|.3,..1,.|
|00004760| 09 33 2e 30 09 09 7d 2c | 0a 09 09 7b 09 33 2c 09 |.3.0..},|...{.3,.|
|00004770| 09 32 2c 09 09 33 2e 30 | 09 09 7d 2c 0a 09 09 7b |.2,..3.0|..},...{|
|00004780| 09 35 2c 09 09 33 2c 09 | 09 31 30 2e 30 09 09 7d |.5,..3,.|.10.0..}|
|00004790| 2c 0a 09 09 7b 09 31 30 | 2c 09 09 34 2c 09 09 32 |,...{.10|,..4,..2|
|000047a0| 31 30 2e 30 09 09 7d 2c | 0a 09 09 7b 09 31 30 2c |10.0..},|...{.10,|
|000047b0| 09 09 35 2c 09 09 32 35 | 32 2e 30 09 09 7d 2c 0a |..5,..25|2.0..},.|
|000047c0| 09 09 7b 09 34 30 2c 09 | 09 36 2c 09 09 33 38 33 |..{.40,.|.6,..383|
|000047d0| 38 33 38 30 2e 30 09 7d | 2c 0a 09 09 7b 09 35 30 |8380.0.}|,...{.50|
|000047e0| 2c 09 09 32 30 2c 09 09 | 34 37 31 32 39 32 31 32 |,..20,..|47129212|
|000047f0| 32 34 33 39 36 30 2e 30 | 7d 2c 0a 09 09 7d 3b 0a |243960.0|},...};.|
|00004800| 09 72 65 67 69 73 74 65 | 72 20 69 6e 74 09 69 3b |.registe|r int.i;|
|00004810| 09 09 2f 2a 20 69 6e 64 | 65 78 65 73 20 74 62 6c |../* ind|exes tbl|
|00004820| 5b 5d 20 74 65 73 74 20 | 76 61 6c 75 65 73 20 2a |[] test |values *|
|00004830| 2f 0a 0a 09 66 6f 72 20 | 28 20 69 20 3d 20 30 3b |/...for |( i = 0;|
|00004840| 20 69 20 3c 20 73 69 7a | 65 6f 66 20 74 62 6c 20 | i < siz|eof tbl |
|00004850| 2f 20 73 69 7a 65 6f 66 | 20 74 62 6c 5b 30 5d 3b |/ sizeof| tbl[0];|
|00004860| 20 2b 2b 69 20 29 0a 09 | 09 52 43 68 65 63 6b 28 | ++i )..|.RCheck(|
|00004870| 20 42 43 6f 65 66 66 28 | 20 74 62 6c 5b 69 5d 2e | BCoeff(| tbl[i].|
|00004880| 6e 2c 20 74 62 6c 5b 69 | 5d 2e 6b 20 29 2c 20 74 |n, tbl[i|].k ), t|
|00004890| 62 6c 5b 69 5d 2e 65 78 | 70 20 29 3b 0a 09 7d 0a |bl[i].ex|p );..}.|
|000048a0| 0a 76 6f 69 64 0a 45 54 | 65 73 74 28 29 0a 09 7b |.void.ET|est()..{|
|000048b0| 0a 09 73 74 61 74 69 63 | 20 73 74 72 75 63 74 0a |..static| struct.|
|000048c0| 09 09 7b 0a 09 09 64 6f | 75 62 6c 65 09 69 6e 3b |..{...do|uble.in;|
|000048d0| 09 09 09 2f 2a 20 69 6e | 70 75 74 20 76 61 6c 75 |.../* in|put valu|
|000048e0| 65 73 20 2a 2f 0a 09 09 | 64 6f 75 62 6c 65 09 65 |es */...|double.e|
|000048f0| 78 70 3b 09 09 09 2f 2a | 20 65 78 70 65 63 74 65 |xp;.../*| expecte|
|00004900| 64 20 6f 75 74 70 75 74 | 20 76 61 6c 75 65 73 20 |d output| values |
|00004910| 2a 2f 0a 09 09 7d 09 74 | 62 6c 5b 5d 20 3d 20 09 |*/...}.t|bl[] = .|
|00004920| 09 2f 2a 20 74 61 62 6c | 65 20 6f 66 20 74 65 73 |./* tabl|e of tes|
|00004930| 74 20 76 61 6c 75 65 73 | 20 2a 2f 0a 09 09 7b 0a |t values| */...{.|
|00004940| 09 09 7b 09 30 2e 30 2c | 09 09 30 2e 30 09 09 09 |..{.0.0,|..0.0...|
|00004950| 7d 2c 0a 09 09 7b 09 30 | 2e 31 2c 09 09 30 2e 31 |},...{.0|.1,..0.1|
|00004960| 31 32 34 36 32 39 31 36 | 30 09 09 7d 2c 0a 09 09 |12462916|0..},...|
|00004970| 7b 09 30 2e 32 2c 09 09 | 30 2e 32 32 32 37 30 32 |{.0.2,..|0.222702|
|00004980| 35 38 39 32 09 09 7d 2c | 0a 09 09 7b 09 30 2e 35 |5892..},|...{.0.5|
|00004990| 2c 09 09 30 2e 35 32 30 | 34 39 39 38 37 37 38 09 |,..0.520|4998778.|
|000049a0| 09 7d 2c 0a 09 09 7b 09 | 30 2e 38 2c 09 09 30 2e |.},...{.|0.8,..0.|
|000049b0| 37 34 32 31 30 30 39 36 | 34 37 09 09 7d 2c 0a 09 |74210096|47..},..|
|000049c0| 09 7b 09 31 2e 30 2c 09 | 09 30 2e 38 34 32 37 30 |.{.1.0,.|.0.84270|
|000049d0| 30 37 39 32 39 09 09 7d | 2c 0a 09 09 7b 09 31 2e |07929..}|,...{.1.|
|000049e0| 35 2c 09 09 30 2e 39 36 | 36 31 30 35 31 34 36 35 |5,..0.96|61051465|
|000049f0| 09 09 7d 2c 0a 09 09 7b | 09 32 2e 30 2c 09 09 30 |..},...{|.2.0,..0|
|00004a00| 2e 39 39 35 33 32 32 32 | 36 35 30 09 09 7d 2c 0a |.9953222|650..},.|
|00004a10| 09 09 7d 3b 0a 09 72 65 | 67 69 73 74 65 72 20 69 |..};..re|gister i|
|00004a20| 6e 74 09 69 3b 09 09 2f | 2a 20 69 6e 64 65 78 65 |nt.i;../|* indexe|
|00004a30| 73 20 74 62 6c 5b 5d 20 | 74 65 73 74 20 76 61 6c |s tbl[] |test val|
|00004a40| 75 65 73 20 2a 2f 0a 0a | 09 66 6f 72 20 28 20 69 |ues */..|.for ( i|
|00004a50| 20 3d 20 30 3b 20 69 20 | 3c 20 73 69 7a 65 6f 66 | = 0; i |< sizeof|
|00004a60| 20 74 62 6c 20 2f 20 73 | 69 7a 65 6f 66 20 74 62 | tbl / s|izeof tb|
|00004a70| 6c 5b 30 5d 3b 20 2b 2b | 69 20 29 0a 09 09 52 43 |l[0]; ++|i )...RC|
|00004a80| 68 65 63 6b 28 20 45 72 | 66 28 20 74 62 6c 5b 69 |heck( Er|f( tbl[i|
|00004a90| 5d 2e 69 6e 20 29 2c 20 | 74 62 6c 5b 69 5d 2e 65 |].in ), |tbl[i].e|
|00004aa0| 78 70 20 29 3b 0a 09 7d | 0a 0a 76 6f 69 64 0a 51 |xp );..}|..void.Q|
|00004ab0| 43 54 65 73 74 28 29 0a | 09 7b 0a 09 73 74 61 74 |CTest().|.{..stat|
|00004ac0| 69 63 20 73 74 72 75 63 | 74 0a 09 09 7b 0a 09 09 |ic struc|t...{...|
|00004ad0| 64 6f 75 62 6c 65 09 63 | 68 69 73 71 3b 09 09 09 |double.c|hisq;...|
|00004ae0| 2f 2a 20 63 68 69 2d 73 | 71 75 61 72 65 20 6c 69 |/* chi-s|quare li|
|00004af0| 6d 69 74 20 2a 2f 0a 09 | 09 69 6e 74 09 64 66 3b |mit */..|.int.df;|
|00004b00| 09 09 09 2f 2a 20 64 65 | 67 72 65 65 73 20 6f 66 |.../* de|grees of|
|00004b10| 20 66 72 65 65 64 6f 6d | 20 2a 2f 0a 09 09 64 6f | freedom| */...do|
|00004b20| 75 62 6c 65 09 65 78 70 | 3b 09 09 09 2f 2a 20 65 |uble.exp|;.../* e|
|00004b30| 78 70 65 63 74 65 64 20 | 6f 75 74 70 75 74 20 76 |xpected |output v|
|00004b40| 61 6c 75 65 73 20 2a 2f | 0a 09 09 7d 09 74 62 6c |alues */|...}.tbl|
|00004b50| 5b 5d 20 3d 20 09 09 2f | 2a 20 74 61 62 6c 65 20 |[] = ../|* table |
|00004b60| 6f 66 20 74 65 73 74 20 | 76 61 6c 75 65 73 20 2a |of test |values *|
|00004b70| 2f 0a 09 09 7b 0a 09 09 | 7b 09 30 2e 30 30 31 2c |/...{...|{.0.001,|
|00004b80| 09 09 31 2c 09 09 30 2e | 39 37 34 37 37 09 09 7d |..1,..0.|97477..}|
|00004b90| 2c 0a 09 09 7b 09 30 2e | 30 31 2c 09 09 31 2c 09 |,...{.0.|01,..1,.|
|00004ba0| 09 30 2e 39 32 30 33 34 | 09 09 7d 2c 0a 09 09 7b |.0.92034|..},...{|
|00004bb0| 09 30 2e 30 31 2c 09 09 | 32 2c 09 09 30 2e 39 39 |.0.01,..|2,..0.99|
|00004bc0| 35 30 31 09 09 7d 2c 0a | 09 09 7b 09 30 2e 30 35 |501..},.|..{.0.05|
|00004bd0| 2c 09 09 31 2c 09 09 30 | 2e 38 32 33 30 36 09 09 |,..1,..0|.82306..|
|00004be0| 7d 2c 0a 09 09 7b 09 30 | 2e 31 2c 09 09 31 2c 09 |},...{.0|.1,..1,.|
|00004bf0| 09 30 2e 37 35 31 38 33 | 09 09 7d 2c 0a 09 09 7b |.0.75183|..},...{|
|00004c00| 09 30 2e 31 2c 09 09 32 | 2c 09 09 30 2e 39 35 31 |.0.1,..2|,..0.951|
|00004c10| 32 33 09 09 7d 2c 0a 09 | 09 7b 09 31 2e 30 2c 09 |23..},..|.{.1.0,.|
|00004c20| 09 31 2c 09 09 30 2e 33 | 31 37 33 31 09 09 7d 2c |.1,..0.3|1731..},|
|00004c30| 0a 09 09 7b 09 31 2e 30 | 2c 09 09 32 2c 09 09 30 |...{.1.0|,..2,..0|
|00004c40| 2e 36 30 36 35 33 09 09 | 7d 2c 0a 09 09 7b 09 31 |.60653..|},...{.1|
|00004c50| 2e 30 2c 09 09 33 2c 09 | 09 30 2e 38 30 31 32 35 |.0,..3,.|.0.80125|
|00004c60| 09 09 7d 2c 0a 09 09 7b | 09 31 2e 30 2c 09 09 34 |..},...{|.1.0,..4|
|00004c70| 2c 09 09 30 2e 39 30 39 | 38 30 09 09 7d 2c 0a 09 |,..0.909|80..},..|
|00004c80| 09 7b 09 31 2e 30 2c 09 | 09 35 2c 09 09 30 2e 39 |.{.1.0,.|.5,..0.9|
|00004c90| 36 32 35 37 09 09 7d 2c | 0a 09 09 7b 09 31 2e 35 |6257..},|...{.1.5|
|00004ca0| 2c 09 09 32 2c 09 09 30 | 2e 34 37 32 33 37 09 09 |,..2,..0|.47237..|
|00004cb0| 7d 2c 0a 09 09 7b 09 32 | 2e 30 2c 09 09 31 2c 09 |},...{.2|.0,..1,.|
|00004cc0| 09 30 2e 31 35 37 33 30 | 09 09 7d 2c 0a 09 09 7b |.0.15730|..},...{|
|00004cd0| 09 32 2e 30 2c 09 09 33 | 2c 09 09 30 2e 35 37 32 |.2.0,..3|,..0.572|
|00004ce0| 34 31 09 09 7d 2c 0a 09 | 09 7b 09 32 2e 30 2c 09 |41..},..|.{.2.0,.|
|00004cf0| 09 35 2c 09 09 30 2e 38 | 34 39 31 35 09 09 7d 2c |.5,..0.8|4915..},|
|00004d00| 0a 09 09 7b 09 34 2e 30 | 2c 09 09 36 2c 09 09 30 |...{.4.0|,..6,..0|
|00004d10| 2e 36 37 36 36 38 09 09 | 7d 2c 0a 09 09 7b 09 35 |.67668..|},...{.5|
|00004d20| 2e 30 2c 09 09 35 2c 09 | 09 30 2e 34 31 35 38 38 |.0,..5,.|.0.41588|
|00004d30| 09 09 7d 2c 0a 09 09 7b | 09 31 30 2e 30 2c 09 09 |..},...{|.10.0,..|
|00004d40| 37 2c 09 09 30 2e 31 38 | 38 35 37 33 09 7d 2c 0a |7,..0.18|8573.},.|
|00004d50| 09 09 7b 09 31 30 2e 30 | 2c 09 09 31 30 2c 09 09 |..{.10.0|,..10,..|
|00004d60| 30 2e 34 34 30 34 39 09 | 09 7d 2c 0a 09 09 7b 09 |0.44049.|.},...{.|
|00004d70| 31 30 2e 30 2c 09 09 32 | 30 2c 09 09 30 2e 39 36 |10.0,..2|0,..0.96|
|00004d80| 38 31 37 09 09 7d 2c 0a | 09 09 7b 09 32 30 2e 30 |817..},.|..{.20.0|
|00004d90| 2c 09 09 33 30 2c 09 09 | 30 2e 39 31 36 35 34 09 |,..30,..|0.91654.|
|00004da0| 09 7d 2c 0a 09 09 7b 09 | 34 30 2e 30 2c 09 09 33 |.},...{.|40.0,..3|
|00004db0| 30 2c 09 09 30 2e 31 30 | 34 38 36 34 09 7d 2c 0a |0,..0.10|4864.},.|
|00004dc0| 09 09 7b 09 38 2e 32 36 | 30 34 30 2c 09 32 30 2c |..{.8.26|040,.20,|
|00004dd0| 09 09 30 2e 39 39 09 09 | 7d 2c 0a 09 09 7b 09 31 |..0.99..|},...{.1|
|00004de0| 30 2e 38 35 30 38 2c 09 | 32 30 2c 09 09 30 2e 39 |0.8508,.|20,..0.9|
|00004df0| 35 09 09 7d 2c 0a 09 09 | 7b 09 31 32 2e 34 34 32 |5..},...|{.12.442|
|00004e00| 36 2c 09 32 30 2c 09 09 | 30 2e 39 30 09 09 7d 2c |6,.20,..|0.90..},|
|00004e10| 0a 09 09 7b 09 31 35 2e | 34 35 31 38 2c 09 32 30 |...{.15.|4518,.20|
|00004e20| 2c 09 09 30 2e 37 35 09 | 09 7d 2c 0a 09 09 7b 09 |,..0.75.|.},...{.|
|00004e30| 31 39 2e 33 33 37 34 2c | 09 32 30 2c 09 09 30 2e |19.3374,|.20,..0.|
|00004e40| 35 30 09 09 7d 2c 0a 09 | 09 7b 09 32 33 2e 38 32 |50..},..|.{.23.82|
|00004e50| 37 37 2c 09 32 30 2c 09 | 09 30 2e 32 35 09 09 7d |77,.20,.|.0.25..}|
|00004e60| 2c 0a 09 09 7b 09 32 38 | 2e 34 31 32 30 2c 09 32 |,...{.28|.4120,.2|
|00004e70| 30 2c 09 09 30 2e 31 30 | 09 09 7d 2c 0a 09 09 7b |0,..0.10|..},...{|
|00004e80| 09 33 31 2e 34 31 30 34 | 2c 09 32 30 2c 09 09 30 |.31.4104|,.20,..0|
|00004e90| 2e 30 35 09 09 7d 2c 0a | 09 09 7b 09 33 37 2e 35 |.05..},.|..{.37.5|
|00004ea0| 36 36 32 2c 09 32 30 2c | 09 09 30 2e 30 31 09 09 |662,.20,|..0.01..|
|00004eb0| 7d 2c 0a 09 09 7d 3b 0a | 09 72 65 67 69 73 74 65 |},...};.|.registe|
|00004ec0| 72 20 69 6e 74 09 69 3b | 09 09 2f 2a 20 69 6e 64 |r int.i;|../* ind|
|00004ed0| 65 78 65 73 20 74 62 6c | 5b 5d 20 74 65 73 74 20 |exes tbl|[] test |
|00004ee0| 76 61 6c 75 65 73 20 2a | 2f 0a 0a 09 66 6f 72 20 |values *|/...for |
|00004ef0| 28 20 69 20 3d 20 30 3b | 20 69 20 3c 20 73 69 7a |( i = 0;| i < siz|
|00004f00| 65 6f 66 20 74 62 6c 20 | 2f 20 73 69 7a 65 6f 66 |eof tbl |/ sizeof|
|00004f10| 20 74 62 6c 5b 30 5d 3b | 20 2b 2b 69 20 29 0a 09 | tbl[0];| ++i )..|
|00004f20| 09 52 43 68 65 63 6b 28 | 20 51 43 68 69 53 71 28 |.RCheck(| QChiSq(|
|00004f30| 20 74 62 6c 5b 69 5d 2e | 63 68 69 73 71 2c 20 74 | tbl[i].|chisq, t|
|00004f40| 62 6c 5b 69 5d 2e 64 66 | 20 29 2c 20 74 62 6c 5b |bl[i].df| ), tbl[|
|00004f50| 69 5d 2e 65 78 70 20 29 | 3b 0a 09 7d 0a 0a 2f 2a |i].exp )|;..}../*|
|00004f60| 41 52 47 53 55 53 45 44 | 2a 2f 0a 6d 61 69 6e 28 |ARGSUSED|*/.main(|
|00004f70| 20 61 72 67 63 2c 20 61 | 72 67 76 20 29 0a 09 69 | argc, a|rgv )..i|
|00004f80| 6e 74 09 09 61 72 67 63 | 3b 0a 09 63 68 61 72 09 |nt..argc|;..char.|
|00004f90| 09 2a 61 72 67 76 5b 5d | 3b 0a 09 7b 0a 09 47 54 |.*argv[]|;..{..GT|
|00004fa0| 65 73 74 28 29 3b 0a 09 | 46 54 65 73 74 28 29 3b |est();..|FTest();|
|00004fb0| 0a 09 42 43 54 65 73 74 | 28 29 3b 0a 09 45 54 65 |..BCTest|();..ETe|
|00004fc0| 73 74 28 29 3b 0a 09 51 | 43 54 65 73 74 28 29 3b |st();..Q|CTest();|
|00004fd0| 0a 0a 09 72 65 74 75 72 | 6e 20 65 72 72 73 3b 0a |...retur|n errs;.|
|00004fe0| 09 7d 0a 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.}......|........|
|00004ff0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005000| 74 6f 74 2d 69 6e 66 6f | 2f 67 61 6d 6d 61 2e 33 |tot-info|/gamma.3|
|00005010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005040| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005050| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005060| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00005070| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00005080| 20 20 20 36 30 31 30 20 | 20 34 34 35 37 32 31 30 |   6010 | 4457210|
|00005090| 34 37 31 20 20 20 36 37 | 36 34 00 20 00 00 00 00 |471   67|64. ....|
|000050a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000050b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000050c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000050d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000050e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000050f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005100| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005110| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005120| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005130| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005140| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005150| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005160| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005170| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005180| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005190| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000051f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005200| 27 5c 22 20 65 0a 2e 54 | 48 20 47 41 4d 4d 41 20 |'\" e..T|H GAMMA |
|00005210| 33 56 20 56 4d 42 0a 27 | 5c 22 09 6c 61 73 74 20 |3V VMB.'|\".last |
|00005220| 65 64 69 74 3a 09 38 38 | 2f 30 39 2f 30 39 09 44 |edit:.88|/09/09.D|
|00005230| 20 41 20 47 77 79 6e 0a | 27 5c 22 09 53 43 43 53 | A Gwyn.|'\".SCCS|
|00005240| 20 49 44 3a 09 40 28 23 | 29 67 61 6d 6d 61 2e 33 | ID:.@(#|)gamma.3|
|00005250| 09 31 2e 32 20 28 65 64 | 69 74 65 64 20 66 6f 72 |.1.2 (ed|ited for|
|00005260| 20 70 75 62 6c 69 63 61 | 74 69 6f 6e 29 0a 27 5c | publica|tion).'\|
|00005270| 22 0a 27 5c 22 09 54 68 | 69 73 20 73 6f 75 72 63 |".'\".Th|is sourc|
|00005280| 65 20 6d 75 73 74 20 62 | 65 20 74 79 70 65 73 65 |e must b|e typese|
|00005290| 74 20 76 69 61 20 22 65 | 71 6e 20 7c 20 74 72 6f |t via "e|qn | tro|
|000052a0| 66 66 20 2d 6d 61 6e 22 | 0a 27 5c 22 0a 2e 45 51 |ff -man"|.'\"..EQ|
|000052b0| 0a 64 65 6c 69 6d 20 40 | 40 0a 2e 45 4e 0a 2e 53 |.delim @|@..EN..S|
|000052c0| 48 20 4e 41 4d 45 0a 67 | 61 6d 6d 61 20 5c 2d 20 |H NAME.g|amma \- |
|000052d0| 67 61 6d 6d 61 20 61 6e | 64 20 72 65 6c 61 74 65 |gamma an|d relate|
|000052e0| 64 20 66 75 6e 63 74 69 | 6f 6e 73 0a 2e 53 48 20 |d functi|ons..SH |
|000052f0| 53 59 4e 4f 50 53 49 53 | 0a 2e 64 73 20 63 57 20 |SYNOPSIS|..ds cW |
|00005300| 28 43 57 5c 22 20 63 68 | 61 6e 67 65 20 74 6f 20 |(CW\" ch|ange to |
|00005310| 42 20 28 77 69 74 68 6f | 75 74 20 74 68 65 20 70 |B (witho|ut the p|
|00005320| 61 72 65 6e 29 20 69 66 | 20 79 6f 75 20 64 6f 6e |aren) if| you don|
|00005330| 27 74 20 68 61 76 65 20 | 61 20 43 57 20 66 6f 6e |'t have |a CW fon|
|00005340| 74 0a 5c 66 5c 2a 28 63 | 57 0a 23 69 6e 63 6c 75 |t.\f\*(c|W.#inclu|
|00005350| 64 65 20 22 67 61 6d 6d | 61 2e 68 22 0a 2e 62 72 |de "gamm|a.h"..br|
|00005360| 0a 2f 2a 20 5c 66 50 41 | 6c 6c 20 74 68 65 20 66 |./* \fPA|ll the f|
|00005370| 6f 6c 6c 6f 77 69 6e 67 | 20 66 75 6e 63 74 69 6f |ollowing| functio|
|00005380| 6e 73 20 61 72 65 20 64 | 65 63 6c 61 72 65 64 20 |ns are d|eclared |
|00005390| 69 6e 20 74 68 69 73 20 | 68 65 61 64 65 72 20 66 |in this |header f|
|000053a0| 69 6c 65 2e 5c 66 5c 2a | 28 63 57 20 2a 2f 0a 2e |ile.\f\*|(cW */..|
|000053b0| 73 70 0a 64 6f 75 62 6c | 65 20 47 61 6d 6d 61 28 |sp.doubl|e Gamma(|
|000053c0| 64 6f 75 62 6c 65 20 7a | 29 3b 0a 2e 62 72 0a 64 |double z|);..br.d|
|000053d0| 6f 75 62 6c 65 20 4c 47 | 61 6d 6d 61 28 64 6f 75 |ouble LG|amma(dou|
|000053e0| 62 6c 65 20 7a 29 3b 0a | 2e 62 72 0a 64 6f 75 62 |ble z);.|.br.doub|
|000053f0| 6c 65 20 46 61 63 74 6f | 72 69 61 6c 28 69 6e 74 |le Facto|rial(int|
|00005400| 20 6e 29 3b 0a 2e 62 72 | 0a 64 6f 75 62 6c 65 20 | n);..br|.double |
|00005410| 4c 46 61 63 74 6f 72 69 | 61 6c 28 69 6e 74 20 6e |LFactori|al(int n|
|00005420| 29 3b 0a 2e 62 72 0a 64 | 6f 75 62 6c 65 20 42 43 |);..br.d|ouble BC|
|00005430| 6f 65 66 66 28 69 6e 74 | 20 6e 2c 20 69 6e 74 20 |oeff(int| n, int |
|00005440| 6b 29 3b 0a 2e 62 72 0a | 64 6f 75 62 6c 65 20 42 |k);..br.|double B|
|00005450| 65 74 61 28 64 6f 75 62 | 6c 65 20 7a 2c 20 64 6f |eta(doub|le z, do|
|00005460| 75 62 6c 65 20 77 29 3b | 0a 2e 62 72 0a 64 6f 75 |uble w);|..br.dou|
|00005470| 62 6c 65 20 50 47 61 6d | 6d 61 28 64 6f 75 62 6c |ble PGam|ma(doubl|
|00005480| 65 20 61 2c 20 64 6f 75 | 62 6c 65 20 78 29 3b 0a |e a, dou|ble x);.|
|00005490| 2e 62 72 0a 64 6f 75 62 | 6c 65 20 51 47 61 6d 6d |.br.doub|le QGamm|
|000054a0| 61 28 64 6f 75 62 6c 65 | 20 61 2c 20 64 6f 75 62 |a(double| a, doub|
|000054b0| 6c 65 20 78 29 3b 0a 2e | 62 72 0a 64 6f 75 62 6c |le x);..|br.doubl|
|000054c0| 65 20 45 72 66 28 64 6f | 75 62 6c 65 20 78 29 3b |e Erf(do|uble x);|
|000054d0| 0a 2e 62 72 0a 64 6f 75 | 62 6c 65 20 45 72 66 63 |..br.dou|ble Erfc|
|000054e0| 28 64 6f 75 62 6c 65 20 | 78 29 3b 0a 2e 62 72 0a |(double |x);..br.|
|000054f0| 64 6f 75 62 6c 65 20 43 | 50 6f 69 73 73 6f 6e 28 |double C|Poisson(|
|00005500| 64 6f 75 62 6c 65 20 78 | 2c 20 69 6e 74 20 6b 29 |double x|, int k)|
|00005510| 3b 0a 2e 62 72 0a 64 6f | 75 62 6c 65 20 50 43 68 |;..br.do|uble PCh|
|00005520| 69 53 71 28 64 6f 75 62 | 6c 65 20 63 68 69 73 71 |iSq(doub|le chisq|
|00005530| 2c 20 69 6e 74 20 6e 75 | 29 3b 0a 2e 62 72 0a 64 |, int nu|);..br.d|
|00005540| 6f 75 62 6c 65 20 51 43 | 68 69 53 71 28 64 6f 75 |ouble QC|hiSq(dou|
|00005550| 62 6c 65 20 63 68 69 73 | 71 2c 20 69 6e 74 20 6e |ble chis|q, int n|
|00005560| 75 29 3b 0a 2e 66 74 0a | 2e 53 48 20 44 45 53 43 |u);..ft.|.SH DESC|
|00005570| 52 49 50 54 49 4f 4e 0a | 5c 66 5c 2a 28 63 57 47 |RIPTION.|\f\*(cWG|
|00005580| 61 6d 6d 61 5c 66 50 0a | 72 65 74 75 72 6e 73 20 |amma\fP.|returns |
|00005590| 74 68 65 20 72 65 61 6c | 20 67 61 6d 6d 61 20 66 |the real| gamma f|
|000055a0| 75 6e 63 74 69 6f 6e 20 | 76 61 6c 75 65 0a 40 47 |unction |value.@G|
|000055b0| 41 4d 4d 41 20 28 7a 29 | 7e 3d 3d 7e 20 69 6e 74 |AMMA (z)|~==~ int|
|000055c0| 20 22 22 20 66 72 6f 6d | 20 30 20 74 6f 20 69 6e | "" from| 0 to in|
|000055d0| 66 20 74 20 73 75 70 20 | 7a 2d 31 20 65 20 73 75 |f t sup |z-1 e su|
|000055e0| 70 20 2d 74 20 64 74 40 | 0a 66 6f 72 20 40 7a 3e |p -t dt@|.for @z>|
|000055f0| 30 40 2e 0a 2e 50 0a 5c | 66 5c 2a 28 63 57 4c 47 |0@...P.\|f\*(cWLG|
|00005600| 61 6d 6d 61 5c 66 50 0a | 72 65 74 75 72 6e 73 20 |amma\fP.|returns |
|00005610| 74 68 65 20 76 61 6c 75 | 65 0a 40 6c 6e 20 7e 20 |the valu|e.@ln ~ |
|00005620| 47 41 4d 4d 41 20 28 7a | 29 40 0a 66 6f 72 20 40 |GAMMA (z|)@.for @|
|00005630| 7a 3e 30 40 2e 0a 2e 50 | 0a 5c 66 5c 2a 28 63 57 |z>0@...P|.\f\*(cW|
|00005640| 46 61 63 74 6f 72 69 61 | 6c 5c 66 50 0a 72 65 74 |Factoria|l\fP.ret|
|00005650| 75 72 6e 73 20 74 68 65 | 20 66 61 63 74 6f 72 69 |urns the| factori|
|00005660| 61 6c 20 66 75 6e 63 74 | 69 6f 6e 20 76 61 6c 75 |al funct|ion valu|
|00005670| 65 0a 40 6e 21 20 7e 3d | 7e 20 47 41 4d 4d 41 20 |e.@n! ~=|~ GAMMA |
|00005680| 28 6e 2b 31 29 40 0a 66 | 6f 72 20 40 6e 3e 3d 30 |(n+1)@.f|or @n>=0|
|00005690| 40 2e 0a 2e 50 0a 5c 66 | 5c 2a 28 63 57 4c 46 61 |@...P.\f|\*(cWLFa|
|000056a0| 63 74 6f 72 69 61 6c 5c | 66 50 0a 72 65 74 75 72 |ctorial\|fP.retur|
|000056b0| 6e 73 20 74 68 65 20 76 | 61 6c 75 65 0a 40 6c 6e |ns the v|alue.@ln|
|000056c0| 20 7e 6e 21 40 5c 26 0a | 66 6f 72 20 40 6e 3e 3d | ~n!@\&.|for @n>=|
|000056d0| 30 40 2e 0a 2e 50 0a 5c | 66 5c 2a 28 63 57 42 43 |0@...P.\|f\*(cWBC|
|000056e0| 6f 65 66 66 5c 66 50 0a | 72 65 74 75 72 6e 73 20 |oeff\fP.|returns |
|000056f0| 74 68 65 20 62 69 6e 6f | 6d 69 61 6c 20 63 6f 65 |the bino|mial coe|
|00005700| 66 66 69 63 69 65 6e 74 | 20 76 61 6c 75 65 0a 40 |fficient| value.@|
|00005710| 73 69 7a 65 20 2b 38 20 | 28 20 70 69 6c 65 20 7b |size +8 |( pile {|
|00005720| 20 6e 20 61 62 6f 76 65 | 20 6b 20 7d 20 73 69 7a | n above| k } siz|
|00005730| 65 20 2b 38 20 29 20 7e | 3d 3d 7e 6e 21 20 6f 76 |e +8 ) ~|==~n! ov|
|00005740| 65 72 20 6b 21 28 6e 2d | 6b 29 21 40 5c 26 0a 66 |er k!(n-|k)!@\&.f|
|00005750| 6f 72 20 40 30 3c 3d 6b | 3c 3d 6e 40 2e 0a 2e 50 |or @0<=k|<=n@...P|
|00005760| 0a 5c 66 5c 2a 28 63 57 | 42 65 74 61 5c 66 50 0a |.\f\*(cW|Beta\fP.|
|00005770| 72 65 74 75 72 6e 73 20 | 74 68 65 20 72 65 61 6c |returns |the real|
|00005780| 20 62 65 74 61 20 66 75 | 6e 63 74 69 6f 6e 20 76 | beta fu|nction v|
|00005790| 61 6c 75 65 0a 40 42 28 | 7a 2c 77 29 7e 3d 3d 7e |alue.@B(|z,w)~==~|
|000057a0| 20 69 6e 74 20 22 22 20 | 66 72 6f 6d 20 30 20 74 | int "" |from 0 t|
|000057b0| 6f 20 31 20 74 20 73 75 | 70 20 7a 2d 31 20 28 31 |o 1 t su|p z-1 (1|
|000057c0| 2d 74 29 20 73 75 70 20 | 77 2d 31 20 64 74 40 0a |-t) sup |w-1 dt@.|
|000057d0| 66 6f 72 20 40 7a 3e 30 | 40 20 61 6e 64 20 40 77 |for @z>0|@ and @w|
|000057e0| 3e 30 40 2e 0a 2e 50 0a | 5c 66 5c 2a 28 63 57 50 |>0@...P.|\f\*(cWP|
|000057f0| 47 61 6d 6d 61 5c 66 50 | 0a 72 65 74 75 72 6e 73 |Gamma\fP|.returns|
|00005800| 20 74 68 65 20 69 6e 63 | 6f 6d 70 6c 65 74 65 20 | the inc|omplete |
|00005810| 67 61 6d 6d 61 20 66 75 | 6e 63 74 69 6f 6e 20 76 |gamma fu|nction v|
|00005820| 61 6c 75 65 0a 40 50 28 | 61 2c 78 29 7e 3d 3d 7e |alue.@P(|a,x)~==~|
|00005830| 20 31 20 6f 76 65 72 20 | 7b 20 47 41 4d 4d 41 20 | 1 over |{ GAMMA |
|00005840| 28 61 29 20 7d 20 69 6e | 74 20 22 22 20 66 72 6f |(a) } in|t "" fro|
|00005850| 6d 20 30 20 74 6f 20 78 | 20 65 20 73 75 70 20 2d |m 0 to x| e sup -|
|00005860| 74 20 74 20 73 75 70 20 | 61 2d 31 20 64 74 40 0a |t t sup |a-1 dt@.|
|00005870| 66 6f 72 20 40 61 3e 30 | 40 20 61 6e 64 20 40 30 |for @a>0|@ and @0|
|00005880| 3c 3d 78 3c 3d 20 69 6e | 66 40 2e 0a 2e 50 0a 5c |<=x<= in|f@...P.\|
|00005890| 66 5c 2a 28 63 57 51 47 | 61 6d 6d 61 5c 66 50 0a |f\*(cWQG|amma\fP.|
|000058a0| 72 65 74 75 72 6e 73 20 | 74 68 65 20 76 61 6c 75 |returns |the valu|
|000058b0| 65 0a 40 31 2d 50 28 61 | 2c 78 29 40 2e 0a 2e 50 |e.@1-P(a|,x)@...P|
|000058c0| 0a 5c 66 5c 2a 28 63 57 | 45 72 66 5c 66 50 0a 72 |.\f\*(cW|Erf\fP.r|
|000058d0| 65 74 75 72 6e 73 20 74 | 68 65 20 65 72 72 6f 72 |eturns t|he error|
|000058e0| 20 66 75 6e 63 74 69 6f | 6e 20 76 61 6c 75 65 0a | functio|n value.|
|000058f0| 40 72 6f 6d 61 6e 20 65 | 72 66 20 28 78 29 7e 3d |@roman e|rf (x)~=|
|00005900| 3d 7e 32 20 6f 76 65 72 | 20 73 71 72 74 20 70 69 |=~2 over| sqrt pi|
|00005910| 20 69 6e 74 20 22 22 20 | 66 72 6f 6d 20 30 20 74 | int "" |from 0 t|
|00005920| 6f 20 78 20 65 20 73 75 | 70 20 7b 2d 74 20 73 75 |o x e su|p {-t su|
|00005930| 70 20 32 20 7d 20 64 74 | 40 2e 0a 2e 50 0a 5c 66 |p 2 } dt|@...P.\f|
|00005940| 5c 2a 28 63 57 45 72 66 | 63 5c 66 50 0a 72 65 74 |\*(cWErf|c\fP.ret|
|00005950| 75 72 6e 73 20 74 68 65 | 20 63 6f 6d 70 6c 65 6d |urns the| complem|
|00005960| 65 6e 74 61 72 79 20 65 | 72 72 6f 72 20 66 75 6e |entary e|rror fun|
|00005970| 63 74 69 6f 6e 20 76 61 | 6c 75 65 0a 40 31 2d 20 |ction va|lue.@1- |
|00005980| 72 6f 6d 61 6e 20 65 72 | 66 20 28 78 29 40 2e 0a |roman er|f (x)@..|
|00005990| 2e 50 0a 5c 66 5c 2a 28 | 63 57 43 50 6f 69 73 73 |.P.\f\*(|cWCPoiss|
|000059a0| 6f 6e 5c 66 50 0a 72 65 | 74 75 72 6e 73 20 74 68 |on\fP.re|turns th|
|000059b0| 65 20 63 75 6d 75 6c 61 | 74 69 76 65 20 50 6f 69 |e cumula|tive Poi|
|000059c0| 73 73 6f 6e 20 70 72 6f | 62 61 62 69 6c 69 74 79 |sson pro|bability|
|000059d0| 20 76 61 6c 75 65 0a 40 | 50 20 73 75 62 20 78 20 | value.@|P sub x |
|000059e0| 28 3c 6b 29 40 2c 0a 77 | 68 69 63 68 20 69 73 20 |(<k)@,.w|hich is |
|000059f0| 74 68 65 20 70 72 6f 62 | 61 62 69 6c 69 74 79 20 |the prob|ability |
|00005a00| 74 68 61 74 20 74 68 65 | 20 6e 75 6d 62 65 72 20 |that the| number |
|00005a10| 6f 66 20 50 6f 69 73 73 | 6f 6e 20 72 61 6e 64 6f |of Poiss|on rando|
|00005a20| 6d 20 65 76 65 6e 74 73 | 20 6f 63 63 75 72 72 69 |m events| occurri|
|00005a30| 6e 67 0a 77 69 6c 6c 20 | 62 65 20 62 65 74 77 65 |ng.will |be betwe|
|00005a40| 65 6e 20 30 20 61 6e 64 | 20 40 6b 2d 31 40 20 69 |en 0 and| @k-1@ i|
|00005a50| 6e 63 6c 75 73 69 76 65 | 2c 0a 77 68 65 72 65 20 |nclusive|,.where |
|00005a60| 74 68 65 20 65 78 70 65 | 63 74 65 64 20 6d 65 61 |the expe|cted mea|
|00005a70| 6e 20 69 73 20 40 78 40 | 2e 0a 2e 50 0a 5c 66 5c |n is @x@|...P.\f\|
|00005a80| 2a 28 63 57 50 43 68 69 | 53 71 5c 66 50 0a 72 65 |*(cWPChi|Sq\fP.re|
|00005a90| 74 75 72 6e 73 20 74 68 | 65 20 76 61 6c 75 65 0a |turns th|e value.|
|00005aa0| 40 50 28 20 63 68 69 20 | 73 75 70 20 32 20 7c 20 |@P( chi |sup 2 | |
|00005ab0| 6e 75 20 29 40 2c 0a 77 | 68 69 63 68 20 69 73 20 |nu )@,.w|hich is |
|00005ac0| 74 68 65 20 70 72 6f 62 | 61 62 69 6c 69 74 79 20 |the prob|ability |
|00005ad0| 74 68 61 74 20 74 68 65 | 20 6f 62 73 65 72 76 65 |that the| observe|
|00005ae0| 64 20 63 68 69 2d 73 71 | 75 61 72 65 20 66 6f 72 |d chi-sq|uare for|
|00005af0| 20 61 20 63 6f 72 72 65 | 63 74 20 6d 6f 64 65 6c | a corre|ct model|
|00005b00| 0a 77 69 6c 6c 20 62 65 | 20 6c 65 73 73 20 74 68 |.will be| less th|
|00005b10| 61 6e 20 40 63 68 69 20 | 73 75 70 20 32 40 2c 0a |an @chi |sup 2@,.|
|00005b20| 77 68 65 72 65 20 40 6e | 75 40 20 69 73 20 74 68 |where @n|u@ is th|
|00005b30| 65 20 6e 75 6d 62 65 72 | 20 6f 66 20 64 65 67 72 |e number| of degr|
|00005b40| 65 65 73 20 6f 66 20 66 | 72 65 65 64 6f 6d 2e 0a |ees of f|reedom..|
|00005b50| 2e 50 0a 5c 66 5c 2a 28 | 63 57 51 43 68 69 53 71 |.P.\f\*(|cWQChiSq|
|00005b60| 5c 66 50 0a 72 65 74 75 | 72 6e 73 20 74 68 65 20 |\fP.retu|rns the |
|00005b70| 76 61 6c 75 65 0a 40 31 | 2d 50 28 20 63 68 69 20 |value.@1|-P( chi |
|00005b80| 73 75 70 20 32 20 7c 20 | 6e 75 20 29 40 2c 0a 77 |sup 2 | |nu )@,.w|
|00005b90| 68 69 63 68 20 69 73 20 | 74 68 65 20 70 72 6f 62 |hich is |the prob|
|00005ba0| 61 62 69 6c 69 74 79 20 | 74 68 61 74 20 74 68 65 |ability |that the|
|00005bb0| 20 6f 62 73 65 72 76 65 | 64 20 63 68 69 2d 73 71 | observe|d chi-sq|
|00005bc0| 75 61 72 65 0a 77 69 6c | 6c 20 65 78 63 65 65 64 |uare.wil|l exceed|
|00005bd0| 20 40 63 68 69 20 73 75 | 70 20 32 40 20 62 79 20 | @chi su|p 2@ by |
|00005be0| 63 68 61 6e 63 65 20 65 | 76 65 6e 20 66 6f 72 20 |chance e|ven for |
|00005bf0| 61 20 63 6f 72 72 65 63 | 74 20 6d 6f 64 65 6c 2c |a correc|t model,|
|00005c00| 0a 77 68 65 72 65 20 40 | 6e 75 40 20 69 73 20 74 |.where @|nu@ is t|
|00005c10| 68 65 20 6e 75 6d 62 65 | 72 20 6f 66 20 64 65 67 |he numbe|r of deg|
|00005c20| 72 65 65 73 20 6f 66 20 | 66 72 65 65 64 6f 6d 2e |rees of |freedom.|
|00005c30| 0a 2e 50 0a 41 6c 6c 20 | 74 68 65 73 65 20 66 75 |..P.All |these fu|
|00005c40| 6e 63 74 69 6f 6e 73 20 | 72 65 74 75 72 6e 0a 2e |nctions |return..|
|00005c50| 49 20 61 70 70 72 6f 78 | 69 6d 61 74 69 6f 6e 73 |I approx|imations|
|00005c60| 0a 74 6f 20 74 68 65 20 | 65 78 61 63 74 20 72 65 |.to the |exact re|
|00005c70| 73 75 6c 74 73 2c 0a 67 | 65 6e 65 72 61 6c 6c 79 |sults,.g|enerally|
|00005c80| 20 61 63 63 75 72 61 74 | 65 20 74 6f 20 61 62 6f | accurat|e to abo|
|00005c90| 75 74 20 73 65 76 65 6e | 20 73 69 67 6e 69 66 69 |ut seven| signifi|
|00005ca0| 63 61 6e 74 20 64 69 67 | 69 74 73 2e 0a 2e 53 48 |cant dig|its...SH|
|00005cb0| 20 43 41 56 45 41 54 53 | 0a 4f 76 65 72 66 6c 6f | CAVEATS|.Overflo|
|00005cc0| 77 20 63 61 6e 20 6f 63 | 63 75 72 2e 0a 2e 50 0a |w can oc|cur...P.|
|00005cd0| 49 6e 76 61 6c 69 64 20 | 70 61 72 61 6d 65 74 65 |Invalid |paramete|
|00005ce0| 72 73 20 63 61 75 73 65 | 20 61 20 64 69 61 67 6e |rs cause| a diagn|
|00005cf0| 6f 73 74 69 63 20 61 6e | 64 20 74 65 72 6d 69 6e |ostic an|d termin|
|00005d00| 61 74 69 6f 6e 20 6f 66 | 20 74 68 65 20 70 72 6f |ation of| the pro|
|00005d10| 63 65 73 73 2e 0a 2e 53 | 48 20 46 49 4c 45 53 0a |cess...S|H FILES.|
|00005d20| 2e 74 61 20 33 69 0a 67 | 61 6d 6d 61 2e 68 09 68 |.ta 3i.g|amma.h.h|
|00005d30| 65 61 64 65 72 20 66 69 | 6c 65 0a 2e 62 72 0a 67 |eader fi|le..br.g|
|00005d40| 61 6d 6d 61 2e 6f 09 6f | 62 6a 65 63 74 20 6d 6f |amma.o.o|bject mo|
|00005d50| 64 75 6c 65 0a 2e 53 48 | 20 41 55 54 48 4f 52 53 |dule..SH| AUTHORS|
|00005d60| 0a 57 69 6c 6c 69 61 6d | 20 48 2e 5c 26 20 50 72 |.William| H.\& Pr|
|00005d70| 65 73 73 2c 0a 42 72 69 | 61 6e 20 50 2e 5c 26 20 |ess,.Bri|an P.\& |
|00005d80| 46 6c 61 6e 6e 65 72 79 | 2c 0a 53 61 75 6c 20 41 |Flannery|,.Saul A|
|00005d90| 2e 5c 26 20 54 65 75 6b | 6f 6c 73 6b 79 2c 0a 61 |.\& Teuk|olsky,.a|
|00005da0| 6e 64 0a 57 69 6c 6c 69 | 61 6d 20 54 2e 5c 26 20 |nd.Willi|am T.\& |
|00005db0| 56 65 74 74 65 72 6c 69 | 6e 67 0a 28 5c 66 49 4e |Vetterli|ng.(\fIN|
|00005dc0| 75 6d 65 72 69 63 61 6c | 20 52 65 63 69 70 65 73 |umerical| Recipes|
|00005dd0| 20 69 6e 20 43 5c 66 50 | 2c 20 43 68 61 70 74 65 | in C\fP|, Chapte|
|00005de0| 72 20 36 29 0a 2e 62 72 | 0a 44 6f 75 67 6c 61 73 |r 6)..br|.Douglas|
|00005df0| 20 41 2e 5c 26 20 47 77 | 79 6e 2c 20 42 52 4c 2f | A.\& Gw|yn, BRL/|
|00005e00| 56 4c 44 2d 56 4d 42 0a | 2c 09 09 30 2e 39 35 31 |VLD-VMB.|,..0.951|
|00005e10| 32 33 09 09 7d 2c 0a 09 | 09 7b 09 31 2e 30 2c 09 |23..},..|.{.1.0,.|
|00005e20| 09 31 2c 09 09 30 2e 33 | 31 37 33 31 09 09 7d 2c |.1,..0.3|1731..},|
|00005e30| 0a 09 09 7b 09 31 2e 30 | 2c 09 09 32 2c 09 09 30 |...{.1.0|,..2,..0|
|00005e40| 2e 36 30 36 35 33 09 09 | 7d 2c 0a 09 09 7b 09 31 |.60653..|},...{.1|
|00005e50| 2e 30 2c 09 09 33 2c 09 | 09 30 2e 38 30 31 32 35 |.0,..3,.|.0.80125|
|00005e60| 09 09 7d 2c 0a 09 09 7b | 09 31 2e 30 2c 09 09 34 |..},...{|.1.0,..4|
|00005e70| 2c 09 09 30 2e 39 30 39 | 38 30 09 09 7d 2c 0a 09 |,..0.909|80..},..|
|00005e80| 09 7b 09 31 2e 30 2c 09 | 09 35 2c 09 09 30 2e 39 |.{.1.0,.|.5,..0.9|
|00005e90| 36 32 35 37 09 09 7d 2c | 0a 09 09 7b 09 31 2e 35 |6257..},|...{.1.5|
|00005ea0| 2c 09 09 32 2c 09 09 30 | 2e 34 37 32 33 37 09 09 |,..2,..0|.47237..|
|00005eb0| 7d 2c 0a 09 09 7b 09 32 | 2e 30 2c 09 09 31 2c 09 |},...{.2|.0,..1,.|
|00005ec0| 09 30 2e 31 35 37 33 30 | 09 09 7d 2c 0a 09 09 7b |.0.15730|..},...{|
|00005ed0| 09 32 2e 30 2c 09 09 33 | 2c 09 09 30 2e 35 37 32 |.2.0,..3|,..0.572|
|00005ee0| 34 31 09 09 7d 2c 0a 09 | 09 7b 09 32 2e 30 2c 09 |41..},..|.{.2.0,.|
|00005ef0| 09 35 2c 09 09 30 2e 38 | 34 39 31 35 09 09 7d 2c |.5,..0.8|4915..},|
|00005f00| 0a 09 09 7b 09 34 2e 30 | 2c 09 09 36 2c 09 09 30 |...{.4.0|,..6,..0|
|00005f10| 2e 36 37 36 36 38 09 09 | 7d 2c 0a 09 09 7b 09 35 |.67668..|},...{.5|
|00005f20| 2e 30 2c 09 09 35 2c 09 | 09 30 2e 34 31 35 38 38 |.0,..5,.|.0.41588|
|00005f30| 09 09 7d 2c 0a 09 09 7b | 09 31 30 2e 30 2c 09 09 |..},...{|.10.0,..|
|00005f40| 37 2c 09 09 30 2e 31 38 | 38 35 37 33 09 7d 2c 0a |7,..0.18|8573.},.|
|00005f50| 09 09 7b 09 31 30 2e 30 | 2c 09 09 31 30 2c 09 09 |..{.10.0|,..10,..|
|00005f60| 30 2e 34 34 30 34 39 09 | 09 7d 2c 0a 09 09 7b 09 |0.44049.|.},...{.|
|00005f70| 31 30 2e 30 2c 09 09 32 | 30 2c 09 09 30 2e 39 36 |10.0,..2|0,..0.96|
|00005f80| 38 31 37 09 09 7d 2c 0a | 09 09 7b 09 32 30 2e 30 |817..},.|..{.20.0|
|00005f90| 2c 09 09 33 30 2c 09 09 | 30 2e 39 31 36 35 34 09 |,..30,..|0.91654.|
|00005fa0| 09 7d 2c 0a 09 09 7b 09 | 34 30 2e 30 2c 09 09 33 |.},...{.|40.0,..3|
|00005fb0| 30 2c 09 09 30 2e 31 30 | 34 38 36 34 09 7d 2c 0a |0,..0.10|4864.},.|
|00005fc0| 09 09 7b 09 38 2e 32 36 | 30 34 30 2c 09 32 30 2c |..{.8.26|040,.20,|
|00005fd0| 09 09 30 2e 39 39 09 09 | 7d 2c 0a 09 09 7b 09 31 |..0.99..|},...{.1|
|00005fe0| 30 2e 38 35 30 38 2c 09 | 32 30 2c 09 09 30 2e 39 |0.8508,.|20,..0.9|
|00005ff0| 35 09 09 7d 2c 0a 09 09 | 7b 09 31 32 2e 34 34 32 |5..},...|{.12.442|
|00006000| 74 6f 74 2d 69 6e 66 6f | 2f 67 61 6d 6d 61 2e 63 |tot-info|/gamma.c|
|00006010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006040| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006050| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006060| 00 00 00 00 20 20 20 36 | 34 34 20 00 20 20 20 32 |....   6|44 .   2|
|00006070| 34 33 20 00 20 20 20 20 | 20 30 20 00 20 20 20 20 |43 .    | 0 .    |
|00006080| 20 20 20 36 37 37 30 20 | 20 34 34 35 37 32 31 30 |   6770 | 4457210|
|00006090| 34 37 32 20 20 20 37 30 | 36 32 00 20 00 00 00 00 |472   70|62. ....|
|000060a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000060b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000060c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000060d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000060e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000060f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006100| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006110| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006120| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006130| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006140| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006150| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006160| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006170| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006180| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006190| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000061f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00006200| 2f 2a 0a 09 47 61 6d 6d | 61 20 2d 2d 20 67 61 6d |/*..Gamm|a -- gam|
|00006210| 6d 61 20 61 6e 64 20 72 | 65 6c 61 74 65 64 20 66 |ma and r|elated f|
|00006220| 75 6e 63 74 69 6f 6e 73 | 0a 0a 09 6c 61 73 74 20 |unctions|...last |
|00006230| 65 64 69 74 3a 09 38 38 | 2f 30 39 2f 30 39 09 44 |edit:.88|/09/09.D|
|00006240| 20 41 20 47 77 79 6e 0a | 0a 09 53 43 43 53 20 49 | A Gwyn.|..SCCS I|
|00006250| 44 3a 09 40 28 23 29 67 | 61 6d 6d 61 2e 63 09 31 |D:.@(#)g|amma.c.1|
|00006260| 2e 31 20 28 65 64 69 74 | 65 64 20 66 6f 72 20 70 |.1 (edit|ed for p|
|00006270| 75 62 6c 69 63 61 74 69 | 6f 6e 29 0a 0a 41 63 6b |ublicati|on)..Ack|
|00006280| 6e 6f 77 6c 65 64 67 65 | 6d 65 6e 74 3a 0a 09 43 |nowledge|ment:..C|
|00006290| 6f 64 65 20 62 61 73 65 | 64 20 6f 6e 20 74 68 61 |ode base|d on tha|
|000062a0| 74 20 66 6f 75 6e 64 20 | 69 6e 20 22 4e 75 6d 65 |t found |in "Nume|
|000062b0| 72 69 63 61 6c 20 4d 65 | 74 68 6f 64 73 20 69 6e |rical Me|thods in|
|000062c0| 20 43 22 2e 0a 2a 2f 0a | 0a 23 69 6e 63 6c 75 64 | C"..*/.|.#includ|
|000062d0| 65 09 3c 61 73 73 65 72 | 74 2e 68 3e 0a 23 69 6e |e.<asser|t.h>.#in|
|000062e0| 63 6c 75 64 65 09 3c 6d | 61 74 68 2e 68 3e 0a 0a |clude.<m|ath.h>..|
|000062f0| 23 69 6e 63 6c 75 64 65 | 09 22 73 74 64 2e 68 22 |#include|."std.h"|
|00006300| 0a 0a 64 6f 75 62 6c 65 | 0a 4c 47 61 6d 6d 61 28 |..double|.LGamma(|
|00006310| 20 78 20 29 0a 09 64 6f | 75 62 6c 65 09 09 09 78 | x )..do|uble...x|
|00006320| 3b 0a 09 7b 0a 09 73 74 | 61 74 69 63 20 63 6f 6e |;..{..st|atic con|
|00006330| 73 74 20 64 6f 75 62 6c | 65 09 63 6f 66 5b 36 5d |st doubl|e.cof[6]|
|00006340| 20 3d 0a 09 09 7b 0a 09 | 09 37 36 2e 31 38 30 30 | =...{..|.76.1800|
|00006350| 39 31 37 33 2c 09 2d 38 | 36 2e 35 30 35 33 32 30 |9173,.-8|6.505320|
|00006360| 33 33 2c 09 32 34 2e 30 | 31 34 30 39 38 32 32 2c |33,.24.0|1409822,|
|00006370| 0a 09 09 2d 31 2e 32 33 | 31 37 33 39 35 31 36 2c |...-1.23|1739516,|
|00006380| 09 30 2e 31 32 30 38 35 | 38 30 30 33 65 2d 32 2c |.0.12085|8003e-2,|
|00006390| 09 2d 30 2e 35 33 36 33 | 38 32 65 2d 35 0a 09 09 |.-0.5363|82e-5...|
|000063a0| 7d 3b 0a 09 64 6f 75 62 | 6c 65 09 09 09 74 6d 70 |};..doub|le...tmp|
|000063b0| 2c 20 73 65 72 3b 0a 09 | 72 65 67 69 73 74 65 72 |, ser;..|register|
|000063c0| 20 69 6e 74 09 09 6a 3b | 0a 0a 09 61 73 73 65 72 | int..j;|...asser|
|000063d0| 74 28 78 20 3e 20 30 2e | 30 29 3b 0a 0a 09 69 66 |t(x > 0.|0);...if|
|000063e0| 20 28 20 2d 2d 78 20 3c | 20 30 2e 30 20 29 09 2f | ( --x <| 0.0 )./|
|000063f0| 2a 20 75 73 65 20 72 65 | 66 6c 65 63 74 69 6f 6e |* use re|flection|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.