diff --git a/.gitignore b/.gitignore index 36cfc713..10aea2bf 100644 --- a/.gitignore +++ b/.gitignore @@ -106,3 +106,4 @@ venv.bak/ # mypy .mypy_cache/ +physical_property/logD/calc_logD2.nb diff --git a/physical_property/logD/calc_logD.nb b/physical_property/logD/calc_logD.nb index 688d9287..5c39d81d 100644 --- a/physical_property/logD/calc_logD.nb +++ b/physical_property/logD/calc_logD.nb @@ -10,76 +10,42 @@ NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] -NotebookDataLength[ 30423, 777] -NotebookOptionsPosition[ 29114, 749] -NotebookOutlinePosition[ 29550, 766] -CellTagsIndexPosition[ 29507, 763] +NotebookDataLength[ 41761, 958] +NotebookOptionsPosition[ 40559, 932] +NotebookOutlinePosition[ 41101, 951] +CellTagsIndexPosition[ 41058, 948] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ -Cell[BoxData[ +Cell[BoxData[{ RowBox[{ RowBox[{"SetDirectory", "[", - RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]], "Input", - CellChangeTimes->{{3.8115299356460733`*^9, 3.811529955884295*^9}, { - 3.822720374371791*^9, 3.822720408954632*^9}, {3.823066797094819*^9, - 3.8230667997537003`*^9}, {3.8230669395343513`*^9, 3.823066946458387*^9}, { - 3.8235820182085996`*^9, 3.8235820192650223`*^9}}, - CellLabel->"In[2]:=",ExpressionUUID->"e07bf09a-837b-4f9f-8ae3-5348be5ddf17"], - -Cell[BoxData[{ + RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"exp", "=", RowBox[{ RowBox[{ - RowBox[{ - "Import", "[", - "\"\<../experimental_data/Experimental_Properties_of_SAMPL7_Compounds.\ + "Import", "[", + "\"\<../experimental_data/Experimental_Properties_of_SAMPL7_Compounds.\ csv\>\"", "]"}], "[", - RowBox[{"[", - RowBox[{"All", ",", - RowBox[{"{", - RowBox[{ - "1", ",", "5", ",", "6", ",", "7", ",", "8", ",", "13", ",", "14"}], - "}"}]}], "]"}], "]"}], "/.", - RowBox[{"\"\<>12\>\"", "\[Rule]", "\"\\""}]}]}], - ";"}], "\[IndentingNewLine]", - RowBox[{ - RowBox[{"logDexp", "=", - RowBox[{"exp", "[", RowBox[{"[", RowBox[{ RowBox[{"2", ";;"}], ",", RowBox[{"{", - RowBox[{"1", ",", "4", ",", "5"}], "}"}]}], "]"}], "]"}]}], + RowBox[{"1", ",", "7", ",", "8"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{"RT", "=", + RowBox[{"(*", RowBox[{ - RowBox[{"(", - RowBox[{"8.3144598", "*", "298.15"}], ")"}], "/", "4184"}]}], - ";"}], "\[IndentingNewLine]", - RowBox[{ - RowBox[{"MSE", "[", "list_", "]"}], ":=", - RowBox[{ - RowBox[{"Total", "[", "list", "]"}], "/", - RowBox[{"Length", "[", "list", "]"}]}]}], "\[IndentingNewLine]", - RowBox[{ - RowBox[{"MAE", "[", "list_", "]"}], ":=", - RowBox[{ - RowBox[{"Total", "[", - RowBox[{"Abs", "[", "list", "]"}], "]"}], "/", - RowBox[{"Length", "[", "list", "]"}]}]}]}], "Input", - CellChangeTimes->CompressedData[" -1:eJxTTMoPSmViYGAQBmIQfaXFZPar1DeOUY9zwXTuMa3lIPoC4/qVINpEifU4 -iBb9UH4GRGvMlX8NobXA9LUP0WKvgfSKl5aSIPr0rK/eEHp6EJiWSo2F0BvB -9Ka1YcUg2s00pQRE/1uQVg2mFdbXgOhPNyV3nc0A0fkHQHSM7MltTtVvHIVD -QvaA6Or9BYdAtJNNP5gubpwd7QykF7xdAKZt/5XedgXSPzjawPTLogdm7kA6 -7WSfBYiu0Iltewek/VxbwLTFqi3/+WrfOHb9PwSmV7z5KyAJpJWncgiC6Llb -CkVBtJN3PZgGANeCpS4= - "], - CellLabel->"In[3]:=",ExpressionUUID->"276a1f97-6d10-40ce-9758-2b5774cf6a58"], + "Set", " ", "active", " ", "directory", " ", "and", " ", "import", " ", + "experimental", " ", "logD", " ", "data"}], "*)"}]}]}], "Input", + CellChangeTimes->{{3.8115299356460733`*^9, 3.811529955884295*^9}, { + 3.822720374371791*^9, 3.822720408954632*^9}, {3.823066797094819*^9, + 3.8230667997537003`*^9}, {3.8230669395343513`*^9, 3.823066946458387*^9}, { + 3.8235820182085996`*^9, 3.8235820192650223`*^9}, {3.8250565297613287`*^9, + 3.82505655441084*^9}}, + CellLabel->"In[1]:=",ExpressionUUID->"e07bf09a-837b-4f9f-8ae3-5348be5ddf17"], Cell[BoxData[{ RowBox[{ @@ -135,7 +101,7 @@ submission_collection.csv\>\"", "]"}], "[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], "=", "\"\<#\>\""}]}], ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{"pkamodels", "=", + RowBox[{"pkamods", "=", RowBox[{"DeleteDuplicates", "[", RowBox[{"Flatten", "@", RowBox[{"pka", "[", @@ -143,7 +109,7 @@ submission_collection.csv\>\"", "]"}], "[", RowBox[{"All", ",", "All", ",", "3"}], "]"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{"logPmodels", "=", + RowBox[{"logPmods", "=", RowBox[{"DeleteDuplicates", "[", RowBox[{"Flatten", "@", RowBox[{"logP", "[", @@ -153,9 +119,17 @@ submission_collection.csv\>\"", "]"}], "[", RowBox[{ RowBox[{"(*", RowBox[{ - "Load", " ", "all", " ", "ranked", " ", "pKa", " ", "and", " ", "logP", - " ", "submissions", " ", "and", " ", "gather", " ", "method", " ", - "names"}], "*)"}]}]}], "Input", + RowBox[{ + "Load", " ", "all", " ", "ranked", " ", "pKa", " ", "and", " ", "logP", + " ", "submissions"}], ",", " ", + RowBox[{ + "group", " ", "them", " ", "by", " ", "their", " ", + "\"\\""}], ",", " ", + RowBox[{ + "and", " ", "gather", " ", "the", " ", "pka", " ", "and", " ", "logP", + " ", "\"\\"", " ", "entries", " ", "in", " ", "seperate", + " ", + RowBox[{"lists", "."}]}]}], "*)"}]}]}], "Input", CellChangeTimes->{{3.82272085317456*^9, 3.822720919245789*^9}, { 3.822720955332499*^9, 3.8227209576312923`*^9}, {3.822720988382762*^9, 3.822721015997744*^9}, {3.822721208958939*^9, 3.8227212298784103`*^9}, { @@ -167,103 +141,43 @@ submission_collection.csv\>\"", "]"}], "[", 3.823072318110059*^9}, {3.823072655259222*^9, 3.82307265610286*^9}, { 3.823072690324131*^9, 3.823072691099353*^9}, {3.823074787092334*^9, 3.823074787720675*^9}, 3.8230791607452707`*^9, {3.8236631253101797`*^9, - 3.823663148058075*^9}, {3.8236937160189795`*^9, 3.8236938249806705`*^9}}, - CellLabel->"In[8]:=",ExpressionUUID->"8416a644-c7b6-49ed-b9a7-00441e0da565"], - -Cell[CellGroupData[{ - -Cell[BoxData[ - RowBox[{"pka", "[", - RowBox[{"[", - RowBox[{"1", ",", "1"}], "]"}], "]"}]], "Input", - CellChangeTimes->{{3.823694284196485*^9, 3.823694291388811*^9}}, - CellLabel->"In[15]:=",ExpressionUUID->"0764d71c-3fcd-41b9-a11f-fab699d2274d"], - -Cell[BoxData[ - RowBox[{"{", - RowBox[{"\<\"#\"\>", ",", "\<\"method_name\"\>", ",", "\<\"file name\"\>", - ",", "\<\"category\"\>", ",", "\<\"Molecule ID\"\>", - ",", "\<\"pKa (calc)\"\>", ",", "\<\"pKa SEM (calc)\"\>", - ",", "\<\"pKa model uncertainty\"\>"}], "}"}]], "Output", - CellChangeTimes->{{3.8236942858842196`*^9, 3.823694291842765*^9}, - 3.823695077738521*^9, 3.8236962674074764`*^9, 3.823859970333976*^9, - 3.8239377158002024`*^9, 3.823939246130381*^9, 3.823940116726141*^9, - 3.8239454938348675`*^9}, - CellLabel->"Out[15]=",ExpressionUUID->"3b525fa6-532a-451d-8731-8a73b2a3dd81"] -}, Open ]], + 3.823663148058075*^9}, {3.8236937160189795`*^9, 3.8236938249806705`*^9}, { + 3.8250570116931624`*^9, 3.8250570303660994`*^9}, {3.8250570615801506`*^9, + 3.8250571382860146`*^9}, {3.8250577383213468`*^9, 3.825057740022767*^9}}, + CellLabel->"In[3]:=",ExpressionUUID->"8416a644-c7b6-49ed-b9a7-00441e0da565"], Cell[CellGroupData[{ -Cell[BoxData[ - RowBox[{"Flatten", "[", - RowBox[{"Join", "@", - RowBox[{"{", - RowBox[{ - RowBox[{"logP", "[", - RowBox[{"[", - RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", - RowBox[{"pka", "[", - RowBox[{"[", - RowBox[{"1", ",", "1", ",", - RowBox[{"{", - RowBox[{"2", ",", "3", ",", "4", ",", "7", ",", "8"}], "}"}]}], - "]"}], "]"}]}], "}"}]}], "]"}]], "Input", - CellChangeTimes->{{3.823694205602987*^9, 3.8236942994160314`*^9}}, - CellLabel->"In[16]:=",ExpressionUUID->"df622103-1cd4-4446-b960-e460cee42968"], - -Cell[BoxData[ - RowBox[{"{", - RowBox[{"\<\"#\"\>", ",", "\<\"method_name\"\>", ",", "\<\"file name\"\>", - ",", "\<\"category\"\>", ",", "\<\"Molecule ID\"\>", - ",", "\<\"logP (calc)\"\>", ",", "\<\"logP SEM (calc)\"\>", - ",", "\<\"logP model uncertainty\"\>", ",", "\<\"method_name\"\>", - ",", "\<\"file name\"\>", ",", "\<\"category\"\>", - ",", "\<\"pKa SEM (calc)\"\>", ",", "\<\"pKa model uncertainty\"\>"}], - "}"}]], "Output", - CellChangeTimes->{ - 3.8236942096625586`*^9, {3.823694244893876*^9, 3.82369430003522*^9}, - 3.8236950777919188`*^9, 3.823696267447751*^9, 3.8238599703885484`*^9, - 3.823937715853605*^9, 3.8239392461929*^9, 3.823940116763912*^9, - 3.8239454938817596`*^9}, - CellLabel->"Out[16]=",ExpressionUUID->"073cb049-5de2-4f1c-9e28-913cfd9a0973"] -}, Open ]], - Cell[BoxData[{ RowBox[{ - RowBox[{ - RowBox[{"assoc", "=", - RowBox[{"{", - RowBox[{ - RowBox[{"{", - RowBox[{"2", ",", "4"}], "}"}], ",", - RowBox[{"{", - RowBox[{"8", ",", "5"}], "}"}], ",", - RowBox[{"{", - RowBox[{"9", ",", "7"}], "}"}], ",", - RowBox[{"{", - RowBox[{"10", ",", "8"}], "}"}], ",", - RowBox[{"{", - RowBox[{"12", ",", "2"}], "}"}], ",", - RowBox[{"{", - RowBox[{"17", ",", "6"}], "}"}]}], "}"}]}], ";"}], - "\[IndentingNewLine]", - RowBox[{"(*", - RowBox[{ - "Manual", " ", "combination", " ", "of", " ", "logP", " ", "and", " ", - "pKa", " ", "submission", " ", "from", " ", "the", " ", "same", " ", - "\"\\""}], "*)"}]}], "\[IndentingNewLine]", + RowBox[{"assoc", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{"2", ",", "4"}], "}"}], ",", + RowBox[{"{", + RowBox[{"8", ",", "5"}], "}"}], ",", + RowBox[{"{", + RowBox[{"9", ",", "7"}], "}"}], ",", + RowBox[{"{", + RowBox[{"10", ",", "8"}], "}"}], ",", + RowBox[{"{", + RowBox[{"12", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"17", ",", "6"}], "}"}]}], "}"}]}], + ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{"logDmodsfull", "=", + RowBox[{"logDmods", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ - RowBox[{"logPmodels", "[", + RowBox[{"logPmods", "[", RowBox[{"[", RowBox[{"assoc", "[", RowBox[{"[", RowBox[{"ii", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "<>", "\"\\"", "<>", - RowBox[{"pkamodels", "[", + RowBox[{"pkamods", "[", RowBox[{"[", RowBox[{"assoc", "[", RowBox[{"[", @@ -272,31 +186,23 @@ Cell[BoxData[{ RowBox[{"ii", ",", RowBox[{"Length", "@", "assoc"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{"logDmods", "//", "TableForm"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{"logDmods", "=", - RowBox[{"logPmodels", "[", - RowBox[{"[", - RowBox[{"assoc", "[", - RowBox[{"[", - RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], - ";"}], "\[IndentingNewLine]", - RowBox[{ - RowBox[{"logDmods", "=", - RowBox[{"Table", "[", + RowBox[{"(*", + RowBox[{ RowBox[{ - RowBox[{"StringReplace", "[", - RowBox[{ - RowBox[{"logDmods", "[", - RowBox[{"[", "ii", "]"}], "]"}], ",", - RowBox[{"{", - RowBox[{ - RowBox[{"\"\\"", "\[Rule]", "\"\\""}], ",", - RowBox[{"\"\\"", "\[Rule]", "\"\\""}]}], "}"}]}], - "]"}], ",", - RowBox[{"{", - RowBox[{"ii", ",", - RowBox[{"Length", "@", "logDmods"}]}], "}"}]}], "]"}]}], - ";"}]}], "Input", + "Manual", " ", "combination", " ", "of", " ", "logP", " ", "and", " ", + "pKa", " ", "submission", " ", "from", " ", "the", " ", "same", " ", + "\"\\""}], ",", " ", + RowBox[{ + RowBox[{"e", ".", "g", ".", " ", "logP"}], " ", "entry", " ", "2", " ", + "is", " ", "associated", " ", "with", " ", "pKa", " ", "entry", " ", + "4"}], ",", " ", + RowBox[{"8", " ", "with", " ", "5", " ", "and", " ", "so", " ", + RowBox[{"on", ".", " ", "List"}], " ", "the", " ", + RowBox[{"logP", "/", "pKa"}], " ", "model", " ", "combinations", " ", + "here", " ", "as", " ", "a", " ", "sanity", " ", + RowBox[{"check", "."}]}]}], "*)"}]}]}], "Input", CellChangeTimes->{{3.8227229928172483`*^9, 3.822723062870565*^9}, { 3.822729800685411*^9, 3.822729879042808*^9}, {3.822734120023357*^9, 3.8227341202943363`*^9}, {3.822738150202759*^9, 3.822738176533668*^9}, { @@ -305,8 +211,41 @@ Cell[BoxData[{ 3.8230730274766817`*^9, 3.823073036130043*^9}, {3.823074382242488*^9, 3.82307438961628*^9}, {3.823074434976383*^9, 3.823074445626651*^9}, 3.823079158462379*^9, {3.823144654829192*^9, 3.823144656136723*^9}, { - 3.82368752851817*^9, 3.8236875342090707`*^9}}, - CellLabel->"In[17]:=",ExpressionUUID->"8379e00e-31fd-45ee-947a-f7bdd9f9030f"], + 3.82368752851817*^9, 3.8236875342090707`*^9}, {3.8250572677204437`*^9, + 3.8250573166948385`*^9}, {3.825057354982473*^9, 3.8250573764850693`*^9}, { + 3.825057420355569*^9, 3.825057485812785*^9}, {3.825061458999426*^9, + 3.8250614795426702`*^9}},ExpressionUUID->"8379e00e-31fd-45ee-947a-\ +f7bdd9f9030f"], + +Cell[BoxData[ + TagBox[ + TagBox[GridBox[{ + {"\<\"logP_prediction_Iorga_Beckstein_CGenFF/pKa_prediction_Iorga_\ +Beckstein_1\"\>"}, + {"\<\"logP-IEFPCMMST-1/pKa-IEFPCMMST-1\"\>"}, + {"\<\"logP-ECRISM-1/pKa-ECRISM-1\"\>"}, + {"\<\"logp-nhlbi-1/pka-nhlbi-1c\"\>"}, + {"\<\"logP_RodriguezPaluch_SMD_2/pKa_RodriguezPaluch_SMD_1\"\>"}, + {"\<\"logP-EvrimArslan-6/pKa-VA-2-charge-correction\"\>"} + }, + GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.5599999999999999]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], + Column], + Function[BoxForm`e$, + TableForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{{3.825057364814217*^9, 3.825057376783749*^9}, { + 3.8250574223433933`*^9, 3.8250574722635593`*^9}, 3.8250589329926434`*^9, + 3.825060262803254*^9}, + CellLabel-> + "Out[12]//TableForm=",ExpressionUUID->"1205d5b1-89f5-45e5-a9c1-\ +06eb476de438"] +}, Open ]], Cell[BoxData[{ RowBox[{ @@ -378,44 +317,127 @@ Cell[BoxData[{ RowBox[{"Length", "@", "logD"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{ - RowBox[{"logD", "=", - RowBox[{"Table", "[", - RowBox[{ - RowBox[{"Select", "[", + RowBox[{"logD", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Select", "[", + RowBox[{ + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "All"}], "]"}], "]"}], ",", RowBox[{ - RowBox[{"logD", "[", - RowBox[{"[", - RowBox[{"ii", ",", "All"}], "]"}], "]"}], ",", RowBox[{ RowBox[{ - RowBox[{ - RowBox[{"#", "[", - RowBox[{"[", "8", "]"}], "]"}], "!=", "\"\\""}], "&&", - RowBox[{ - RowBox[{"#", "[", - RowBox[{"[", "8", "]"}], "]"}], "!=", "\"\\""}]}], "&"}]}], - "]"}], ",", - RowBox[{"{", - RowBox[{"ii", ",", - RowBox[{"Length", "@", "logD"}]}], "}"}]}], "]"}]}], ";"}], "\n", + RowBox[{"#", "[", + RowBox[{"[", "8", "]"}], "]"}], "!=", "\"\\""}], "&&", + RowBox[{ + RowBox[{"#", "[", + RowBox[{"[", "8", "]"}], "]"}], "!=", "\"\\""}]}], "&"}]}], + "]"}], ",", + RowBox[{"{", + RowBox[{"ii", ",", + RowBox[{"Length", "@", "logD"}]}], "}"}]}], "]"}]}], ";"}], "\n", + RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ - "Skip", " ", "SM28", " ", "and", " ", "SM33", " ", "because", " ", "for", - " ", "these", " ", "molecules", " ", "pKa"}], " ", ">", " ", + "Initiate", " ", "the", " ", "logD", " ", "list", " ", "to", " ", "be", + " ", "exported", " ", "later", " ", "so", " ", "that", " ", "it", " ", + "already", " ", "contains", " ", "the", " ", "pka", " ", "and", " ", + "logP", " ", "data"}], ",", " ", "errors", ",", " ", RowBox[{ - "12", " ", "and", " ", "no", " ", "predicted", " ", "pKa", " ", "is", " ", - RowBox[{"given", ".", " ", "Use"}], " ", "tmp", " ", "to", " ", - "remerge", " ", "SM28", " ", "and", " ", - RowBox[{"SM33", "."}]}]}], "*)"}]}], "\n", + RowBox[{"method", " ", "names", " ", + RowBox[{"etc", ".", "\[IndentingNewLine]", "Skip"}], " ", "SM28", " ", + "and", " ", "SM33", " ", "because", " ", "for", " ", "these", " ", + "molecules", " ", "pKa"}], " ", ">", " ", + RowBox[{ + "12", " ", "and", " ", "no", " ", "predicted", " ", "pKa", " ", "is", + " ", + RowBox[{"given", ".", " ", "Use"}], " ", "the", " ", "list", " ", + "\"\\"", " ", "to", " ", "remerge", " ", "SM28", " ", "and", " ", + "SM33", " ", "entries", " ", "after", " ", "calculating", " ", "logD", + " ", "for", " ", "all", " ", "other", " ", + RowBox[{"compounds", "."}]}]}]}], "*)"}]}]}], "Input", + CellChangeTimes->{{3.813313473421191*^9, 3.813313589514941*^9}, + 3.813313671483691*^9, {3.813313853358691*^9, 3.813313869561816*^9}, { + 3.813313905468066*^9, 3.813313908952441*^9}, {3.813313940671191*^9, + 3.813313943264941*^9}, {3.822721346762309*^9, 3.822721347792636*^9}, { + 3.822721766034278*^9, 3.822721936410342*^9}, {3.82272204662801*^9, + 3.822722089131222*^9}, {3.8227221384138412`*^9, 3.822722175380028*^9}, { + 3.8227231080243397`*^9, 3.822723267841703*^9}, {3.822723544460298*^9, + 3.8227235565425377`*^9}, {3.8227235938108397`*^9, 3.8227236757791*^9}, { + 3.822724012182016*^9, 3.822724017467703*^9}, {3.822725423161537*^9, + 3.82272545924118*^9}, {3.822725508763267*^9, 3.8227255165931873`*^9}, { + 3.822725556074925*^9, 3.822725595568293*^9}, {3.8227259171263857`*^9, + 3.822725929141288*^9}, {3.822726603845252*^9, 3.822726639092874*^9}, { + 3.8227266695655107`*^9, 3.8227266833945723`*^9}, {3.822726725335083*^9, + 3.8227267318018312`*^9}, {3.822729944451618*^9, 3.8227301142522507`*^9}, { + 3.822730192548806*^9, 3.822730291181621*^9}, {3.8227335102063026`*^9, + 3.822733515000246*^9}, {3.822735544535851*^9, 3.822735581979568*^9}, { + 3.822735641213106*^9, 3.822735646052837*^9}, {3.822736709859652*^9, + 3.822736782259831*^9}, {3.8230786544705687`*^9, 3.823078676471931*^9}, { + 3.82307874767752*^9, 3.823078749413064*^9}, {3.823678121353013*^9, + 3.823678125399952*^9}, {3.8236782247689877`*^9, 3.8236783322398105`*^9}, { + 3.8236796365113897`*^9, 3.823679651769524*^9}, {3.823679779120246*^9, + 3.8236797964284253`*^9}, {3.8236798991641235`*^9, 3.823679975381322*^9}, { + 3.823680011512501*^9, 3.8236800403628726`*^9}, 3.823680083640499*^9, { + 3.8236802054423437`*^9, 3.8236802148253746`*^9}, 3.8236803474147067`*^9, { + 3.8236803929106684`*^9, 3.8236803930337095`*^9}, 3.8236807122253532`*^9, { + 3.8236809045583506`*^9, 3.823680947906972*^9}, {3.823681278161251*^9, + 3.8236813472678576`*^9}, 3.823683782587288*^9, {3.8236838676414857`*^9, + 3.823683908554928*^9}, {3.8236839387300186`*^9, 3.8236840904295044`*^9}, { + 3.823684349393591*^9, 3.8236844171909504`*^9}, {3.8236852431813107`*^9, + 3.8236852633212123`*^9}, {3.823687499923328*^9, 3.823687511297411*^9}, { + 3.8236941703738494`*^9, 3.823694194440398*^9}, {3.8236943172685785`*^9, + 3.823694431892218*^9}, {3.823696254253566*^9, 3.8236962543160844`*^9}, + 3.8238600363712325`*^9, 3.8238691492203074`*^9, {3.8238692459658766`*^9, + 3.8238692932472935`*^9}, {3.8238710891500287`*^9, + 3.8238711149814725`*^9}, {3.8238712254815006`*^9, + 3.8238712286898937`*^9}, {3.8239370851849346`*^9, 3.823937101705805*^9}, { + 3.8239376673934155`*^9, 3.823937669917285*^9}, {3.8239401003548517`*^9, + 3.8239401123256245`*^9}, {3.8239415383293805`*^9, + 3.8239415413394456`*^9}, {3.825057597050869*^9, 3.8250576724434853`*^9}, { + 3.8250577829726515`*^9, 3.825057784093907*^9}, {3.8250578468296676`*^9, + 3.8250578568350163`*^9}, 3.8250579676140537`*^9, {3.825058059146756*^9, + 3.8250580849243584`*^9}, {3.8250586839387293`*^9, 3.825058691180912*^9}}, + CellLabel->"In[13]:=",ExpressionUUID->"9f3db9fc-4bd0-4901-a278-800e3ef14180"], + +Cell[BoxData[{ RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", - RowBox[{"True", ",", "\[IndentingNewLine]", + RowBox[{"True", ",", + RowBox[{"(*", + RowBox[{ + RowBox[{"Only", " ", "this", " ", "branch", " ", "is", " ", + RowBox[{"evaluated", ".", " ", "Calculates"}], " ", "the", " ", + "log", " ", "D7", ".4", " ", "according", " ", "to", " ", + RowBox[{"e", ".", "g", ".", " ", "Scherrer"}], " ", "RA"}], ",", + " ", + RowBox[{ + RowBox[{"Howard", " ", "SM", " ", + RowBox[{"(", "1977", ")"}], " ", "J", " ", "Med", " ", "Chem", + " ", "20"}], ":", + RowBox[{"53", "-", "58"}]}], ",", " ", + RowBox[{ + RowBox[{ + "and", " ", "the", " ", "combined", " ", "SEM", " ", "and", " ", + "model", " ", + RowBox[{"uncertainties", ".", "\[IndentingNewLine]", "logD"}]}], + "=", + RowBox[{"logP", "+", + RowBox[{"Log10", + RowBox[{"(", + RowBox[{"1", "/", + RowBox[{"(", + RowBox[{"1", "+", + RowBox[{"10", "^", + RowBox[{"(", + RowBox[{"7.4", "-", "pKa"}], ")"}]}]}], ")"}]}]}]}]}]}]}], + "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"logD", "[", @@ -522,6 +544,16 @@ Cell[BoxData[{ RowBox[{"[", RowBox[{"ii", ",", "jj", ",", "13"}], "]"}], "]"}]}]}]}], ";"}], ",", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + RowBox[{ + "This", " ", "branch", " ", "is", " ", "never", " ", "evaluated"}], + ",", " ", + RowBox[{ + "but", " ", "would", " ", "do", " ", "the", " ", "same", " ", "as", + " ", "the", " ", "other", " ", "branch", " ", "only", " ", "for", + " ", + RowBox[{"\"\\"", "."}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"logD", "[", @@ -547,8 +579,89 @@ Cell[BoxData[{ RowBox[{"[", RowBox[{"ii", ",", "2"}], "]"}], "]"}], ",", "jj", ",", "6"}], "]"}], "]"}], "-", "7.4"}], ")"}]}]}], ")"}]}], - "]"}]}], ",", "0.01"}], "]"}]}], ";"}]}], "\[IndentingNewLine]", - "]"}], "\[IndentingNewLine]", ",", + "]"}]}], ",", "0.01"}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "10"}], "]"}], "]"}], "=", + RowBox[{ + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "10"}], "]"}], "]"}], "+", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{"-", + RowBox[{"(", + RowBox[{"10", "^", + RowBox[{"(", + RowBox[{ + RowBox[{"pka", "[", + RowBox[{"[", + RowBox[{ + RowBox[{"assoc", "[", + RowBox[{"[", + RowBox[{"ii", ",", "2"}], "]"}], "]"}], ",", "jj", ",", + "6"}], "]"}], "]"}], "-", "7.4"}], ")"}]}], ")"}]}], "/", + RowBox[{"(", + RowBox[{"1", "+", + RowBox[{"10", "^", + RowBox[{"(", + RowBox[{ + RowBox[{"pka", "[", + RowBox[{"[", + RowBox[{ + RowBox[{"assoc", "[", + RowBox[{"[", + RowBox[{"ii", ",", "2"}], "]"}], "]"}], ",", "jj", ",", + "6"}], "]"}], "]"}], "-", "7.4"}], ")"}]}]}], ")"}]}], + ")"}], "^", "2"}], "*", + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "11"}], "]"}], "]"}]}]}]}], ";", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "12"}], "]"}], "]"}], "=", + RowBox[{ + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "12"}], "]"}], "]"}], "+", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{"-", + RowBox[{"(", + RowBox[{"10", "^", + RowBox[{"(", + RowBox[{ + RowBox[{"pka", "[", + RowBox[{"[", + RowBox[{ + RowBox[{"assoc", "[", + RowBox[{"[", + RowBox[{"ii", ",", "2"}], "]"}], "]"}], ",", "jj", ",", + "6"}], "]"}], "]"}], "-", "7.4"}], ")"}]}], ")"}]}], "/", + RowBox[{"(", + RowBox[{"1", "+", + RowBox[{"10", "^", + RowBox[{"(", + RowBox[{ + RowBox[{"pka", "[", + RowBox[{"[", + RowBox[{ + RowBox[{"assoc", "[", + RowBox[{"[", + RowBox[{"ii", ",", "2"}], "]"}], "]"}], ",", "jj", ",", + "6"}], "]"}], "]"}], "-", "7.4"}], ")"}]}]}], ")"}]}], + ")"}], "^", "2"}], "*", + RowBox[{"logD", "[", + RowBox[{"[", + RowBox[{"ii", ",", "jj", ",", "13"}], "]"}], "]"}]}]}]}], + ";"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"jj", ",", RowBox[{"Length", "@", @@ -574,46 +687,95 @@ Cell[BoxData[{ RowBox[{"Length", "@", "logD"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ - RowBox[{ - RowBox[{"Clear", "[", "tmp", "]"}], ";"}], "\[IndentingNewLine]", + RowBox[{"Clear", "[", "tmp", "]"}], ";"}], "\[IndentingNewLine]", + RowBox[{ RowBox[{"(*", RowBox[{ - RowBox[{"log", " ", "D7", ".4", " ", "according", " ", "to", " ", - RowBox[{"e", ".", "g", ".", " ", "Scherrer"}], " ", "RA"}], ",", " ", RowBox[{ - RowBox[{"Howard", " ", "SM", " ", - RowBox[{"(", "1977", ")"}], " ", "J", " ", "Med", " ", "Chem", " ", - "20"}], ":", - RowBox[{"53", "-", + "Technically", " ", "a", " ", "list", " ", "of", " ", "the", " ", + "ionization", " ", "states", " ", "for", " ", "the", " ", "most", " ", + "relevant", " ", "pKa", " ", "transitions", " ", + RowBox[{"(", RowBox[{ - "58.", "\[IndentingNewLine]", " ", "Technically", " ", "a", " ", "list", - " ", "of", " ", "the", " ", "ionization", " ", "states", " ", "for", - " ", "the", " ", "most", " ", "relevant", " ", "pKa", " ", - "transitions", " ", - RowBox[{"(", + RowBox[{"i", ".", "e", ".", " ", "0"}], "\[Rule]", + RowBox[{ RowBox[{ - RowBox[{"i", ".", "e", ".", " ", "0"}], "\[Rule]", - RowBox[{ - RowBox[{ - RowBox[{"-", "1"}], " ", "or", " ", "1"}], "\[Rule]", - RowBox[{"0", " ", - RowBox[{"etc", "."}]}]}]}], ")"}], " ", "for", " ", "each", " ", - "submission", " ", "is", " ", "needed", " ", "to", " ", - "differentiate", " ", "\"\\"", " ", "and", " ", - "\"\\""}]}]}], ",", " ", + RowBox[{"-", "1"}], " ", "or", " ", "1"}], "\[Rule]", + RowBox[{"0", " ", + RowBox[{"etc", "."}]}]}]}], ")"}], " ", "for", " ", "each", " ", + "compound", " ", "and", " ", "each", " ", "submission", " ", "is", " ", + "needed", " ", "to", " ", "differentiate", " ", "\"\\"", " ", + "and", " ", "\"\\""}], ",", " ", RowBox[{ + "especially", " ", "if", " ", "there", " ", "is", " ", "disagreement", + " ", "between", " ", "different", " ", + RowBox[{"methods", ".", " ", "Submissions"}], " ", "that", " ", + "included", " ", "macrostate", " ", "pKa", " ", "values", " ", "agreed", + " ", "that", " ", "every", " ", "relevant", " ", "transition", " ", "is", + " ", "\"\<0\[Rule]-1\>\"", " ", "so", " ", "only", " ", "this", " ", + "equation", " ", "is", " ", "used", " ", "here", " ", RowBox[{ - "especially", " ", "if", " ", "there", " ", "is", " ", "disagreement", - " ", "between", " ", "different", " ", - RowBox[{"methods", ".", "\[IndentingNewLine]", "Submissions"}], " ", - "that", " ", "included", " ", "macrostate", " ", "pKa", " ", "values", - " ", "agreed", " ", "that", " ", "every", " ", "relevant", " ", - "transition", " ", "is", " ", "0"}], "\[Rule]", - RowBox[{"-", "1"}]}], ",", " ", - RowBox[{ - "so", " ", "only", " ", "this", " ", "equation", " ", "is", " ", "used", - " ", - RowBox[{"here", "."}]}]}], "*)"}]}], "\[IndentingNewLine]", + RowBox[{"(", + RowBox[{"If", "[", + RowBox[{"True", ",", "..."}], "]"}], ")"}], ".", + "\[IndentingNewLine]", "After"}], " ", "calculation", " ", "of", " ", + "the", " ", "logD", " ", "values", " ", "the", " ", "\"\\"", " ", + "list", " ", "containing", " ", "the", " ", "logP", " ", "data", " ", + "for", " ", "SM28", " ", "and", " ", "SM33", " ", "is", " ", "merged", + " ", "with", " ", "the", " ", "logD", " ", "list", " ", + RowBox[{"again", "."}]}]}], "*)"}]}]}], "Input", + CellChangeTimes->{{3.813313473421191*^9, 3.813313589514941*^9}, + 3.813313671483691*^9, {3.813313853358691*^9, 3.813313869561816*^9}, { + 3.813313905468066*^9, 3.813313908952441*^9}, {3.813313940671191*^9, + 3.813313943264941*^9}, {3.822721346762309*^9, 3.822721347792636*^9}, { + 3.822721766034278*^9, 3.822721936410342*^9}, {3.82272204662801*^9, + 3.822722089131222*^9}, {3.8227221384138412`*^9, 3.822722175380028*^9}, { + 3.8227231080243397`*^9, 3.822723267841703*^9}, {3.822723544460298*^9, + 3.8227235565425377`*^9}, {3.8227235938108397`*^9, 3.8227236757791*^9}, { + 3.822724012182016*^9, 3.822724017467703*^9}, {3.822725423161537*^9, + 3.82272545924118*^9}, {3.822725508763267*^9, 3.8227255165931873`*^9}, { + 3.822725556074925*^9, 3.822725595568293*^9}, {3.8227259171263857`*^9, + 3.822725929141288*^9}, {3.822726603845252*^9, 3.822726639092874*^9}, { + 3.8227266695655107`*^9, 3.8227266833945723`*^9}, {3.822726725335083*^9, + 3.8227267318018312`*^9}, {3.822729944451618*^9, 3.8227301142522507`*^9}, { + 3.822730192548806*^9, 3.822730291181621*^9}, {3.8227335102063026`*^9, + 3.822733515000246*^9}, {3.822735544535851*^9, 3.822735581979568*^9}, { + 3.822735641213106*^9, 3.822735646052837*^9}, {3.822736709859652*^9, + 3.822736782259831*^9}, {3.8230786544705687`*^9, 3.823078676471931*^9}, { + 3.82307874767752*^9, 3.823078749413064*^9}, {3.823678121353013*^9, + 3.823678125399952*^9}, {3.8236782247689877`*^9, 3.8236783322398105`*^9}, { + 3.8236796365113897`*^9, 3.823679651769524*^9}, {3.823679779120246*^9, + 3.8236797964284253`*^9}, {3.8236798991641235`*^9, 3.823679975381322*^9}, { + 3.823680011512501*^9, 3.8236800403628726`*^9}, 3.823680083640499*^9, { + 3.8236802054423437`*^9, 3.8236802148253746`*^9}, 3.8236803474147067`*^9, { + 3.8236803929106684`*^9, 3.8236803930337095`*^9}, 3.8236807122253532`*^9, { + 3.8236809045583506`*^9, 3.823680947906972*^9}, {3.823681278161251*^9, + 3.8236813472678576`*^9}, 3.823683782587288*^9, {3.8236838676414857`*^9, + 3.823683908554928*^9}, {3.8236839387300186`*^9, 3.8236840904295044`*^9}, { + 3.823684349393591*^9, 3.8236844171909504`*^9}, {3.8236852431813107`*^9, + 3.8236852633212123`*^9}, {3.823687499923328*^9, 3.823687511297411*^9}, { + 3.8236941703738494`*^9, 3.823694194440398*^9}, {3.8236943172685785`*^9, + 3.823694431892218*^9}, {3.823696254253566*^9, 3.8236962543160844`*^9}, + 3.8238600363712325`*^9, 3.8238691492203074`*^9, {3.8238692459658766`*^9, + 3.8238692932472935`*^9}, {3.8238710891500287`*^9, + 3.8238711149814725`*^9}, {3.8238712254815006`*^9, + 3.8238712286898937`*^9}, {3.8239370851849346`*^9, 3.823937101705805*^9}, { + 3.8239376673934155`*^9, 3.823937669917285*^9}, {3.8239401003548517`*^9, + 3.8239401123256245`*^9}, {3.8239415383293805`*^9, + 3.8239415413394456`*^9}, {3.825057597050869*^9, 3.8250576724434853`*^9}, { + 3.8250577829726515`*^9, 3.825057784093907*^9}, {3.8250578468296676`*^9, + 3.8250578568350163`*^9}, 3.8250579676140537`*^9, {3.825058133521*^9, + 3.8250581727535563`*^9}, {3.825058214540615*^9, 3.8250582835920258`*^9}, { + 3.825058473490491*^9, 3.8250586033852425`*^9}, {3.8250586762239475`*^9, + 3.8250587265442257`*^9}, {3.8250588359787335`*^9, + 3.8250588568068843`*^9}, {3.8250589120544124`*^9, 3.825059012864532*^9}, { + 3.8250591321885853`*^9, 3.825059141005411*^9}, {3.825061624401218*^9, + 3.8250617547800684`*^9}, {3.8250636546437287`*^9, + 3.8250636672786117`*^9}, {3.825063698811824*^9, + 3.8250637007755604`*^9}},ExpressionUUID->"8dcda01e-2c9e-4dc5-8077-\ +5ebe246709ad"], + +Cell[BoxData[{ RowBox[{ RowBox[{"logD", "=", RowBox[{"Partition", "[", @@ -641,7 +803,7 @@ Cell[BoxData[{ RowBox[{"Flatten", "[", RowBox[{"Table", "[", RowBox[{ - RowBox[{"logDexp", "[", + RowBox[{"exp", "[", RowBox[{"[", RowBox[{"All", ",", RowBox[{"{", @@ -675,7 +837,17 @@ Cell[BoxData[{ "12"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", - RowBox[{"Prepare", " ", "logD", " ", "Table", " ", "for", " ", "export"}], + RowBox[{ + RowBox[{"Prepare", " ", "logD", " ", "Table", " ", "for", " ", + RowBox[{"export", ".", " ", "Experimental"}], " ", "data", " ", "is", + " ", "added", " ", "to", " ", "the", " ", "list"}], ",", " ", + RowBox[{ + "the", " ", "difference", " ", "between", " ", "predicted", " ", "and", + " ", "experimental", " ", "values", " ", "calculated"}], ",", " ", + RowBox[{ + "and", " ", "the", " ", "list", " ", "sorted", " ", "to", " ", "have", + " ", "the", " ", "same", " ", "shape", " ", "as", " ", "the", " ", + "other", " ", "\"\\"", " ", "files"}]}], "*)"}]}]}], "Input", CellChangeTimes->{{3.813313473421191*^9, 3.813313589514941*^9}, 3.813313671483691*^9, {3.813313853358691*^9, 3.813313869561816*^9}, { @@ -714,8 +886,18 @@ Cell[BoxData[{ 3.8238711149814725`*^9}, {3.8238712254815006`*^9, 3.8238712286898937`*^9}, {3.8239370851849346`*^9, 3.823937101705805*^9}, { 3.8239376673934155`*^9, 3.823937669917285*^9}, {3.8239401003548517`*^9, - 3.8239401123256245`*^9}, {3.8239415383293805`*^9, 3.8239415413394456`*^9}}, - CellLabel->"In[21]:=",ExpressionUUID->"9f3db9fc-4bd0-4901-a278-800e3ef14180"], + 3.8239401123256245`*^9}, {3.8239415383293805`*^9, + 3.8239415413394456`*^9}, {3.825057597050869*^9, 3.8250576724434853`*^9}, { + 3.8250577829726515`*^9, 3.825057784093907*^9}, {3.8250578468296676`*^9, + 3.8250578568350163`*^9}, 3.8250579676140537`*^9, {3.825058133521*^9, + 3.8250581727535563`*^9}, {3.825058214540615*^9, 3.8250582835920258`*^9}, { + 3.825058473490491*^9, 3.8250586033852425`*^9}, {3.8250586762239475`*^9, + 3.8250587265442257`*^9}, {3.8250588359787335`*^9, + 3.8250588568068843`*^9}, {3.8250589120544124`*^9, 3.825059012864532*^9}, { + 3.8250591321885853`*^9, 3.825059141005411*^9}, {3.8250599290648427`*^9, + 3.8250600019288597`*^9}, {3.825060669366186*^9, + 3.825060676597318*^9}},ExpressionUUID->"b8094c35-5837-436f-8f85-\ +46a68c4f4313"], Cell[BoxData[ RowBox[{ @@ -744,11 +926,14 @@ Cell[BoxData[ 3.823694487436466*^9, 3.8236945549583187`*^9}, {3.823695019671069*^9, 3.8236950204267473`*^9}, {3.823937660332825*^9, 3.8239376630256453`*^9}, { 3.8239380771485143`*^9, 3.8239380793477545`*^9}, 3.8239392401103687`*^9, { - 3.8239454404906864`*^9, 3.8239454533620615`*^9}}, - CellLabel->"In[32]:=",ExpressionUUID->"a980ffde-5f1c-41eb-859a-645f6f681bab"] + 3.8239454404906864`*^9, 3.8239454533620615`*^9}, {3.825058923544533*^9, + 3.825058923791282*^9}, {3.8250600067264957`*^9, 3.8250600068515368`*^9}}, + CellLabel->"In[24]:=",ExpressionUUID->"a980ffde-5f1c-41eb-859a-645f6f681bab"] }, -WindowSize->{1140.6, 574.8}, +WindowSize->{1141.2, 580.1999999999999}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, +TaggingRules->{ + "WelcomeScreenSettings" -> {"FEStarting" -> False}, "TryRealOnly" -> False}, Magnification:>0.75 Inherited, FrontEndVersion->"12.2 for Microsoft Windows (64-bit) (December 12, 2020)", StyleDefinitions->"Default.nb", @@ -765,20 +950,16 @@ CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ -Cell[558, 20, 470, 8, 21, "Input",ExpressionUUID->"e07bf09a-837b-4f9f-8ae3-5348be5ddf17"], -Cell[1031, 30, 1738, 51, 81, "Input",ExpressionUUID->"276a1f97-6d10-40ce-9758-2b5774cf6a58"], -Cell[2772, 83, 3306, 87, 123, "Input",ExpressionUUID->"8416a644-c7b6-49ed-b9a7-00441e0da565"], -Cell[CellGroupData[{ -Cell[6103, 174, 248, 5, 21, "Input",ExpressionUUID->"0764d71c-3fcd-41b9-a11f-fab699d2274d"], -Cell[6354, 181, 605, 10, 24, "Output",ExpressionUUID->"3b525fa6-532a-451d-8731-8a73b2a3dd81"] -}, Open ]], +Cell[558, 20, 1066, 27, 51, "Input",ExpressionUUID->"e07bf09a-837b-4f9f-8ae3-5348be5ddf17"], +Cell[1627, 49, 3745, 97, 139, "Input",ExpressionUUID->"8416a644-c7b6-49ed-b9a7-00441e0da565"], Cell[CellGroupData[{ -Cell[6996, 196, 566, 15, 21, "Input",ExpressionUUID->"df622103-1cd4-4446-b960-e460cee42968"], -Cell[7565, 213, 789, 14, 24, "Output",ExpressionUUID->"073cb049-5de2-4f1c-9e28-913cfd9a0973"] +Cell[5397, 150, 2828, 67, 79, "Input",ExpressionUUID->"8379e00e-31fd-45ee-947a-f7bdd9f9030f"], +Cell[8228, 219, 1063, 27, 88, "Output",ExpressionUUID->"1205d5b1-89f5-45e5-a9c1-06eb476de438"] }, Open ]], -Cell[8369, 230, 2881, 78, 79, "Input",ExpressionUUID->"8379e00e-31fd-45ee-947a-f7bdd9f9030f"], -Cell[11253, 310, 16219, 407, 364, "Input",ExpressionUUID->"9f3db9fc-4bd0-4901-a278-800e3ef14180"], -Cell[27475, 719, 1635, 28, 109, "Input",ExpressionUUID->"a980ffde-5f1c-41eb-859a-645f6f681bab"] +Cell[9306, 249, 7140, 153, 93, "Input",ExpressionUUID->"9f3db9fc-4bd0-4901-a278-800e3ef14180"], +Cell[16449, 404, 16316, 371, 264, "Input",ExpressionUUID->"8dcda01e-2c9e-4dc5-8077-5ebe246709ad"], +Cell[32768, 777, 6049, 122, 108, "Input",ExpressionUUID->"b8094c35-5837-436f-8f85-46a68c4f4313"], +Cell[38820, 901, 1735, 29, 52, "Input",ExpressionUUID->"a980ffde-5f1c-41eb-859a-645f6f681bab"] } ] *)