June 24 2002 3:47:45.060 PM BDMLIB_PRB Test BDMLIB, Bayes Dirichlet Mixture estimate weights of a Dirichlet mixture. MIXTURE_READ - Note: Number of records read was 38 Number of goofy records was 0 Numeric key for amino acid abbreviations. I Amino Acid Symbol 1 A Alanine 2 C Cysteine 3 D Aspartic acid 4 E Glutamic acid 5 F Phenylalanine 6 G Glycine 7 H Histidine 8 I Isoleucine 9 K Lysine 10 L Leucine 11 M Methionine 12 N Asparagine 13 P Proline 14 Q Glutamine 15 R Arginine 16 S Serine 17 T Threonine 18 V Valine 19 W Tryptophan 20 Y Tyrosine Number of components = 9 Compon: 1 2 3 4 5 6 7 8 9 Weight: 0.1781 0.0566 0.0960 0.0781 0.0835 0.0904 0.1145 0.0682 0.2346 1 A 0.2707 0.0215 0.5615 0.0701 0.0411 0.1156 0.0935 0.4522 0.0052 2 C 0.0398 0.0103 0.0454 0.0111 0.0148 0.0374 0.0047 0.1146 0.0040 3 D 0.0176 0.0117 0.4384 0.0195 0.0056 0.0124 0.3873 0.0625 0.0067 4 E 0.0164 0.0109 0.7642 0.0947 0.0102 0.0182 0.3478 0.1157 0.0061 5 F 0.0143 0.3857 0.0874 0.0132 0.1536 0.0518 0.0108 0.2842 0.0035 6 G 0.1319 0.0164 0.2591 0.0480 0.0078 0.0173 0.1059 0.1402 0.0169 7 H 0.0124 0.0762 0.2149 0.0770 0.0072 0.0049 0.0498 0.1004 0.0036 8 I 0.0226 0.0353 0.1459 0.0329 0.2996 0.7969 0.0150 0.5502 0.0022 9 K 0.0204 0.0139 0.7622 0.5766 0.0108 0.0171 0.0943 0.1440 0.0050 10 L 0.0307 0.0935 0.2473 0.0723 0.9994 0.2859 0.0278 0.7006 0.0060 11 M 0.0153 0.0220 0.1187 0.0282 0.2102 0.0758 0.0100 0.2766 0.0015 12 N 0.0483 0.0286 0.4416 0.0804 0.0061 0.0145 0.1879 0.1186 0.0042 13 P 0.0538 0.0131 0.1748 0.0377 0.0130 0.0151 0.0500 0.0975 0.0091 14 Q 0.0207 0.0230 0.5308 0.1850 0.0198 0.0114 0.1100 0.1267 0.0036 15 R 0.0236 0.0189 0.4655 0.5068 0.0145 0.0127 0.0387 0.1436 0.0066 16 S 0.2161 0.0292 0.5834 0.0737 0.0120 0.0275 0.1195 0.2790 0.0032 17 T 0.1472 0.0182 0.4456 0.0716 0.0358 0.0883 0.0658 0.3585 0.0037 18 V 0.0654 0.0361 0.2271 0.0425 0.1801 0.9443 0.0254 0.6618 0.0030 19 W 0.0038 0.0718 0.0295 0.0113 0.0127 0.0044 0.0032 0.0615 0.0028 20 Y 0.0096 0.4196 0.1211 0.0287 0.0265 0.0167 0.0187 0.1994 0.0027 Sum 1.1806 1.3558 6.6644 2.0814 2.0810 2.5682 1.7661 4.9877 0.0995 TEST01 Generate a (nonrandom) sequence of 10 isoleucine results in a row. Exact component weights: 1 0.178091 2 0.565910E-01 3 0.960191E-01 4 0.781233E-01 5 0.834977E-01 6 0.904123E-01 7 0.114468 8 0.682132E-01 9 0.234585 Initial ALPHA: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 Initial estimated component weights: 1 0.111111 2 0.111111 3 0.111111 4 0.111111 5 0.111111 6 0.111111 7 0.111111 8 0.111111 9 0.111111 Current ALPHA: 1 1.02823 2 1.03844 3 1.03230 4 1.02334 5 1.21239 6 1.45770 7 1.01250 8 1.16273 9 1.03238 Estimated component weights: 1 0.102823 2 0.103844 3 0.103230 4 0.102334 5 0.121239 6 0.145770 7 0.101250 8 0.116273 9 0.103238 Current ALPHA: 1 1.05081 2 1.06947 3 1.05822 4 1.04192 5 1.41260 6 1.97645 7 1.02234 8 1.30984 9 1.05837 Estimated component weights: 1 0.955279E-01 2 0.972245E-01 3 0.962021E-01 4 0.947197E-01 5 0.128418 6 0.179677 7 0.929396E-01 8 0.119076 9 0.962151E-01 Current ALPHA: 1 1.06942 2 1.09526 3 1.07966 4 1.05717 5 1.60080 6 2.54392 7 1.03035 8 1.44355 9 1.07986 Estimated component weights: 1 0.891182E-01 2 0.912713E-01 3 0.899721E-01 4 0.880979E-01 5 0.133400 6 0.211994 7 0.858626E-01 8 0.120296 9 0.899886E-01 Current ALPHA: 1 1.08515 2 1.11719 3 1.09784 4 1.07003 5 1.77798 6 3.15069 7 1.03706 8 1.56596 9 1.09808 Estimated component weights: 1 0.834734E-01 2 0.859380E-01 3 0.844491E-01 4 0.823104E-01 5 0.136768 6 0.242361 7 0.797740E-01 8 0.120459 9 0.844679E-01 Current ALPHA: 1 1.09873 2 1.13622 3 1.11355 4 1.08110 5 1.94529 6 3.78961 7 1.04280 8 1.67886 9 1.11383 Estimated component weights: 1 0.784806E-01 2 0.811585E-01 3 0.795392E-01 4 0.772216E-01 5 0.138949 6 0.270687 7 0.744860E-01 8 0.119919 9 0.795596E-01 Current ALPHA: 1 1.11063 2 1.15298 3 1.12735 4 1.09079 5 2.10384 6 4.45523 7 1.04781 8 1.78370 9 1.12767 Estimated component weights: 1 0.740422E-01 2 0.768652E-01 3 0.751567E-01 4 0.727191E-01 5 0.140256 6 0.297016 7 0.698537E-01 8 0.118913 9 0.751783E-01 Current ALPHA: 1 1.12122 2 1.16793 3 1.13964 4 1.09938 5 2.25462 6 5.14335 7 1.05222 8 1.88165 9 1.13999 Estimated component weights: 1 0.700760E-01 2 0.729958E-01 3 0.712274E-01 4 0.687112E-01 5 0.140914 6 0.321459 7 0.657640E-01 8 0.117603 9 0.712497E-01 Current ALPHA: 1 1.13073 2 1.18142 3 1.15070 4 1.10709 5 2.39850 6 5.85066 7 1.05618 8 1.97365 9 1.15109 Estimated component weights: 1 0.665134E-01 2 0.694953E-01 3 0.676881E-01 4 0.651229E-01 5 0.141088 6 0.344157 7 0.621279E-01 8 0.116097 9 0.677109E-01 Current ALPHA: 1 1.13936 2 1.19370 3 1.16074 4 1.11408 5 2.53620 6 6.57455 7 1.05974 8 2.06047 9 1.16116 Estimated component weights: 1 0.632976E-01 2 0.663164E-01 3 0.644858E-01 4 0.618931E-01 5 0.140900 6 0.365253 7 0.588746E-01 8 0.114470 9 0.645089E-01 Current ALPHA: 1 1.14725 2 1.20495 3 1.16994 4 1.12046 5 2.66837 6 7.31290 7 1.06299 8 2.14274 9 1.17038 Estimated component weights: 1 0.603816E-01 2 0.634186E-01 3 0.615760E-01 4 0.589714E-01 5 0.140441 6 0.384890 7 0.559470E-01 8 0.112776 9 0.615992E-01 TEST02 Generate many random events. We should be able to approximate the exact component weights. Exact component weights: 1 0.178091 2 0.565910E-01 3 0.960191E-01 4 0.781233E-01 5 0.834977E-01 6 0.904123E-01 7 0.114468 8 0.682132E-01 9 0.234585 Initial ALPHA: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 Initial estimated component weights: 1 0.111111 2 0.111111 3 0.111111 4 0.111111 5 0.111111 6 0.111111 7 0.111111 8 0.111111 9 0.111111 Event 1 Current ALPHA: 1 1.01674 2 1.04149 3 1.44309 4 1.00270 5 1.00268 6 1.00083 7 1.22195 8 1.27049 9 1.00004 Estimated component weights: 1 0.101674 2 0.104149 3 0.144309 4 0.100270 5 0.100268 6 0.100083 7 0.122195 8 0.127049 9 0.100004 Event 2 Current ALPHA: 1 1.13158 2 1.46980 3 1.44480 4 1.01577 5 1.00444 6 1.00184 7 1.27219 8 1.27693 9 1.38265 Estimated component weights: 1 0.102871 2 0.133619 3 0.131346 4 0.923429E-01 5 0.913129E-01 6 0.910762E-01 7 0.115653 8 0.116085 9 0.125695 Event 3 Current ALPHA: 1 1.14529 2 1.47139 3 1.73046 4 1.62874 5 1.00717 6 1.00235 7 1.28294 8 1.34882 9 1.38285 Estimated component weights: 1 0.954408E-01 2 0.122616 3 0.144205 4 0.135728 5 0.839308E-01 6 0.835289E-01 7 0.106912 8 0.112402 9 0.115238 Event 4 Current ALPHA: 1 1.15060 2 1.47160 3 1.78379 4 1.63095 5 1.00717 6 1.00237 7 2.22002 8 1.34928 9 1.38421 Estimated component weights: 1 0.885080E-01 2 0.113200 3 0.137214 4 0.125458 5 0.774750E-01 6 0.771051E-01 7 0.170770 8 0.103791 9 0.106478 Event 5 Current ALPHA: 1 1.15746 2 2.16491 3 1.78384 4 1.63488 5 1.00939 6 1.00292 7 2.22651 8 1.34975 9 1.67034 Estimated component weights: 1 0.826757E-01 2 0.154637 3 0.127417 4 0.116777 5 0.720996E-01 6 0.716369E-01 7 0.159037 8 0.964108E-01 9 0.119310 Event 6 Current ALPHA: 1 1.16665 2 2.30520 3 1.79027 4 2.15978 5 1.01120 6 1.00332 7 2.24743 8 1.35263 9 1.96353 Estimated component weights: 1 0.777767E-01 2 0.153680 3 0.119351 4 0.143985 5 0.674137E-01 6 0.668877E-01 7 0.149829 8 0.901752E-01 9 0.130902 Event 7 Current ALPHA: 1 1.21236 2 2.32762 3 1.79369 4 2.17108 5 1.27295 6 1.44374 7 2.25385 8 1.54940 9 1.97532 Estimated component weights: 1 0.757724E-01 2 0.145476 3 0.112106 4 0.135692 5 0.795592E-01 6 0.902338E-01 7 0.140866 8 0.968373E-01 9 0.123457 Event 8 Current ALPHA: 1 1.25622 2 2.33104 3 1.79598 4 2.17448 5 1.27303 6 1.44393 7 2.86109 8 1.54979 9 2.31444 Estimated component weights: 1 0.738953E-01 2 0.137120 3 0.105646 4 0.127911 5 0.748841E-01 6 0.849373E-01 7 0.168300 8 0.911641E-01 9 0.136143 Event 9 Current ALPHA: 1 1.40060 2 2.33755 3 1.82438 4 2.75631 5 1.27412 6 1.44643 7 2.98325 8 1.55562 9 2.42173 Estimated component weights: 1 0.778113E-01 2 0.129864 3 0.101354 4 0.153128 5 0.707846E-01 6 0.803575E-01 7 0.165736 8 0.864232E-01 9 0.134541 Event 10 Current ALPHA: 1 1.42010 2 2.39690 3 1.82636 4 3.03346 5 1.65044 6 1.45633 7 2.99824 8 1.57919 9 2.63897 Estimated component weights: 1 0.747422E-01 2 0.126153 3 0.961244E-01 4 0.159656 5 0.868654E-01 6 0.766490E-01 7 0.157802 8 0.831153E-01 9 0.138893 Event 100 Current ALPHA: 1 10.2818 2 12.3647 3 14.0043 4 10.9758 5 6.32278 6 9.68196 7 13.5657 8 10.8210 9 20.9820 Estimated component weights: 1 0.943287E-01 2 0.113438 3 0.128480 4 0.100695 5 0.580072E-01 6 0.888253E-01 7 0.124456 8 0.992752E-01 9 0.192495 Event 200 Current ALPHA: 1 32.7519 2 22.7969 3 21.0525 4 16.2798 5 14.1537 6 17.7913 7 23.9892 8 17.1589 9 43.0257 Estimated component weights: 1 0.156708 2 0.109076 3 0.100730 4 0.778939E-01 5 0.677211E-01 6 0.851258E-01 7 0.114781 8 0.821002E-01 9 0.205864 Event 300 Current ALPHA: 1 49.9672 2 28.3626 3 31.4624 4 23.7526 5 21.3108 6 29.2743 7 37.7077 8 24.0932 9 63.0692 Estimated component weights: 1 0.161706 2 0.917882E-01 3 0.101820 4 0.768692E-01 5 0.689669E-01 6 0.947389E-01 7 0.122031 8 0.779716E-01 9 0.204107 RVEC_UNIT_SUM - Fatal error! The vector entries sum to 0.