January 19 2008 3:39:36.783 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.00909 2 1.13330 3 1.00002 4 1.00356 5 1.00405 6 1.00055 7 1.00191 8 1.00029 9 1.84722 Estimated component weights: 1 0.100909 2 0.113330 3 0.100002 4 0.100356 5 0.100405 6 0.100055 7 0.100191 8 0.100029 9 0.184722 Event 2 Current ALPHA: 1 1.03689 2 1.13761 3 1.00006 4 1.00593 5 1.00482 6 1.00091 7 1.00815 8 1.00039 9 2.80524 Estimated component weights: 1 0.942623E-01 2 0.103419 3 0.909143E-01 4 0.914486E-01 5 0.913474E-01 6 0.909917E-01 7 0.916499E-01 8 0.909444E-01 9 0.255022 Event 3 Current ALPHA: 1 1.03699 2 2.09690 3 1.00009 4 1.00603 5 1.03850 6 1.00134 7 1.00818 8 1.00626 9 2.80570 Estimated component weights: 1 0.864161E-01 2 0.174742 3 0.833406E-01 4 0.838356E-01 5 0.865416E-01 6 0.834449E-01 7 0.840153E-01 8 0.838553E-01 9 0.233809 Event 4 Current ALPHA: 1 1.06836 2 2.60555 3 1.02478 4 1.11397 5 1.09847 6 1.00433 7 1.12726 8 1.10365 9 2.85363 Estimated component weights: 1 0.821818E-01 2 0.200427 3 0.788291E-01 4 0.856899E-01 5 0.844980E-01 6 0.772562E-01 7 0.867120E-01 8 0.848962E-01 9 0.219510 Event 5 Current ALPHA: 1 1.18109 2 2.65531 3 1.02567 4 1.12422 5 1.70612 6 1.04570 7 1.13630 8 1.13846 9 2.98713 Estimated component weights: 1 0.843636E-01 2 0.189665 3 0.732619E-01 4 0.803013E-01 5 0.121866 6 0.746925E-01 7 0.811646E-01 8 0.813185E-01 9 0.213367 Event 6 Current ALPHA: 1 1.43616 2 2.68523 3 1.02607 4 1.13152 5 1.70765 6 1.04660 7 1.16223 8 1.13908 9 3.66546 Estimated component weights: 1 0.957439E-01 2 0.179015 3 0.684048E-01 4 0.754344E-01 5 0.113843 6 0.697732E-01 7 0.774821E-01 8 0.759389E-01 9 0.244364 Event 7 Current ALPHA: 1 2.01231 2 2.72692 3 1.02968 4 1.14839 5 1.81059 6 1.07800 7 1.17896 8 1.17417 9 3.84099 Estimated component weights: 1 0.125769 2 0.170432 3 0.643550E-01 4 0.717741E-01 5 0.113162 6 0.673751E-01 7 0.736851E-01 8 0.733856E-01 9 0.240062 Event 8 Current ALPHA: 1 2.16620 2 2.73860 3 1.04518 4 1.17560 5 1.81119 6 1.07861 7 1.76960 8 1.17595 9 4.03908 Estimated component weights: 1 0.127423 2 0.161094 3 0.614811E-01 4 0.691527E-01 5 0.106541 6 0.634475E-01 7 0.104094 8 0.691736E-01 9 0.237593 Event 9 Current ALPHA: 1 2.17635 2 2.73891 3 1.46931 4 1.44056 5 1.81120 6 1.07861 7 2.06686 8 1.17846 9 4.03974 Estimated component weights: 1 0.120908 2 0.152161 3 0.816282E-01 4 0.800311E-01 5 0.100622 6 0.599231E-01 7 0.114826 8 0.654701E-01 9 0.224430 Event 10 Current ALPHA: 1 2.18157 2 3.00618 3 1.47620 4 1.49938 5 2.38469 6 1.08504 7 2.07113 8 1.24279 9 4.05301 Estimated component weights: 1 0.114820 2 0.158220 3 0.776948E-01 4 0.789148E-01 5 0.125510 6 0.571075E-01 7 0.109007 8 0.654102E-01 9 0.213316 Event 100 Current ALPHA: 1 16.4615 2 10.1909 3 6.61620 4 5.66023 5 10.5664 6 9.59402 7 16.5402 8 7.35014 9 26.0203 Estimated component weights: 1 0.151023 2 0.934948E-01 3 0.606991E-01 4 0.519287E-01 5 0.969395E-01 6 0.880186E-01 7 0.151745 8 0.674324E-01 9 0.238719 Event 200 Current ALPHA: 1 32.1841 2 15.3834 3 16.2055 4 15.9438 5 23.2464 6 17.7950 7 27.2256 8 14.5857 9 46.4305 Estimated component weights: 1 0.153991 2 0.736046E-01 3 0.775381E-01 4 0.762862E-01 5 0.111227 6 0.851436E-01 7 0.130266 8 0.697881E-01 9 0.222156 Event 300 Current ALPHA: 1 49.7323 2 23.2261 3 23.0306 4 22.1708 5 36.7206 6 25.0658 7 39.9202 8 18.4870 9 70.6465 Estimated component weights: 1 0.160946 2 0.751654E-01 3 0.745328E-01 4 0.717502E-01 5 0.118837 6 0.811190E-01 7 0.129192 8 0.598285E-01 9 0.228630 Event 400 Current ALPHA: 1 71.2324 2 28.9674 3 31.7188 4 32.7040 5 44.4116 6 32.2053 7 49.2361 8 20.7530 9 97.7715 Estimated component weights: 1 0.174162 2 0.708249E-01 3 0.775520E-01 4 0.799608E-01 5 0.108586 6 0.787416E-01 7 0.120382 8 0.507407E-01 9 0.239050 Event 500 Current ALPHA: 1 87.3646 2 36.6847 3 42.5902 4 43.4131 5 51.5568 6 37.2513 7 60.9027 8 25.9483 9 123.288 Estimated component weights: 1 0.171640 2 0.720722E-01 3 0.836742E-01 4 0.852910E-01 5 0.101290 6 0.731853E-01 7 0.119652 8 0.509789E-01 9 0.242217 Event 600 Current ALPHA: 1 106.768 2 44.9311 3 50.5822 4 53.2200 5 57.0394 6 44.0028 7 69.9109 8 30.8151 9 151.731 Estimated component weights: 1 0.175317 2 0.737784E-01 3 0.830578E-01 4 0.873891E-01 5 0.936608E-01 6 0.722542E-01 7 0.114796 8 0.505995E-01 9 0.249147 Event 700 Current ALPHA: 1 126.510 2 47.9069 3 58.0739 4 61.9788 5 66.7178 6 56.0331 7 78.0092 8 37.6520 9 176.118 Estimated component weights: 1 0.178435 2 0.675697E-01 3 0.819095E-01 4 0.874172E-01 5 0.941013E-01 6 0.790312E-01 7 0.110027 8 0.531058E-01 9 0.248403 Event 800 Current ALPHA: 1 146.237 2 51.4601 3 65.7884 4 72.0842 5 73.9994 6 67.0908 7 91.8561 8 43.7551 9 196.729 Estimated component weights: 1 0.180762 2 0.636095E-01 3 0.813206E-01 4 0.891028E-01 5 0.914702E-01 6 0.829306E-01 7 0.113543 8 0.540854E-01 9 0.243176 Event 900 Current ALPHA: 1 165.409 2 58.1584 3 74.4548 4 78.0365 5 82.6897 6 75.7129 7 102.031 8 47.8523 9 224.656 Estimated component weights: 1 0.181968 2 0.639806E-01 3 0.819085E-01 4 0.858488E-01 5 0.909678E-01 6 0.832926E-01 7 0.112245 8 0.526428E-01 9 0.247146 Event 1000 Current ALPHA: 1 180.999 2 65.4581 3 78.4854 4 85.8808 5 93.1849 6 86.9078 7 111.528 8 55.8917 9 250.664 Estimated component weights: 1 0.179385 2 0.648743E-01 3 0.777853E-01 4 0.851148E-01 5 0.923537E-01 6 0.861327E-01 7 0.110533 8 0.553932E-01 9 0.248428 BDMLIB_PRB Normal end of execution. January 19 2008 3:39:36.990 PM