In [1]:

    

import sys
sys.path.append('/Users/simon/git/lda/code/')
sys.path.append('/Users/simon/git/efcompute/')
%load_ext autoreload
%autoreload 2

import numpy as np

import plotly as plotly
from plotly.graph_objs import *
plotly.offline.init_notebook_mode()


    

In [5]:

    

datadir = '/Users/simon/Dropbox/BioResearch/Meta_clustering/MS2LDA/fingerid-104-traindata/spectra_massbank/'
import glob
in_files = glob.glob(datadir + '*.ms')
print len(in_files)


    

    




2132

In [6]:

    

import pickle
with open('formulas.dict','r') as f:
    formulas = pickle.load(f)
print [f for f in formulas if f == 'C4H8N']


    

    




['C4H8N']

In [7]:

    

from feature_extraction import CorpusMaker
c = CorpusMaker('gnps',in_files,min_intensity = 500,min_loss = 10.0,max_loss = 200.0,algorithm='kde',replace='max',formulas = formulas)


    

    




Processed 100 spectra
Processed 200 spectra
Processed 300 spectra
Processed 400 spectra
Processed 500 spectra
Processed 600 spectra
Processed 700 spectra
Processed 800 spectra
Processed 900 spectra
Processed 1000 spectra
Processed 1100 spectra
Processed 1200 spectra
Processed 1300 spectra
Processed 1400 spectra
Processed 1500 spectra
Processed 1600 spectra
Processed 1700 spectra
Processed 1800 spectra
Processed 1900 spectra
Processed 2000 spectra
Processed 2100 spectra
Done kde for 0 of 65793
Done kde for 1000 of 65793
Done kde for 2000 of 65793
Done kde for 3000 of 65793
Done kde for 4000 of 65793
Done kde for 5000 of 65793
Done kde for 6000 of 65793
Done kde for 7000 of 65793
Done kde for 8000 of 65793
Done kde for 9000 of 65793
Done kde for 10000 of 65793
Done kde for 11000 of 65793
Done kde for 12000 of 65793
Done kde for 13000 of 65793
Done kde for 14000 of 65793
Done kde for 15000 of 65793
Done kde for 16000 of 65793
Done kde for 17000 of 65793
Done kde for 18000 of 65793
Done kde for 19000 of 65793
Done kde for 20000 of 65793
Done kde for 21000 of 65793
Done kde for 22000 of 65793
Done kde for 23000 of 65793
Done kde for 24000 of 65793
Done kde for 25000 of 65793
Done kde for 26000 of 65793
Done kde for 27000 of 65793
Done kde for 28000 of 65793
Done kde for 29000 of 65793
Done kde for 30000 of 65793
Done kde for 31000 of 65793
Done kde for 32000 of 65793
Done kde for 33000 of 65793
Done kde for 34000 of 65793
Done kde for 35000 of 65793
Done kde for 36000 of 65793
Done kde for 37000 of 65793
Done kde for 38000 of 65793
Done kde for 39000 of 65793
Done kde for 40000 of 65793
Done kde for 41000 of 65793
Done kde for 42000 of 65793
Done kde for 43000 of 65793
Done kde for 44000 of 65793
Done kde for 45000 of 65793
Done kde for 46000 of 65793
Done kde for 47000 of 65793
Done kde for 48000 of 65793
Done kde for 49000 of 65793
Done kde for 50000 of 65793
Done kde for 51000 of 65793
Done kde for 52000 of 65793
Done kde for 53000 of 65793
Done kde for 54000 of 65793
Done kde for 55000 of 65793
Done kde for 56000 of 65793
Done kde for 57000 of 65793
Done kde for 58000 of 65793
Done kde for 59000 of 65793
Done kde for 60000 of 65793
Done kde for 61000 of 65793
Done kde for 62000 of 65793
Done kde for 63000 of 65793
Done kde for 64000 of 65793
Done kde for 65000 of 65793
Made kde
Performing kernel density estimation
Atoms being considered = ['H', 'C', 'N', 'O']
Atoms being considered = ['H', 'C', 'N', 'O', 'P']
Atoms being considered = ['H', 'C', 'N', 'O', 'P']
Found 100 groups
Found 200 groups
Found 300 groups
Found 400 groups
Found 500 groups
Found 600 groups
Found 700 groups
Found 800 groups
Found 900 groups
Found 1000 groups
Found 1100 groups
Found 1200 groups
Found 1300 groups
Found 1400 groups
Found 1500 groups
Found 1600 groups
Found 1700 groups
Found 1800 groups
Found 1900 groups
Found 2000 groups
Found 2100 groups
Found 2200 groups
Found 2300 groups
Found 2400 groups
Found 2500 groups
Found 2600 groups
Found 2700 groups
Found 2800 groups
Found 2900 groups
Found 3000 groups
Found 3100 groups
Found 3200 groups
Found 3300 groups
Found 3400 groups
Found 3500 groups
Found 3600 groups
Found 3700 groups
Found 3800 groups
Found 3900 groups
Found 4000 groups
Found 4100 groups
Found 4200 groups
Found 4300 groups
Found 4400 groups
Found 4500 groups
Found 4600 groups
Found 4700 groups
Found 4800 groups
Found 4900 groups
Found 5000 groups
Found 5100 groups
Found 5200 groups
Found 5300 groups
Found 5400 groups
Found 5500 groups
Found 5600 groups
Found 5700 groups
Found 5800 groups
Found 5900 groups
Found 6000 groups
Found 6100 groups
Found 6200 groups
Found 6300 groups
Found 6400 groups
Found 6500 groups
Found 6600 groups
Found 6700 groups
Found 6800 groups
Found 6900 groups
Found 7000 groups
Found 7100 groups
Found 7200 groups
Found 7300 groups
Found 7400 groups
Found 7500 groups
Found 7600 groups
Found 7700 groups
Found 7800 groups
Finished losses, total words now: 7870
Done kde for 0 of 109045
Done kde for 1000 of 109045
Done kde for 2000 of 109045
Done kde for 3000 of 109045
Done kde for 4000 of 109045
Done kde for 5000 of 109045
Done kde for 6000 of 109045
Done kde for 7000 of 109045
Done kde for 8000 of 109045
Done kde for 9000 of 109045
Done kde for 10000 of 109045
Done kde for 11000 of 109045
Done kde for 12000 of 109045
Done kde for 13000 of 109045
Done kde for 14000 of 109045
Done kde for 15000 of 109045
Done kde for 16000 of 109045
Done kde for 17000 of 109045
Done kde for 18000 of 109045
Done kde for 19000 of 109045
Done kde for 20000 of 109045
Done kde for 21000 of 109045
Done kde for 22000 of 109045
Done kde for 23000 of 109045
Done kde for 24000 of 109045
Done kde for 25000 of 109045
Done kde for 26000 of 109045
Done kde for 27000 of 109045
Done kde for 28000 of 109045
Done kde for 29000 of 109045
Done kde for 30000 of 109045
Done kde for 31000 of 109045
Done kde for 32000 of 109045
Done kde for 33000 of 109045
Done kde for 34000 of 109045
Done kde for 35000 of 109045
Done kde for 36000 of 109045
Done kde for 37000 of 109045
Done kde for 38000 of 109045
Done kde for 39000 of 109045
Done kde for 40000 of 109045
Done kde for 41000 of 109045
Done kde for 42000 of 109045
Done kde for 43000 of 109045
Done kde for 44000 of 109045
Done kde for 45000 of 109045
Done kde for 46000 of 109045
Done kde for 47000 of 109045
Done kde for 48000 of 109045
Done kde for 49000 of 109045
Done kde for 50000 of 109045
Done kde for 51000 of 109045
Done kde for 52000 of 109045
Done kde for 53000 of 109045
Done kde for 54000 of 109045
Done kde for 55000 of 109045
Done kde for 56000 of 109045
Done kde for 57000 of 109045
Done kde for 58000 of 109045
Done kde for 59000 of 109045
Done kde for 60000 of 109045
Done kde for 61000 of 109045
Done kde for 62000 of 109045
Done kde for 63000 of 109045
Done kde for 64000 of 109045
Done kde for 65000 of 109045
Done kde for 66000 of 109045
Done kde for 67000 of 109045
Done kde for 68000 of 109045
Done kde for 69000 of 109045
Done kde for 70000 of 109045
Done kde for 71000 of 109045
Done kde for 72000 of 109045
Done kde for 73000 of 109045
Done kde for 74000 of 109045
Done kde for 75000 of 109045
Done kde for 76000 of 109045
Done kde for 77000 of 109045
Done kde for 78000 of 109045
Done kde for 79000 of 109045
Done kde for 80000 of 109045
Done kde for 81000 of 109045
Done kde for 82000 of 109045
Done kde for 83000 of 109045
Done kde for 84000 of 109045
Done kde for 85000 of 109045
Done kde for 86000 of 109045
Done kde for 87000 of 109045
Done kde for 88000 of 109045
Done kde for 89000 of 109045
Done kde for 90000 of 109045
Done kde for 91000 of 109045
Done kde for 92000 of 109045
Done kde for 93000 of 109045
Done kde for 94000 of 109045
Done kde for 95000 of 109045
Done kde for 96000 of 109045
Done kde for 97000 of 109045
Done kde for 98000 of 109045
Done kde for 99000 of 109045
Done kde for 100000 of 109045
Done kde for 101000 of 109045
Done kde for 102000 of 109045
Done kde for 103000 of 109045
Done kde for 104000 of 109045
Done kde for 105000 of 109045
Done kde for 106000 of 109045
Done kde for 107000 of 109045
Done kde for 108000 of 109045
Done kde for 109000 of 109045
Made kde
Performing kernel density estimation
Atoms being considered = ['H', 'C', 'N', 'O']
Atoms being considered = ['H', 'C', 'N', 'O', 'P']
Atoms being considered = ['H', 'C', 'N', 'O', 'P']
Found 100 groups
Found 200 groups
Found 300 groups
Found 400 groups
Found 500 groups
Found 600 groups
Found 700 groups
Found 800 groups
Found 900 groups
Found 1000 groups
Found 1100 groups
Found 1200 groups
Found 1300 groups
Found 1400 groups
Found 1500 groups
Found 1600 groups
Found 1700 groups
Found 1800 groups
Found 1900 groups
Found 2000 groups
Found 2100 groups
Found 2200 groups
Found 2300 groups
Found 2400 groups
Found 2500 groups
Found 2600 groups
Found 2700 groups
Found 2800 groups
Found 2900 groups
Found 3000 groups
Found 3100 groups
Found 3200 groups
Found 3300 groups
Found 3400 groups
Found 3500 groups
Found 3600 groups
Found 3700 groups
Found 3800 groups
Found 3900 groups
Found 4000 groups
Found 4100 groups
Found 4200 groups
Found 4300 groups
Found 4400 groups
Found 4500 groups
Found 4600 groups
Found 4700 groups
Found 4800 groups
Found 4900 groups
Found 5000 groups
Found 5100 groups
Found 5200 groups
Found 5300 groups
Found 5400 groups
Found 5500 groups
Found 5600 groups
Found 5700 groups
Found 5800 groups
Found 5900 groups
Found 6000 groups
Found 6100 groups
Found 6200 groups
Found 6300 groups
Found 6400 groups
Found 6500 groups
Found 6600 groups
Found 6700 groups
Found 6800 groups
Found 6900 groups
Found 7000 groups
Found 7100 groups
Found 7200 groups
Found 7300 groups
Found 7400 groups
Found 7500 groups
Found 7600 groups
Found 7700 groups
Found 7800 groups
Found 7900 groups
Found 8000 groups
Found 8100 groups
Found 8200 groups
Found 8300 groups
Found 8400 groups
Found 8500 groups
Found 8600 groups
Found 8700 groups
Found 8800 groups
Found 8900 groups
Found 9000 groups
Found 9100 groups
Found 9200 groups
Found 9300 groups
Found 9400 groups
Found 9500 groups
Found 9600 groups
Found 9700 groups
Found 9800 groups
Found 9900 groups
Found 10000 groups
Found 10100 groups
Found 10200 groups
Found 10300 groups
Found 10400 groups
Found 10500 groups
Found 10600 groups
Found 10700 groups
Found 10800 groups
Found 10900 groups
Found 11000 groups
Found 11100 groups
Found 11200 groups
Found 11300 groups
Found 11400 groups
Found 11500 groups
Found 11600 groups
Found 11700 groups
Found 11800 groups
Found 11900 groups
Found 12000 groups
Found 12100 groups
Found 12200 groups
Found 12300 groups
Found 12400 groups
Found 12500 groups
Found 12600 groups
Found 12700 groups
Found 12800 groups
Found 12900 groups
Found 13000 groups
Found 13100 groups
Found 13200 groups
Found 13300 groups
Found 13400 groups
Found 13500 groups
Found 13600 groups
Found 13700 groups
Found 13800 groups
Found 13900 groups
Found 14000 groups
Found 14100 groups
Found 14200 groups
Found 14300 groups
Found 14400 groups
Found 14500 groups
Finished fragments, total words now: 22415

In [8]:

    

with open('annotated_beer3_topics.dict','r') as f:
    beer_topics = pickle.load(f)


    

In [9]:

    

from lda import VariationalLDA
vlda = VariationalLDA(c.get_first_file(),K = 500,normalise = 1000.0)
vlda.add_fixed_topics_formulas(beer_topics)


    

    




Found 22415 unique words
Object created with 1953 documents
Normalising intensities
Matching topics based on formulas
Topic: motif_19, 0.617635834566 probability matched (Fragments indicative for ferulic acid based substructure (MzCloud))
Match accepted, storing as motif_0
Topic: motif_23, 0.372373516156 probability matched (Unclear yet what these losses relate to.)
Topic: motif_53, 0.0 probability matched (Losses indicative of a hexose with NH2 group – EF fits)
Topic: motif_68, 0.215117619926 probability matched (Possibly suggests iminosugar like substructure. Fragments have losses (H2O CO) related to sugar (polyhydroxylated structure))
Topic: motif_69, 0.823195681326 probability matched (Fragment ions possibly indicative for ribose substructure (MzCloud))
Match accepted, storing as motif_1
Topic: motif_40, 0.636447126738 probability matched (Imidazole group linked to a carboxylgroup through one CH2 group i.e. like in imidazole acetic acid - Prevalent in Beer1)
Match accepted, storing as motif_2
Topic: motif_41, 0.754906735499 probability matched (Unclear yet what these fragments relate to.)
Match accepted, storing as motif_3
Topic: motif_42, 0.604603544447 probability matched (Fragments indicative for asparagine substructure (MzCloud) – prevalent in Beer 3)
Match accepted, storing as motif_4
Topic: motif_45, 0.830947583677 probability matched (Fragment ions indicative for pipecolic acid (pipecolate) (MzCloud) - Quite prevalent especially in Beer2)
Match accepted, storing as motif_5
Topic: motif_60, 0.912849221206 probability matched (Double water loss i.e. 2*H20 – Generic feature for metabolites containing several free OH groups attached to a aliphatic chain like sugars.)
Match accepted, storing as motif_6
Topic: motif_162, 0.537748774156 probability matched (Fragment ions indicative for guanine (C5H5N5O) based substructure)
Match accepted, storing as motif_7
Topic: motif_249, 0.344138219021 probability matched (Unclear yet what these fragments relate to.)
Topic: motif_284, 0.70321297041 probability matched (Fragments indicative for dihydroxylated benzene ring substructure (MzCloud) – C6H5O2 fragment corresponds to positively charged fragment with two hydroxyl groups.)
Match accepted, storing as motif_8
Topic: motif_139, 0.495712943387 probability matched (Unclear yet what these fragments relate to.)
Topic: motif_115, 0.6874516628 probability matched (Fragments indicative for [phenylalanine-CHOOH] based substructure)
Match accepted, storing as motif_9
Topic: motif_280, 0.741039520706 probability matched (Loss possibly indicative of carboxylic acid group with 1-carbon attached.)
Match accepted, storing as motif_10
Topic: motif_243, 0.736018716879 probability matched (Fragment ion indicative for conjugation of a phosphate group  (H4O4P))
Match accepted, storing as motif_11
Topic: motif_272, 0.590570443678 probability matched (Unclear if fragments relate to a specific substructure.)
Match accepted, storing as motif_12
Topic: motif_241, 0.515158360819 probability matched (Fragments indicative for histidine (C6H10N3O2) substructure (MzCloud))
Match accepted, storing as motif_13
Topic: motif_136, 0.556527723162 probability matched (Unclear what these fragments relate to. Possibly small ring structure (CHN loss))
Match accepted, storing as motif_14
Topic: motif_131, 0.780011746548 probability matched (Fragments indicative of a [hexose-H2O] substructure – i.e. indicative for a hexose (like glucose) conjugation (MzCloud))
Match accepted, storing as motif_15
Topic: motif_130, 0.347917279953 probability matched (Fragment ions indicative for cytosine and a loss of conjugated deoxyribose – possibly combined due to many spectra that combine these two substructures. Loss of NH2 group is likely from remaining fragment after loss of deoxyribose. – Quite prevalent in Beer1)
Topic: motif_156, 0.68298797801 probability matched ([Ribose (pentose C5-sugar)-H2O] related loss –  indicative for conjugated ribose sugar - EF fits)
Match accepted, storing as motif_16
Topic: motif_154, 0.0 probability matched (C13 isotope peaks of proline (abundant ions taken for fragmentation))
Topic: motif_260, 0.941792063677 probability matched (Water loss  - indicative of a free hydroxyl group) – often seen in sugary structures)
Match accepted, storing as motif_17
Topic: motif_209, 0.0298432273405 probability matched (Possibly suggests 2-pyyrolidine substructure – Mass2Motif not consistent over the four beers.)
Topic: motif_176, 0.889775721552 probability matched (Small abundant and generic aromatic fragment found across several mass patterns.)
Match accepted, storing as motif_18
Topic: motif_97, 0.852554253889 probability matched (Oxygen-rich losses and fragments also occurring in hexose spectra – related to M2M 211 (hexose [glucose] conjugatation) – possibly hydrated-hexose loss?)
Match accepted, storing as motif_19
Topic: motif_188, 0.361525086935 probability matched (Unclear yet – possibly related to methylsuccinic acid….)
Topic: motif_262, 0.856459026116 probability matched (Combined loss of H2O and CO – indicative for free carboxylic acid group (COOH) – generic substructure in amino acids and organic acids)
Match accepted, storing as motif_20
Topic: motif_222, 0.170204617573 probability matched (Fragments (and loss) indicative for serine substructure (MzCloud) - Present in Beer 2 & Beer 3)
Topic: motif_67, 0.153752853881 probability matched (Fragment ions indicative for aconitic acid substructure (C3HO2 fragment is quite specific))
Topic: motif_202, 0.406436487639 probability matched (Fragments indicative of [tryptophan-NH3] related substructure (C8H8N is the basic indole skeleton a fused benzene and 5 membered N-containing ring))
Topic: motif_227, 0.554095045501 probability matched (Possibly related to 14-Cyclohex-2-enedione substructure – double CO loss between fragments. Could be alkaloid fragments as well.)
Match accepted, storing as motif_21
Topic: motif_158, 0.842882860513 probability matched (Leucine related substructure (mzCloud) – prevalent in Beer 2)
Match accepted, storing as motif_22
Topic: motif_151, 0.573826369881 probability matched (Fragment and loss of [proline-H2O] - indicative for conjugated proline – EF fits)
Match accepted, storing as motif_23
Topic: motif_2, 0.979260119981 probability matched (Small nitrogen containing fragment ion – often proline or ornithine derived – most abundant fragment in all four beers)
Match accepted, storing as motif_24
Topic: motif_35, 0.625997207818 probability matched (Unclear yet what these fragments relate to.)
Match accepted, storing as motif_25
Topic: motif_34, 0.692840948727 probability matched (Combination of small losses – free carboxylgroup and acetyl group loss + loss of NH2 group in some cases)
Match accepted, storing as motif_26
Topic: motif_98, 0.837618169241 probability matched (Fragment ions possibly indicative for N-Methyl-oxo-pyrrolidinecarboxylic acid like structure (loss of free carboxyl group))
Match accepted, storing as motif_27
Topic: motif_290, 0.812259068894 probability matched (Two small fragments unclear if they represent a substructure)
Match accepted, storing as motif_28
Topic: motif_17, 0.490827038163 probability matched ( Fragments indicative for tyrosine substructure (MzCloud))
Topic: motif_221, 0.836582863007 probability matched (Unclear yet what these fragments relate to.)
Match accepted, storing as motif_29
Topic: motif_58, 0.749255271738 probability matched (Fragments indicative for ethylphenol substructure (i.e. resulting from Tyramine – MzCloud))
Match accepted, storing as motif_30
Topic: motif_13, 0.0 probability matched (Fragment ions indicative for glutamine (C5H11N2O3) substructure)
Topic: motif_12, 0.834044885479 probability matched (Combination of small losses (CO2 H2O etc.) – Unclear if they relate to a substructure loss.)
Match accepted, storing as motif_31
Topic: motif_220, 0.727018288858 probability matched (Fragments indicative adenine (C5H6N5) substructure – most prevalent in Beer3)
Match accepted, storing as motif_32
Topic: motif_79, 0.399986624758 probability matched (Fragment ions indicative for trimethylated amine connected to a carboxylic acid group i.e. like in betaine (MzCloud))
Topic: motif_291, 0.311336540473 probability matched (Unclear if fragments relate to a specific substructure. H2O loss between fragments.)
Topic: motif_91, 0.45862793319 probability matched (Fragments possibly suggest benzene ring substituted with one hydroxyl and one NH2 group (fragments point to orientation from 3-hydroxyanthranilic acid – i.e. MzCloud) – prevalent in Beer 2 and 4)
Topic: motif_294, 0.0810948765393 probability matched (Unclear what these fragments relate to.)
Topic: motif_166, 0.430940524288 probability matched (Fragment ions indicative for 5-methylcytosine substructure (MzCloud) – prevalent in Beer 2)
Topic: motif_55, 0.689722307613 probability matched (Fragments indicative for cinnamic acid (cinnamate) substructure (MzCloud))
Match accepted, storing as motif_33
Topic: motif_104, 0.440421892816 probability matched (Unclear what these fragments relate to.)
Topic: motif_149, 0.453024526262 probability matched (Possibly suggests 2-pyridone/ol substructure)
Topic: motif_7, 0.532447426139 probability matched (Possibly suggests 5-aminopentanol substructure)
Match accepted, storing as motif_34
Topic: motif_102, 0.461427445467 probability matched (Possibly suggests alkylbenzene substructure.)
Topic: motif_128, 0.734015213816 probability matched (N/O-Acetylation (Beer 3))
Match accepted, storing as motif_35
Topic: motif_4, 0.0 probability matched (Unclear yet what these fragments relate to.)
Topic: motif_174, 0.955353233105 probability matched (Fragment ions indicative for pyroglutamic acid (pyroglutamate) or glutamine (both in MzCloud) – structure can be formed from glutamic acid (glutamate) in the mass spectrometer as well.)
Match accepted, storing as motif_36
Topic: motif_207, 0.528437815976 probability matched (Unclear what these fragments relate to.)
Match accepted, storing as motif_37
Topic: motif_250, 0.749514110239 probability matched (Combination of small losses (i.e. NH3 and CH2O2))
Match accepted, storing as motif_38
Topic: motif_121, 0.274479448874 probability matched (Possibly suggests 4-aminooxane-4-carboxylic acid like substructure?)
Topic: motif_276, 0.740255663984 probability matched (Alkyl aromatic substructure – indicative for aromatic ring with 2-carbon alkyl chain attached i.e. phenylethene fragment from ethylbenzene as a result of the fragmentation process.)
Match accepted, storing as motif_39
Topic: motif_211, 0.731347980318 probability matched (Fragment ions indicative for N-acetylputrescine substructure (MzCloud))
Match accepted, storing as motif_40
Topic: motif_168, 0.334416991764 probability matched (Unclear yet what these fragments relate to.)
Topic: motif_72, 0.366283059705 probability matched (Possibly suggests methoxy-1H-indole-23-dione (methoxy-isatin) substructure)
Topic: motif_214, 0.882578554434 probability matched (Amine loss - Indicative for free NH2 group in fragmented molecule)
Match accepted, storing as motif_41
Topic: motif_195, 0.874175848988 probability matched (Aliphatic amine (NH3 loss indicates free NH2 group coupled to aliphatic chain))
Match accepted, storing as motif_42
Topic: motif_270, 0.0428663014337 probability matched (Unclear yet what these fragments relate to.)
Topic: motif_233, 0.383066818961 probability matched (Combination of small losses (CO2 H2O etc.) – Unclear if they relate to a substructure loss.)
Topic: motif_230, 0.788598574822 probability matched (Fragments indicative for y-aminobutyric acid (amimobutyrate) substructure (MzCloud) – present in Beer 2 & Beer 3 – in beer 3 mainly based on C4H7O2 fragment.)
Match accepted, storing as motif_43
Topic: motif_226, 0.797411605481 probability matched (Loss of [hexose-H2O] – indication of hexose conjugation (for example glucose))
Match accepted, storing as motif_44
Topic: motif_82, 0.302499850129 probability matched (Possibly suggests agmatine based substructure (C5H15N4) with unknown conjugation….)
Topic: motif_90, 0.363529613931 probability matched (Combinations of small generic losses like CH2O2 + CH4O = C2H6O3)

In [10]:

    

vlda.run_vb(n_its = 1000)


    

    




Initialising
Starting iterations
Iteration 0 (change = 568.742672439) (8.523331 seconds, I think I'll finish in 142.055516667 minutes)
Iteration 1 (change = 44.7966636041) (7.99451 seconds, I think I'll finish in 133.1085915 minutes)
Iteration 2 (change = 59.0949779968) (6.935181 seconds, I think I'll finish in 115.3551773 minutes)
Iteration 3 (change = 85.0372323891) (6.895802 seconds, I think I'll finish in 114.585243233 minutes)
Iteration 4 (change = 131.379432109) (6.927623 seconds, I think I'll finish in 114.9985418 minutes)
Iteration 5 (change = 201.31674332) (7.231017 seconds, I think I'll finish in 119.91436525 minutes)
Iteration 6 (change = 253.524984811) (7.281971 seconds, I think I'll finish in 120.637986233 minutes)
Iteration 7 (change = 224.653527426) (7.720328 seconds, I think I'll finish in 127.7714284 minutes)
Iteration 8 (change = 155.075221082) (7.411892 seconds, I think I'll finish in 122.543281067 minutes)
Iteration 9 (change = 99.9469279427) (7.627754 seconds, I think I'll finish in 125.985070233 minutes)
Iteration 10 (change = 65.3990553115) (7.885334 seconds, I think I'll finish in 130.108011 minutes)
Iteration 11 (change = 45.4961395477) (7.736636 seconds, I think I'll finish in 127.525550067 minutes)
Iteration 12 (change = 33.0501402398) (7.711667 seconds, I think I'll finish in 126.985449933 minutes)
Iteration 13 (change = 24.7344051558) (8.081112 seconds, I think I'll finish in 132.9342924 minutes)
Iteration 14 (change = 19.4245604766) (7.815317 seconds, I think I'll finish in 128.431709367 minutes)
Iteration 15 (change = 15.3848851756) (7.578243 seconds, I think I'll finish in 124.40948925 minutes)
Iteration 16 (change = 12.2151229649) (7.576874 seconds, I think I'll finish in 124.2607336 minutes)
Iteration 17 (change = 9.72158882041) (7.651091 seconds, I think I'll finish in 125.350374217 minutes)
Iteration 18 (change = 7.98258883429) (7.562184 seconds, I think I'll finish in 123.7677448 minutes)
Iteration 19 (change = 7.01272695779) (8.011844 seconds, I think I'll finish in 130.9936494 minutes)
Iteration 20 (change = 6.40529765458) (7.645175 seconds, I think I'll finish in 124.871191667 minutes)
Iteration 21 (change = 5.68209344981) (7.582463 seconds, I think I'll finish in 123.720521283 minutes)
Iteration 22 (change = 5.11495891989) (7.686282 seconds, I think I'll finish in 125.2863966 minutes)
Iteration 23 (change = 4.61724205242) (7.685362 seconds, I think I'll finish in 125.143311233 minutes)
Iteration 24 (change = 4.0422730283) (7.494076 seconds, I think I'll finish in 121.903636267 minutes)
Iteration 25 (change = 3.59915108833) (7.611305 seconds, I think I'll finish in 123.68370625 minutes)
Iteration 26 (change = 3.27492182938) (7.507082 seconds, I think I'll finish in 121.864964467 minutes)
Iteration 27 (change = 2.96662494154) (7.561596 seconds, I think I'll finish in 122.6238818 minutes)
Iteration 28 (change = 2.67001750896) (7.527933 seconds, I think I'll finish in 121.9525146 minutes)
Iteration 29 (change = 2.46502306481) (7.693579 seconds, I think I'll finish in 124.507753483 minutes)
Iteration 30 (change = 2.33189249205) (7.688314 seconds, I think I'll finish in 124.294409667 minutes)
Iteration 31 (change = 2.16157642284) (7.597374 seconds, I think I'll finish in 122.6975901 minutes)
Iteration 32 (change = 2.06956062894) (7.533102 seconds, I think I'll finish in 121.5340456 minutes)
Iteration 33 (change = 2.03851294973) (7.523398 seconds, I think I'll finish in 121.252097767 minutes)
Iteration 34 (change = 2.00450589571) (7.697989 seconds, I think I'll finish in 123.9376229 minutes)
Iteration 35 (change = 1.99475084821) (7.573346 seconds, I think I'll finish in 121.804648167 minutes)
Iteration 36 (change = 2.05967277353) (7.669309 seconds, I think I'll finish in 123.220231267 minutes)
Iteration 37 (change = 1.98743865837) (7.574826 seconds, I think I'll finish in 121.5759573 minutes)
Iteration 38 (change = 1.8328573211) (7.546532 seconds, I think I'll finish in 120.996063067 minutes)
Iteration 39 (change = 1.57229742319) (7.786263 seconds, I think I'll finish in 124.70997905 minutes)
Iteration 40 (change = 1.30960496564) (7.563227 seconds, I think I'll finish in 121.011632 minutes)
Iteration 41 (change = 1.2244078664) (7.591329 seconds, I think I'll finish in 121.33474185 minutes)
Iteration 42 (change = 1.18460988859) (7.573065 seconds, I think I'll finish in 120.9166045 minutes)
Iteration 43 (change = 1.19143353023) (7.516205 seconds, I think I'll finish in 119.88346975 minutes)
Iteration 44 (change = 1.25533425499) (7.623199 seconds, I think I'll finish in 121.462970733 minutes)
Iteration 45 (change = 1.27735959498) (7.637552 seconds, I think I'll finish in 121.564369333 minutes)
Iteration 46 (change = 1.23958676298) (7.897739 seconds, I think I'll finish in 125.5740501 minutes)
Iteration 47 (change = 1.11534903758) (8.001067 seconds, I think I'll finish in 127.083614183 minutes)
Iteration 48 (change = 0.94366571387) (7.992958 seconds, I think I'll finish in 126.821600267 minutes)
Iteration 49 (change = 0.792556624695) (8.572418 seconds, I think I'll finish in 135.8728253 minutes)
Iteration 50 (change = 0.697883624207) (8.037858 seconds, I think I'll finish in 127.266085 minutes)
Iteration 51 (change = 0.650610905123) (8.624029 seconds, I think I'll finish in 136.403392017 minutes)
Iteration 52 (change = 0.66055781521) (7.4908 seconds, I think I'll finish in 118.35464 minutes)
Iteration 53 (change = 0.69414558744) (7.67099 seconds, I think I'll finish in 121.073792167 minutes)
Iteration 54 (change = 0.731983501321) (7.893651 seconds, I think I'll finish in 124.4565641 minutes)
Iteration 55 (change = 0.745769633527) (7.742607 seconds, I think I'll finish in 121.94606025 minutes)
Iteration 56 (change = 0.696372752013) (7.676616 seconds, I think I'll finish in 120.7787584 minutes)
Iteration 57 (change = 0.633007783022) (7.805183 seconds, I think I'll finish in 122.671459483 minutes)
Iteration 58 (change = 0.619643750559) (7.663074 seconds, I think I'll finish in 120.3102618 minutes)
Iteration 59 (change = 0.627457642898) (7.754679 seconds, I think I'll finish in 121.61921565 minutes)
Iteration 60 (change = 0.648934911928) (7.817396 seconds, I think I'll finish in 122.472537333 minutes)
Iteration 61 (change = 0.649699959575) (7.916766 seconds, I think I'll finish in 123.8973879 minutes)
Iteration 62 (change = 0.560210805182) (7.920188 seconds, I think I'll finish in 123.818939067 minutes)
Iteration 63 (change = 0.462805902228) (7.673511 seconds, I think I'll finish in 119.83466345 minutes)
Iteration 64 (change = 0.43360138444) (7.670503 seconds, I think I'll finish in 119.6598468 minutes)
Iteration 65 (change = 0.377165567657) (7.71829 seconds, I think I'll finish in 120.276685833 minutes)
Iteration 66 (change = 0.324188862394) (7.695904 seconds, I think I'll finish in 119.799572267 minutes)
Iteration 67 (change = 0.305232587995) (7.646352 seconds, I think I'll finish in 118.9007736 minutes)
Iteration 68 (change = 0.299443926181) (7.685413 seconds, I think I'll finish in 119.380081933 minutes)
Iteration 69 (change = 0.306491347931) (7.780985 seconds, I think I'll finish in 120.734950583 minutes)
Iteration 70 (change = 0.328201044445) (7.81942 seconds, I think I'll finish in 121.20101 minutes)
Iteration 71 (change = 0.369390429178) (7.727798 seconds, I think I'll finish in 119.652072367 minutes)
Iteration 72 (change = 0.422620832368) (7.710243 seconds, I think I'll finish in 119.2517584 minutes)
Iteration 73 (change = 0.421817551179) (7.718725 seconds, I think I'll finish in 119.25430125 minutes)
Iteration 74 (change = 0.374240316133) (7.639341 seconds, I think I'll finish in 117.9004961 minutes)
Iteration 75 (change = 0.38655333581) (7.659931 seconds, I think I'll finish in 118.090602917 minutes)
Iteration 76 (change = 0.422439143345) (7.669457 seconds, I think I'll finish in 118.1096378 minutes)
Iteration 77 (change = 0.472836354022) (7.99441 seconds, I think I'll finish in 122.980673833 minutes)
Iteration 78 (change = 0.52195635918) (7.795138 seconds, I think I'll finish in 119.785287267 minutes)
Iteration 79 (change = 0.550329461981) (7.678993 seconds, I think I'll finish in 117.87254255 minutes)
Iteration 80 (change = 0.531804761339) (7.558955 seconds, I think I'll finish in 115.903976667 minutes)
Iteration 81 (change = 0.401126121078) (7.770324 seconds, I think I'll finish in 119.0154626 minutes)
Iteration 82 (change = 0.27833782187) (7.699896 seconds, I think I'll finish in 117.8084088 minutes)
Iteration 83 (change = 0.248834164284) (7.844701 seconds, I think I'll finish in 119.893180283 minutes)
Iteration 84 (change = 0.234815089917) (7.85852 seconds, I think I'll finish in 119.973405333 minutes)
Iteration 85 (change = 0.224051466555) (7.679686 seconds, I think I'll finish in 117.1152115 minutes)
Iteration 86 (change = 0.212100630445) (7.731519 seconds, I think I'll finish in 117.7768061 minutes)
Iteration 87 (change = 0.200318818546) (7.658877 seconds, I think I'll finish in 116.54257835 minutes)
Iteration 88 (change = 0.192410362004) (7.610183 seconds, I think I'll finish in 115.6747816 minutes)
Iteration 89 (change = 0.187019989296) (7.849564 seconds, I think I'll finish in 119.182546733 minutes)
Iteration 90 (change = 0.184226936284) (8.022594 seconds, I think I'll finish in 121.676009 minutes)
Iteration 91 (change = 0.184488414569) (9.964705 seconds, I think I'll finish in 150.96528075 minutes)
Iteration 92 (change = 0.190069162977) (7.823315 seconds, I think I'll finish in 118.392833667 minutes)
Iteration 93 (change = 0.205507061714) (7.936868 seconds, I think I'll finish in 119.978987933 minutes)
Iteration 94 (change = 0.236541561182) (7.728191 seconds, I think I'll finish in 116.6956841 minutes)
Iteration 95 (change = 0.267027150347) (8.2312 seconds, I think I'll finish in 124.153933333 minutes)
Iteration 96 (change = 0.227057434233) (7.686874 seconds, I think I'll finish in 115.815568267 minutes)
Iteration 97 (change = 0.176843281158) (7.89551 seconds, I think I'll finish in 118.8274255 minutes)
Iteration 98 (change = 0.198587468808) (9.096054 seconds, I think I'll finish in 136.7440118 minutes)
Iteration 99 (change = 0.225161116003) (8.236661 seconds, I think I'll finish in 123.687192683 minutes)
Iteration 100 (change = 0.185477083494) (7.645539 seconds, I think I'll finish in 114.683085 minutes)
Iteration 101 (change = 0.143829238527) (7.6557 seconds, I think I'll finish in 114.707905 minutes)
Iteration 102 (change = 0.135609415771) (7.62559 seconds, I think I'll finish in 114.129663667 minutes)
Iteration 103 (change = 0.137851308598) (7.625466 seconds, I think I'll finish in 114.0007167 minutes)
Iteration 104 (change = 0.142490660422) (7.694181 seconds, I think I'll finish in 114.8997696 minutes)
Iteration 105 (change = 0.148352184334) (8.046696 seconds, I think I'll finish in 120.029882 minutes)
Iteration 106 (change = 0.155726594636) (8.025332 seconds, I think I'll finish in 119.5774468 minutes)
Iteration 107 (change = 0.157454823712) (8.465775 seconds, I think I'll finish in 125.99895125 minutes)
Iteration 108 (change = 0.153487696283) (8.601755 seconds, I think I'll finish in 127.879424333 minutes)
Iteration 109 (change = 0.152366332803) (9.40379 seconds, I think I'll finish in 139.6462815 minutes)
Iteration 110 (change = 0.152223526773) (9.467647 seconds, I think I'll finish in 140.436763833 minutes)
Iteration 111 (change = 0.149459449202) (10.047117 seconds, I think I'll finish in 148.86478355 minutes)
Iteration 112 (change = 0.140184364045) (10.431666 seconds, I think I'll finish in 154.3886568 minutes)
Iteration 113 (change = 0.128835769007) (10.850441 seconds, I think I'll finish in 160.405686117 minutes)
Iteration 114 (change = 0.11834160095) (11.080603 seconds, I think I'll finish in 163.623570967 minutes)
Iteration 115 (change = 0.108725620001) (11.264465 seconds, I think I'll finish in 166.15085875 minutes)
Iteration 116 (change = 0.0993621602424) (11.130906 seconds, I think I'll finish in 163.9953484 minutes)
Iteration 117 (change = 0.0900744475912) (10.953414 seconds, I think I'll finish in 161.1977427 minutes)
Iteration 118 (change = 0.0815855844793) (10.743175 seconds, I think I'll finish in 157.9246725 minutes)
Iteration 119 (change = 0.074762678041) (11.203126 seconds, I think I'll finish in 164.499233433 minutes)
Iteration 120 (change = 0.0700377345257) (10.837401 seconds, I think I'll finish in 158.948548 minutes)
Iteration 121 (change = 0.0670432948912) (10.764615 seconds, I think I'll finish in 157.70160975 minutes)
Iteration 122 (change = 0.0650858789942) (10.852542 seconds, I think I'll finish in 158.8088646 minutes)
Iteration 123 (change = 0.0634178315213) (11.182634 seconds, I think I'll finish in 163.452833633 minutes)
Iteration 124 (change = 0.0613108965662) (11.021292 seconds, I think I'll finish in 160.9108632 minutes)
Iteration 125 (change = 0.0587274027368) (10.472634 seconds, I think I'll finish in 152.7259125 minutes)
Iteration 126 (change = 0.0560955302129) (10.051652 seconds, I think I'll finish in 146.419064133 minutes)
Iteration 127 (change = 0.053676024258) (10.089024 seconds, I think I'll finish in 146.7952992 minutes)
Iteration 128 (change = 0.0516420409627) (9.372124 seconds, I think I'll finish in 136.208202133 minutes)
Iteration 129 (change = 0.0501578828294) (9.356147 seconds, I think I'll finish in 135.820067283 minutes)
Iteration 130 (change = 0.0494337643194) (9.307945 seconds, I think I'll finish in 134.9652025 minutes)
Iteration 131 (change = 0.0497611462086) (9.204735 seconds, I think I'll finish in 133.31524525 minutes)
Iteration 132 (change = 0.0509672005275) (9.424375 seconds, I think I'll finish in 136.339291667 minutes)
Iteration 133 (change = 0.0526468437896) (9.278066 seconds, I think I'll finish in 134.0680537 minutes)
Iteration 134 (change = 0.0545570879254) (9.160397 seconds, I think I'll finish in 132.215063367 minutes)
Iteration 135 (change = 0.0565878672336) (9.59934 seconds, I think I'll finish in 138.390485 minutes)
Iteration 136 (change = 0.0587138644516) (9.389468 seconds, I think I'll finish in 135.2083392 minutes)
Iteration 137 (change = 0.0608934753596) (9.253261 seconds, I think I'll finish in 133.092737383 minutes)
Iteration 138 (change = 0.0629407662269) (9.230015 seconds, I think I'll finish in 132.604548833 minutes)
Iteration 139 (change = 0.0647771658113) (8.856319 seconds, I think I'll finish in 127.08817765 minutes)
Iteration 140 (change = 0.0664184024573) (8.709448 seconds, I think I'll finish in 124.835421333 minutes)
Iteration 141 (change = 0.0678606600672) (8.954001 seconds, I think I'll finish in 128.19144765 minutes)
Iteration 142 (change = 0.0688337793828) (8.941816 seconds, I think I'll finish in 127.8679688 minutes)
Iteration 143 (change = 0.0687990819605) (8.981734 seconds, I think I'll finish in 128.289100633 minutes)
Iteration 144 (change = 0.0670604792455) (8.929174 seconds, I think I'll finish in 127.389549067 minutes)
Iteration 145 (change = 0.0638059025694) (9.000134 seconds, I think I'll finish in 128.2519095 minutes)
Iteration 146 (change = 0.0603676655691) (9.024708 seconds, I think I'll finish in 128.4516772 minutes)
Iteration 147 (change = 0.0573431557238) (8.882936 seconds, I think I'll finish in 126.285740133 minutes)
Iteration 148 (change = 0.0547868324098) (8.806784 seconds, I think I'll finish in 125.0563328 minutes)
Iteration 149 (change = 0.0525477854473) (8.83302 seconds, I think I'll finish in 125.281667 minutes)
Iteration 150 (change = 0.0503260548779) (8.736806 seconds, I think I'll finish in 123.771418333 minutes)
Iteration 151 (change = 0.0481243891055) (8.881476 seconds, I think I'll finish in 125.6728854 minutes)
Iteration 152 (change = 0.0459661746071) (8.408954 seconds, I think I'll finish in 118.846549867 minutes)
Iteration 153 (change = 0.0438310269431) (8.493588 seconds, I think I'll finish in 119.9011506 minutes)
Iteration 154 (change = 0.0419369914061) (8.328192 seconds, I think I'll finish in 117.4275072 minutes)
Iteration 155 (change = 0.0402888811987) (8.277849 seconds, I think I'll finish in 116.57970675 minutes)
Iteration 156 (change = 0.0388143171563) (8.374708 seconds, I think I'll finish in 117.804225867 minutes)
Iteration 157 (change = 0.0374157414553) (8.910047 seconds, I think I'll finish in 125.18616035 minutes)
Iteration 158 (change = 0.0361778868755) (8.292485 seconds, I think I'll finish in 116.371206167 minutes)
Iteration 159 (change = 0.0353088181275) (8.190958 seconds, I think I'll finish in 114.809927967 minutes)
Iteration 160 (change = 0.0349797884551) (8.409045 seconds, I think I'll finish in 117.72663 minutes)
Iteration 161 (change = 0.0353214205305) (8.16439 seconds, I think I'll finish in 114.165386833 minutes)
Iteration 162 (change = 0.0368256970305) (8.162301 seconds, I think I'll finish in 114.0001373 minutes)
Iteration 163 (change = 0.039723380469) (8.182328 seconds, I think I'll finish in 114.1434756 minutes)
Iteration 164 (change = 0.043811586478) (8.181812 seconds, I think I'll finish in 113.999913867 minutes)
Iteration 165 (change = 0.0477392129297) (8.265491 seconds, I think I'll finish in 115.028083083 minutes)
Iteration 166 (change = 0.0495065697885) (8.140461 seconds, I think I'll finish in 113.1524079 minutes)
Iteration 167 (change = 0.0479111675656) (8.138621 seconds, I think I'll finish in 112.991188217 minutes)
Iteration 168 (change = 0.0436814112284) (8.325559 seconds, I think I'll finish in 115.447751467 minutes)
Iteration 169 (change = 0.0385998341344) (8.180119 seconds, I think I'll finish in 113.29464815 minutes)
Iteration 170 (change = 0.0340805256095) (8.01012 seconds, I think I'll finish in 110.80666 minutes)
Iteration 171 (change = 0.0306925666187) (8.341125 seconds, I think I'll finish in 115.24654375 minutes)
Iteration 172 (change = 0.0284277051154) (8.589849 seconds, I think I'll finish in 118.5399162 minutes)
Iteration 173 (change = 0.0271025040929) (8.680925 seconds, I think I'll finish in 119.652082917 minutes)
Iteration 174 (change = 0.0264834423611) (8.724483 seconds, I think I'll finish in 120.1070493 minutes)
Iteration 175 (change = 0.0262702710285) (8.941139 seconds, I think I'll finish in 122.94066125 minutes)
Iteration 176 (change = 0.0261826041388) (8.669423 seconds, I think I'll finish in 119.060075867 minutes)
Iteration 177 (change = 0.0260798468978) (8.665156 seconds, I think I'll finish in 118.857056467 minutes)
Iteration 178 (change = 0.0259047095038) (8.57334 seconds, I think I'll finish in 117.454758 minutes)
Iteration 179 (change = 0.0257627939341) (8.569795 seconds, I think I'll finish in 117.263361583 minutes)
Iteration 180 (change = 0.0258628783146) (8.36355 seconds, I think I'll finish in 114.30185 minutes)
Iteration 181 (change = 0.0264978354111) (8.304935 seconds, I think I'll finish in 113.36236275 minutes)
Iteration 182 (change = 0.0279313653067) (8.23064 seconds, I think I'll finish in 112.211058667 minutes)
Iteration 183 (change = 0.0304596227364) (8.222462 seconds, I think I'll finish in 111.962524233 minutes)
Iteration 184 (change = 0.0341254265705) (7.97635 seconds, I think I'll finish in 108.47836 minutes)
Iteration 185 (change = 0.0380067801568) (7.966171 seconds, I think I'll finish in 108.207156083 minutes)
Iteration 186 (change = 0.0392652192073) (8.064847 seconds, I think I'll finish in 109.413090967 minutes)
Iteration 187 (change = 0.0342159255043) (8.198616 seconds, I think I'll finish in 111.0912468 minutes)
Iteration 188 (change = 0.0280177987219) (8.781405 seconds, I think I'll finish in 118.841681 minutes)
Iteration 189 (change = 0.0284723300849) (8.307494 seconds, I think I'll finish in 112.289627233 minutes)
Iteration 190 (change = 0.0304185924416) (8.536856 seconds, I think I'll finish in 115.247556 minutes)
Iteration 191 (change = 0.0291111524118) (8.217342 seconds, I think I'll finish in 110.7971613 minutes)
Iteration 192 (change = 0.0273787010223) (8.204364 seconds, I think I'll finish in 110.4854352 minutes)
Iteration 193 (change = 0.0264122712243) (8.283309 seconds, I think I'll finish in 111.41050605 minutes)
Iteration 194 (change = 0.0257128355395) (9.924893 seconds, I think I'll finish in 133.324395967 minutes)
Iteration 195 (change = 0.0247955923635) (9.619848 seconds, I think I'll finish in 129.066294 minutes)
Iteration 196 (change = 0.0236285596815) (8.103635 seconds, I think I'll finish in 108.588709 minutes)
Iteration 197 (change = 0.0225073575752) (8.042392 seconds, I think I'll finish in 107.634012933 minutes)
Iteration 198 (change = 0.0216028834261) (8.107381 seconds, I think I'll finish in 108.368659367 minutes)
Iteration 199 (change = 0.0209751392896) (8.01683 seconds, I think I'll finish in 107.0246805 minutes)
Iteration 200 (change = 0.0206644974448) (8.114702 seconds, I think I'll finish in 108.196026667 minutes)
Iteration 201 (change = 0.0207120793529) (8.07995 seconds, I think I'll finish in 107.598000833 minutes)
Iteration 202 (change = 0.021224865158) (7.97635 seconds, I think I'll finish in 106.085455 minutes)
Iteration 203 (change = 0.022247691235) (8.079443 seconds, I think I'll finish in 107.321934517 minutes)
Iteration 204 (change = 0.0240510076124) (8.147608 seconds, I think I'll finish in 108.091599467 minutes)
Iteration 205 (change = 0.0271395044129) (8.030238 seconds, I think I'll finish in 106.4006535 minutes)
Iteration 206 (change = 0.0325334927878) (7.994761 seconds, I think I'll finish in 105.797337233 minutes)
Iteration 207 (change = 0.0422911501053) (7.971553 seconds, I think I'll finish in 105.357358817 minutes)
Iteration 208 (change = 0.0592566835294) (8.178358 seconds, I think I'll finish in 107.9543256 minutes)
Iteration 209 (change = 0.080849705768) (7.995829 seconds, I think I'll finish in 105.411678983 minutes)
Iteration 210 (change = 0.0840941565198) (7.977234 seconds, I think I'll finish in 105.033581 minutes)
Iteration 211 (change = 0.0595420589569) (8.016425 seconds, I think I'll finish in 105.41598875 minutes)
Iteration 212 (change = 0.0282428781842) (8.027997 seconds, I think I'll finish in 105.4343606 minutes)
Iteration 213 (change = 0.0142403265714) (8.481913 seconds, I think I'll finish in 111.254425517 minutes)
Iteration 214 (change = 0.0130840590779) (8.406818 seconds, I think I'll finish in 110.1293158 minutes)
Iteration 215 (change = 0.0133338982804) (8.38962 seconds, I think I'll finish in 109.764195 minutes)
Iteration 216 (change = 0.0141708896539) (8.393455 seconds, I think I'll finish in 109.674478667 minutes)
Iteration 217 (change = 0.0156728450678) (8.281247 seconds, I think I'll finish in 108.07027335 minutes)
Iteration 218 (change = 0.0179125229065) (8.341585 seconds, I think I'll finish in 108.718657833 minutes)
Iteration 219 (change = 0.0213130888707) (8.449195 seconds, I think I'll finish in 109.980354917 minutes)
Iteration 220 (change = 0.0264459856197) (8.470195 seconds, I think I'll finish in 110.112535 minutes)
Iteration 221 (change = 0.0337688805953) (8.48131 seconds, I think I'll finish in 110.115674833 minutes)
Iteration 222 (change = 0.0419844058452) (8.389671 seconds, I think I'll finish in 108.7860673 minutes)
Iteration 223 (change = 0.0455204424799) (8.256053 seconds, I think I'll finish in 106.91588635 minutes)
Iteration 224 (change = 0.0409808882921) (8.290912 seconds, I think I'll finish in 107.229128533 minutes)
Iteration 225 (change = 0.043232605306) (8.217086 seconds, I think I'll finish in 106.137360833 minutes)
Iteration 226 (change = 0.0438643144951) (8.282977 seconds, I think I'll finish in 106.8504033 minutes)
Iteration 227 (change = 0.0165455911489) (8.377599 seconds, I think I'll finish in 107.93140045 minutes)
Iteration 228 (change = 0.00632479280955) (8.311256 seconds, I think I'll finish in 106.938160533 minutes)
Iteration 229 (change = 0.00611809958495) (8.284744 seconds, I think I'll finish in 106.4589604 minutes)
Iteration 230 (change = 0.00604880498911) (8.2185 seconds, I think I'll finish in 105.47075 minutes)
Iteration 231 (change = 0.00596345997136) (9.007741 seconds, I think I'll finish in 115.449213817 minutes)
Iteration 232 (change = 0.00593585833522) (8.34543 seconds, I think I'll finish in 106.821504 minutes)
Iteration 233 (change = 0.00592856554574) (8.260802 seconds, I think I'll finish in 105.600585567 minutes)
Iteration 234 (change = 0.00592281606852) (8.480084 seconds, I think I'll finish in 108.262405733 minutes)
Iteration 235 (change = 0.00591921143468) (8.227148 seconds, I think I'll finish in 104.896137 minutes)
Iteration 236 (change = 0.00591801904929) (8.299494 seconds, I think I'll finish in 105.6802236 minutes)
Iteration 237 (change = 0.0059185842402) (8.282039 seconds, I think I'll finish in 105.319929283 minutes)
Iteration 238 (change = 0.00590796968507) (8.300754 seconds, I think I'll finish in 105.4195758 minutes)
Iteration 239 (change = 0.00589925018165) (8.385884 seconds, I think I'll finish in 106.360962067 minutes)
Iteration 240 (change = 0.00589969156829) (8.140867 seconds, I think I'll finish in 103.117648667 minutes)
Iteration 241 (change = 0.00590210027947) (8.150482 seconds, I think I'll finish in 103.1035973 minutes)
Iteration 242 (change = 0.00593216105353) (8.221701 seconds, I think I'll finish in 103.8674893 minutes)
Iteration 243 (change = 0.00598956618617) (8.123863 seconds, I think I'll finish in 102.496071517 minutes)
Iteration 244 (change = 0.00605615741017) (8.27333 seconds, I think I'll finish in 104.243958 minutes)
Iteration 245 (change = 0.00614462093546) (8.087549 seconds, I think I'll finish in 101.768324917 minutes)
Iteration 246 (change = 0.00625417563657) (8.146064 seconds, I think I'll finish in 102.368870933 minutes)
Iteration 247 (change = 0.00637075069283) (8.099263 seconds, I think I'll finish in 101.64575065 minutes)
Iteration 248 (change = 0.00649093061782) (8.056954 seconds, I think I'll finish in 100.980490133 minutes)
Iteration 249 (change = 0.00660966985615) (8.194004 seconds, I think I'll finish in 102.561616733 minutes)
Iteration 250 (change = 0.00671225438999) (8.417356 seconds, I think I'll finish in 105.21695 minutes)
Iteration 251 (change = 0.00679824388935) (8.137043 seconds, I think I'll finish in 101.577420117 minutes)
Iteration 252 (change = 0.006899389191) (7.817282 seconds, I think I'll finish in 97.4554489333 minutes)
Iteration 253 (change = 0.00697050558819) (7.860239 seconds, I think I'll finish in 97.85997555 minutes)
Iteration 254 (change = 0.00703689852377) (7.744525 seconds, I think I'll finish in 96.2902608333 minutes)
Iteration 255 (change = 0.0070948715612) (7.875548 seconds, I think I'll finish in 97.7880543333 minutes)
Iteration 256 (change = 0.00713055375403) (7.628033 seconds, I think I'll finish in 94.5876092 minutes)
Iteration 257 (change = 0.00712725502317) (7.737874 seconds, I think I'll finish in 95.8206730333 minutes)
Iteration 258 (change = 0.00708844971016) (7.670888 seconds, I think I'll finish in 94.8633149333 minutes)
Iteration 259 (change = 0.00700065524436) (7.878668 seconds, I think I'll finish in 97.3015498 minutes)
Iteration 260 (change = 0.00688556458066) (7.654028 seconds, I think I'll finish in 94.3996786667 minutes)
Iteration 261 (change = 0.00677828930336) (7.729201 seconds, I think I'll finish in 95.1979923167 minutes)
Iteration 262 (change = 0.00666912284668) (7.678396 seconds, I think I'll finish in 94.4442708 minutes)
Iteration 263 (change = 0.00654876419373) (7.694953 seconds, I think I'll finish in 94.5196726833 minutes)
Iteration 264 (change = 0.00641849610622) (7.645033 seconds, I think I'll finish in 93.7790714667 minutes)
Iteration 265 (change = 0.00627607124508) (7.998849 seconds, I think I'll finish in 97.98590025 minutes)
Iteration 266 (change = 0.00610186496867) (7.846491 seconds, I think I'll finish in 95.9887399 minutes)
Iteration 267 (change = 0.00589094587594) (7.774831 seconds, I think I'll finish in 94.9825187167 minutes)
Iteration 268 (change = 0.00574222145802) (7.607608 seconds, I think I'll finish in 92.8128176 minutes)
Iteration 269 (change = 0.00562514784943) (7.738673 seconds, I think I'll finish in 94.2828327167 minutes)
Iteration 270 (change = 0.00548800824308) (8.285124 seconds, I think I'll finish in 100.802342 minutes)
Iteration 271 (change = 0.00533257803694) (8.407379 seconds, I think I'll finish in 102.14965485 minutes)
Iteration 272 (change = 0.00520743840732) (7.793164 seconds, I think I'll finish in 94.5570565333 minutes)
Iteration 273 (change = 0.0051426911195) (7.677272 seconds, I think I'll finish in 93.0229457333 minutes)
Iteration 274 (change = 0.00509422950841) (7.679768 seconds, I think I'll finish in 92.9251928 minutes)
Iteration 275 (change = 0.00505609786127) (7.607434 seconds, I think I'll finish in 91.9231608333 minutes)
Iteration 276 (change = 0.00502330126167) (7.63536 seconds, I think I'll finish in 92.133344 minutes)
Iteration 277 (change = 0.00498156587812) (7.641047 seconds, I think I'll finish in 92.07461635 minutes)
Iteration 278 (change = 0.0049505479858) (7.633905 seconds, I think I'll finish in 91.8613235 minutes)
Iteration 279 (change = 0.00496708682747) (7.666722 seconds, I think I'll finish in 92.1284427 minutes)
Iteration 280 (change = 0.00499450594662) (7.732897 seconds, I think I'll finish in 92.794764 minutes)
Iteration 281 (change = 0.00503062645446) (7.939909 seconds, I think I'll finish in 95.1465761833 minutes)
Iteration 282 (change = 0.00507638236712) (7.677858 seconds, I think I'll finish in 91.8783674 minutes)
Iteration 283 (change = 0.00513245511081) (7.696402 seconds, I think I'll finish in 91.9720039 minutes)
Iteration 284 (change = 0.00520015924312) (7.663021 seconds, I think I'll finish in 91.4453839333 minutes)
Iteration 285 (change = 0.00528139616341) (7.703436 seconds, I think I'll finish in 91.799279 minutes)
Iteration 286 (change = 0.00537813020271) (7.815679 seconds, I think I'll finish in 93.0065801 minutes)
Iteration 287 (change = 0.00549115025932) (7.841412 seconds, I think I'll finish in 93.1821126 minutes)
Iteration 288 (change = 0.00561250720326) (7.723032 seconds, I think I'll finish in 91.6466464 minutes)
Iteration 289 (change = 0.00574271359925) (7.826788 seconds, I think I'll finish in 92.7474378 minutes)
Iteration 290 (change = 0.00596344358518) (7.795174 seconds, I think I'll finish in 92.2428923333 minutes)
Iteration 291 (change = 0.00624417640947) (7.825721 seconds, I think I'll finish in 92.4739364833 minutes)
Iteration 292 (change = 0.00659525496707) (7.688284 seconds, I think I'll finish in 90.7217512 minutes)
Iteration 293 (change = 0.00703839219235) (7.738314 seconds, I think I'll finish in 91.1831333 minutes)
Iteration 294 (change = 0.00760252985876) (7.780108 seconds, I think I'll finish in 91.5459374667 minutes)
Iteration 295 (change = 0.00832291215446) (7.818009 seconds, I think I'll finish in 91.86160575 minutes)
Iteration 296 (change = 0.0092281694987) (7.869873 seconds, I think I'll finish in 92.3398432 minutes)
Iteration 297 (change = 0.0103900879477) (7.96139 seconds, I think I'll finish in 93.2809528333 minutes)
Iteration 298 (change = 0.0119037309158) (7.832468 seconds, I think I'll finish in 91.6398756 minutes)
Iteration 299 (change = 0.0137050894503) (7.837631 seconds, I think I'll finish in 91.5696555167 minutes)
Iteration 300 (change = 0.0156601647772) (7.740063 seconds, I think I'll finish in 90.300735 minutes)
Iteration 301 (change = 0.0174816714131) (7.73632 seconds, I think I'll finish in 90.128128 minutes)
Iteration 302 (change = 0.0187173180104) (7.845417 seconds, I think I'll finish in 91.2683511 minutes)
Iteration 303 (change = 0.0188763377603) (7.819161 seconds, I think I'll finish in 90.83258695 minutes)
Iteration 304 (change = 0.0177873782816) (7.750515 seconds, I think I'll finish in 89.905974 minutes)
Iteration 305 (change = 0.0159341748652) (7.869016 seconds, I think I'll finish in 91.1494353333 minutes)
Iteration 306 (change = 0.014387632779) (7.772115 seconds, I think I'll finish in 89.8974635 minutes)
Iteration 307 (change = 0.0140419696844) (7.818983 seconds, I think I'll finish in 90.30925365 minutes)
Iteration 308 (change = 0.0152508858929) (7.8465 seconds, I think I'll finish in 90.4963 minutes)
Iteration 309 (change = 0.0168959124952) (8.569735 seconds, I think I'll finish in 98.6947814167 minutes)
Iteration 310 (change = 0.0167269369761) (9.888475 seconds, I think I'll finish in 113.7174625 minutes)
Iteration 311 (change = 0.0147030664724) (7.917466 seconds, I think I'll finish in 90.9189012333 minutes)
Iteration 312 (change = 0.0126889512532) (8.065775 seconds, I think I'll finish in 92.4875533333 minutes)
Iteration 313 (change = 0.0111745089977) (7.893196 seconds, I think I'll finish in 90.3770942 minutes)
Iteration 314 (change = 0.00986611838877) (8.417223 seconds, I think I'll finish in 96.2369163 minutes)
Iteration 315 (change = 0.00857961609661) (7.81062 seconds, I think I'll finish in 89.171245 minutes)
Iteration 316 (change = 0.00730812665269) (7.69616 seconds, I think I'll finish in 87.736224 minutes)
Iteration 317 (change = 0.00618602627446) (8.35431 seconds, I think I'll finish in 95.0998955 minutes)
Iteration 318 (change = 0.00524612614646) (7.951555 seconds, I think I'll finish in 90.3826751667 minutes)
Iteration 319 (change = 0.0044978864665) (7.782954 seconds, I think I'll finish in 88.3365279 minutes)
Iteration 320 (change = 0.00393577384842) (7.810076 seconds, I think I'll finish in 88.5141946667 minutes)
Iteration 321 (change = 0.00353304208722) (7.738912 seconds, I think I'll finish in 87.5786874667 minutes)
Iteration 322 (change = 0.00325002742039) (7.799133 seconds, I think I'll finish in 88.1302029 minutes)
Iteration 323 (change = 0.00306148716765) (7.807757 seconds, I think I'll finish in 88.0975248167 minutes)
Iteration 324 (change = 0.00292542160591) (7.782296 seconds, I think I'll finish in 87.6805349333 minutes)
Iteration 325 (change = 0.00282471586692) (7.823173 seconds, I think I'll finish in 88.01069625 minutes)
Iteration 326 (change = 0.00282490454775) (7.868146 seconds, I think I'll finish in 88.3855067333 minutes)
Iteration 327 (change = 0.00287264494263) (7.982352 seconds, I think I'll finish in 89.5353816 minutes)
Iteration 328 (change = 0.00295087110787) (7.935129 seconds, I think I'll finish in 88.8734448 minutes)
Iteration 329 (change = 0.0030554002998) (7.811105 seconds, I think I'll finish in 87.3541909167 minutes)
Iteration 330 (change = 0.00318409420062) (7.829614 seconds, I think I'll finish in 87.4306896667 minutes)
Iteration 331 (change = 0.00333657083861) (7.892359 seconds, I think I'll finish in 87.99980285 minutes)
Iteration 332 (change = 0.00351262962365) (7.863574 seconds, I think I'll finish in 87.5477905333 minutes)
Iteration 333 (change = 0.00371165210629) (7.880554 seconds, I think I'll finish in 87.6054919667 minutes)
Iteration 334 (change = 0.00394184565095) (7.7884 seconds, I think I'll finish in 86.45124 minutes)
Iteration 335 (change = 0.00420872616739) (7.821142 seconds, I think I'll finish in 86.6843238333 minutes)
Iteration 336 (change = 0.00451276946504) (7.915805 seconds, I think I'll finish in 87.6015753333 minutes)
Iteration 337 (change = 0.00485895312462) (7.891456 seconds, I think I'll finish in 87.2005888 minutes)
Iteration 338 (change = 0.00525287795431) (7.783275 seconds, I think I'll finish in 85.8754675 minutes)
Iteration 339 (change = 0.00570086074297) (7.862578 seconds, I think I'll finish in 86.6194009667 minutes)
Iteration 340 (change = 0.00620963016955) (7.851423 seconds, I think I'll finish in 86.365653 minutes)
Iteration 341 (change = 0.00678461656039) (7.859015 seconds, I think I'll finish in 86.3181814167 minutes)
Iteration 342 (change = 0.00742950017982) (7.84971 seconds, I think I'll finish in 86.085153 minutes)
Iteration 343 (change = 0.00813921501542) (8.00278 seconds, I think I'll finish in 87.630441 minutes)
Iteration 344 (change = 0.00891170786275) (8.000647 seconds, I think I'll finish in 87.4737405333 minutes)
Iteration 345 (change = 0.00976466333193) (7.968973 seconds, I think I'll finish in 86.9946219167 minutes)
Iteration 346 (change = 0.0106477279737) (7.838611 seconds, I think I'll finish in 85.4408599 minutes)
Iteration 347 (change = 0.0115507609101) (7.980206 seconds, I think I'll finish in 86.8512419667 minutes)
Iteration 348 (change = 0.0124293598816) (7.817402 seconds, I think I'll finish in 84.9491017333 minutes)
Iteration 349 (change = 0.0132060990076) (7.983644 seconds, I think I'll finish in 86.6225374 minutes)
Iteration 350 (change = 0.0138041047732) (8.485968 seconds, I think I'll finish in 91.93132 minutes)
Iteration 351 (change = 0.0141981545354) (7.878722 seconds, I think I'll finish in 85.2215096333 minutes)
Iteration 352 (change = 0.0143893905944) (7.955528 seconds, I think I'll finish in 85.9197024 minutes)
Iteration 353 (change = 0.0143877814803) (8.003106 seconds, I think I'll finish in 86.3001597 minutes)
Iteration 354 (change = 0.0142930721622) (7.949785 seconds, I think I'll finish in 85.5926851667 minutes)
Iteration 355 (change = 0.0142746316719) (7.870867 seconds, I think I'll finish in 84.61182025 minutes)
Iteration 356 (change = 0.0145694383607) (7.983388 seconds, I think I'll finish in 85.6883645333 minutes)
Iteration 357 (change = 0.0154452826912) (7.916303 seconds, I think I'll finish in 84.8363804833 minutes)
Iteration 358 (change = 0.0171621082429) (7.792617 seconds, I think I'll finish in 83.3810019 minutes)
Iteration 359 (change = 0.0199059293134) (8.075133 seconds, I think I'll finish in 86.26933755 minutes)
Iteration 360 (change = 0.0236380503339) (7.90624 seconds, I think I'll finish in 84.3332266667 minutes)
Iteration 361 (change = 0.0281476212342) (7.887465 seconds, I think I'll finish in 84.00150225 minutes)
Iteration 362 (change = 0.0321270438223) (7.855421 seconds, I think I'll finish in 83.5293099667 minutes)
Iteration 363 (change = 0.034095039201) (7.862287 seconds, I think I'll finish in 83.4712803167 minutes)
Iteration 364 (change = 0.0345197739906) (8.439829 seconds, I think I'll finish in 89.4621874 minutes)
Iteration 365 (change = 0.034548623083) (7.841518 seconds, I think I'll finish in 82.9893988333 minutes)
Iteration 366 (change = 0.0342306680337) (7.746141 seconds, I think I'll finish in 81.8508899 minutes)
Iteration 367 (change = 0.0341861627552) (7.896984 seconds, I think I'll finish in 83.3131812 minutes)
Iteration 368 (change = 0.0345187656659) (7.816126 seconds, I think I'll finish in 82.3298605333 minutes)
Iteration 369 (change = 0.0350729385549) (7.917356 seconds, I think I'll finish in 83.2641939333 minutes)
Iteration 370 (change = 0.0356288288459) (7.871048 seconds, I think I'll finish in 82.646004 minutes)
Iteration 371 (change = 0.0363677060841) (7.936361 seconds, I think I'll finish in 83.1995178167 minutes)
Iteration 372 (change = 0.0366863390418) (7.809811 seconds, I think I'll finish in 81.7426884667 minutes)
Iteration 373 (change = 0.036294662025) (7.858819 seconds, I think I'll finish in 82.12465855 minutes)
Iteration 374 (change = 0.0352207720708) (8.016887 seconds, I think I'll finish in 83.6428543667 minutes)
Iteration 375 (change = 0.0336139911737) (7.928332 seconds, I think I'll finish in 82.5867916667 minutes)
Iteration 376 (change = 0.0315743681379) (7.765664 seconds, I think I'll finish in 80.7629056 minutes)
Iteration 377 (change = 0.0297275246329) (7.765116 seconds, I think I'll finish in 80.6277878 minutes)
Iteration 378 (change = 0.0279453935994) (7.849581 seconds, I think I'll finish in 81.3739897 minutes)
Iteration 379 (change = 0.0258517494995) (7.798854 seconds, I think I'll finish in 80.7181389 minutes)
Iteration 380 (change = 0.0234100039912) (7.876035 seconds, I think I'll finish in 81.385695 minutes)
Iteration 381 (change = 0.0210611071351) (7.781261 seconds, I think I'll finish in 80.2766759833 minutes)
Iteration 382 (change = 0.0189699650208) (7.738006 seconds, I think I'll finish in 79.7014618 minutes)
Iteration 383 (change = 0.0171919287549) (7.855806 seconds, I think I'll finish in 80.7838717 minutes)
Iteration 384 (change = 0.0157673597923) (7.708914 seconds, I think I'll finish in 79.1448504 minutes)
Iteration 385 (change = 0.0147280400632) (7.772313 seconds, I think I'll finish in 79.66620825 minutes)
Iteration 386 (change = 0.0141147150257) (7.980038 seconds, I think I'll finish in 81.6623888667 minutes)
Iteration 387 (change = 0.0140244440235) (7.942253 seconds, I think I'll finish in 81.1433514833 minutes)
Iteration 388 (change = 0.0148872311222) (7.723744 seconds, I think I'll finish in 78.7821888 minutes)
Iteration 389 (change = 0.0181686701475) (7.755049 seconds, I think I'll finish in 78.9722489833 minutes)
Iteration 390 (change = 0.0242059807546) (7.86992 seconds, I think I'll finish in 80.0108533333 minutes)
Iteration 391 (change = 0.0340001416263) (8.190984 seconds, I think I'll finish in 83.1384876 minutes)
Iteration 392 (change = 0.0419665609) (7.80125 seconds, I think I'll finish in 79.0526666667 minutes)
Iteration 393 (change = 0.0309517062367) (7.713142 seconds, I think I'll finish in 78.0312865667 minutes)
Iteration 394 (change = 0.0154477914227) (7.744417 seconds, I think I'll finish in 78.2186117 minutes)
Iteration 395 (change = 0.0147314888648) (7.736145 seconds, I think I'll finish in 78.00612875 minutes)
Iteration 396 (change = 0.0148693904793) (7.662183 seconds, I think I'll finish in 77.1326422 minutes)
Iteration 397 (change = 0.0148427488542) (7.8715 seconds, I think I'll finish in 79.108575 minutes)
Iteration 398 (change = 0.0146036937076) (7.736144 seconds, I think I'll finish in 77.6193114667 minutes)
Iteration 399 (change = 0.0142330047357) (7.934486 seconds, I think I'll finish in 79.4771014333 minutes)
Iteration 400 (change = 0.0138352773206) (7.810089 seconds, I think I'll finish in 78.10089 minutes)
Iteration 401 (change = 0.0134618384442) (7.904717 seconds, I think I'll finish in 78.9154247167 minutes)
Iteration 402 (change = 0.0130702814944) (7.713279 seconds, I think I'll finish in 76.8756807 minutes)
Iteration 403 (change = 0.0126042205073) (7.887036 seconds, I think I'll finish in 78.4760082 minutes)
Iteration 404 (change = 0.0119155697532) (7.928655 seconds, I think I'll finish in 78.757973 minutes)
Iteration 405 (change = 0.0109983845271) (7.688348 seconds, I think I'll finish in 76.2427843333 minutes)
Iteration 406 (change = 0.00990794605675) (8.004518 seconds, I think I'll finish in 79.2447282 minutes)
Iteration 407 (change = 0.00873961018397) (7.952507 seconds, I think I'll finish in 78.5972775167 minutes)
Iteration 408 (change = 0.00759551339896) (7.807386 seconds, I think I'll finish in 77.0328752 minutes)
Iteration 409 (change = 0.00655720551327) (7.833016 seconds, I think I'll finish in 77.1552076 minutes)
Iteration 410 (change = 0.00567594141035) (7.767025 seconds, I think I'll finish in 76.3757458333 minutes)
Iteration 411 (change = 0.00496612612997) (8.485212 seconds, I think I'll finish in 83.2964978 minutes)
Iteration 412 (change = 0.00440903282313) (7.798444 seconds, I think I'll finish in 76.4247512 minutes)
Iteration 413 (change = 0.00397230942276) (7.725916 seconds, I think I'll finish in 75.5852115333 minutes)
Iteration 414 (change = 0.00362620169657) (7.889937 seconds, I think I'll finish in 77.0583847 minutes)
Iteration 415 (change = 0.00334554523038) (7.933755 seconds, I think I'll finish in 77.35411125 minutes)
Iteration 416 (change = 0.00310920510067) (7.77886 seconds, I think I'll finish in 75.7142373333 minutes)
Iteration 417 (change = 0.00290115174361) (7.795076 seconds, I think I'll finish in 75.7421551333 minutes)
Iteration 418 (change = 0.0027118489673) (7.715507 seconds, I think I'll finish in 74.8404179 minutes)
Iteration 419 (change = 0.00253496597644) (7.724484 seconds, I think I'll finish in 74.7987534 minutes)
Iteration 420 (change = 0.00236085660151) (7.660545 seconds, I think I'll finish in 74.051935 minutes)
Iteration 421 (change = 0.0021848343649) (7.723193 seconds, I think I'll finish in 74.52881245 minutes)
Iteration 422 (change = 0.0020574734977) (8.043451 seconds, I think I'll finish in 77.4852446333 minutes)
Iteration 423 (change = 0.00194831236647) (8.110867 seconds, I think I'll finish in 77.9995043167 minutes)
Iteration 424 (change = 0.00184672421589) (7.918528 seconds, I think I'll finish in 76.0178688 minutes)
Iteration 425 (change = 0.00174259779268) (7.869691 seconds, I think I'll finish in 75.4178720833 minutes)
Iteration 426 (change = 0.00159268336954) (8.085069 seconds, I think I'll finish in 77.3471601 minutes)
Iteration 427 (change = 0.00144858891625) (10.927125 seconds, I think I'll finish in 104.35404375 minutes)
Iteration 428 (change = 0.0014159438165) (8.646899 seconds, I think I'll finish in 82.4337704667 minutes)
Iteration 429 (change = 0.00139715259597) (8.072179 seconds, I think I'll finish in 76.8202368167 minutes)
Iteration 430 (change = 0.00137418344795) (7.940795 seconds, I think I'll finish in 75.4375525 minutes)
Iteration 431 (change = 0.00136499868079) (7.933381 seconds, I think I'll finish in 75.2348964833 minutes)
Iteration 432 (change = 0.00136203390203) (7.898034 seconds, I think I'll finish in 74.7680552 minutes)
Iteration 433 (change = 0.00136172857573) (7.895138 seconds, I think I'll finish in 74.6090541 minutes)
Iteration 434 (change = 0.00136303887905) (7.815069 seconds, I think I'll finish in 73.7221509 minutes)
Iteration 435 (change = 0.00136551306282) (7.864428 seconds, I think I'll finish in 74.056697 minutes)
Iteration 436 (change = 0.00136897252464) (7.885797 seconds, I think I'll finish in 74.1264918 minutes)
Iteration 437 (change = 0.00137337358143) (8.248349 seconds, I think I'll finish in 77.3970081167 minutes)
Iteration 438 (change = 0.0013787286542) (7.874182 seconds, I think I'll finish in 73.7548380667 minutes)
Iteration 439 (change = 0.00138504615611) (7.908825 seconds, I think I'll finish in 73.94751375 minutes)
Iteration 440 (change = 0.0013924755076) (7.782273 seconds, I think I'll finish in 72.634548 minutes)
Iteration 441 (change = 0.00140120633155) (7.809191 seconds, I think I'll finish in 72.7556294833 minutes)
Iteration 442 (change = 0.00141129485309) (7.831189 seconds, I think I'll finish in 72.8300577 minutes)
Iteration 443 (change = 0.00142326901198) (7.896361 seconds, I think I'll finish in 73.3045512833 minutes)
Iteration 444 (change = 0.00143732463378) (7.888546 seconds, I think I'll finish in 73.1005262667 minutes)
Iteration 445 (change = 0.00145332427475) (7.907324 seconds, I think I'll finish in 73.142747 minutes)
Iteration 446 (change = 0.00147128714602) (7.754038 seconds, I think I'll finish in 71.5956175333 minutes)
Iteration 447 (change = 0.00149144470288) (7.800019 seconds, I think I'll finish in 71.8901751167 minutes)
Iteration 448 (change = 0.00151407960468) (7.746991 seconds, I think I'll finish in 71.2723172 minutes)
Iteration 449 (change = 0.00153952879458) (7.938977 seconds, I think I'll finish in 72.9062721167 minutes)
Iteration 450 (change = 0.0015681912009) (7.849358 seconds, I think I'll finish in 71.9524483333 minutes)
Iteration 451 (change = 0.00160047308145) (7.970934 seconds, I think I'll finish in 72.9340461 minutes)
Iteration 452 (change = 0.0016368091781) (8.119606 seconds, I think I'll finish in 74.1590681333 minutes)
Iteration 453 (change = 0.00167779495397) (8.183057 seconds, I think I'll finish in 74.6022029833 minutes)
Iteration 454 (change = 0.00172406472547) (7.960929 seconds, I think I'll finish in 72.4444539 minutes)
Iteration 455 (change = 0.00177635770671) (8.012344 seconds, I think I'll finish in 72.7787913333 minutes)
Iteration 456 (change = 0.00183556441737) (7.995467 seconds, I think I'll finish in 72.4922341333 minutes)
Iteration 457 (change = 0.00190276765398) (8.51574 seconds, I think I'll finish in 77.067447 minutes)
Iteration 458 (change = 0.00197908762731) (7.95843 seconds, I think I'll finish in 71.891151 minutes)
Iteration 459 (change = 0.0020659692297) (7.97147 seconds, I think I'll finish in 71.8760878333 minutes)
Iteration 460 (change = 0.00216509233778) (7.931403 seconds, I think I'll finish in 71.382627 minutes)
Iteration 461 (change = 0.00227845597) (8.0839 seconds, I think I'll finish in 72.6203683333 minutes)
Iteration 462 (change = 0.00240848375781) (7.929883 seconds, I think I'll finish in 71.1046175667 minutes)
Iteration 463 (change = 0.00255814028766) (7.963444 seconds, I think I'll finish in 71.2728238 minutes)
Iteration 464 (change = 0.0027310457668) (7.975523 seconds, I think I'll finish in 71.2480054667 minutes)
Iteration 465 (change = 0.00293164031945) (7.975728 seconds, I think I'll finish in 71.116908 minutes)
Iteration 466 (change = 0.00316529854763) (8.287932 seconds, I think I'll finish in 73.7625948 minutes)
Iteration 467 (change = 0.00343799684796) (8.116016 seconds, I think I'll finish in 72.0972754667 minutes)
Iteration 468 (change = 0.00375341961537) (8.095387 seconds, I think I'll finish in 71.7790980667 minutes)
Iteration 469 (change = 0.00411507449951) (8.11232 seconds, I think I'll finish in 71.794032 minutes)
Iteration 470 (change = 0.00458067638809) (7.990865 seconds, I think I'll finish in 70.5859741667 minutes)
Iteration 471 (change = 0.00516224474971) (8.03588 seconds, I think I'll finish in 70.8496753333 minutes)
Iteration 472 (change = 0.00588841720144) (7.953257 seconds, I think I'll finish in 69.9886616 minutes)
Iteration 473 (change = 0.00681570659961) (8.007067 seconds, I think I'll finish in 70.3287384833 minutes)
Iteration 474 (change = 0.00803527551588) (7.827694 seconds, I think I'll finish in 68.6227840667 minutes)
Iteration 475 (change = 0.00970153401537) (8.197847 seconds, I think I'll finish in 71.73116125 minutes)
Iteration 476 (change = 0.012088098989) (8.06672 seconds, I think I'll finish in 70.4493546667 minutes)
Iteration 477 (change = 0.0157137657783) (8.028345 seconds, I think I'll finish in 69.98040725 minutes)
Iteration 478 (change = 0.0217135730237) (8.022831 seconds, I think I'll finish in 69.7986297 minutes)
Iteration 479 (change = 0.0322614832274) (8.064784 seconds, I think I'll finish in 70.0292077333 minutes)
Iteration 480 (change = 0.050759648305) (7.910249 seconds, I think I'll finish in 68.5554913333 minutes)
Iteration 481 (change = 0.0768455327199) (8.01953 seconds, I think I'll finish in 69.3689345 minutes)
Iteration 482 (change = 0.082956424635) (7.941255 seconds, I think I'll finish in 68.5595015 minutes)
Iteration 483 (change = 0.0483173893085) (8.086156 seconds, I think I'll finish in 69.6757108667 minutes)
Iteration 484 (change = 0.0281186025059) (8.156525 seconds, I think I'll finish in 70.146115 minutes)
Iteration 485 (change = 0.0216716022885) (7.983602 seconds, I think I'll finish in 68.5259171667 minutes)
Iteration 486 (change = 0.0163997662173) (7.924149 seconds, I think I'll finish in 67.8835431 minutes)
Iteration 487 (change = 0.0169026213765) (7.96565 seconds, I think I'll finish in 68.1063075 minutes)
Iteration 488 (change = 0.0181870006092) (7.869807 seconds, I think I'll finish in 67.1556864 minutes)
Iteration 489 (change = 0.0180630378485) (8.012465 seconds, I think I'll finish in 68.2394935833 minutes)
Iteration 490 (change = 0.0154090753093) (8.047748 seconds, I think I'll finish in 68.405858 minutes)
Iteration 491 (change = 0.012228588912) (7.963046 seconds, I think I'll finish in 67.5531735667 minutes)
Iteration 492 (change = 0.00963824337059) (8.063031 seconds, I think I'll finish in 68.2669958 minutes)
Iteration 493 (change = 0.00797997561882) (7.968008 seconds, I think I'll finish in 67.3296676 minutes)
Iteration 494 (change = 0.00708355727119) (7.984611 seconds, I think I'll finish in 67.3368861 minutes)
Iteration 495 (change = 0.00651608024565) (7.945865 seconds, I think I'll finish in 66.8776970833 minutes)
Iteration 496 (change = 0.00614205817907) (8.158572 seconds, I think I'll finish in 68.5320048 minutes)
Iteration 497 (change = 0.00587153805014) (8.492781 seconds, I think I'll finish in 71.19781405 minutes)
Iteration 498 (change = 0.0056397373811) (7.98548 seconds, I think I'll finish in 66.8118493333 minutes)
Iteration 499 (change = 0.00540223263822) (8.417087 seconds, I think I'll finish in 70.28267645 minutes)
Iteration 500 (change = 0.00515108478834) (8.152088 seconds, I think I'll finish in 67.9340666667 minutes)
Iteration 501 (change = 0.00490296606239) (8.173328 seconds, I think I'll finish in 67.9748445333 minutes)
Iteration 502 (change = 0.00463317019067) (7.977276 seconds, I think I'll finish in 66.2113908 minutes)
Iteration 503 (change = 0.00434389097733) (8.315036 seconds, I think I'll finish in 68.8762148667 minutes)
Iteration 504 (change = 0.00404570832989) (8.687808 seconds, I think I'll finish in 71.8192128 minutes)
Iteration 505 (change = 0.0037456072915) (8.020725 seconds, I think I'll finish in 66.17098125 minutes)
Iteration 506 (change = 0.00344469387014) (8.641881 seconds, I think I'll finish in 71.1514869 minutes)
Iteration 507 (change = 0.00316453504697) (8.170956 seconds, I think I'll finish in 67.1380218 minutes)
Iteration 508 (change = 0.00291733899899) (7.936742 seconds, I think I'll finish in 65.0812844 minutes)
Iteration 509 (change = 0.00269606663043) (8.088065 seconds, I think I'll finish in 66.1873319167 minutes)
Iteration 510 (change = 0.00249820096452) (7.933858 seconds, I think I'll finish in 64.7931736667 minutes)
Iteration 511 (change = 0.00232280135786) (8.111862 seconds, I think I'll finish in 66.1116753 minutes)
Iteration 512 (change = 0.00216736715879) (8.015586 seconds, I think I'll finish in 65.1934328 minutes)
Iteration 513 (change = 0.00202925568731) (8.062805 seconds, I think I'll finish in 65.4431005833 minutes)
Iteration 514 (change = 0.00190629553053) (8.16178 seconds, I think I'll finish in 66.110418 minutes)
Iteration 515 (change = 0.00179607393617) (8.120735 seconds, I think I'll finish in 65.6426079167 minutes)
Iteration 516 (change = 0.0016963229468) (8.102064 seconds, I think I'll finish in 65.3566496 minutes)
Iteration 517 (change = 0.00160528994773) (8.101007 seconds, I think I'll finish in 65.21310635 minutes)
Iteration 518 (change = 0.0015215478568) (7.980245 seconds, I think I'll finish in 64.1079681667 minutes)
Iteration 519 (change = 0.00144397696717) (8.068191 seconds, I think I'll finish in 64.67999785 minutes)
Iteration 520 (change = 0.00137170303217) (7.944941 seconds, I think I'll finish in 63.559528 minutes)
Iteration 521 (change = 0.00130407467938) (8.031381 seconds, I think I'll finish in 64.11719165 minutes)
Iteration 522 (change = 0.00124058669808) (8.001739 seconds, I think I'll finish in 63.7471873667 minutes)
Iteration 523 (change = 0.00118084082548) (8.173451 seconds, I think I'll finish in 64.97893545 minutes)
Iteration 524 (change = 0.0011244211669) (8.001566 seconds, I think I'll finish in 63.4790902667 minutes)
Iteration 525 (change = 0.0010707729968) (8.019252 seconds, I think I'll finish in 63.485745 minutes)
Iteration 526 (change = 0.00101859975588) (7.939205 seconds, I think I'll finish in 62.7197195 minutes)
Iteration 527 (change = 0.000963461899545) (8.074375 seconds, I think I'll finish in 63.6529895833 minutes)
Iteration 528 (change = 0.000891237237277) (7.993208 seconds, I think I'll finish in 62.8799029333 minutes)
Iteration 529 (change = 0.00083036315651) (8.018958 seconds, I think I'll finish in 62.9488203 minutes)
Iteration 530 (change = 0.000793791138481) (8.159148 seconds, I think I'll finish in 63.913326 minutes)
Iteration 531 (change = 0.000761238469521) (8.054171 seconds, I think I'll finish in 62.9567699833 minutes)
Iteration 532 (change = 0.000725583607819) (7.890352 seconds, I think I'll finish in 61.5447456 minutes)
Iteration 533 (change = 0.000685661012739) (7.969749 seconds, I think I'll finish in 62.03121305 minutes)
Iteration 534 (change = 0.000653748944033) (7.938209 seconds, I think I'll finish in 61.6534232333 minutes)
Iteration 535 (change = 0.000625529550425) (8.012179 seconds, I think I'll finish in 62.09438725 minutes)
Iteration 536 (change = 0.000599876577947) (7.974136 seconds, I think I'll finish in 61.6666517333 minutes)
Iteration 537 (change = 0.000576280589569) (8.006573 seconds, I think I'll finish in 61.7840549833 minutes)
Iteration 538 (change = 0.000554526774728) (7.888116 seconds, I think I'll finish in 60.7384932 minutes)
Iteration 539 (change = 0.000534372091572) (8.058035 seconds, I think I'll finish in 61.9125689167 minutes)
Iteration 540 (change = 0.000515694060129) (8.033178 seconds, I think I'll finish in 61.587698 minutes)
Iteration 541 (change = 0.000498395314334) (7.922537 seconds, I think I'll finish in 60.60740805 minutes)
Iteration 542 (change = 0.000482335580925) (8.14556 seconds, I think I'll finish in 62.1777746667 minutes)
Iteration 543 (change = 0.000467706229699) (9.971799 seconds, I think I'll finish in 75.95186905 minutes)
Iteration 544 (change = 0.000454156376163) (8.062633 seconds, I think I'll finish in 61.2760108 minutes)
Iteration 545 (change = 0.00044156094523) (8.263368 seconds, I think I'll finish in 62.663874 minutes)
Iteration 546 (change = 0.000429844301484) (8.185862 seconds, I think I'll finish in 61.9396891333 minutes)
Iteration 547 (change = 0.000418942956821) (8.074301 seconds, I think I'll finish in 60.96097255 minutes)
Iteration 548 (change = 0.000408802607367) (8.453322 seconds, I think I'll finish in 63.6816924 minutes)
Iteration 549 (change = 0.000399374139373) (8.00919 seconds, I think I'll finish in 60.2024115 minutes)
Iteration 550 (change = 0.000390637404593) (8.555347 seconds, I think I'll finish in 64.1651025 minutes)
Iteration 551 (change = 0.000382534356321) (8.07714 seconds, I think I'll finish in 60.443931 minutes)
Iteration 552 (change = 0.000375022497374) (7.900283 seconds, I think I'll finish in 58.9887797333 minutes)
Iteration 553 (change = 0.000368062252122) (8.031949 seconds, I think I'll finish in 59.83802005 minutes)
Iteration 554 (change = 0.00036161557893) (7.998563 seconds, I think I'll finish in 59.4559849667 minutes)
Iteration 555 (change = 0.000355652075663) (8.034823 seconds, I think I'll finish in 59.5916039167 minutes)
Iteration 556 (change = 0.000350155240089) (7.945032 seconds, I think I'll finish in 58.7932368 minutes)
Iteration 557 (change = 0.000345107044686) (7.973617 seconds, I think I'll finish in 58.8718721833 minutes)
Iteration 558 (change = 0.000340476581767) (7.925181 seconds, I think I'll finish in 58.3821667 minutes)
Iteration 559 (change = 0.000336242337217) (7.975554 seconds, I think I'll finish in 58.6203219 minutes)
Iteration 560 (change = 0.000332410222209) (8.087001 seconds, I think I'll finish in 59.304674 minutes)
Iteration 561 (change = 0.000328937348714) (8.061362 seconds, I think I'll finish in 58.9822986333 minutes)
Iteration 562 (change = 0.000325820267092) (8.04466 seconds, I think I'll finish in 58.726018 minutes)
Iteration 563 (change = 0.000323093558838) (8.063882 seconds, I think I'll finish in 58.7319405667 minutes)
Iteration 564 (change = 0.000320674927178) (7.93594 seconds, I think I'll finish in 57.6678306667 minutes)
Iteration 565 (change = 0.000318556064105) (8.094425 seconds, I think I'll finish in 58.68458125 minutes)
Iteration 566 (change = 0.000316731303494) (8.064322 seconds, I think I'll finish in 58.3319291333 minutes)
Iteration 567 (change = 0.00031520137003) (7.94812 seconds, I think I'll finish in 57.3589326667 minutes)
Iteration 568 (change = 0.000313961082154) (7.991616 seconds, I think I'll finish in 57.5396352 minutes)
Iteration 569 (change = 0.00031299583733) (8.028991 seconds, I think I'll finish in 57.6749186833 minutes)
Iteration 570 (change = 0.000312305641921) (7.969892 seconds, I think I'll finish in 57.1175593333 minutes)
Iteration 571 (change = 0.000311892994513) (8.022983 seconds, I think I'll finish in 57.36432845 minutes)
Iteration 572 (change = 0.000311788891209) (7.876079 seconds, I think I'll finish in 56.1826968667 minutes)
Iteration 573 (change = 0.00031196176082) (7.99406 seconds, I think I'll finish in 56.8910603333 minutes)
Iteration 574 (change = 0.000312415143712) (8.025907 seconds, I think I'll finish in 56.9839397 minutes)
Iteration 575 (change = 0.000313224436329) (8.058059 seconds, I think I'll finish in 57.0779179167 minutes)
Iteration 576 (change = 0.00031431736378) (8.279065 seconds, I think I'll finish in 58.5053926667 minutes)
Iteration 577 (change = 0.00031569936788) (8.18756 seconds, I think I'll finish in 57.722298 minutes)
Iteration 578 (change = 0.000317370519821) (8.096948 seconds, I think I'll finish in 56.9485342667 minutes)
Iteration 579 (change = 0.000319348199991) (8.059836 seconds, I think I'll finish in 56.5531826 minutes)
Iteration 580 (change = 0.000321635370958) (8.026464 seconds, I think I'll finish in 56.185248 minutes)
Iteration 581 (change = 0.000324244880979) (8.024599 seconds, I think I'll finish in 56.0384496833 minutes)
Iteration 582 (change = 0.000327189370193) (8.026705 seconds, I think I'll finish in 55.9193781667 minutes)
Iteration 583 (change = 0.000330483949787) (8.381802 seconds, I think I'll finish in 58.2535239 minutes)
Iteration 584 (change = 0.00033414727093) (8.185736 seconds, I think I'll finish in 56.7544362667 minutes)
Iteration 585 (change = 0.000338201233852) (8.102634 seconds, I think I'll finish in 56.0432185 minutes)
Iteration 586 (change = 0.000342669454788) (8.047381 seconds, I think I'll finish in 55.5269289 minutes)
Iteration 587 (change = 0.0003475812565) (8.011192 seconds, I think I'll finish in 55.1437049333 minutes)
Iteration 588 (change = 0.000352964650912) (7.950235 seconds, I think I'll finish in 54.5916136667 minutes)
Iteration 589 (change = 0.000358852029974) (8.006648 seconds, I think I'll finish in 54.8455388 minutes)
Iteration 590 (change = 0.00036528057261) (7.892068 seconds, I think I'll finish in 53.9291313333 minutes)
Iteration 591 (change = 0.000372293089307) (8.21321 seconds, I think I'll finish in 55.9867148333 minutes)
Iteration 592 (change = 0.000379939149817) (8.131025 seconds, I think I'll finish in 55.29097 minutes)
Iteration 593 (change = 0.0003882767062) (8.105806 seconds, I think I'll finish in 54.9843840333 minutes)
Iteration 594 (change = 0.000397384124602) (8.107894 seconds, I think I'll finish in 54.8634160667 minutes)
Iteration 595 (change = 0.000407328216711) (8.07531 seconds, I think I'll finish in 54.5083425 minutes)
Iteration 596 (change = 0.000418119669187) (8.696074 seconds, I think I'll finish in 58.5535649333 minutes)
Iteration 597 (change = 0.000429906432253) (8.105641 seconds, I think I'll finish in 54.4428887167 minutes)
Iteration 598 (change = 0.000442820820793) (7.983102 seconds, I think I'll finish in 53.4867834 minutes)
Iteration 599 (change = 0.000456945991319) (8.151048 seconds, I think I'll finish in 54.4761708 minutes)
Iteration 600 (change = 0.000472425388618) (7.962168 seconds, I think I'll finish in 53.08112 minutes)
Iteration 601 (change = 0.000489405995345) (8.004302 seconds, I think I'll finish in 53.2286083 minutes)
Iteration 602 (change = 0.000508069407646) (7.922091 seconds, I think I'll finish in 52.5498703 minutes)
Iteration 603 (change = 0.000528611508776) (7.974709 seconds, I think I'll finish in 52.7659912167 minutes)
Iteration 604 (change = 0.000551273564713) (8.007779 seconds, I think I'll finish in 52.8513414 minutes)
Iteration 605 (change = 0.000576340705076) (7.995809 seconds, I think I'll finish in 52.6390759167 minutes)
Iteration 606 (change = 0.000604138574459) (8.003886 seconds, I think I'll finish in 52.5588514 minutes)
Iteration 607 (change = 0.000635019523676) (8.217133 seconds, I think I'll finish in 53.82222115 minutes)
Iteration 608 (change = 0.000669455310023) (8.094722 seconds, I think I'll finish in 52.8855170667 minutes)
Iteration 609 (change = 0.000708003303172) (8.029017 seconds, I think I'll finish in 52.32242745 minutes)
Iteration 610 (change = 0.000751307528978) (7.94924 seconds, I think I'll finish in 51.67006 minutes)
Iteration 611 (change = 0.000800148696618) (7.999561 seconds, I think I'll finish in 51.8638204833 minutes)
Iteration 612 (change = 0.000855481468278) (7.937408 seconds, I think I'll finish in 51.3285717333 minutes)
Iteration 613 (change = 0.000918462107371) (8.05519 seconds, I think I'll finish in 51.9559755 minutes)
Iteration 614 (change = 0.000990508132701) (7.978619 seconds, I think I'll finish in 51.3291155667 minutes)
Iteration 615 (change = 0.00107337988038) (8.072905 seconds, I think I'll finish in 51.8011404167 minutes)
Iteration 616 (change = 0.00116927571503) (8.135859 seconds, I think I'll finish in 52.0694976 minutes)
Iteration 617 (change = 0.00128096202104) (8.009016 seconds, I think I'll finish in 51.1242188 minutes)
Iteration 618 (change = 0.00141194714097) (8.262886 seconds, I think I'll finish in 52.6070408667 minutes)
Iteration 619 (change = 0.00156671856253) (8.107025 seconds, I think I'll finish in 51.47960875 minutes)
Iteration 620 (change = 0.00175105046376) (7.843322 seconds, I think I'll finish in 49.6743726667 minutes)
Iteration 621 (change = 0.00197235148831) (7.99945 seconds, I think I'll finish in 50.5298591667 minutes)
Iteration 622 (change = 0.00223982829627) (8.001106 seconds, I think I'll finish in 50.4069678 minutes)
Iteration 623 (change = 0.00256297999813) (8.135555 seconds, I think I'll finish in 51.1184039167 minutes)
Iteration 624 (change = 0.00294853036771) (7.918212 seconds, I think I'll finish in 49.6207952 minutes)
Iteration 625 (change = 0.00344883095112) (7.975594 seconds, I think I'll finish in 49.8474625 minutes)
Iteration 626 (change = 0.00410816078008) (7.91178 seconds, I think I'll finish in 49.316762 minutes)
Iteration 627 (change = 0.00497141333716) (8.008876 seconds, I think I'll finish in 49.7885124667 minutes)
Iteration 628 (change = 0.00611895898095) (7.960181 seconds, I think I'll finish in 49.3531222 minutes)
Iteration 629 (change = 0.00768286952539) (8.014265 seconds, I think I'll finish in 49.5548719167 minutes)
Iteration 630 (change = 0.00986405975361) (7.944128 seconds, I think I'll finish in 48.9887893333 minutes)
Iteration 631 (change = 0.0129801088396) (8.133275 seconds, I think I'll finish in 50.01964125 minutes)
Iteration 632 (change = 0.0175746368387) (7.925601 seconds, I think I'll finish in 48.6103528 minutes)
Iteration 633 (change = 0.0246179049849) (7.899039 seconds, I think I'll finish in 48.31578855 minutes)
Iteration 634 (change = 0.035866361362) (7.922599 seconds, I think I'll finish in 48.3278539 minutes)
Iteration 635 (change = 0.0532600042278) (8.022662 seconds, I think I'll finish in 48.8045271667 minutes)
Iteration 636 (change = 0.0673724167175) (7.820912 seconds, I think I'll finish in 47.4468661333 minutes)
Iteration 637 (change = 0.0420538150905) (7.9861 seconds, I think I'll finish in 48.315905 minutes)
Iteration 638 (change = 0.00641859655981) (8.110482 seconds, I think I'll finish in 48.9332414 minutes)
Iteration 639 (change = 0.000512221921047) (8.038354 seconds, I think I'll finish in 48.3640965667 minutes)
Iteration 640 (change = 0.000418933618599) (8.014171 seconds, I think I'll finish in 48.085026 minutes)
Iteration 641 (change = 0.000374770377384) (8.020469 seconds, I think I'll finish in 47.9891395167 minutes)
Iteration 642 (change = 0.000341527655249) (8.159839 seconds, I think I'll finish in 48.6870393667 minutes)
Iteration 643 (change = 0.000315852361022) (8.000037 seconds, I think I'll finish in 47.60022015 minutes)
Iteration 644 (change = 0.00029515473229) (7.934287 seconds, I think I'll finish in 47.0767695333 minutes)
Iteration 645 (change = 0.000278560330111) (7.868031 seconds, I think I'll finish in 46.55251675 minutes)
Iteration 646 (change = 0.000265230015015) (7.912899 seconds, I think I'll finish in 46.6861041 minutes)
Iteration 647 (change = 0.000254234846706) (8.011198 seconds, I think I'll finish in 47.1325482333 minutes)
Iteration 648 (change = 0.000245071497911) (7.965128 seconds, I think I'll finish in 46.7287509333 minutes)
Iteration 649 (change = 0.00023737876156) (7.914263 seconds, I think I'll finish in 46.29843855 minutes)
Iteration 650 (change = 0.000230895380271) (7.884028 seconds, I think I'll finish in 45.9901633333 minutes)
Iteration 651 (change = 0.00022543574572) (7.953472 seconds, I think I'll finish in 46.2626954667 minutes)
Iteration 652 (change = 0.000220827997628) (7.904672 seconds, I think I'll finish in 45.8470976 minutes)
Iteration 653 (change = 0.00021696427813) (7.966796 seconds, I think I'll finish in 46.0746368667 minutes)
Iteration 654 (change = 0.00021375314689) (8.15512 seconds, I think I'll finish in 47.0278586667 minutes)
Iteration 655 (change = 0.000211113904475) (8.225386 seconds, I think I'll finish in 47.2959695 minutes)
Iteration 656 (change = 0.000208981356472) (7.992716 seconds, I think I'll finish in 45.8249050667 minutes)
Iteration 657 (change = 0.000207304851533) (8.031428 seconds, I think I'll finish in 45.9129967333 minutes)
Iteration 658 (change = 0.000206047177906) (8.384191 seconds, I think I'll finish in 47.7898887 minutes)
Iteration 659 (change = 0.000205165389164) (8.170305 seconds, I think I'll finish in 46.43456675 minutes)
Iteration 660 (change = 0.0002046289058) (7.954522 seconds, I think I'll finish in 45.0756246667 minutes)
Iteration 661 (change = 0.000204417981216) (8.1263 seconds, I think I'll finish in 45.913595 minutes)
Iteration 662 (change = 0.000204507294738) (8.625778 seconds, I think I'll finish in 48.5918827333 minutes)
Iteration 663 (change = 0.000204876986084) (8.027213 seconds, I think I'll finish in 45.0861796833 minutes)
Iteration 664 (change = 0.000205521934099) (7.984055 seconds, I think I'll finish in 44.710708 minutes)
Iteration 665 (change = 0.000206429413051) (8.012102 seconds, I think I'll finish in 44.7342361667 minutes)
Iteration 666 (change = 0.000207587673015) (8.004432 seconds, I think I'll finish in 44.5580048 minutes)
Iteration 667 (change = 0.000208992201645) (8.001197 seconds, I think I'll finish in 44.40664335 minutes)
Iteration 668 (change = 0.000210642078084) (8.003431 seconds, I think I'll finish in 44.2856515333 minutes)
Iteration 669 (change = 0.000212506232433) (8.315402 seconds, I think I'll finish in 45.8733010333 minutes)
Iteration 670 (change = 0.000214610699669) (8.072653 seconds, I think I'll finish in 44.3995915 minutes)
Iteration 671 (change = 0.000216961561354) (8.025285 seconds, I think I'll finish in 44.00531275 minutes)
Iteration 672 (change = 0.000219561904254) (7.945969 seconds, I think I'll finish in 43.4379638667 minutes)
Iteration 673 (change = 0.000222408208647) (7.995756 seconds, I think I'll finish in 43.5768702 minutes)
Iteration 674 (change = 0.00022549710407) (7.971635 seconds, I think I'll finish in 43.3125501667 minutes)
Iteration 675 (change = 0.00022880807116) (8.115 seconds, I think I'll finish in 43.95625 minutes)
Iteration 676 (change = 0.000232294336739) (8.024423 seconds, I think I'll finish in 43.3318842 minutes)
Iteration 677 (change = 0.000235729187773) (8.015415 seconds, I think I'll finish in 43.14965075 minutes)
Iteration 678 (change = 0.000238356798945) (8.005214 seconds, I think I'll finish in 42.9613151333 minutes)
Iteration 679 (change = 0.000237179223739) (7.984052 seconds, I think I'll finish in 42.7146782 minutes)
Iteration 680 (change = 0.000231628038771) (7.924977 seconds, I think I'll finish in 42.266544 minutes)
Iteration 681 (change = 0.00023601843892) (8.071058 seconds, I think I'll finish in 42.9111250333 minutes)
Iteration 682 (change = 0.000240818031385) (7.849645 seconds, I think I'll finish in 41.6031185 minutes)
Iteration 683 (change = 0.0002430339327) (7.88033 seconds, I think I'll finish in 41.6344101667 minutes)
Iteration 684 (change = 0.000234762852178) (7.896252 seconds, I think I'll finish in 41.5869272 minutes)
Iteration 685 (change = 0.000215704708273) (8.176025 seconds, I think I'll finish in 42.92413125 minutes)
Iteration 686 (change = 0.000200990432461) (8.087276 seconds, I think I'll finish in 42.3234110667 minutes)
Iteration 687 (change = 0.00018827266561) (7.964587 seconds, I think I'll finish in 41.5485955167 minutes)
Iteration 688 (change = 0.000177145064703) (7.931552 seconds, I think I'll finish in 41.2440704 minutes)
Iteration 689 (change = 0.000167404798937) (7.908729 seconds, I think I'll finish in 40.99357865 minutes)
Iteration 690 (change = 0.000158873007326) (7.929768 seconds, I think I'll finish in 40.970468 minutes)
Iteration 691 (change = 0.000151393561131) (8.110038 seconds, I think I'll finish in 41.7666957 minutes)
Iteration 692 (change = 0.000144827294858) (7.975377 seconds, I think I'll finish in 40.9402686 minutes)
Iteration 693 (change = 0.000139047721055) (7.994861 seconds, I think I'll finish in 40.9070387833 minutes)
Iteration 694 (change = 0.0001339298394) (7.85498 seconds, I think I'll finish in 40.060398 minutes)
Iteration 695 (change = 0.000129300220821) (8.060176 seconds, I think I'll finish in 40.9725613333 minutes)
Iteration 696 (change = 0.00012485255908) (7.968211 seconds, I think I'll finish in 40.3722690667 minutes)
Iteration 697 (change = 0.000119671325724) (7.928845 seconds, I think I'll finish in 40.04066725 minutes)
Iteration 698 (change = 0.000111134187496) (7.913121 seconds, I think I'll finish in 39.8293757 minutes)
Iteration 699 (change = 0.000104412705368) (7.939599 seconds, I think I'll finish in 39.83032165 minutes)
Iteration 700 (change = 0.00010222024527) (7.896079 seconds, I think I'll finish in 39.480395 minutes)
Iteration 701 (change = 0.000100338753839) (10.099036 seconds, I think I'll finish in 50.3268627333 minutes)
Iteration 702 (change = 9.87044658632e-05) (9.186152 seconds, I think I'll finish in 45.6245549333 minutes)
Iteration 703 (change = 9.72946612354e-05) (8.063241 seconds, I think I'll finish in 39.91304295 minutes)
Iteration 704 (change = 9.60865355942e-05) (7.952406 seconds, I think I'll finish in 39.2318696 minutes)
Iteration 705 (change = 9.50625040235e-05) (8.083019 seconds, I think I'll finish in 39.7415100833 minutes)
Iteration 706 (change = 9.42063433943e-05) (7.905901 seconds, I think I'll finish in 38.7389149 minutes)
Iteration 707 (change = 9.35032399152e-05) (7.870681 seconds, I think I'll finish in 38.4351588833 minutes)
Iteration 708 (change = 9.29392765718e-05) (7.818416 seconds, I think I'll finish in 38.0496245333 minutes)
Iteration 709 (change = 9.25019876898e-05) (7.867489 seconds, I think I'll finish in 38.15732165 minutes)
Iteration 710 (change = 9.22223151337e-05) (7.910671 seconds, I think I'll finish in 38.2349098333 minutes)
Iteration 711 (change = 9.205828254e-05) (7.990704 seconds, I think I'll finish in 38.4885576 minutes)
Iteration 712 (change = 9.19808211593e-05) (7.791272 seconds, I think I'll finish in 37.3981056 minutes)
Iteration 713 (change = 9.19836720106e-05) (7.866377 seconds, I think I'll finish in 37.6275033167 minutes)
Iteration 714 (change = 9.20604387459e-05) (7.934333 seconds, I think I'll finish in 37.8203206333 minutes)
Iteration 715 (change = 9.22063325821e-05) (7.84925 seconds, I think I'll finish in 37.2839375 minutes)
Iteration 716 (change = 9.24172086477e-05) (7.989522 seconds, I think I'll finish in 37.8170708 minutes)
Iteration 717 (change = 9.26930462386e-05) (8.021733 seconds, I think I'll finish in 37.83584065 minutes)
Iteration 718 (change = 9.30360917799e-05) (7.829614 seconds, I think I'll finish in 36.7991858 minutes)
Iteration 719 (change = 9.34019113868e-05) (7.911673 seconds, I think I'll finish in 37.0530018833 minutes)
Iteration 720 (change = 9.38185098615e-05) (7.872722 seconds, I think I'll finish in 36.7393693333 minutes)
Iteration 721 (change = 9.42927849652e-05) (7.913376 seconds, I think I'll finish in 36.7971984 minutes)
Iteration 722 (change = 9.48357416731e-05) (7.832271 seconds, I think I'll finish in 36.2895223 minutes)
Iteration 723 (change = 9.54379540216e-05) (7.905473 seconds, I think I'll finish in 36.4969336833 minutes)
Iteration 724 (change = 9.59998460826e-05) (7.861754 seconds, I think I'll finish in 36.1640684 minutes)
Iteration 725 (change = 9.660002887e-05) (7.933968 seconds, I think I'll finish in 36.36402 minutes)
Iteration 726 (change = 9.72558852661e-05) (7.744184 seconds, I think I'll finish in 35.3651069333 minutes)
Iteration 727 (change = 9.79492880113e-05) (7.962694 seconds, I think I'll finish in 36.2302577 minutes)
Iteration 728 (change = 9.86772381784e-05) (7.846017 seconds, I think I'll finish in 35.5686104 minutes)
Iteration 729 (change = 9.94401982645e-05) (8.082524 seconds, I think I'll finish in 36.5060667333 minutes)
Iteration 730 (change = 0.000100231147064) (7.913269 seconds, I think I'll finish in 35.6097105 minutes)
Iteration 731 (change = 0.000101049485958) (7.900765 seconds, I think I'll finish in 35.4217630833 minutes)
Iteration 732 (change = 0.000101897811784) (8.308754 seconds, I think I'll finish in 37.1124345333 minutes)
Iteration 733 (change = 0.000102774516529) (7.969395 seconds, I think I'll finish in 35.46380775 minutes)
Iteration 734 (change = 0.000103677860841) (8.048153 seconds, I think I'll finish in 35.6801449667 minutes)
Iteration 735 (change = 0.000104607228824) (8.319822 seconds, I think I'll finish in 36.7458805 minutes)
Iteration 736 (change = 0.000105562149546) (8.26015 seconds, I think I'll finish in 36.34466 minutes)
Iteration 737 (change = 0.000106541731931) (8.004691 seconds, I think I'll finish in 35.0872288833 minutes)
Iteration 738 (change = 0.000107545460254) (8.14207 seconds, I think I'll finish in 35.5537056667 minutes)
Iteration 739 (change = 0.000108573780665) (8.538085 seconds, I think I'll finish in 37.14066975 minutes)
Iteration 740 (change = 0.000109625216253) (8.601305 seconds, I think I'll finish in 37.2723216667 minutes)
Iteration 741 (change = 0.000110699275395) (8.06236 seconds, I think I'll finish in 34.8025206667 minutes)
Iteration 742 (change = 0.000111795525669) (7.81637 seconds, I think I'll finish in 33.610391 minutes)
Iteration 743 (change = 0.000112913795708) (7.945549 seconds, I think I'll finish in 34.0334348833 minutes)
Iteration 744 (change = 0.000114054195198) (7.878449 seconds, I think I'll finish in 33.6147157333 minutes)
Iteration 745 (change = 0.000115215934003) (7.903674 seconds, I think I'll finish in 33.5906145 minutes)
Iteration 746 (change = 0.000116398693945) (7.855231 seconds, I think I'll finish in 33.2538112333 minutes)
Iteration 747 (change = 0.000117601925777) (8.19007 seconds, I think I'll finish in 34.5347951667 minutes)
Iteration 748 (change = 0.000118825286337) (7.958825 seconds, I think I'll finish in 33.427065 minutes)
Iteration 749 (change = 0.00012006846434) (8.061642 seconds, I think I'll finish in 33.7245357 minutes)
Iteration 750 (change = 0.000121331114325) (7.793464 seconds, I think I'll finish in 32.4727666667 minutes)
Iteration 751 (change = 0.000122613501043) (7.94067 seconds, I think I'll finish in 32.9537805 minutes)
Iteration 752 (change = 0.000123915168477) (7.918421 seconds, I think I'll finish in 32.7294734667 minutes)
Iteration 753 (change = 0.000125235332183) (7.924406 seconds, I think I'll finish in 32.6221380333 minutes)
Iteration 754 (change = 0.000126573855007) (7.900473 seconds, I think I'll finish in 32.3919393 minutes)
Iteration 755 (change = 0.000127930695747) (7.952752 seconds, I think I'll finish in 32.4737373333 minutes)
Iteration 756 (change = 0.000129304966894) (7.973835 seconds, I think I'll finish in 32.426929 minutes)
Iteration 757 (change = 0.000130696358139) (7.917854 seconds, I think I'll finish in 32.0673087 minutes)
Iteration 758 (change = 0.000132104636265) (7.755183 seconds, I think I'll finish in 31.2792381 minutes)
Iteration 759 (change = 0.00013352914552) (7.719232 seconds, I think I'll finish in 31.0055818667 minutes)
Iteration 760 (change = 0.000134970007891) (7.785202 seconds, I think I'll finish in 31.140808 minutes)
Iteration 761 (change = 0.000136426227321) (7.810557 seconds, I think I'll finish in 31.11205205 minutes)
Iteration 762 (change = 0.000137897304822) (7.903898 seconds, I think I'll finish in 31.3521287333 minutes)
Iteration 763 (change = 0.000139382656175) (8.29819 seconds, I think I'll finish in 32.7778505 minutes)
Iteration 764 (change = 0.000140884154347) (7.821698 seconds, I think I'll finish in 30.7653454667 minutes)
Iteration 765 (change = 0.000142399447866) (7.923167 seconds, I think I'll finish in 31.0324040833 minutes)
Iteration 766 (change = 0.000143927968575) (7.876825 seconds, I think I'll finish in 30.7196175 minutes)
Iteration 767 (change = 0.000145481960781) (7.928035 seconds, I think I'll finish in 30.7872025833 minutes)
Iteration 768 (change = 0.000147034808807) (7.843514 seconds, I think I'll finish in 30.3282541333 minutes)
Iteration 769 (change = 0.00014858180853) (7.860945 seconds, I think I'll finish in 30.26463825 minutes)
Iteration 770 (change = 0.000150146439711) (7.822397 seconds, I think I'll finish in 29.9858551667 minutes)
Iteration 771 (change = 0.000151723371284) (7.916519 seconds, I think I'll finish in 30.2147141833 minutes)
Iteration 772 (change = 0.00015330868345) (7.842073 seconds, I think I'll finish in 29.7998774 minutes)
Iteration 773 (change = 0.000154901460907) (7.896881 seconds, I think I'll finish in 29.8765331167 minutes)
Iteration 774 (change = 0.000156500167323) (7.83952 seconds, I think I'll finish in 29.5288586667 minutes)
Iteration 775 (change = 0.000158103699959) (8.018973 seconds, I think I'll finish in 30.07114875 minutes)
Iteration 776 (change = 0.000159710810786) (7.858861 seconds, I think I'll finish in 29.3397477333 minutes)
Iteration 777 (change = 0.000161320306466) (7.850446 seconds, I think I'll finish in 29.1774909667 minutes)
Iteration 778 (change = 0.000162931018679) (7.937508 seconds, I think I'll finish in 29.3687796 minutes)
Iteration 779 (change = 0.000164541665807) (7.914404 seconds, I think I'll finish in 29.1513880667 minutes)
Iteration 780 (change = 0.000166150972466) (7.818167 seconds, I think I'll finish in 28.6666123333 minutes)
Iteration 781 (change = 0.000167757628205) (7.886677 seconds, I think I'll finish in 28.78637105 minutes)
Iteration 782 (change = 0.000169360218925) (8.414399 seconds, I think I'll finish in 30.5723163667 minutes)
Iteration 783 (change = 0.000170957375378) (8.246841 seconds, I think I'll finish in 29.82607495 minutes)
Iteration 784 (change = 0.000172547626049) (7.855157 seconds, I think I'll finish in 28.2785652 minutes)
Iteration 785 (change = 0.000174129414015) (7.89144 seconds, I think I'll finish in 28.27766 minutes)
Iteration 786 (change = 0.000175701224114) (7.870798 seconds, I think I'll finish in 28.0725128667 minutes)
Iteration 787 (change = 0.000177261600843) (8.282315 seconds, I think I'll finish in 29.40221825 minutes)
Iteration 788 (change = 0.00017880899112) (7.878563 seconds, I think I'll finish in 27.8375892667 minutes)
Iteration 789 (change = 0.000180342388815) (7.866604 seconds, I think I'll finish in 27.6642240667 minutes)
Iteration 790 (change = 0.000181859861493) (7.822846 seconds, I think I'll finish in 27.379961 minutes)
Iteration 791 (change = 0.000183355149419) (7.91642 seconds, I think I'll finish in 27.5755296667 minutes)
Iteration 792 (change = 0.000184832146169) (7.759579 seconds, I think I'll finish in 26.8998738667 minutes)
Iteration 793 (change = 0.000186288490309) (7.853158 seconds, I think I'll finish in 27.0933951 minutes)
Iteration 794 (change = 0.000187721817438) (7.950404 seconds, I think I'll finish in 27.2963870667 minutes)
Iteration 795 (change = 0.000189130380137) (7.988454 seconds, I think I'll finish in 27.2938845 minutes)
Iteration 796 (change = 0.000190512598509) (7.913741 seconds, I think I'll finish in 26.9067194 minutes)
Iteration 797 (change = 0.000191866924325) (7.941713 seconds, I think I'll finish in 26.8694623167 minutes)
Iteration 798 (change = 0.000193191798967) (7.811156 seconds, I think I'll finish in 26.2975585333 minutes)
Iteration 799 (change = 0.000194485756401) (7.788624 seconds, I think I'll finish in 26.0918904 minutes)
Iteration 800 (change = 0.000195747426664) (7.710539 seconds, I think I'll finish in 25.7017966667 minutes)
Iteration 801 (change = 0.000196975500683) (7.807191 seconds, I think I'll finish in 25.89385015 minutes)
Iteration 802 (change = 0.000198168735911) (7.688612 seconds, I think I'll finish in 25.3724196 minutes)
Iteration 803 (change = 0.000199326080932) (7.972651 seconds, I think I'll finish in 26.1768707833 minutes)
Iteration 804 (change = 0.000200446699603) (7.73311 seconds, I think I'll finish in 25.2614926667 minutes)
Iteration 805 (change = 0.000201529496123) (7.880376 seconds, I think I'll finish in 25.611222 minutes)
Iteration 806 (change = 0.000202573439595) (7.889316 seconds, I think I'll finish in 25.5087884 minutes)
Iteration 807 (change = 0.000203577536584) (7.873845 seconds, I think I'll finish in 25.32753475 minutes)
Iteration 808 (change = 0.000204540584622) (7.788039 seconds, I think I'll finish in 24.9217248 minutes)
Iteration 809 (change = 0.000205460845794) (8.149936 seconds, I think I'll finish in 25.9439629333 minutes)
Iteration 810 (change = 0.000206335487568) (8.103308 seconds, I think I'll finish in 25.6604753333 minutes)
Iteration 811 (change = 0.000207159352701) (8.106886 seconds, I think I'll finish in 25.5366909 minutes)
Iteration 812 (change = 0.00020792241464) (7.845554 seconds, I think I'll finish in 24.5827358667 minutes)
Iteration 813 (change = 0.000208603875264) (7.923084 seconds, I think I'll finish in 24.6936118 minutes)
Iteration 814 (change = 0.00020915665563) (8.216936 seconds, I think I'll finish in 25.4725016 minutes)
Iteration 815 (change = 0.0002094629442) (8.018063 seconds, I think I'll finish in 24.7223609167 minutes)
Iteration 816 (change = 0.000209152848664) (7.944743 seconds, I think I'll finish in 24.3638785333 minutes)
Iteration 817 (change = 0.000206836997511) (7.917756 seconds, I think I'll finish in 24.1491558 minutes)
Iteration 818 (change = 0.000200234152009) (8.847295 seconds, I think I'll finish in 26.8367948333 minutes)
Iteration 819 (change = 0.000197141239174) (10.543142 seconds, I think I'll finish in 31.8051450333 minutes)
Iteration 820 (change = 0.000195959208012) (7.925015 seconds, I think I'll finish in 23.775045 minutes)
Iteration 821 (change = 0.000194831834011) (7.976783 seconds, I think I'll finish in 23.7974026167 minutes)
Iteration 822 (change = 0.000193709593499) (7.884478 seconds, I think I'll finish in 23.3906180667 minutes)
Iteration 823 (change = 0.000192565889471) (7.894632 seconds, I think I'll finish in 23.2891644 minutes)
Iteration 824 (change = 0.000191375633916) (7.919501 seconds, I think I'll finish in 23.2305362667 minutes)
Iteration 825 (change = 0.000190136037045) (8.071489 seconds, I think I'll finish in 23.5418429167 minutes)
Iteration 826 (change = 0.000188833086813) (8.018418 seconds, I think I'll finish in 23.2534122 minutes)
Iteration 827 (change = 0.000187456487631) (8.046445 seconds, I think I'll finish in 23.2005830833 minutes)
Iteration 828 (change = 0.000185997345628) (7.899561 seconds, I think I'll finish in 22.6454082 minutes)
Iteration 829 (change = 0.000184446934684) (7.939942 seconds, I think I'll finish in 22.6288347 minutes)
Iteration 830 (change = 0.00018279482234) (8.413889 seconds, I think I'll finish in 23.8393521667 minutes)
Iteration 831 (change = 0.000181028085176) (8.017904 seconds, I think I'll finish in 22.5837629333 minutes)
Iteration 832 (change = 0.000179128438655) (7.90473 seconds, I think I'll finish in 22.133244 minutes)
Iteration 833 (change = 0.000177067779704) (8.031747 seconds, I think I'll finish in 22.35502915 minutes)
Iteration 834 (change = 0.000174799302973) (7.958786 seconds, I think I'll finish in 22.0193079333 minutes)
Iteration 835 (change = 0.000172238505896) (8.014205 seconds, I think I'll finish in 22.03906375 minutes)
Iteration 836 (change = 0.000169217638588) (7.849692 seconds, I think I'll finish in 21.4558248 minutes)
Iteration 837 (change = 0.000165359951067) (7.832882 seconds, I think I'll finish in 21.2793294333 minutes)
Iteration 838 (change = 0.000159666074206) (7.847425 seconds, I think I'll finish in 21.1880475 minutes)
Iteration 839 (change = 0.000148815733603) (7.949845 seconds, I think I'll finish in 21.3320840833 minutes)
Iteration 840 (change = 0.000119286052236) (8.147927 seconds, I think I'll finish in 21.7278053333 minutes)
Iteration 841 (change = 6.54632465092e-05) (8.083667 seconds, I think I'll finish in 21.42171755 minutes)
Iteration 842 (change = 5.88355851013e-05) (7.872644 seconds, I think I'll finish in 20.7312958667 minutes)
Iteration 843 (change = 5.71356410201e-05) (7.843144 seconds, I think I'll finish in 20.5228934667 minutes)
Iteration 844 (change = 5.55991964689e-05) (7.755361 seconds, I think I'll finish in 20.1639386 minutes)
Iteration 845 (change = 5.41601187711e-05) (7.871586 seconds, I think I'll finish in 20.3349305 minutes)
Iteration 846 (change = 5.27864227269e-05) (7.912747 seconds, I think I'll finish in 20.3093839667 minutes)
Iteration 847 (change = 5.14688261382e-05) (7.857053 seconds, I think I'll finish in 20.03548515 minutes)
Iteration 848 (change = 5.01958363525e-05) (7.925282 seconds, I think I'll finish in 20.0773810667 minutes)
Iteration 849 (change = 4.89619172449e-05) (7.938199 seconds, I think I'll finish in 19.9778008167 minutes)
Iteration 850 (change = 4.77635572428e-05) (7.863178 seconds, I think I'll finish in 19.657945 minutes)
Iteration 851 (change = 4.6599265926e-05) (7.83868 seconds, I think I'll finish in 19.4660553333 minutes)
Iteration 852 (change = 4.54667484283e-05) (7.922846 seconds, I think I'll finish in 19.5430201333 minutes)
Iteration 853 (change = 4.43651062483e-05) (7.876182 seconds, I think I'll finish in 19.2966459 minutes)
Iteration 854 (change = 4.32943177718e-05) (7.702139 seconds, I think I'll finish in 18.7418715667 minutes)
Iteration 855 (change = 4.2255216777e-05) (7.919063 seconds, I think I'll finish in 19.1377355833 minutes)
Iteration 856 (change = 4.12425029504e-05) (8.008349 seconds, I think I'll finish in 19.2200376 minutes)
Iteration 857 (change = 4.02559976679e-05) (7.959664 seconds, I think I'll finish in 18.9705325333 minutes)
Iteration 858 (change = 3.93008477027e-05) (7.74277 seconds, I think I'll finish in 18.3245556667 minutes)
Iteration 859 (change = 3.83758097602e-05) (7.953959 seconds, I think I'll finish in 18.69180365 minutes)
Iteration 860 (change = 3.74790019588e-05) (7.799074 seconds, I think I'll finish in 18.1978393333 minutes)
Iteration 861 (change = 3.66101449211e-05) (7.911592 seconds, I think I'll finish in 18.3285214667 minutes)
Iteration 862 (change = 3.57692456673e-05) (7.856487 seconds, I think I'll finish in 18.0699201 minutes)
Iteration 863 (change = 3.49559208004e-05) (7.878364 seconds, I think I'll finish in 17.9889311333 minutes)
Iteration 864 (change = 3.41700020427e-05) (7.7421 seconds, I think I'll finish in 17.54876 minutes)
Iteration 865 (change = 3.34108913881e-05) (7.905561 seconds, I think I'll finish in 17.78751225 minutes)
Iteration 866 (change = 3.26782337747e-05) (7.786984 seconds, I think I'll finish in 17.3909309333 minutes)
Iteration 867 (change = 3.19716616717e-05) (7.924972 seconds, I think I'll finish in 17.5670212667 minutes)
Iteration 868 (change = 3.12907615072e-05) (7.880259 seconds, I think I'll finish in 17.3365698 minutes)
Iteration 869 (change = 3.06351666327e-05) (7.88417 seconds, I think I'll finish in 17.2137711667 minutes)
Iteration 870 (change = 3.00044277612e-05) (7.787005 seconds, I think I'll finish in 16.8718441667 minutes)
Iteration 871 (change = 2.93990008201e-05) (8.157189 seconds, I think I'll finish in 17.53795635 minutes)
Iteration 872 (change = 2.88189632763e-05) (8.113105 seconds, I think I'll finish in 17.3079573333 minutes)
Iteration 873 (change = 2.82624361216e-05) (8.002446 seconds, I think I'll finish in 16.9385107 minutes)
Iteration 874 (change = 2.77287840114e-05) (7.769162 seconds, I think I'll finish in 16.3152402 minutes)
Iteration 875 (change = 2.72174659107e-05) (7.918061 seconds, I think I'll finish in 16.4959604167 minutes)
Iteration 876 (change = 2.67279869331e-05) (7.867392 seconds, I think I'll finish in 16.2592768 minutes)
Iteration 877 (change = 2.62597602218e-05) (8.364219 seconds, I think I'll finish in 17.14664895 minutes)
Iteration 878 (change = 2.58126024871e-05) (7.733817 seconds, I think I'll finish in 15.7254279 minutes)
Iteration 879 (change = 2.53856483947e-05) (7.856087 seconds, I think I'll finish in 15.8431087833 minutes)
Iteration 880 (change = 2.49783265624e-05) (7.788848 seconds, I think I'll finish in 15.577696 minutes)
Iteration 881 (change = 2.45901301973e-05) (7.958504 seconds, I think I'll finish in 15.7843662667 minutes)
Iteration 882 (change = 2.42205225437e-05) (7.80384 seconds, I think I'll finish in 15.347552 minutes)
Iteration 883 (change = 2.38689783123e-05) (8.201296 seconds, I think I'll finish in 15.9925272 minutes)
Iteration 884 (change = 2.35349707538e-05) (7.898509 seconds, I think I'll finish in 15.2704507333 minutes)
Iteration 885 (change = 2.32179812852e-05) (7.897741 seconds, I think I'll finish in 15.1373369167 minutes)
Iteration 886 (change = 2.29174970213e-05) (7.828428 seconds, I think I'll finish in 14.8740132 minutes)
Iteration 887 (change = 2.2633040199e-05) (7.820292 seconds, I think I'll finish in 14.7282166 minutes)
Iteration 888 (change = 2.23641333358e-05) (7.971523 seconds, I think I'll finish in 14.8801762667 minutes)
Iteration 889 (change = 2.21103180218e-05) (8.076995 seconds, I think I'll finish in 14.94244075 minutes)
Iteration 890 (change = 2.18710914064e-05) (7.989964 seconds, I think I'll finish in 14.6482673333 minutes)
Iteration 891 (change = 2.16461923222e-05) (7.911238 seconds, I think I'll finish in 14.3720823667 minutes)
Iteration 892 (change = 2.14358406592e-05) (7.970091 seconds, I think I'll finish in 14.3461638 minutes)
Iteration 893 (change = 2.12386806532e-05) (8.022223 seconds, I think I'll finish in 14.3062976833 minutes)
Iteration 894 (change = 2.10542911821e-05) (7.827774 seconds, I think I'll finish in 13.8290674 minutes)
Iteration 895 (change = 2.08822657145e-05) (7.951164 seconds, I think I'll finish in 13.914537 minutes)
Iteration 896 (change = 2.07222316252e-05) (7.809842 seconds, I think I'll finish in 13.5370594667 minutes)
Iteration 897 (change = 2.05738012738e-05) (8.571586 seconds, I think I'll finish in 14.7145559667 minutes)
Iteration 898 (change = 2.04366203926e-05) (7.865682 seconds, I think I'll finish in 13.3716594 minutes)
Iteration 899 (change = 2.03103671209e-05) (7.961275 seconds, I think I'll finish in 13.4014795833 minutes)
Iteration 900 (change = 2.01946954599e-05) (7.938451 seconds, I think I'll finish in 13.2307516667 minutes)
Iteration 901 (change = 2.00892392272e-05) (7.930898 seconds, I think I'll finish in 13.0859817 minutes)
Iteration 902 (change = 1.99936626642e-05) (7.858527 seconds, I think I'll finish in 12.8355941 minutes)
Iteration 903 (change = 1.99076613739e-05) (7.958368 seconds, I think I'll finish in 12.8660282667 minutes)
Iteration 904 (change = 1.98309560966e-05) (8.0419 seconds, I think I'll finish in 12.86704 minutes)
Iteration 905 (change = 1.97632597331e-05) (7.946864 seconds, I think I'll finish in 12.5825346667 minutes)
Iteration 906 (change = 1.97042957792e-05) (7.922653 seconds, I think I'll finish in 12.4121563667 minutes)
Iteration 907 (change = 1.96537890249e-05) (7.983749 seconds, I think I'll finish in 12.37481095 minutes)
Iteration 908 (change = 1.96114896365e-05) (7.956953 seconds, I think I'll finish in 12.2006612667 minutes)
Iteration 909 (change = 1.95771641166e-05) (7.9414 seconds, I think I'll finish in 12.0444566667 minutes)
Iteration 910 (change = 1.95505819631e-05) (7.790311 seconds, I think I'll finish in 11.6854665 minutes)
Iteration 911 (change = 1.95315273453e-05) (7.853723 seconds, I think I'll finish in 11.6496891167 minutes)
Iteration 912 (change = 1.95197856938e-05) (7.882589 seconds, I think I'll finish in 11.5611305333 minutes)
Iteration 913 (change = 1.95151585513e-05) (7.910962 seconds, I think I'll finish in 11.4708949 minutes)
Iteration 914 (change = 1.95174610446e-05) (7.864721 seconds, I think I'll finish in 11.2727667667 minutes)
Iteration 915 (change = 1.95265306446e-05) (7.920484 seconds, I think I'll finish in 11.2206856667 minutes)
Iteration 916 (change = 1.95422931358e-05) (7.926813 seconds, I think I'll finish in 11.0975382 minutes)
Iteration 917 (change = 1.95644724495e-05) (7.977759 seconds, I think I'll finish in 11.03589995 minutes)
Iteration 918 (change = 1.95929178173e-05) (7.668556 seconds, I think I'll finish in 10.4803598667 minutes)
Iteration 919 (change = 1.96274879932e-05) (8.159613 seconds, I think I'll finish in 11.01547755 minutes)
Iteration 920 (change = 1.96680436944e-05) (8.117922 seconds, I think I'll finish in 10.823896 minutes)
Iteration 921 (change = 1.9714452612e-05) (7.902894 seconds, I think I'll finish in 10.4054771 minutes)
Iteration 922 (change = 1.97665900518e-05) (7.854557 seconds, I think I'll finish in 10.2109241 minutes)
Iteration 923 (change = 1.98243452566e-05) (7.877448 seconds, I think I'll finish in 10.1093916 minutes)
Iteration 924 (change = 1.98876146596e-05) (8.222416 seconds, I think I'll finish in 10.4150602667 minutes)
Iteration 925 (change = 1.995630814e-05) (7.916737 seconds, I think I'll finish in 9.89592125 minutes)
Iteration 926 (change = 2.00303352314e-05) (7.912768 seconds, I think I'll finish in 9.75908053333 minutes)
Iteration 927 (change = 2.01096054628e-05) (7.933727 seconds, I think I'll finish in 9.65270118333 minutes)
Iteration 928 (change = 2.01940382874e-05) (7.929746 seconds, I think I'll finish in 9.5156952 minutes)
Iteration 929 (change = 2.0283571488e-05) (7.910654 seconds, I think I'll finish in 9.36094056667 minutes)
Iteration 930 (change = 2.03781406947e-05) (7.716122 seconds, I think I'll finish in 9.00214233333 minutes)
Iteration 931 (change = 2.04776982308e-05) (7.973131 seconds, I think I'll finish in 9.16910065 minutes)
Iteration 932 (change = 2.05822262054e-05) (7.926373 seconds, I think I'll finish in 8.98322273333 minutes)
Iteration 933 (change = 2.06917727858e-05) (7.794501 seconds, I think I'll finish in 8.70385945 minutes)
Iteration 934 (change = 2.08064794058e-05) (7.777919 seconds, I think I'll finish in 8.5557109 minutes)
Iteration 935 (change = 2.09271198577e-05) (8.106485 seconds, I think I'll finish in 8.78202541667 minutes)
Iteration 936 (change = 2.10566481353e-05) (10.56805 seconds, I think I'll finish in 11.2725866667 minutes)
Iteration 937 (change = 2.11872080779e-05) (8.096793 seconds, I think I'll finish in 8.50163265 minutes)
Iteration 938 (change = 2.13131200131e-05) (7.819002 seconds, I think I'll finish in 8.0796354 minutes)
Iteration 939 (change = 2.14485560597e-05) (7.979723 seconds, I think I'll finish in 8.11271838333 minutes)
Iteration 940 (change = 2.15903833237e-05) (7.854919 seconds, I think I'll finish in 7.854919 minutes)
Iteration 941 (change = 2.17374080329e-05) (7.887417 seconds, I think I'll finish in 7.75596005 minutes)
Iteration 942 (change = 2.18892718997e-05) (7.799271 seconds, I think I'll finish in 7.5392953 minutes)
Iteration 943 (change = 2.20458725153e-05) (8.004132 seconds, I think I'll finish in 7.6039254 minutes)
Iteration 944 (change = 2.22071853772e-05) (7.861115 seconds, I think I'll finish in 7.33704066667 minutes)
Iteration 945 (change = 2.23732198645e-05) (7.916184 seconds, I think I'll finish in 7.256502 minutes)
Iteration 946 (change = 2.25439952879e-05) (7.92106 seconds, I think I'll finish in 7.128954 minutes)
Iteration 947 (change = 2.27195309878e-05) (7.900315 seconds, I think I'll finish in 6.97861158333 minutes)
Iteration 948 (change = 2.28998615032e-05) (7.827117 seconds, I think I'll finish in 6.7835014 minutes)
Iteration 949 (change = 2.30850141565e-05) (7.854927 seconds, I think I'll finish in 6.67668795 minutes)
Iteration 950 (change = 2.32750275337e-05) (8.040301 seconds, I think I'll finish in 6.70025083333 minutes)
Iteration 951 (change = 2.34699449544e-05) (8.019027 seconds, I think I'll finish in 6.54887205 minutes)
Iteration 952 (change = 2.36698157578e-05) (7.808348 seconds, I think I'll finish in 6.2466784 minutes)
Iteration 953 (change = 2.38747008324e-05) (7.855073 seconds, I think I'll finish in 6.15314051667 minutes)
Iteration 954 (change = 2.40846526539e-05) (7.873948 seconds, I think I'll finish in 6.03669346667 minutes)
Iteration 955 (change = 2.42997292836e-05) (7.890738 seconds, I think I'll finish in 5.9180535 minutes)
Iteration 956 (change = 2.45199981211e-05) (7.976751 seconds, I think I'll finish in 5.8496174 minutes)
Iteration 957 (change = 2.47455330591e-05) (7.881537 seconds, I think I'll finish in 5.64843485 minutes)
Iteration 958 (change = 2.49764111437e-05) (7.887751 seconds, I think I'll finish in 5.5214257 minutes)
Iteration 959 (change = 2.52127133305e-05) (7.984716 seconds, I think I'll finish in 5.4562226 minutes)
Iteration 960 (change = 2.54545312261e-05) (7.728053 seconds, I think I'll finish in 5.15203533333 minutes)
Iteration 961 (change = 2.570195036e-05) (7.931796 seconds, I think I'll finish in 5.1556674 minutes)
Iteration 962 (change = 2.59550740268e-05) (7.771704 seconds, I think I'll finish in 4.9220792 minutes)
Iteration 963 (change = 2.62139977038e-05) (7.888064 seconds, I think I'll finish in 4.86430613333 minutes)
Iteration 964 (change = 2.64788212835e-05) (7.717902 seconds, I think I'll finish in 4.6307412 minutes)
Iteration 965 (change = 2.67496609689e-05) (7.779655 seconds, I think I'll finish in 4.53813208333 minutes)
Iteration 966 (change = 2.7026627596e-05) (8.054627 seconds, I think I'll finish in 4.56428863333 minutes)
Iteration 967 (change = 2.73098394652e-05) (8.063611 seconds, I think I'll finish in 4.43498605 minutes)
Iteration 968 (change = 2.75994260004e-05) (7.851332 seconds, I think I'll finish in 4.18737706667 minutes)
Iteration 969 (change = 2.78955261658e-05) (7.872856 seconds, I think I'll finish in 4.06764226667 minutes)
Iteration 970 (change = 2.81982891227e-05) (8.007278 seconds, I think I'll finish in 4.003639 minutes)
Iteration 971 (change = 2.85078429645e-05) (8.62191 seconds, I think I'll finish in 4.1672565 minutes)
Iteration 972 (change = 2.88243383866e-05) (7.924996 seconds, I think I'll finish in 3.69833146667 minutes)
Iteration 973 (change = 2.91479364086e-05) (7.885022 seconds, I think I'll finish in 3.5482599 minutes)
Iteration 974 (change = 2.94788104104e-05) (7.89679 seconds, I think I'll finish in 3.42194233333 minutes)
Iteration 975 (change = 2.98171598371e-05) (8.356461 seconds, I think I'll finish in 3.48185875 minutes)
Iteration 976 (change = 3.01631258227e-05) (8.049073 seconds, I think I'll finish in 3.2196292 minutes)
Iteration 977 (change = 3.0516889822e-05) (7.904269 seconds, I think I'll finish in 3.02996978333 minutes)
Iteration 978 (change = 3.08786442914e-05) (7.943137 seconds, I think I'll finish in 2.91248356667 minutes)
Iteration 979 (change = 3.1248589002e-05) (7.844137 seconds, I think I'll finish in 2.74544795 minutes)
Iteration 980 (change = 3.16269304643e-05) (8.011509 seconds, I think I'll finish in 2.670503 minutes)
Iteration 981 (change = 3.20138828677e-05) (7.915534 seconds, I think I'll finish in 2.50658576667 minutes)
Iteration 982 (change = 3.24096698153e-05) (8.111078 seconds, I think I'll finish in 2.4333234 minutes)
Iteration 983 (change = 3.28145381962e-05) (7.967826 seconds, I think I'll finish in 2.2575507 minutes)
Iteration 984 (change = 3.32287383364e-05) (7.943647 seconds, I think I'll finish in 2.11830586667 minutes)
Iteration 985 (change = 3.36525396535e-05) (7.898419 seconds, I think I'll finish in 1.97460475 minutes)
Iteration 986 (change = 3.408618004e-05) (7.90317 seconds, I think I'll finish in 1.844073 minutes)
Iteration 987 (change = 3.45299201469e-05) (7.917384 seconds, I think I'll finish in 1.7154332 minutes)
Iteration 988 (change = 3.498405741e-05) (7.843314 seconds, I think I'll finish in 1.5686628 minutes)
Iteration 989 (change = 3.54488663874e-05) (7.944069 seconds, I think I'll finish in 1.45641265 minutes)
Iteration 990 (change = 3.59246502593e-05) (7.919394 seconds, I think I'll finish in 1.319899 minutes)
Iteration 991 (change = 3.64117248865e-05) (8.04609 seconds, I think I'll finish in 1.2069135 minutes)
Iteration 992 (change = 3.69104047264e-05) (7.944153 seconds, I think I'll finish in 1.0592204 minutes)
Iteration 993 (change = 3.74210283157e-05) (7.855713 seconds, I think I'll finish in 0.91649985 minutes)
Iteration 994 (change = 3.79439560301e-05) (7.810677 seconds, I think I'll finish in 0.7810677 minutes)
Iteration 995 (change = 3.84795710084e-05) (7.894401 seconds, I think I'll finish in 0.65786675 minutes)
Iteration 996 (change = 3.9028261851e-05) (7.820488 seconds, I think I'll finish in 0.521365866667 minutes)
Iteration 997 (change = 3.95904086398e-05) (7.906998 seconds, I think I'll finish in 0.3953499 minutes)
Iteration 998 (change = 4.01664066806e-05) (7.988731 seconds, I think I'll finish in 0.266291033333 minutes)
Iteration 999 (change = 4.07566843393e-05) (7.952874 seconds, I think I'll finish in 0.1325479 minutes)

In [519]:

    

out_name = 'massbank_wider_fixed_kde.dict'
mbd = vlda.make_dictionary(filename=out_name,metadata=c.metadata)


    

    




Done 500
Done 1000
Done 1500

In [520]:

    

# Merge the chemspider stuff from the other dict
with open('massbank_to_chemspider.dict') as f:
    mb2cs = pickle.load(f)


    

In [521]:

    

for doc in mbd['doc_metadata']:
    if doc in mb2cs:
        if len(mb2cs[doc])>0:
            mbd['doc_metadata'][doc]['csid'] = mb2cs[doc][0]


    

In [522]:

    

with open(out_name,'w') as f:
    pickle.dump(mbd,f)


    

In [458]:

    

# Alternative matching


    

In [523]:

    

with open(out_name,'r') as f:
    massbank_kde = pickle.load(f)


    

In [524]:

    

beta = vlda.beta_matrix.copy()
# Post analysis topic matching
prob_thresh = 0.5
topic_matches = {}
for topic in beer_topics['beta']:
    word_matches = {}
    topic_matches[topic] = []
    for word in beer_topics['beta'][topic]:
        split_word = word.split('_')
        if len(split_word) == 3:
            word_matches[word] = [w for w in massbank_kde['word_index'] if w.startswith(split_word[0]) and w.endswith(split_word[2])]

    for t in range(vlda.K):
        match_probability = 0
        other_probability = 0
        for word,matches in word_matches.items():
            for match in matches:
                match_probability += beta[t,massbank_kde['word_index'][match]]
                other_probability += beer_topics['beta'][topic][word]
        if other_probability > prob_thresh and match_probability > 0.2:
            topic_matches[topic].append(t)


    

In [525]:

    

for topic in beer_topics['beta']:
    print beer_topics['topic_metadata'][topic]['annotation']
    print "\t",topic_matches[topic]


    

    




Fragments indicative for ferulic acid based substructure (MzCloud)
	[0]
Unclear yet what these losses relate to.
	[]
Losses indicative of a hexose with NH2 group – EF fits
	[]
Possibly suggests iminosugar like substructure. Fragments have losses (H2O CO) related to sugar (polyhydroxylated structure)
	[]
Fragment ions possibly indicative for ribose substructure (MzCloud)
	[1]
Imidazole group linked to a carboxylgroup through one CH2 group i.e. like in imidazole acetic acid - Prevalent in Beer1
	[2, 12]
Unclear yet what these fragments relate to.
	[3]
Fragments indicative for asparagine substructure (MzCloud) – prevalent in Beer 3
	[4]
Fragment ions indicative for pipecolic acid (pipecolate) (MzCloud) - Quite prevalent especially in Beer2
	[5, 12]
Double water loss i.e. 2*H20 – Generic feature for metabolites containing several free OH groups attached to a aliphatic chain like sugars.
	[6, 10]
Fragment ions indicative for guanine (C5H5N5O) based substructure
	[7]
Unclear yet what these fragments relate to.
	[]
Fragments indicative for dihydroxylated benzene ring substructure (MzCloud) – C6H5O2 fragment corresponds to positively charged fragment with two hydroxyl groups.
	[8]
Unclear yet what these fragments relate to.
	[]
Fragments indicative for [phenylalanine-CHOOH] based substructure
	[9, 18, 30]
Loss possibly indicative of carboxylic acid group with 1-carbon attached.
	[10]
Fragment ion indicative for conjugation of a phosphate group  (H4O4P)
	[11]
Unclear if fragments relate to a specific substructure.
	[12]
Fragments indicative for histidine (C6H10N3O2) substructure (MzCloud)
	[13, 14, 373]
Unclear what these fragments relate to. Possibly small ring structure (CHN loss)
	[14]
Fragments indicative of a [hexose-H2O] substructure – i.e. indicative for a hexose (like glucose) conjugation (MzCloud)
	[15]
Fragment ions indicative for cytosine and a loss of conjugated deoxyribose – possibly combined due to many spectra that combine these two substructures. Loss of NH2 group is likely from remaining fragment after loss of deoxyribose. – Quite prevalent in Beer1
	[]
[Ribose (pentose C5-sugar)-H2O] related loss –  indicative for conjugated ribose sugar - EF fits
	[16]
C13 isotope peaks of proline (abundant ions taken for fragmentation)
	[]
Water loss  - indicative of a free hydroxyl group) – often seen in sugary structures
	[17]
Possibly suggests 2-pyyrolidine substructure – Mass2Motif not consistent over the four beers.
	[]
Small abundant and generic aromatic fragment found across several mass patterns.
	[18]
Oxygen-rich losses and fragments also occurring in hexose spectra – related to M2M 211 (hexose [glucose] conjugatation) – possibly hydrated-hexose loss?
	[1, 8, 19]
Unclear yet – possibly related to methylsuccinic acid….
	[]
Combined loss of H2O and CO – indicative for free carboxylic acid group (COOH) – generic substructure in amino acids and organic acids
	[20, 35]
Fragments (and loss) indicative for serine substructure (MzCloud) - Present in Beer 2 & Beer 3
	[]
Fragment ions indicative for aconitic acid substructure (C3HO2 fragment is quite specific)
	[]
Fragments indicative of [tryptophan-NH3] related substructure (C8H8N is the basic indole skeleton a fused benzene and 5 membered N-containing ring)
	[]
Possibly related to 14-Cyclohex-2-enedione substructure – double CO loss between fragments. Could be alkaloid fragments as well.
	[21]
Leucine related substructure (mzCloud) – prevalent in Beer 2
	[22]
Fragment and loss of [proline-H2O] - indicative for conjugated proline – EF fits
	[23]
Small nitrogen containing fragment ion – often proline or ornithine derived – most abundant fragment in all four beers
	[24]
Unclear yet what these fragments relate to.
	[25]
Combination of small losses – free carboxylgroup and acetyl group loss + loss of NH2 group in some cases
	[26]
Fragment ions possibly indicative for N-Methyl-oxo-pyrrolidinecarboxylic acid like structure (loss of free carboxyl group)
	[27]
Two small fragments unclear if they represent a substructure
	[28]
 Fragments indicative for tyrosine substructure (MzCloud)
	[]
Unclear yet what these fragments relate to.
	[29]
Fragments indicative for ethylphenol substructure (i.e. resulting from Tyramine – MzCloud)
	[8, 18, 30]
Fragment ions indicative for glutamine (C5H11N2O3) substructure
	[]
Combination of small losses (CO2 H2O etc.) – Unclear if they relate to a substructure loss.
	[31]
Fragments indicative adenine (C5H6N5) substructure – most prevalent in Beer3
	[32, 314]
Fragment ions indicative for trimethylated amine connected to a carboxylic acid group i.e. like in betaine (MzCloud)
	[]
Unclear if fragments relate to a specific substructure. H2O loss between fragments.
	[]
Fragments possibly suggest benzene ring substituted with one hydroxyl and one NH2 group (fragments point to orientation from 3-hydroxyanthranilic acid – i.e. MzCloud) – prevalent in Beer 2 and 4
	[]
Unclear what these fragments relate to.
	[]
Fragment ions indicative for 5-methylcytosine substructure (MzCloud) – prevalent in Beer 2
	[]
Fragments indicative for cinnamic acid (cinnamate) substructure (MzCloud)
	[18, 33]
Unclear what these fragments relate to.
	[]
Possibly suggests 2-pyridone/ol substructure
	[]
Possibly suggests 5-aminopentanol substructure
	[34]
Possibly suggests alkylbenzene substructure.
	[]
N/O-Acetylation (Beer 3)
	[20, 35]
Unclear yet what these fragments relate to.
	[]
Fragment ions indicative for pyroglutamic acid (pyroglutamate) or glutamine (both in MzCloud) – structure can be formed from glutamic acid (glutamate) in the mass spectrometer as well.
	[12, 36]
Unclear what these fragments relate to.
	[37]
Combination of small losses (i.e. NH3 and CH2O2)
	[38]
Possibly suggests 4-aminooxane-4-carboxylic acid like substructure?
	[]
Alkyl aromatic substructure – indicative for aromatic ring with 2-carbon alkyl chain attached i.e. phenylethene fragment from ethylbenzene as a result of the fragmentation process.
	[8, 39]
Fragment ions indicative for N-acetylputrescine substructure (MzCloud)
	[40, 232]
Unclear yet what these fragments relate to.
	[]
Possibly suggests methoxy-1H-indole-23-dione (methoxy-isatin) substructure
	[]
Amine loss - Indicative for free NH2 group in fragmented molecule
	[]
Aliphatic amine (NH3 loss indicates free NH2 group coupled to aliphatic chain)
	[42]
Unclear yet what these fragments relate to.
	[]
Combination of small losses (CO2 H2O etc.) – Unclear if they relate to a substructure loss.
	[]
Fragments indicative for y-aminobutyric acid (amimobutyrate) substructure (MzCloud) – present in Beer 2 & Beer 3 – in beer 3 mainly based on C4H7O2 fragment.
	[43]
Loss of [hexose-H2O] – indication of hexose conjugation (for example glucose)
	[44]
Possibly suggests agmatine based substructure (C5H15N4) with unknown conjugation….
	[]
Combinations of small generic losses like CH2O2 + CH4O = C2H6O3
	[]

In [526]:

    

ti = [(topic,vlda.topic_index[topic]) for topic in vlda.topic_index]
ti = sorted(ti,key = lambda x:x[1])
reverse,_ = zip(*ti)
for topic in topic_matches:
    for match in topic_matches[topic]:
        tname = reverse[match]
        if not 'annotation' in vlda.topic_metadata[tname]:
            vlda.topic_metadata[tname]['annotation'] = beer_topics['topic_metadata'][topic]['annotation']


    

In [527]:

    

mbd = vlda.make_dictionary(filename = out_name,metadata = c.metadata)


    

    




Done 500
Done 1000
Done 1500

In [ ]:

    




    

In [ ]:

    




    

In [ ]:

    




    

In [ ]:

    




    

In [515]:

    

formula = 'fragment_C8H5O'

tokens = formula.split('_')

mz = Formula(tokens[1]).compute_exact_mass()

ppm = 50
delta = ppm * mz/1e6
min_mz = mz - delta
max_mz = mz + delta

if tokens[0] == 'loss':
    is_loss = True
else:
    is_loss = False    


if is_loss:
    kde = c.loss_kde
    groups = c.loss_groups
    array = c.loss_array
    group_list = c.loss_group_list
else:
    kde = c.frag_kde
    groups = c.frag_groups
    array = c.frag_array
    group_list = c.frag_group_list


pos = np.where((array > min_mz) & (array < max_mz))[0]
print len(pos)
data = []

unique_groups = list(set(groups[pos]))

data.append(
    Scatter(
        x = array[pos],
        y = np.zeros_like(array[pos]),
        mode = 'markers',
    )
)
for group in unique_groups:
    if group == -1:
        continue
    gpos = np.where((groups == group) & (array > min_mz) & (array < max_mz))[0]
    group_info = [g for g in group_list if g[0] == group][0]
    data.append(
        Scatter(
            x = array[gpos],
            y = kde[gpos],
            mode = 'markers',
            name = group_info[2],
        )
    )
plotly.offline.iplot({'data':data})


    

    




42






    

In [511]:

    

# Search for multiple formulas
form_found = {}
form_counts = {}
form_ppm = {}
for word in c.word_counts:
    split_word = word.split('_')
    if len(split_word) == 3:
        formula = Formula(split_word[2])
        newword = split_word[0] + '_' + split_word[2]
        mz = float(split_word[1])
        ppm = 1e6*np.abs(mz - formula.compute_exact_mass())/mz
        if newword in form_found:
            form_found[newword].append(word)
            form_counts[newword] += 1
            form_ppm[newword].append(ppm)
        else:
            form_found[newword] = [word]
            form_counts[newword] = 1
            form_ppm[newword] = [ppm]


    

In [514]:

    

form_counts


    

    
Out[514]:





{'fragment_C10H12NO2': 1,
 'fragment_C10H9O3': 1,
 'fragment_C11H10NO2': 1,
 'fragment_C11H17': 1,
 'fragment_C11H20NO4': 1,
 'fragment_C12H22NO6': 1,
 'fragment_C15H11O': 1,
 'fragment_C17H14O4': 1,
 'fragment_C18H29O2': 1,
 'fragment_C2H3N2': 1,
 'fragment_C2H4NO2': 1,
 'fragment_C2H6NO': 1,
 'fragment_C3H10N': 1,
 'fragment_C3H4NO': 1,
 'fragment_C3H5O': 1,
 'fragment_C3H5O2': 1,
 'fragment_C3H6N': 1,
 'fragment_C3H6NO': 1,
 'fragment_C3H6NO2': 1,
 'fragment_C3H7N2O': 1,
 'fragment_C3H7O': 1,
 'fragment_C3H8N': 2,
 'fragment_C3H8NO': 1,
 'fragment_C3HO2': 1,
 'fragment_C4H10N': 1,
 'fragment_C4H3N2O': 1,
 'fragment_C4H4N3O': 2,
 'fragment_C4H4O': 1,
 'fragment_C4H5': 1,
 'fragment_C4H5O': 1,
 'fragment_C4H5O2': 2,
 'fragment_C4H6N': 1,
 'fragment_C4H6N3O': 1,
 'fragment_C4H6NO': 1,
 'fragment_C4H7N2': 1,
 'fragment_C4H7O2': 1,
 'fragment_C4H8N': 1,
 'fragment_C4H9': 1,
 'fragment_C4H9O': 2,
 'fragment_C4H9O2': 1,
 'fragment_C5H10N': 1,
 'fragment_C5H10NO2': 1,
 'fragment_C5H10O5': 1,
 'fragment_C5H12N': 1,
 'fragment_C5H12N3': 1,
 'fragment_C5H12NO2': 1,
 'fragment_C5H14NO': 1,
 'fragment_C5H3N4': 2,
 'fragment_C5H3N4O': 1,
 'fragment_C5H3O3': 1,
 'fragment_C5H5': 1,
 'fragment_C5H5N2': 1,
 'fragment_C5H5N2O': 2,
 'fragment_C5H5N4O2': 1,
 'fragment_C5H5O': 1,
 'fragment_C5H5O2': 1,
 'fragment_C5H5O4': 1,
 'fragment_C5H6N': 1,
 'fragment_C5H6N3': 2,
 'fragment_C5H6N5': 2,
 'fragment_C5H6N5O': 1,
 'fragment_C5H6NO': 1,
 'fragment_C5H7': 1,
 'fragment_C5H7N2': 1,
 'fragment_C5H7O': 1,
 'fragment_C5H7O3': 1,
 'fragment_C5H8N3': 2,
 'fragment_C5H8N3O': 1,
 'fragment_C5H8NO': 1,
 'fragment_C5H8NO3': 1,
 'fragment_C5H8O4': 1,
 'fragment_C5H9': 1,
 'fragment_C5H9N2O2': 1,
 'fragment_C5H9O': 1,
 'fragment_C6H10N3O2': 1,
 'fragment_C6H10NO2': 1,
 'fragment_C6H10NO3': 1,
 'fragment_C6H12NO': 2,
 'fragment_C6H12NO3': 1,
 'fragment_C6H12O6': 1,
 'fragment_C6H13N4O': 1,
 'fragment_C6H14NO2': 1,
 'fragment_C6H5O2': 1,
 'fragment_C6H5O5': 1,
 'fragment_C6H6NO': 1,
 'fragment_C6H6NO2': 1,
 'fragment_C6H7': 1,
 'fragment_C6H7O': 1,
 'fragment_C6H7O2': 1,
 'fragment_C6H7O3': 1,
 'fragment_C6H9': 1,
 'fragment_C6H9N2O': 1,
 'fragment_C7H10NO2': 1,
 'fragment_C7H10O6': 1,
 'fragment_C7H5NO2': 1,
 'fragment_C7H6NO2': 1,
 'fragment_C7H7': 1,
 'fragment_C7H7O2': 1,
 'fragment_C7H8NO2': 1,
 'fragment_C7H9': 1,
 'fragment_C8H10N': 1,
 'fragment_C8H10NO': 1,
 'fragment_C8H10NO2': 1,
 'fragment_C8H10NO3': 1,
 'fragment_C8H12NO4': 1,
 'fragment_C8H5O': 2,
 'fragment_C8H6NO3': 1,
 'fragment_C8H7': 1,
 'fragment_C8H7N': 1,
 'fragment_C8H7O': 1,
 'fragment_C8H8N': 1,
 'fragment_C8H8NO': 1,
 'fragment_C8H8NO2': 1,
 'fragment_C8H9': 1,
 'fragment_C8H9O': 1,
 'fragment_C8H9O2': 1,
 'fragment_C9H12NO3': 1,
 'fragment_C9H13N2O2': 1,
 'fragment_C9H5O2': 1,
 'fragment_C9H7': 1,
 'fragment_C9H7O2': 1,
 'fragment_C9H8N': 1,
 'fragment_C9H8NO3': 1,
 'fragment_C9H9O2': 1,
 'fragment_CH6N3': 1,
 'fragment_CH6O3': 1,
 'fragment_H4O4P': 1,
 'loss_C10H14NO3': 1,
 'loss_C10H9O3': 1,
 'loss_C2H3N2': 1,
 'loss_C2H4NO2': 1,
 'loss_C2H4O2': 1,
 'loss_C2H4O4': 1,
 'loss_C2H5NO': 1,
 'loss_C2H5NO3': 1,
 'loss_C2H6NO': 1,
 'loss_C2H6O3': 1,
 'loss_C2O2': 1,
 'loss_C3H10N': 1,
 'loss_C3H2O': 1,
 'loss_C3H4NO': 1,
 'loss_C3H4O3': 1,
 'loss_C3H5O': 1,
 'loss_C3H5O2': 1,
 'loss_C3H6N': 1,
 'loss_C3H6NO': 1,
 'loss_C3H6NO2': 1,
 'loss_C3H7N2O': 1,
 'loss_C3H7O': 1,
 'loss_C3H8N': 1,
 'loss_C3H8NO': 1,
 'loss_C3H8NO3': 1,
 'loss_C3H8O4': 1,
 'loss_C3H9N': 1,
 'loss_C4H10N': 1,
 'loss_C4H3N2O': 1,
 'loss_C4H4O': 1,
 'loss_C4H5': 1,
 'loss_C4H5O2': 1,
 'loss_C4H6N': 1,
 'loss_C4H6NO': 1,
 'loss_C4H6NO3': 1,
 'loss_C4H6O2': 1,
 'loss_C4H7': 1,
 'loss_C4H7N2': 1,
 'loss_C4H7O2': 1,
 'loss_C4H8N': 1,
 'loss_C4H9': 1,
 'loss_C4H9N2O3': 1,
 'loss_C4H9O': 1,
 'loss_C4H9O2': 1,
 'loss_C5H10N': 1,
 'loss_C5H10NO2': 1,
 'loss_C5H10O5': 1,
 'loss_C5H11N2O3': 1,
 'loss_C5H11NO4': 1,
 'loss_C5H12N': 1,
 'loss_C5H12NO2': 1,
 'loss_C5H14NO': 1,
 'loss_C5H3O3': 1,
 'loss_C5H5N2': 1,
 'loss_C5H5N4O2': 1,
 'loss_C5H5O': 1,
 'loss_C5H5O4': 1,
 'loss_C5H6N': 1,
 'loss_C5H6N5': 1,
 'loss_C5H6NO': 1,
 'loss_C5H7N2': 1,
 'loss_C5H7O3': 1,
 'loss_C5H8NO': 1,
 'loss_C5H8NO3': 1,
 'loss_C5H8O4': 1,
 'loss_C5H9': 1,
 'loss_C5H9N2O2': 1,
 'loss_C5H9NO2': 1,
 'loss_C5H9O4': 1,
 'loss_C6H10N3O2': 1,
 'loss_C6H10NO2': 1,
 'loss_C6H10NO3': 1,
 'loss_C6H12NO': 1,
 'loss_C6H12NO3': 1,
 'loss_C6H12O6': 1,
 'loss_C6H13N4O': 1,
 'loss_C6H13O3': 1,
 'loss_C6H14NO2': 1,
 'loss_C6H14O7': 1,
 'loss_C6H5O2': 1,
 'loss_C6H5O5': 1,
 'loss_C6H6NO': 1,
 'loss_C6H6NO2': 1,
 'loss_C6H8O4': 1,
 'loss_C6H9N2O': 1,
 'loss_C7H10NO2': 1,
 'loss_C7H10O6': 1,
 'loss_C7H5NO2': 1,
 'loss_C7H6NO2': 1,
 'loss_C7H7': 1,
 'loss_C7H7O2': 1,
 'loss_C7H8NO2': 1,
 'loss_C7H9': 1,
 'loss_C8H10N': 1,
 'loss_C8H10NO': 1,
 'loss_C8H10NO2': 1,
 'loss_C8H10NO3': 1,
 'loss_C8H12NO4': 1,
 'loss_C8H6NO3': 1,
 'loss_C8H7': 1,
 'loss_C8H7N': 1,
 'loss_C8H7O': 1,
 'loss_C8H8N': 1,
 'loss_C8H8NO': 1,
 'loss_C8H8NO2': 1,
 'loss_C8H9': 1,
 'loss_C8H9O': 1,
 'loss_C8H9O2': 1,
 'loss_C9H12NO3': 1,
 'loss_C9H13N2O2': 1,
 'loss_C9H5O2': 1,
 'loss_C9H7': 1,
 'loss_C9H7O2': 1,
 'loss_C9H8N': 1,
 'loss_C9H8NO3': 1,
 'loss_C9H9O2': 1,
 'loss_CH2O2': 1,
 'loss_CH4O3': 1,
 'loss_CH5NO2': 1,
 'loss_CH6O3': 1,
 'loss_CO': 1,
 'loss_H2O': 1,
 'loss_H3N': 1,
 'loss_H4O2': 1,
 'loss_H4O4P': 1}

In [497]:

    

print form_ppm['loss_C4H5O']
print form_found[formula]


    

    




[0.57627835541356298, 19.407944350838701, 13.909135945577521]
['loss_69.034_C4H5O', 'loss_69.0327_C4H5O', 'loss_69.035_C4H5O']

In [512]:

    

print len(form_counts)


    

    




252

In [ ]: