September 20 2010 12:49:47.641 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.6617 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.07967 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.15070 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.772215E-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.15108 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.621280E-01 8 0.116097 9 0.677109E-01 Current ALPHA: 1 1.13936 2 1.19370 3 1.16075 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.31291 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.00889 2 1.95747 3 1.00080 4 1.00060 5 1.00438 6 1.00103 7 1.00084 8 1.02507 9 1.00092 Estimated component weights: 1 0.100889 2 0.195747 3 0.100080 4 0.100060 5 0.100438 6 0.100103 7 0.100084 8 0.102507 9 0.100092 Event 2 Current ALPHA: 1 1.03644 2 1.95887 3 1.01036 4 1.00270 5 1.00439 6 1.00107 7 1.94756 8 1.02556 9 1.01305 Estimated component weights: 1 0.942220E-01 2 0.178079 3 0.918514E-01 4 0.911548E-01 5 0.913078E-01 6 0.910061E-01 7 0.177051 8 0.932325E-01 9 0.920958E-01 Event 3 Current ALPHA: 1 1.03703 2 1.96060 3 1.01041 4 1.00286 5 1.11958 6 1.86300 7 1.94761 8 1.04580 9 1.01311 Estimated component weights: 1 0.864195E-01 2 0.163383 3 0.842007E-01 4 0.835714E-01 5 0.932980E-01 6 0.155250 7 0.162301 8 0.871500E-01 9 0.844256E-01 Event 4 Current ALPHA: 1 1.11322 2 1.96892 3 1.01047 4 1.00560 5 1.12002 6 1.86368 7 1.97395 8 1.04594 9 1.89820 Estimated component weights: 1 0.856320E-01 2 0.151456 3 0.777282E-01 4 0.773537E-01 5 0.861554E-01 6 0.143360 7 0.151842 8 0.804569E-01 9 0.146016 Event 5 Current ALPHA: 1 1.38131 2 1.96920 3 1.18714 4 1.01424 5 1.12003 6 1.86385 7 2.50704 8 1.05833 9 1.89884 Estimated component weights: 1 0.986652E-01 2 0.140657 3 0.847960E-01 4 0.724455E-01 5 0.800023E-01 6 0.133132 7 0.179075 8 0.755952E-01 9 0.135632 Event 6 Current ALPHA: 1 1.93038 2 1.98114 3 1.18938 4 1.03000 5 1.12277 6 1.86941 7 2.60070 8 1.06245 9 2.21377 Estimated component weights: 1 0.128692 2 0.132076 3 0.792919E-01 4 0.686670E-01 5 0.748511E-01 6 0.124627 7 0.173380 8 0.708300E-01 9 0.147585 Event 7 Current ALPHA: 1 2.00946 2 2.07186 3 1.18967 4 1.03475 5 1.12362 6 1.86986 7 2.61753 8 1.06375 9 3.01950 Estimated component weights: 1 0.125591 2 0.129491 3 0.743544E-01 4 0.646719E-01 5 0.702263E-01 6 0.116866 7 0.163595 8 0.664844E-01 9 0.188719 Event 8 Current ALPHA: 1 2.48874 2 2.11929 3 1.19297 4 1.04700 5 1.15246 6 2.13253 7 2.66365 8 1.10184 9 3.10152 Estimated component weights: 1 0.146396 2 0.124664 3 0.701747E-01 4 0.615884E-01 5 0.677917E-01 6 0.125443 7 0.156685 8 0.648140E-01 9 0.182443 Event 9 Current ALPHA: 1 2.48929 2 3.00345 3 1.19319 4 1.04711 5 1.19616 6 2.14428 7 2.66370 8 1.16127 9 3.10155 Estimated component weights: 1 0.138294 2 0.166858 3 0.662881E-01 4 0.581726E-01 5 0.664533E-01 6 0.119127 7 0.147984 8 0.645149E-01 9 0.172308 Event 10 Current ALPHA: 1 3.41143 2 3.00679 3 1.19533 4 1.05082 5 1.19725 6 2.15135 7 2.68668 8 1.16709 9 3.13325 Estimated component weights: 1 0.179549 2 0.158252 3 0.629121E-01 4 0.553064E-01 5 0.630132E-01 6 0.113229 7 0.141404 8 0.614260E-01 9 0.164908 Event 100 Current ALPHA: 1 15.2750 2 13.4623 3 11.9291 4 8.55082 5 11.9915 6 12.5468 7 15.0678 8 6.66970 9 13.5070 Estimated component weights: 1 0.140138 2 0.123508 3 0.109441 4 0.784479E-01 5 0.110013 6 0.115108 7 0.138236 8 0.611899E-01 9 0.123918 Event 200 Current ALPHA: 1 30.9418 2 23.4248 3 22.5598 4 14.9013 5 18.6640 6 23.1332 7 25.8691 8 13.3409 9 36.1651 Estimated component weights: 1 0.148047 2 0.112080 3 0.107942 4 0.712981E-01 5 0.893017E-01 6 0.110685 7 0.123776 8 0.638320E-01 9 0.173039 Event 300 Current ALPHA: 1 44.3985 2 32.0377 3 32.5453 4 22.0773 5 26.7846 6 35.1134 7 37.5417 8 17.8121 9 60.6895 Estimated component weights: 1 0.143685 2 0.103682 3 0.105324 4 0.714475E-01 5 0.866814E-01 6 0.113636 7 0.121494 8 0.576442E-01 9 0.196406 Event 400 Current ALPHA: 1 57.0163 2 38.6226 3 42.4387 4 28.1958 5 37.2932 6 48.3884 7 48.6541 8 23.4961 9 84.8949 Estimated component weights: 1 0.139404 2 0.944317E-01 3 0.103762 4 0.689384E-01 5 0.911815E-01 6 0.118309 7 0.118959 8 0.574476E-01 9 0.207567 Event 500 Current ALPHA: 1 72.8374 2 45.4876 3 54.7690 4 34.8399 5 49.7125 6 53.1704 7 65.9538 8 28.9174 9 103.312 Estimated component weights: 1 0.143099 2 0.893666E-01 3 0.107601 4 0.684477E-01 5 0.976670E-01 6 0.104461 7 0.129575 8 0.568122E-01 9 0.202970 Event 600 Current ALPHA: 1 86.0096 2 52.6888 3 65.6349 4 40.0380 5 64.6481 6 63.3811 7 80.5213 8 34.1726 9 121.905 Estimated component weights: 1 0.141231 2 0.865169E-01 3 0.107775 4 0.657438E-01 5 0.106155 6 0.104074 7 0.132219 8 0.561127E-01 9 0.200173 Event 700 Current ALPHA: 1 102.805 2 59.8921 3 75.4902 4 44.1746 5 73.0185 6 70.4547 7 97.7034 8 43.5189 9 141.943 Estimated component weights: 1 0.144999 2 0.844740E-01 3 0.106474 4 0.623055E-01 5 0.102988 6 0.993719E-01 7 0.137805 8 0.613807E-01 9 0.200202 Event 800 Current ALPHA: 1 122.624 2 64.5964 3 89.7826 4 50.3531 5 77.1577 6 77.7877 7 113.562 8 48.2484 9 164.888 Estimated component weights: 1 0.151575 2 0.798472E-01 3 0.110980 4 0.622411E-01 5 0.953742E-01 6 0.961529E-01 7 0.140373 8 0.596395E-01 9 0.203817 Event 900 Current ALPHA: 1 134.940 2 71.5775 3 100.358 4 56.2442 5 87.0154 6 87.4436 7 128.885 8 53.4260 9 189.110 Estimated component weights: 1 0.148449 2 0.787432E-01 3 0.110405 4 0.618748E-01 5 0.957265E-01 6 0.961975E-01 7 0.141788 8 0.587745E-01 9 0.208042 Event 1000 Current ALPHA: 1 149.355 2 76.9761 3 113.793 4 62.9874 5 92.6268 6 95.1730 7 142.255 8 57.6758 9 218.159 Estimated component weights: 1 0.148022 2 0.762895E-01 3 0.112778 4 0.624255E-01 5 0.918006E-01 6 0.943241E-01 7 0.140986 8 0.571614E-01 9 0.216213 BDMLIB_PRB Normal end of execution. September 20 2010 12:49:47.852 PM