9 September 2023 7:17:34.553 PM bdmlib_test(): FORTRAN77 version Test bdmlib(). 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.17809100 2: 0.56591000E-01 3: 0.96019100E-01 4: 0.78123300E-01 5: 0.83497700E-01 6: 0.90412300E-01 7: 0.11446800 8: 0.68213200E-01 9: 0.23458500 Initial ALPHA: 1: 1.0000000 2: 1.0000000 3: 1.0000000 4: 1.0000000 5: 1.0000000 6: 1.0000000 7: 1.0000000 8: 1.0000000 9: 1.0000000 Initial estimated component weights: 1: 0.11111111 2: 0.11111111 3: 0.11111111 4: 0.11111111 5: 0.11111111 6: 0.11111111 7: 0.11111111 8: 0.11111111 9: 0.11111111 Current ALPHA: 1: 1.0282345 2: 1.0384360 3: 1.0322992 4: 1.0233434 5: 1.2123878 6: 1.4576978 7: 1.0124976 8: 1.1627264 9: 1.0323774 Estimated component weights: 1: 0.10282345 2: 0.10384360 3: 0.10322992 4: 0.10233434 5: 0.12123878 6: 0.14576978 7: 0.10124976 8: 0.11627264 9: 0.10323774 Current ALPHA: 1: 1.0508072 2: 1.0694694 3: 1.0582236 4: 1.0419170 5: 1.4125959 6: 1.9764465 7: 1.0223361 8: 1.3098378 9: 1.0583664 Estimated component weights: 1: 0.95527928E-01 2: 0.97224489E-01 3: 0.96202143E-01 4: 0.94719731E-01 5: 0.12841781 6: 0.17967696 7: 0.92939644E-01 8: 0.11907617 9: 0.96215130E-01 Current ALPHA: 1: 1.0694190 2: 1.0952558 3: 1.0796650 4: 1.0571745 5: 1.6008013 6: 2.5439239 7: 1.0303511 8: 1.4435467 9: 1.0798627 Estimated component weights: 1: 0.89118249E-01 2: 0.91271318E-01 3: 0.89972084E-01 4: 0.88097876E-01 5: 0.13340011 6: 0.21199366 7: 0.85862590E-01 8: 0.12029556 9: 0.89988555E-01 Current ALPHA: 1: 1.0851541 2: 1.1171938 3: 1.0978379 4: 1.0700349 5: 1.7779790 6: 3.1506953 7: 1.0370615 8: 1.5659607 9: 1.0980828 Estimated component weights: 1: 0.83473394E-01 2: 0.85937983E-01 3: 0.84449066E-01 4: 0.82310375E-01 5: 0.13676762 6: 0.24236117 7: 0.79773965E-01 8: 0.12045852 9: 0.84467909E-01 Current ALPHA: 1: 1.0987289 2: 1.1362189 3: 1.1135484 4: 1.0811017 5: 1.9452870 6: 3.7896142 7: 1.0428039 8: 1.6788621 9: 1.1138349 Estimated component weights: 1: 0.78480635E-01 2: 0.81158492E-01 3: 0.79539169E-01 4: 0.77221550E-01 5: 0.13894907 6: 0.27068673 7: 0.74485993E-01 8: 0.11991872 9: 0.79559634E-01 Current ALPHA: 1: 1.1106337 2: 1.1529781 3: 1.1273507 4: 1.0907863 5: 2.1038367 6: 4.4552335 7: 1.0478052 8: 1.7837017 9: 1.1276742 Estimated component weights: 1: 0.74042248E-01 2: 0.76865204E-01 3: 0.75156715E-01 4: 0.72719089E-01 5: 0.14025578 6: 0.29701556 7: 0.69853678E-01 8: 0.11891345 9: 0.75178279E-01 Current ALPHA: 1: 1.1212156 2: 1.1679325 3: 1.1396382 4: 1.0993787 5: 2.2546197 6: 5.1433475 7: 1.0522241 8: 1.8816488 9: 1.1399949 Estimated component weights: 1: 0.70075974E-01 2: 0.72995782E-01 3: 0.71227386E-01 4: 0.68711172E-01 5: 0.14091373 6: 0.32145922 7: 0.65764006E-01 8: 0.11760305 9: 0.71249682E-01 Current ALPHA: 1: 1.1307272 2: 1.1814204 3: 1.1506979 4: 1.1070895 5: 2.3984962 6: 5.8506607 7: 1.0561752 8: 1.9736479 9: 1.1510849 Estimated component weights: 1: 0.66513368E-01 2: 0.69495316E-01 3: 0.67688114E-01 4: 0.65122914E-01 5: 0.14108801 6: 0.34415651 7: 0.62127954E-01 8: 0.11609693 9: 0.67710877E-01 Current ALPHA: 1: 1.1393576 2: 1.1936957 3: 1.1607451 4: 1.1140757 5: 2.5362039 6: 6.5745516 7: 1.0597434 8: 2.0604674 9: 1.1611597 Estimated component weights: 1: 0.63297644E-01 2: 0.66316427E-01 3: 0.64485838E-01 4: 0.61893093E-01 5: 0.14090021 6: 0.36525287 7: 0.58874635E-01 8: 0.11447041 9: 0.64508875E-01 Current ALPHA: 1: 1.1472509 2: 1.2049534 3: 1.1699442 4: 1.1204568 5: 2.6683739 6: 7.3129054 7: 1.0629931 8: 2.1427377 9: 1.1703845 Estimated component weights: 1: 0.60381627E-01 2: 0.63418602E-01 3: 0.61576013E-01 4: 0.58971411E-01 5: 0.14044073 6: 0.38488976 7: 0.55947008E-01 8: 0.11277567 9: 0.61599182E-01 TEST02 Generate many random events. We should be able to approximate the exact component weights. Exact component weights: 1: 0.17809100 2: 0.56591000E-01 3: 0.96019100E-01 4: 0.78123300E-01 5: 0.83497700E-01 6: 0.90412300E-01 7: 0.11446800 8: 0.68213200E-01 9: 0.23458500 Initial ALPHA: 1: 1.0000000 2: 1.0000000 3: 1.0000000 4: 1.0000000 5: 1.0000000 6: 1.0000000 7: 1.0000000 8: 1.0000000 9: 1.0000000 Initial estimated component weights: 1: 0.11111111 2: 0.11111111 3: 0.11111111 4: 0.11111111 5: 0.11111111 6: 0.11111111 7: 0.11111111 8: 0.11111111 9: 0.11111111 Event 1 Current ALPHA: 1: 1.0088928 2: 1.9574664 3: 1.0007951 4: 1.0006020 5: 1.0043787 6: 1.0010339 7: 1.0008439 8: 1.0250703 9: 1.0009170 Estimated component weights: 1: 0.10088928 2: 0.19574664 3: 0.10007951 4: 0.10006020 5: 0.10043787 6: 0.10010339 7: 0.10008439 8: 0.10250703 9: 0.10009170 Event 2 Current ALPHA: 1: 1.0364424 2: 1.9588685 3: 1.0103649 4: 1.0027025 5: 1.0043855 6: 1.0010671 7: 1.9475574 8: 1.0255578 9: 1.0130540 Estimated component weights: 1: 0.94222034E-01 2: 0.17807895 3: 0.91851354E-01 4: 0.91154775E-01 5: 0.91307775E-01 6: 0.91006096E-01 7: 0.17705067 8: 0.93232528E-01 9: 0.92095814E-01 Event 3 Current ALPHA: 1: 1.0370338 2: 1.9606013 3: 1.0104085 4: 1.0028568 5: 1.1195763 6: 1.8630013 7: 1.9476148 8: 1.0457999 9: 1.0131073 Estimated component weights: 1: 0.86419485E-01 2: 0.16338344 3: 0.84200706E-01 4: 0.83571401E-01 5: 0.93298024E-01 6: 0.15525010 7: 0.16230123 8: 0.87149996E-01 9: 0.84425611E-01 Event 4 Current ALPHA: 1: 1.1132161 2: 1.9689239 3: 1.0104660 4: 1.0055980 5: 1.1200206 6: 1.8636803 7: 1.9739506 8: 1.0459402 9: 1.8982043 Estimated component weights: 1: 0.85632005E-01 2: 0.15145568 3: 0.77728153E-01 4: 0.77353691E-01 5: 0.86155434E-01 6: 0.14336002 7: 0.15184236 8: 0.80456941E-01 9: 0.14601571 Event 5 Current ALPHA: 1: 1.3813125 2: 1.9692006 3: 1.1871440 4: 1.0142367 5: 1.1200317 6: 1.8638531 7: 2.5070443 8: 1.0583329 9: 1.8988443 Estimated component weights: 1: 0.98665177E-01 2: 0.14065718 3: 0.84795997E-01 4: 0.72445477E-01 5: 0.80002267E-01 6: 0.13313237 7: 0.17907459 8: 0.75595208E-01 9: 0.13563173 Event 6 Current ALPHA: 1: 1.9303785 2: 1.9811445 3: 1.1893783 4: 1.0300047 5: 1.1227664 6: 1.8694087 7: 2.6006997 8: 1.0624499 9: 2.2137694 Estimated component weights: 1: 0.12869190 2: 0.13207630 3: 0.79291888E-01 4: 0.68666981E-01 5: 0.74851092E-01 6: 0.12462724 7: 0.17337998 8: 0.70829994E-01 9: 0.14758463 Event 7 Current ALPHA: 1: 2.0094623 2: 2.0718603 3: 1.1896707 4: 1.0347502 5: 1.1236207 6: 1.8698565 7: 2.6175258 8: 1.0637501 9: 3.0195034 Estimated component weights: 1: 0.12559139 2: 0.12949127 3: 0.74354421E-01 4: 0.64671886E-01 5: 0.70226296E-01 6: 0.11686603 7: 0.16359536 8: 0.66484384E-01 9: 0.18871896 Event 8 Current ALPHA: 1: 2.4887371 2: 2.1192906 3: 1.1929704 4: 1.0470026 5: 1.1524597 6: 2.1325304 7: 2.6636488 8: 1.1018379 9: 3.1015226 Estimated component weights: 1: 0.14639630 2: 0.12466415 3: 0.70174728E-01 4: 0.61588388E-01 5: 0.67791745E-01 6: 0.12544296 7: 0.15668522 8: 0.64813996E-01 9: 0.18244250 Event 9 Current ALPHA: 1: 2.4892939 2: 3.0034500 3: 1.1931858 4: 1.0471059 5: 1.1961586 6: 2.1442847 7: 2.6637049 8: 1.1612683 9: 3.1015479 Estimated component weights: 1: 0.13829411 2: 0.16685833 3: 0.66288101E-01 4: 0.58172552E-01 5: 0.66453255E-01 6: 0.11912693 7: 0.14798360 8: 0.64514904E-01 9: 0.17230822 Event 10 Current ALPHA: 1: 3.4114298 2: 3.0067878 3: 1.1953301 4: 1.0508208 5: 1.1972511 6: 2.1513549 7: 2.6866848 8: 1.1670947 9: 3.1332460 Estimated component weights: 1: 0.17954894 2: 0.15825199 3: 0.62912110E-01 4: 0.55306356E-01 5: 0.63013217E-01 6: 0.11322920 7: 0.14140446 8: 0.61426039E-01 9: 0.16490768 Event 100 Current ALPHA: 1: 15.275029 2: 13.462331 3: 11.929075 4: 8.5508245 5: 11.991471 6: 12.546767 7: 15.067756 8: 6.6697042 9: 13.507043 Estimated component weights: 1: 0.14013788 2: 0.12350762 3: 0.10944106 4: 0.78447931E-01 5: 0.11001350 6: 0.11510796 7: 0.13823629 8: 0.61189947E-01 9: 0.12391782 Event 200 Current ALPHA: 1: 30.941837 2: 23.424756 3: 22.559812 4: 14.901298 5: 18.664049 6: 23.133166 7: 25.869103 8: 13.340897 9: 36.165082 Estimated component weights: 1: 0.14804707 2: 0.11208017 3: 0.10794169 4: 0.71298078E-01 5: 0.89301671E-01 6: 0.11068501 7: 0.12377561 8: 0.63832042E-01 9: 0.17303867 Event 300 Current ALPHA: 1: 44.398514 2: 32.037732 3: 32.545268 4: 22.077280 5: 26.784561 6: 35.113381 7: 37.541745 8: 17.812057 9: 60.689461 Estimated component weights: 1: 0.14368451 2: 0.10368198 3: 0.10532449 4: 0.71447509E-01 5: 0.86681428E-01 6: 0.11363554 7: 0.12149432 8: 0.57644196E-01 9: 0.19640602 Event 400 Current ALPHA: 1: 57.016255 2: 38.622573 3: 42.438690 4: 28.195825 5: 37.293232 6: 48.388407 7: 48.654053 8: 23.496054 9: 84.894911 Estimated component weights: 1: 0.13940405 2: 0.94431719E-01 3: 0.10376208 4: 0.68938447E-01 5: 0.91181497E-01 6: 0.11830906 7: 0.11895856 8: 0.57447565E-01 9: 0.20756702 Event 500 Current ALPHA: 1: 72.837412 2: 45.487617 3: 54.768992 4: 34.839886 5: 49.712488 6: 53.170395 7: 65.953828 8: 28.917414 9: 103.31197 Estimated component weights: 1: 0.14309904 2: 0.89366634E-01 3: 0.10760116 4: 0.68447713E-01 5: 0.97666970E-01 6: 0.10446050 7: 0.12957530 8: 0.56812208E-01 9: 0.20297047 Event 600 Current ALPHA: 1: 86.009584 2: 52.688814 3: 65.634880 4: 40.038004 5: 64.648128 6: 63.381137 7: 80.521346 8: 34.172644 9: 121.90546 Estimated component weights: 1: 0.14123084 2: 0.86516936E-01 3: 0.10777484 4: 0.65743849E-01 5: 0.10615456 6: 0.10407412 7: 0.13221896 8: 0.56112716E-01 9: 0.20017317 Event 700 Current ALPHA: 1: 102.80462 2: 59.892095 3: 75.490221 4: 44.174589 5: 73.018452 6: 70.454704 7: 97.703423 8: 43.518927 9: 141.94296 Estimated component weights: 1: 0.14499947 2: 0.84474041E-01 3: 0.10647422 4: 0.62305486E-01 5: 0.10298794 6: 0.99371938E-01 7: 0.13780455 8: 0.61380715E-01 9: 0.20020164 Event 800 Current ALPHA: 1: 122.62435 2: 64.596389 3: 89.782553 4: 50.353083 5: 77.157742 6: 77.787730 7: 113.56167 8: 48.248383 9: 164.88810 Estimated component weights: 1: 0.15157521 2: 0.79847205E-01 3: 0.11097967 4: 0.62241141E-01 5: 0.95374218E-01 6: 0.96152942E-01 7: 0.14037290 8: 0.59639533E-01 9: 0.20381718 Event 900 Current ALPHA: 1: 134.94014 2: 71.577526 3: 100.35815 4: 56.244201 5: 87.015372 6: 87.443555 7: 128.88497 8: 53.426046 9: 189.11003 Estimated component weights: 1: 0.14844900 2: 0.78743153E-01 3: 0.11040501 4: 0.61874809E-01 5: 0.95726482E-01 6: 0.96197530E-01 7: 0.14178765 8: 0.58774528E-01 9: 0.20804184 Event 1000 Current ALPHA: 1: 149.35463 2: 76.976087 3: 113.79300 4: 62.987360 5: 92.626764 6: 95.172967 7: 142.25487 8: 57.675815 9: 218.15850 Estimated component weights: 1: 0.14802243 2: 0.76289482E-01 3: 0.11277800 4: 0.62425530E-01 5: 0.91800559E-01 6: 0.94324051E-01 7: 0.14098599 8: 0.57161363E-01 9: 0.21621259 bdmlib_test(): Normal end of execution. 9 September 2023 7:17:34.588 PM