Main function: simulated_annealing()
This is the main function for the phytoclass package.
It takes in the inputs (listed below) and returns the updated pigment to Chl a ratios, the Chl a biomass of each phytoplankton group, error associated with each group, and a graph displaying the Chl a concentration for each group.
It is important that samples are clustered appropriately before using the function (see the Cluster function).
Arguments:
S = Sample matrix – a matrix of pigment samples. Ensure that Chl a is the final column
F = Pigment to Chl a matrix. If left blank default values will be used. Ensure that pigment columns are in the same order as S and column naming conventions match.
user_defined_min_max = If blank default values are used. To create different min_max values, follow the same structure as the phytoclass::min_max file. See the example below.
do_matrix_checks = this should only be set to true when using the default values. This will remove pigment columns that have column sums of 0. Set to FALSE if using customised names for pigments and phytoplankton groups.
niter = number of iterations. Default value is 500.
step = step ratio used. Default value is 0.009.
weight.upper.bound = the upper limit of the weights applied. Default value is 30.
When using the default values, the only argument required is the sample matrix. However, make sure that the pigment names match those in the built-in pigment to Chl matrix Fm.
For the examples that follow the argument niter equals one for processing speed, but should be set much higher to obtain convergence.
Cluster function
Prior to analysis using simulated annealing, pigment samples require clustering.
The Cluster function divides all pigment concentrations by the total Chl a concentration. Following this the data undergoes BoxCox transformation, and the data is hierarchically clustered using the Ward method based on the Manhattan distances between pigment samples. The DynamicTreeCut method of (Langfelder, Zhang, and Horvath 2008) is then used to prune the dendogram into reasonable clusters of specified size(s).
The function returns a list of the clusters and the cluster dendrogram.
An example, using the built-in sample data set Sm:
# list of clusters
Cluster.result$cluster.list
#> $`1`
#> Per X19but Fuco Neox Pra Viol X19hex
#> 1 0.0000 0.03024000 0.0622500 0.005570000 0.01407000 0.007590000 0.08224000
#> 2 0.0000 0.01084000 0.0286400 0.001110000 0.00351000 0.001440000 0.01497000
#> 3 0.0000 0.01560000 0.2172000 0.006400000 0.00920000 0.000000000 0.01740000
#> 4 0.0000 0.01770000 0.2347000 0.007000000 0.01150000 0.000000000 0.01890000
#> 5 0.0000 0.02520000 0.2952000 0.009900000 0.01300000 0.000000000 0.02110000
#> 6 0.0102 0.02220000 0.2275000 0.007600000 0.01070000 0.000000000 0.01900000
#> 7 0.0000 0.01510000 0.1417000 0.004000000 0.00800000 0.002600000 0.02550000
#> 8 0.0000 0.04520000 0.4484000 0.003700000 0.00650000 0.000000000 0.13170000
#> 9 0.0000 0.05180000 0.6228000 0.004900000 0.00750000 0.004000000 0.13210000
#> 10 0.0000 0.02200000 0.1090000 0.004000000 0.00800000 0.003000000 0.04200000
#> 11 0.0000 0.04906269 0.5266418 0.017334346 0.03446357 0.000000000 0.12044287
#> 12 0.0000 0.02052662 0.2714934 0.017692288 0.01994758 0.000000000 0.03011548
#> 13 0.0000 0.02216177 0.1654257 0.012242414 0.01956286 0.000000000 0.02721049
#> 14 0.0000 0.01438545 0.2555594 0.020671983 0.04736366 0.011507712 0.02167651
#> 15 0.0000 0.02619187 0.6025477 0.038948613 0.04929322 0.034893958 0.10082727
#> 16 0.0000 0.02662521 0.2205439 0.024444776 0.06623835 0.027029275 0.04015983
#> 17 0.0000 0.06542364 0.6603459 0.044215705 0.11285620 0.064560275 0.13896777
#> 18 0.0000 0.05315834 0.2584586 0.016123905 0.04701666 0.011387548 0.08065818
#> 19 0.0000 0.06040911 0.4494852 0.027862628 0.07757631 0.025258248 0.10999693
#> 20 0.0000 0.07577285 0.3920489 0.017239170 0.05136883 0.020899111 0.12859331
#> 21 0.0000 0.06067605 0.4112795 0.021768249 0.07620416 0.022592896 0.08973832
#> 22 0.0000 0.04388275 0.6257830 0.012373717 0.03347165 0.013495185 0.12822965
#> 23 0.0000 0.03971444 0.4332551 0.014402749 0.02657611 0.010560487 0.09220387
#> 24 0.0000 0.03092146 0.1639481 0.013383240 0.01262176 0.008793182 0.05881285
#> 25 0.0000 0.02725515 0.1966991 0.009916299 0.01832477 0.009112193 0.05488214
#> 26 0.0000 0.01678833 0.1755048 0.008821737 0.02676049 0.006637964 0.02189435
#> 27 0.0000 0.02444961 0.1273596 0.007049435 0.01323978 0.004447137 0.02218985
#> 28 0.0000 0.01567523 0.1052905 0.006144340 0.01372236 0.006294480 0.02208074
#> 29 0.0000 0.03108003 0.0874373 0.011728700 0.02004430 0.005772257 0.04855466
#> Allo Zea Lut ChlcMGDG18 ChlcMGDG14 Chl_b Tchla
#> 1 0.001880000 0.002010000 0 0 0 0.08661000 0.4585100
#> 2 0.001440000 0.001910000 0 0 0 0.01473000 0.1457100
#> 3 0.005800000 0.003600000 0 0 0 0.05700000 0.6127000
#> 4 0.005400000 0.003800000 0 0 0 0.06190000 0.6207000
#> 5 0.001400000 0.007600000 0 0 0 0.05780000 0.5302000
#> 6 0.000000000 0.003000000 0 0 0 0.04530000 0.4057000
#> 7 0.000800000 0.001000000 0 0 0 0.03120000 0.2117000
#> 8 0.004400000 0.002300000 0 0 0 0.00930000 0.6885000
#> 9 0.059700000 0.003500000 0 0 0 0.01210000 1.0849000
#> 10 0.003000000 0.003000000 0 0 0 0.03200000 0.2060000
#> 11 0.029990012 0.039966194 0 0 0 0.15325050 1.0396607
#> 12 0.000000000 0.030828916 0 0 0 0.09243326 0.4812043
#> 13 0.000000000 0.017661542 0 0 0 0.08720595 0.3409079
#> 14 0.009032764 0.018938989 0 0 0 0.14213442 0.4892902
#> 15 0.031610292 0.036313122 0 0 0 0.16461333 1.3811642
#> 16 0.007170750 0.013507988 0 0 0 0.18532190 0.5711149
#> 17 0.016001053 0.028670688 0 0 0 0.41729467 1.7413395
#> 18 0.010093841 0.016835362 0 0 0 0.13198111 0.5776707
#> 19 0.019249598 0.020581908 0 0 0 0.26015615 1.0739310
#> 20 0.024550418 0.020066576 0 0 0 0.22391720 0.9909010
#> 21 0.016719215 0.019212773 0 0 0 0.23686078 0.9520484
#> 22 0.018189726 0.017601736 0 0 0 0.15579381 1.2494685
#> 23 0.014131906 0.012984478 0 0 0 0.12090732 0.9190665
#> 24 0.011139374 0.012782671 0 0 0 0.05213182 0.3486558
#> 25 0.013767971 0.012392348 0 0 0 0.08734698 0.4556423
#> 26 0.008302439 0.010946550 0 0 0 0.07005955 0.3497736
#> 27 0.005887899 0.010324368 0 0 0 0.03627212 0.2367185
#> 28 0.004908419 0.009914413 0 0 0 0.03735579 0.2117071
#> 29 0.009256850 0.009871839 0 0 0 0.09414383 0.3014690
# plot of clusters
plot(Cluster.result$cluster.plot)Example without clustering
The example here uses the built-in sample matrix Sm.
Results$`condition number`
#> [1] 490.8151
Results$RMSE
#> [1] 0.0266215
Results$MAE
#> Per X19but Fuco Neox Pra Viol
#> 4.025195e-05 3.796674e-04 4.577666e-03 3.465166e-03 2.550379e-04 3.888018e-03
#> X19hex Allo Zea Chl_b Tchla
#> 4.820585e-04 1.304592e-04 2.194573e-04 2.768145e-03 5.920934e-02
Results$Error
#> Per X19but Fuco Neox Pra
#> [1,] -2.446502e-04 -1.796359e-03 -0.0152477051 -0.0029656500 7.994156e-04
#> [2,] -3.381334e-04 5.794041e-05 -0.0210740043 -0.0014928139 3.660001e-04
#> [3,] -2.401804e-04 4.115579e-05 -0.0149691266 0.0014102150 1.155096e-04
#> [4,] -2.069773e-04 3.546632e-05 -0.0128997597 0.0014613457 8.625066e-05
#> [5,] -2.736381e-05 4.688891e-06 -0.0017054368 0.0047142303 -2.625719e-04
#> [6,] -7.452407e-06 1.276998e-06 -0.0004644677 0.0044995046 -2.649827e-04
#> [7,] 0.000000e+00 -2.163368e-05 0.0078685721 0.0023230746 1.659747e-04
#> [8,] 0.000000e+00 -3.198299e-03 0.0029961332 0.0018850780 -1.555917e-04
#> [9,] 0.000000e+00 -2.880282e-06 0.0010476122 0.0015983826 -1.163054e-04
#> [10,] 0.000000e+00 -1.441037e-05 0.0052413210 0.0021209882 -1.984937e-04
#> [11,] 0.000000e+00 -6.918837e-06 0.0025165095 0.0017931617 -1.300810e-04
#> [12,] 0.000000e+00 -1.158810e-05 0.0042148058 0.0094906748 -6.318108e-04
#> [13,] 0.000000e+00 -9.387938e-06 0.0034145677 0.0062413943 -4.114213e-04
#> [14,] 0.000000e+00 -1.814101e-05 0.0065982231 0.0060149219 -4.556140e-04
#> [15,] -5.653298e-05 -3.837021e-03 -0.0035233913 0.0083017906 -4.886041e-04
#> [16,] 0.000000e+00 -7.428384e-06 0.0027018413 0.0039927502 -1.314450e-05
#> [17,] -1.762040e-05 3.019320e-06 -0.0010981829 0.0006283880 -3.550086e-05
#> [18,] 0.000000e+00 -7.449477e-06 0.0027095134 0.0021322291 -1.641408e-04
#> [19,] 0.000000e+00 -4.066618e-06 0.0014791047 0.0009540358 -7.612857e-05
#> [20,] 0.000000e+00 -3.315335e-06 0.0012058490 -0.0019175524 1.084244e-04
#> [21,] 0.000000e+00 -4.845374e-06 0.0017623528 -0.0009965889 4.324982e-05
#> [22,] -6.768279e-06 -1.874696e-03 -0.0004218298 -0.0008999813 6.479423e-05
#> [23,] -2.162752e-05 3.705955e-06 -0.0013479248 0.0017358908 -8.888653e-05
#> [24,] 0.000000e+00 -1.068377e-05 0.0038858870 0.0110104280 -7.495480e-04
#> [25,] 0.000000e+00 -5.300569e-06 0.0019279156 0.0018672132 -1.411382e-04
#> [26,] 0.000000e+00 -6.148142e-06 0.0022361937 0.0021137566 -1.587922e-04
#> [27,] 0.000000e+00 -8.626152e-06 0.0031374917 0.0063984913 -4.440787e-04
#> [28,] 0.000000e+00 -8.299461e-06 0.0030186684 0.0047253060 -3.447244e-04
#> [29,] 0.000000e+00 -5.603026e-06 0.0020379248 0.0048039661 -3.149219e-04
#> Viol X19hex Allo Zea Chl_b
#> [1,] -0.0024247490 1.281137e-03 -4.610417e-04 -2.222705e-03 0.0011790511
#> [2,] -0.0010095750 -7.448523e-04 -6.372104e-04 -1.228118e-04 -0.0019230658
#> [3,] -0.0064433816 -5.290778e-04 -4.526184e-04 -8.723472e-05 0.0016377465
#> [4,] -0.0066909523 -4.559369e-04 -3.900474e-04 -7.517519e-05 0.0021483036
#> [5,] -0.0063948977 -6.027799e-05 -5.156695e-05 -9.938676e-06 0.0019510032
#> [6,] -0.0063811539 -1.641643e-05 -1.404402e-05 -2.706752e-06 0.0022722005
#> [7,] -0.0017164175 2.781116e-04 2.379204e-04 -2.073154e-03 0.0015108612
#> [8,] -0.0007672260 3.264340e-03 9.059347e-05 1.746039e-05 0.0000206543
#> [9,] 0.0011627367 3.702744e-05 3.167644e-05 6.105109e-06 -0.0016318626
#> [10,] -0.0011406832 1.852524e-04 1.584808e-04 3.054455e-05 0.0005524276
#> [11,] -0.0077855985 8.894503e-05 7.609119e-05 1.466532e-05 0.0054019854
#> [12,] -0.0101910523 1.489706e-04 0.000000e+00 2.456238e-05 0.0031491494
#> [13,] -0.0130683875 1.206865e-04 0.000000e+00 1.989888e-05 0.0070865663
#> [14,] -0.0039435857 2.332116e-04 1.995091e-04 3.845209e-05 0.0007192470
#> [15,] 0.0057337675 3.684050e-03 -1.065361e-04 -2.053306e-05 -0.0097518295
#> [16,] 0.0049873653 9.549551e-05 8.169502e-05 -1.445227e-03 -0.0055385806
#> [17,] 0.0051753164 -3.881484e-05 -3.320553e-05 -6.399817e-06 -0.0045329281
#> [18,] -0.0024927034 9.576667e-05 8.192701e-05 1.579008e-05 0.0011612606
#> [19,] -0.0014940030 5.227837e-05 4.472339e-05 8.619694e-06 0.0008564534
#> [20,] -0.0016293908 4.262025e-05 3.646101e-05 7.027257e-06 0.0025536745
#> [21,] -0.0015978223 6.228966e-05 5.328790e-05 1.027036e-05 0.0020975397
#> [22,] -0.0014098820 1.842354e-03 -1.275478e-05 -2.458273e-06 0.0015292662
#> [23,] -0.0016169078 -4.764187e-05 -4.075693e-05 -7.855224e-06 0.0000282889
#> [24,] 0.0033890633 1.373451e-04 1.174968e-04 2.264556e-05 -0.0082412069
#> [25,] -0.0005809676 6.814141e-05 5.829399e-05 1.123520e-05 -0.0002928274
#> [26,] -0.0015209495 7.903738e-05 6.761533e-05 1.303174e-05 0.0003429731
#> [27,] 0.0005354803 1.108934e-04 9.486770e-05 1.828418e-05 -0.0035356308
#> [28,] 0.0044168148 1.066936e-04 9.127486e-05 1.759172e-05 -0.0055924275
#> [29,] -0.0070516941 7.202964e-05 6.162032e-05 1.187630e-05 0.0030372062
#> Tchla
#> [1,] 0.197219829
#> [2,] 0.272579481
#> [3,] 0.193616586
#> [4,] 0.166850578
#> [5,] 0.022058792
#> [6,] 0.006007609
#> [7,] -0.101775214
#> [8,] -0.038753168
#> [9,] -0.013550229
#> [10,] -0.067793313
#> [11,] -0.032549526
#> [12,] -0.054515960
#> [13,] -0.044165365
#> [14,] -0.085344019
#> [15,] 0.045572932
#> [16,] -0.034946681
#> [17,] 0.014204331
#> [18,] -0.035045915
#> [19,] -0.019131324
#> [20,] -0.015596926
#> [21,] -0.022794965
#> [22,] 0.005456113
#> [23,] 0.017434591
#> [24,] -0.050261594
#> [25,] -0.024936420
#> [26,] -0.028923811
#> [27,] -0.040581555
#> [28,] -0.039044648
#> [29,] -0.026359323
Results$`F matrix`
#> Per X19but Fuco Neox Pra Viol X19hex Allo
#> Prasinophytes 0.0000 0.0000 0.0000 0.1214 0.7119 0.1048 0.000 0.0000
#> Chlorophytes 0.0000 0.0000 0.0000 0.0485 0.0000 0.0655 0.000 0.0000
#> Cryptophytes 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.4753
#> Diatoms-2 0.0000 0.0000 0.7314 0.0000 0.0000 0.0000 0.000 0.0000
#> Dinoflagellates-1 0.8957 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.0000
#> Haptophytes 0.0000 0.2946 0.2374 0.0000 0.0000 0.0000 0.758 0.0000
#> Pelagophytes 0.0000 1.2557 0.9071 0.0000 0.0000 0.0000 0.000 0.0000
#> Syn 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.0000
#> Zea Chl_b Tchla
#> Prasinophytes 0.1495 0.9447 1
#> Chlorophytes 0.0075 0.5983 1
#> Cryptophytes 0.0000 0.0000 1
#> Diatoms-2 0.0000 0.0000 1
#> Dinoflagellates-1 0.0000 0.0000 1
#> Haptophytes 0.0000 0.0000 1
#> Pelagophytes 0.0000 0.0000 1
#> Syn 2.4661 0.0000 1
Results$`Class abundances`
#> Prasinophytes Chlorophytes Cryptophytes Diatoms-2 Dinoflagellates-1
#> 1 0.027948563 0.1675149708 6.918687e-03 0.09743578 3.029893e-04
#> 2 0.008296230 0.0305610550 5.738042e-03 0.05902072 1.458900e-04
#> 3 0.018206062 0.1033925158 1.868051e-02 0.42791216 3.612360e-04
#> 4 0.021782704 0.1013621487 1.652790e-02 0.43456461 3.081319e-04
#> 5 0.019391211 0.0667483376 3.176885e-03 0.39586200 3.059471e-05
#> 6 0.015482030 0.0492355402 2.244591e-05 0.29200398 1.152216e-02
#> 7 0.009184671 0.0275933267 1.206051e-03 0.14265321 0.000000e+00
#> 8 0.008762466 0.0005751572 8.370564e-03 0.51427120 0.000000e+00
#> 9 0.010596713 0.0082818648 1.224397e-01 0.77251874 0.000000e+00
#> 10 0.009944229 0.0307765763 5.400044e-03 0.10711162 0.000000e+00
#> 11 0.045889198 0.1514349929 5.905721e-02 0.62027395 0.000000e+00
#> 12 0.026032220 0.0935987880 0.000000e+00 0.31014629 0.000000e+00
#> 13 0.025584969 0.0858344125 0.000000e+00 0.18339593 0.000000e+00
#> 14 0.056953277 0.1103892512 1.574205e-02 0.27542758 0.000000e+00
#> 15 0.077218741 0.2213451801 7.299963e-02 0.86698127 1.694604e-04
#> 16 0.086788327 0.1620215948 1.387960e-02 0.25093616 0.000000e+00
#> 17 0.163068035 0.4848302327 3.482962e-02 0.86009562 6.649965e-05
#> 18 0.061808916 0.1058131052 1.959853e-02 0.27290927 0.000000e+00
#> 19 0.105215917 0.2499055799 3.883039e-02 0.52581191 0.000000e+00
#> 20 0.069726059 0.2449842800 4.997305e-02 0.43883654 0.000000e+00
#> 21 0.102258188 0.2107449798 3.343559e-02 0.47224765 0.000000e+00
#> 22 0.047282700 0.1824417778 3.872026e-02 0.81197460 1.754289e-05
#> 23 0.038780178 0.1474375403 3.086279e-02 0.57037467 4.199850e-05
#> 24 0.016757691 0.0614580007 2.109172e-02 0.16976219 0.000000e+00
#> 25 0.024718492 0.1006285481 2.752467e-02 0.22628238 0.000000e+00
#> 26 0.035692609 0.0539947186 1.642425e-02 0.20829723 0.000000e+00
#> 27 0.017443079 0.0310642149 1.134164e-02 0.13556430 0.000000e+00
#> 28 0.018044131 0.0330759076 9.486123e-03 0.11640912 0.000000e+00
#> 29 0.026969533 0.1037018832 1.839911e-02 0.08060536 0.000000e+00
#> Haptophytes Pelagophytes Syn
#> 1 0.15838901 0.0000000000 0.0000000000
#> 2 0.03440674 0.0067835342 0.0007577908
#> 3 0.03366761 0.0097691392 0.0007107633
#> 4 0.03467359 0.0109758268 0.0005050903
#> 5 0.02907586 0.0140794142 0.0018356963
#> 6 0.02536438 0.0119267415 0.0001427245
#> 7 0.02761621 0.0034465286 0.0000000000
#> 8 0.15619515 0.0000000000 0.0003254667
#> 9 0.16996725 0.0003831301 0.0007125658
#> 10 0.04841964 0.0039839370 0.0003639526
#> 11 0.14926293 0.0017497417 0.0119927125
#> 12 0.03564665 0.0063818571 0.0093985150
#> 13 0.03284164 0.0084959218 0.0047550255
#> 14 0.02397100 0.0040993762 0.0027076371
#> 15 0.13175246 0.0000000000 0.0106974511
#> 16 0.04926846 0.0082207993 0.0000000000
#> 17 0.18856320 0.0092910086 0.0005952704
#> 18 0.09902752 0.0162270733 0.0022862649
#> 19 0.13968168 0.0135894723 0.0008960068
#> 20 0.16450171 0.0199615038 0.0029178445
#> 21 0.11306864 0.0196904277 0.0006029252
#> 22 0.16524066 0.0000000000 0.0037909285
#> 23 0.12576392 0.0031602580 0.0026451288
#> 24 0.07023888 0.0058559128 0.0034913921
#> 25 0.06898191 0.0045217743 0.0029844838
#> 26 0.02724494 0.0062538769 0.0018659637
#> 27 0.02694860 0.0116485794 0.0027081034
#> 28 0.02691915 0.0052470630 0.0025255900
#> 29 0.06070972 0.0092391362 0.0018442123Example with clustering
Clust1 <- Cluster(Sm, minSamplesPerCluster = 14)$cluster.list[[1]]
#> ..done.
set.seed("7683")
Results <- simulated_annealing(Clust1, niter = 1)Results$`condition number`
#> [1] 490.8151
Results$RMSE
#> [1] 0.0266215
Results$MAE
#> Per X19but Fuco Neox Pra Viol
#> 4.025195e-05 3.796674e-04 4.577666e-03 3.465166e-03 2.550379e-04 3.888018e-03
#> X19hex Allo Zea Chl_b Tchla
#> 4.820585e-04 1.304592e-04 2.194573e-04 2.768145e-03 5.920934e-02
Results$Error
#> Per X19but Fuco Neox Pra
#> 1 -2.446502e-04 -1.796359e-03 -0.0152477051 -0.0029656500 7.994156e-04
#> 2 -3.381334e-04 5.794041e-05 -0.0210740043 -0.0014928139 3.660001e-04
#> 3 -2.401804e-04 4.115579e-05 -0.0149691266 0.0014102150 1.155096e-04
#> 4 -2.069773e-04 3.546632e-05 -0.0128997597 0.0014613457 8.625066e-05
#> 5 -2.736381e-05 4.688891e-06 -0.0017054368 0.0047142303 -2.625719e-04
#> 6 -7.452407e-06 1.276998e-06 -0.0004644677 0.0044995046 -2.649827e-04
#> 7 0.000000e+00 -2.163368e-05 0.0078685721 0.0023230746 1.659747e-04
#> 8 0.000000e+00 -3.198299e-03 0.0029961332 0.0018850780 -1.555917e-04
#> 9 0.000000e+00 -2.880282e-06 0.0010476122 0.0015983826 -1.163054e-04
#> 10 0.000000e+00 -1.441037e-05 0.0052413210 0.0021209882 -1.984937e-04
#> 11 0.000000e+00 -6.918837e-06 0.0025165095 0.0017931617 -1.300810e-04
#> 12 0.000000e+00 -1.158810e-05 0.0042148058 0.0094906748 -6.318108e-04
#> 13 0.000000e+00 -9.387938e-06 0.0034145677 0.0062413943 -4.114213e-04
#> 14 0.000000e+00 -1.814101e-05 0.0065982231 0.0060149219 -4.556140e-04
#> 15 -5.653298e-05 -3.837021e-03 -0.0035233913 0.0083017906 -4.886041e-04
#> 16 0.000000e+00 -7.428384e-06 0.0027018413 0.0039927502 -1.314450e-05
#> 17 -1.762040e-05 3.019320e-06 -0.0010981829 0.0006283880 -3.550086e-05
#> 18 0.000000e+00 -7.449477e-06 0.0027095134 0.0021322291 -1.641408e-04
#> 19 0.000000e+00 -4.066618e-06 0.0014791047 0.0009540358 -7.612857e-05
#> 20 0.000000e+00 -3.315335e-06 0.0012058490 -0.0019175524 1.084244e-04
#> 21 0.000000e+00 -4.845374e-06 0.0017623528 -0.0009965889 4.324982e-05
#> 22 -6.768279e-06 -1.874696e-03 -0.0004218298 -0.0008999813 6.479423e-05
#> 23 -2.162752e-05 3.705955e-06 -0.0013479248 0.0017358908 -8.888653e-05
#> 24 0.000000e+00 -1.068377e-05 0.0038858870 0.0110104280 -7.495480e-04
#> 25 0.000000e+00 -5.300569e-06 0.0019279156 0.0018672132 -1.411382e-04
#> 26 0.000000e+00 -6.148142e-06 0.0022361937 0.0021137566 -1.587922e-04
#> 27 0.000000e+00 -8.626152e-06 0.0031374917 0.0063984913 -4.440787e-04
#> 28 0.000000e+00 -8.299461e-06 0.0030186684 0.0047253060 -3.447244e-04
#> 29 0.000000e+00 -5.603026e-06 0.0020379248 0.0048039661 -3.149219e-04
#> Viol X19hex Allo Zea Chl_b
#> 1 -0.0024247490 1.281137e-03 -4.610417e-04 -2.222705e-03 0.0011790511
#> 2 -0.0010095750 -7.448523e-04 -6.372104e-04 -1.228118e-04 -0.0019230658
#> 3 -0.0064433816 -5.290778e-04 -4.526184e-04 -8.723472e-05 0.0016377465
#> 4 -0.0066909523 -4.559369e-04 -3.900474e-04 -7.517519e-05 0.0021483036
#> 5 -0.0063948977 -6.027799e-05 -5.156695e-05 -9.938676e-06 0.0019510032
#> 6 -0.0063811539 -1.641643e-05 -1.404402e-05 -2.706752e-06 0.0022722005
#> 7 -0.0017164175 2.781116e-04 2.379204e-04 -2.073154e-03 0.0015108612
#> 8 -0.0007672260 3.264340e-03 9.059347e-05 1.746039e-05 0.0000206543
#> 9 0.0011627367 3.702744e-05 3.167644e-05 6.105109e-06 -0.0016318626
#> 10 -0.0011406832 1.852524e-04 1.584808e-04 3.054455e-05 0.0005524276
#> 11 -0.0077855985 8.894503e-05 7.609119e-05 1.466532e-05 0.0054019854
#> 12 -0.0101910523 1.489706e-04 0.000000e+00 2.456238e-05 0.0031491494
#> 13 -0.0130683875 1.206865e-04 0.000000e+00 1.989888e-05 0.0070865663
#> 14 -0.0039435857 2.332116e-04 1.995091e-04 3.845209e-05 0.0007192470
#> 15 0.0057337675 3.684050e-03 -1.065361e-04 -2.053306e-05 -0.0097518295
#> 16 0.0049873653 9.549551e-05 8.169502e-05 -1.445227e-03 -0.0055385806
#> 17 0.0051753164 -3.881484e-05 -3.320553e-05 -6.399817e-06 -0.0045329281
#> 18 -0.0024927034 9.576667e-05 8.192701e-05 1.579008e-05 0.0011612606
#> 19 -0.0014940030 5.227837e-05 4.472339e-05 8.619694e-06 0.0008564534
#> 20 -0.0016293908 4.262025e-05 3.646101e-05 7.027257e-06 0.0025536745
#> 21 -0.0015978223 6.228966e-05 5.328790e-05 1.027036e-05 0.0020975397
#> 22 -0.0014098820 1.842354e-03 -1.275478e-05 -2.458273e-06 0.0015292662
#> 23 -0.0016169078 -4.764187e-05 -4.075693e-05 -7.855224e-06 0.0000282889
#> 24 0.0033890633 1.373451e-04 1.174968e-04 2.264556e-05 -0.0082412069
#> 25 -0.0005809676 6.814141e-05 5.829399e-05 1.123520e-05 -0.0002928274
#> 26 -0.0015209495 7.903738e-05 6.761533e-05 1.303174e-05 0.0003429731
#> 27 0.0005354803 1.108934e-04 9.486770e-05 1.828418e-05 -0.0035356308
#> 28 0.0044168148 1.066936e-04 9.127486e-05 1.759172e-05 -0.0055924275
#> 29 -0.0070516941 7.202964e-05 6.162032e-05 1.187630e-05 0.0030372062
#> Tchla
#> 1 0.197219829
#> 2 0.272579481
#> 3 0.193616586
#> 4 0.166850578
#> 5 0.022058792
#> 6 0.006007609
#> 7 -0.101775214
#> 8 -0.038753168
#> 9 -0.013550229
#> 10 -0.067793313
#> 11 -0.032549526
#> 12 -0.054515960
#> 13 -0.044165365
#> 14 -0.085344019
#> 15 0.045572932
#> 16 -0.034946681
#> 17 0.014204331
#> 18 -0.035045915
#> 19 -0.019131324
#> 20 -0.015596926
#> 21 -0.022794965
#> 22 0.005456113
#> 23 0.017434591
#> 24 -0.050261594
#> 25 -0.024936420
#> 26 -0.028923811
#> 27 -0.040581555
#> 28 -0.039044648
#> 29 -0.026359323
Results$`F matrix`
#> Per X19but Fuco Neox Pra Viol X19hex Allo
#> Prasinophytes 0.0000 0.0000 0.0000 0.1214 0.7119 0.1048 0.000 0.0000
#> Chlorophytes 0.0000 0.0000 0.0000 0.0485 0.0000 0.0655 0.000 0.0000
#> Cryptophytes 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.4753
#> Diatoms-2 0.0000 0.0000 0.7314 0.0000 0.0000 0.0000 0.000 0.0000
#> Dinoflagellates-1 0.8957 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.0000
#> Haptophytes 0.0000 0.2946 0.2374 0.0000 0.0000 0.0000 0.758 0.0000
#> Pelagophytes 0.0000 1.2557 0.9071 0.0000 0.0000 0.0000 0.000 0.0000
#> Syn 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.000 0.0000
#> Zea Chl_b Tchla
#> Prasinophytes 0.1495 0.9447 1
#> Chlorophytes 0.0075 0.5983 1
#> Cryptophytes 0.0000 0.0000 1
#> Diatoms-2 0.0000 0.0000 1
#> Dinoflagellates-1 0.0000 0.0000 1
#> Haptophytes 0.0000 0.0000 1
#> Pelagophytes 0.0000 0.0000 1
#> Syn 2.4661 0.0000 1
Results$`Class abundances`
#> Prasinophytes Chlorophytes Cryptophytes Diatoms-2 Dinoflagellates-1
#> 1 0.027948563 0.1675149708 6.918687e-03 0.09743578 3.029893e-04
#> 2 0.008296230 0.0305610550 5.738042e-03 0.05902072 1.458900e-04
#> 3 0.018206062 0.1033925158 1.868051e-02 0.42791216 3.612360e-04
#> 4 0.021782704 0.1013621487 1.652790e-02 0.43456461 3.081319e-04
#> 5 0.019391211 0.0667483376 3.176885e-03 0.39586200 3.059471e-05
#> 6 0.015482030 0.0492355402 2.244591e-05 0.29200398 1.152216e-02
#> 7 0.009184671 0.0275933267 1.206051e-03 0.14265321 0.000000e+00
#> 8 0.008762466 0.0005751572 8.370564e-03 0.51427120 0.000000e+00
#> 9 0.010596713 0.0082818648 1.224397e-01 0.77251874 0.000000e+00
#> 10 0.009944229 0.0307765763 5.400044e-03 0.10711162 0.000000e+00
#> 11 0.045889198 0.1514349929 5.905721e-02 0.62027395 0.000000e+00
#> 12 0.026032220 0.0935987880 0.000000e+00 0.31014629 0.000000e+00
#> 13 0.025584969 0.0858344125 0.000000e+00 0.18339593 0.000000e+00
#> 14 0.056953277 0.1103892512 1.574205e-02 0.27542758 0.000000e+00
#> 15 0.077218741 0.2213451801 7.299963e-02 0.86698127 1.694604e-04
#> 16 0.086788327 0.1620215948 1.387960e-02 0.25093616 0.000000e+00
#> 17 0.163068035 0.4848302327 3.482962e-02 0.86009562 6.649965e-05
#> 18 0.061808916 0.1058131052 1.959853e-02 0.27290927 0.000000e+00
#> 19 0.105215917 0.2499055799 3.883039e-02 0.52581191 0.000000e+00
#> 20 0.069726059 0.2449842800 4.997305e-02 0.43883654 0.000000e+00
#> 21 0.102258188 0.2107449798 3.343559e-02 0.47224765 0.000000e+00
#> 22 0.047282700 0.1824417778 3.872026e-02 0.81197460 1.754289e-05
#> 23 0.038780178 0.1474375403 3.086279e-02 0.57037467 4.199850e-05
#> 24 0.016757691 0.0614580007 2.109172e-02 0.16976219 0.000000e+00
#> 25 0.024718492 0.1006285481 2.752467e-02 0.22628238 0.000000e+00
#> 26 0.035692609 0.0539947186 1.642425e-02 0.20829723 0.000000e+00
#> 27 0.017443079 0.0310642149 1.134164e-02 0.13556430 0.000000e+00
#> 28 0.018044131 0.0330759076 9.486123e-03 0.11640912 0.000000e+00
#> 29 0.026969533 0.1037018832 1.839911e-02 0.08060536 0.000000e+00
#> Haptophytes Pelagophytes Syn
#> 1 0.15838901 0.0000000000 0.0000000000
#> 2 0.03440674 0.0067835342 0.0007577908
#> 3 0.03366761 0.0097691392 0.0007107633
#> 4 0.03467359 0.0109758268 0.0005050903
#> 5 0.02907586 0.0140794142 0.0018356963
#> 6 0.02536438 0.0119267415 0.0001427245
#> 7 0.02761621 0.0034465286 0.0000000000
#> 8 0.15619515 0.0000000000 0.0003254667
#> 9 0.16996725 0.0003831301 0.0007125658
#> 10 0.04841964 0.0039839370 0.0003639526
#> 11 0.14926293 0.0017497417 0.0119927125
#> 12 0.03564665 0.0063818571 0.0093985150
#> 13 0.03284164 0.0084959218 0.0047550255
#> 14 0.02397100 0.0040993762 0.0027076371
#> 15 0.13175246 0.0000000000 0.0106974511
#> 16 0.04926846 0.0082207993 0.0000000000
#> 17 0.18856320 0.0092910086 0.0005952704
#> 18 0.09902752 0.0162270733 0.0022862649
#> 19 0.13968168 0.0135894723 0.0008960068
#> 20 0.16450171 0.0199615038 0.0029178445
#> 21 0.11306864 0.0196904277 0.0006029252
#> 22 0.16524066 0.0000000000 0.0037909285
#> 23 0.12576392 0.0031602580 0.0026451288
#> 24 0.07023888 0.0058559128 0.0034913921
#> 25 0.06898191 0.0045217743 0.0029844838
#> 26 0.02724494 0.0062538769 0.0018659637
#> 27 0.02694860 0.0116485794 0.0027081034
#> 28 0.02691915 0.0052470630 0.0025255900
#> 29 0.06070972 0.0092391362 0.0018442123Example using non-default values
#Create Fm (F matrix). Alternatively, a .csv file can be uploaded.
#Create Fm (F matrix). Alternatively, a .csv file can be uploaded.
Fu <- data.frame(
Per = c(0, 0, 0, 0, 1, 0, 0, 0),
X19but = c(0, 0, 0, 0, 0, 1, 1, 0),
Fuco = c(0, 0, 0, 1, 0, 1, 1, 0),
Pra = c(1, 0, 0, 0, 0, 0, 0, 0),
X19hex = c(0, 0, 0, 0, 0, 1, 0, 0),
Allo = c(0, 0, 1, 0, 0, 0, 0, 0),
Zea = c(1, 1, 0, 0, 0, 0, 0, 1),
Chl_b = c(1, 1, 0, 0, 0, 0, 0, 0),
Tchla = c(1, 1, 1, 1, 1, 1, 1, 1)
)
rownames(Fu) <- c(
"Prasinophytes", "Chlorophytes", "Cryptophytes"
, "Diatoms-2", "Dinoflagellates-1",
"Haptophytes", "Pelagophytes", "Syn"
)
Min_max <- data.frame(
Class = c(
"Syn", "Chlorophytes", "Chlorophytes", "Prasinophytes", "Prasinophytes",
"Prasinophytes", "Cryptophytes", "Diatoms-2", "Diatoms-2", "Pelagophytes",
"Pelagophytes", "Pelagophytes", "Dinoflagellates-1", "Haptophytes",
"Haptophytes", "Haptophytes", "Haptophytes", "Diatoms-2", "Cryptophytes",
"Prasinophytes", "Chlorophytes", "Syn", "Dinoflagellates-1", "Pelagophytes"
),
Pig_Abbrev = c(
"Zea", "Zea", "Chl_b", "Pra", "Zea", "Chl_b", "Allo", "Chl_c3",
"Fuco", "Chl_c3", "X19but", "Fuco", "Per", "X19but", "X19hex",
"Fuco", "Tchla", "Tchla", "Tchla", "Tchla", "Tchla", "Tchla", "Tchla",
"Tchla"
),
min = as.numeric(c(
0.0800, 0.0063, 0.1666, 0.0642, 0.0151, 0.4993, 0.2118, 0.0189,
0.3315, 0.1471, 0.2457, 0.3092, 0.3421, 0.0819, 0.2107, 0.0090,
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000
)),
max = as.numeric(c(
1.2123, 0.0722, 0.9254, 0.4369, 0.1396, 0.9072, 0.5479, 0.1840,
0.9332, 0.2967, 1.0339, 1.2366, 0.8650, 0.2872, 1.3766, 0.4689,
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000
))
)set.seed("7683")
Results <- simulated_annealing(
S = Sm,
F = Fu,
user_defined_min_max = Min_max,
do_matrix_checks = TRUE,
niter = 1,
step = 0.01,
weight.upper.bound = 30
)set.seed("7683")
Results <- simulated_annealing(
S = Sm,
F = Fu,
user_defined_min_max = Min_max,
do_matrix_checks = TRUE,
niter = 1,
step = 0.01,
weight.upper.bound = 30
)Results$`condition number`
#> [1] 1042.712
Results$RMSE
#> [1] 0.02472112
Results$MAE
#> Per X19but Fuco Pra X19hex Allo
#> 5.063212e-05 5.343446e-05 3.888053e-03 1.397568e-03 1.523488e-04 1.355428e-04
#> Zea Chl_b Tchla
#> 7.020114e-04 3.272465e-03 5.520775e-02
Results$Error
#> Per X19but Fuco Pra X19hex
#> [1,] -1.823283e-04 -9.338014e-05 -7.129245e-03 1.057403e-03 -2.825419e-04
#> [2,] -3.868116e-04 -1.981070e-04 -1.512477e-02 1.873648e-03 -5.994157e-04
#> [3,] -2.809690e-04 -1.438993e-04 -1.098621e-02 1.387444e-03 -4.353986e-04
#> [4,] -2.503103e-04 -1.281974e-04 -9.787419e-03 1.244128e-03 -3.878889e-04
#> [5,] -4.327118e-05 -2.216150e-05 -1.691952e-03 2.095980e-04 -6.705441e-05
#> [6,] -1.460049e-05 -7.477695e-06 -5.708957e-04 9.494176e-05 -2.262539e-05
#> [7,] 0.000000e+00 6.748269e-05 5.152066e-03 -5.442885e-04 2.041835e-04
#> [8,] 0.000000e+00 2.524220e-05 1.927153e-03 1.695728e-03 7.637574e-05
#> [9,] -9.541554e-07 -4.886743e-07 -3.730856e-05 1.167201e-03 -1.478590e-06
#> [10,] 0.000000e+00 6.746479e-05 5.150699e-03 -6.380654e-04 2.041293e-04
#> [11,] 0.000000e+00 4.494477e-05 3.431375e-03 -4.250766e-04 1.359901e-04
#> [12,] 0.000000e+00 6.894816e-05 5.263949e-03 -6.520948e-04 2.086176e-04
#> [13,] 0.000000e+00 7.191380e-05 5.490365e-03 -6.801431e-04 2.175908e-04
#> [14,] 0.000000e+00 3.798773e-05 2.900230e-03 2.757486e-03 1.149401e-04
#> [15,] -1.467133e-04 -1.478730e-04 -5.736656e-03 7.106534e-04 -1.768873e-04
#> [16,] -1.445952e-05 -7.405499e-06 -5.653838e-04 5.754373e-03 -2.240694e-05
#> [17,] -7.650486e-05 -3.918225e-05 -2.991427e-03 4.175204e-04 -1.185544e-04
#> [18,] 0.000000e+00 2.966029e-05 2.264459e-03 3.719941e-03 8.974364e-05
#> [19,] -2.780644e-05 -1.424118e-05 -1.087264e-03 1.475015e-04 -4.308976e-05
#> [20,] 0.000000e+00 3.881392e-05 2.963307e-03 -3.670926e-04 1.174399e-04
#> [21,] -1.214088e-05 -6.217995e-06 -4.747221e-04 1.579300e-03 -1.881389e-05
#> [22,] -6.822903e-06 -3.494375e-06 -2.667833e-04 3.304895e-05 -1.057299e-05
#> [23,] -2.463853e-05 -1.261872e-05 -9.633945e-04 1.193447e-04 -3.818066e-05
#> [24,] 0.000000e+00 5.426024e-05 4.142578e-03 -5.131801e-04 1.641761e-04
#> [25,] 0.000000e+00 4.718996e-05 3.602788e-03 -4.463111e-04 1.427835e-04
#> [26,] 0.000000e+00 9.973554e-06 7.614457e-04 4.974753e-03 3.017715e-05
#> [27,] 0.000000e+00 4.506843e-05 3.440816e-03 2.962516e-03 1.363643e-04
#> [28,] 0.000000e+00 2.997056e-05 2.288147e-03 3.543958e-03 9.068242e-05
#> [29,] 0.000000e+00 8.593361e-05 6.560729e-03 -8.127390e-04 2.600108e-04
#> Allo Zea Chl_b Tchla
#> [1,] -2.747395e-04 -6.492457e-03 0.0007073838 0.1012305130
#> [2,] -5.828629e-04 -3.492230e-04 -0.0058887019 0.2147616655
#> [3,] -4.233751e-04 -1.215321e-03 -0.0037480536 0.1559967872
#> [4,] -3.771774e-04 -1.376115e-03 -0.0031775723 0.1389747970
#> [5,] -6.520271e-05 -3.906628e-05 -0.0006587473 0.0240245914
#> [6,] -2.200059e-05 -8.927764e-04 0.0002618900 0.0081063365
#> [7,] 1.985450e-04 -3.292934e-03 0.0038839522 -0.0731558841
#> [8,] 7.426664e-05 4.449694e-05 -0.0035463265 -0.0273642847
#> [9,] -1.437759e-06 -8.614349e-07 -0.0025967380 0.0005297566
#> [10,] 1.984923e-04 1.189269e-04 0.0020053810 -0.0731364856
#> [11,] 1.322348e-04 7.922863e-05 0.0013359766 -0.0487232279
#> [12,] 0.000000e+00 1.215418e-04 0.0020494740 -0.0747445606
#> [13,] 0.000000e+00 1.267696e-04 0.0021376270 -0.0779595134
#> [14,] 1.117660e-04 6.696476e-05 -0.0057934855 -0.0411813186
#> [15,] -2.210733e-04 -1.324563e-04 -0.0022335183 0.0814566810
#> [16,] -2.178818e-05 -2.599874e-03 -0.0112635451 0.0080280716
#> [17,] -1.152805e-04 -1.773971e-03 -0.0002262428 0.0424762629
#> [18,] 8.726538e-05 5.228515e-05 -0.0080038014 -0.0321538016
#> [19,] -4.189984e-05 -4.904047e-04 -0.0001671973 0.0154384158
#> [20,] 1.141968e-04 6.842118e-05 0.0011537382 -0.0420769701
#> [21,] -1.829435e-05 -5.156516e-04 -0.0032533364 0.0067407356
#> [22,] -1.028102e-05 -6.159884e-06 -0.0001038698 0.0037881439
#> [23,] -3.712630e-05 -2.224427e-05 -0.0003750895 0.0136795574
#> [24,] 1.596424e-04 9.564994e-05 0.0016128777 -0.0588218429
#> [25,] 1.388405e-04 8.318646e-05 0.0014027147 -0.0511571745
#> [26,] 2.934381e-05 1.758138e-05 -0.0109625031 -0.0108120197
#> [27,] 1.325986e-04 7.944662e-05 -0.0061871490 -0.0488572877
#> [28,] 8.817824e-05 5.283210e-05 -0.0076102186 -0.0324901538
#> [29,] 2.528306e-04 1.514837e-04 0.0025543640 -0.0931579607
Results$`F matrix`
#> Per X19but Fuco Pra X19hex Allo Zea Chl_b
#> Prasinophytes 0.0000 0.0000 0.0000 0.2731 0.0000 0.0000 0.0752 0.9055362
#> Chlorophytes 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0216 0.2890000
#> Cryptophytes 0.0000 0.0000 0.0000 0.0000 0.0000 0.4094 0.0000 0.0000000
#> Diatoms-2 0.0000 0.0000 0.8332 0.0000 0.0000 0.0000 0.0000 0.0000000
#> Dinoflagellates-1 0.6169 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000000
#> Haptophytes 0.0000 0.1690 0.2065 0.0000 0.6446 0.0000 0.0000 0.0000000
#> Pelagophytes 0.0000 0.7156 0.3382 0.0000 0.0000 0.0000 0.0000 0.0000000
#> Syn 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.6833 0.0000000
#> Tchla
#> Prasinophytes 1
#> Chlorophytes 1
#> Cryptophytes 1
#> Diatoms-2 1
#> Dinoflagellates-1 1
#> Haptophytes 1
#> Pelagophytes 1
#> Syn 1
Results$`Class abundances`
#> Prasinophytes Chlorophytes Cryptophytes Diatoms-2 Dinoflagellates-1
#> 1 0.05813235 0.173699887 6.077150e-03 0.05303277 2.604977e-04
#> 2 0.01683260 0.029685365 5.694085e-03 0.04258162 2.065513e-04
#> 3 0.03799919 0.156085001 1.989036e-02 0.34154307 5.616466e-04
#> 4 0.04818353 0.136702123 1.802016e-02 0.35620642 5.058270e-04
#> 5 0.04898325 0.057796684 3.732086e-03 0.35211433 6.974818e-05
#> 6 0.03950829 0.034634426 4.056271e-05 0.26034219 1.680152e-02
#> 7 0.02621901 0.006615512 1.515334e-03 0.13286536 0.000000e+00
#> 8 0.01472159 0.000000000 9.975232e-03 0.45430045 0.000000e+00
#> 9 0.01904245 0.000000000 1.459708e-01 0.68770208 3.056739e-06
#> 10 0.02631555 0.011191616 6.188038e-03 0.09197158 0.000000e+00
#> 11 0.11824951 0.106026572 6.640582e-02 0.51924046 0.000000e+00
#> 12 0.06564930 0.067287451 0.000000e+00 0.26259584 0.000000e+00
#> 13 0.06345669 0.057949902 0.000000e+00 0.15156387 0.000000e+00
#> 14 0.15068699 0.000000000 2.010173e-02 0.26754191 0.000000e+00
#> 15 0.20289955 0.048948517 9.140034e-02 0.81583561 6.624879e-04
#> 16 0.22225077 0.000000000 1.785929e-02 0.24477086 2.693010e-05
#> 17 0.44273030 0.180860112 4.334287e-02 0.79566465 4.276621e-04
#> 18 0.14655057 0.000000000 2.290540e-02 0.24203776 0.000000e+00
#> 19 0.29162329 0.015301751 4.868063e-02 0.49672528 9.623219e-05
#> 20 0.17637489 0.157065855 5.498201e-02 0.36128725 0.000000e+00
#> 21 0.27184723 0.000000000 4.146849e-02 0.44495457 3.714952e-05
#> 22 0.12313524 0.157816110 4.479589e-02 0.70148577 2.530729e-05
#> 23 0.09903301 0.120858642 3.554667e-02 0.48970550 6.793005e-05
#> 24 0.04255409 0.024762495 2.412683e-02 0.14506404 0.000000e+00
#> 25 0.06244151 0.075964548 3.038085e-02 0.18590192 0.000000e+00
#> 26 0.08383819 0.000000000 1.981428e-02 0.19132020 0.000000e+00
#> 27 0.03943359 0.000000000 1.295367e-02 0.12001427 0.000000e+00
#> 28 0.04207187 0.000000000 1.117714e-02 0.10439305 0.000000e+00
#> 29 0.06339606 0.071545928 1.875162e-02 0.06002747 0.000000e+00
#> Haptophytes Pelagophytes Syn
#> 1 0.15279608 0.01451126 0.0000000000
#> 2 0.03480569 0.01437420 0.0015298918
#> 3 0.03630156 0.02031919 0.0000000000
#> 4 0.03825464 0.02282731 0.0000000000
#> 5 0.03431180 0.02872944 0.0044626560
#> 6 0.02994669 0.02442632 0.0000000000
#> 7 0.03427013 0.01021466 0.0000000000
#> 8 0.19385841 0.01415085 0.0014934567
#> 9 0.20513559 0.02401207 0.0030339356
#> 10 0.05649731 0.01333526 0.0005006533
#> 11 0.17049977 0.02232296 0.0369156295
#> 12 0.04046404 0.01537199 0.0298356944
#> 13 0.03633243 0.01816138 0.0134436168
#> 14 0.03085693 0.01120814 0.0088944616
#> 15 0.18288037 0.00000000 0.0385373286
#> 16 0.06334808 0.02285901 0.0000000000
#> 17 0.23435666 0.04395723 0.0000000000
#> 18 0.11728874 0.04197639 0.0069118432
#> 19 0.17602370 0.04548007 0.0000000000
#> 20 0.18422583 0.05435068 0.0026144716
#> 21 0.14113121 0.05260976 0.0000000000
#> 22 0.20034659 0.01444466 0.0074189117
#> 23 0.14676234 0.02227286 0.0048195406
#> 24 0.08155547 0.01939145 0.0112013936
#> 25 0.07741931 0.01637514 0.0071589919
#> 26 0.03323533 0.01511917 0.0064464230
#> 27 0.03125042 0.02370004 0.0093665315
#> 28 0.03212207 0.01297398 0.0089689786
#> 29 0.06330828 0.02160906 0.0028305532