(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 74168, 1923] NotebookOptionsPosition[ 70841, 1855] NotebookOutlinePosition[ 71208, 1871] CellTagsIndexPosition[ 71165, 1868] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Mahalanobis Indexes for optimal Mahalanobis partition", FontColor->RGBColor[1, 0, 0]], "\n", StyleBox["Loads K partitions and for each of them calculatse indices", FontSize->14, FontColor->RGBColor[1, 0.5, 0]] }], "Section", CellChangeTimes->{{3.59226736583011*^9, 3.5922674485562344`*^9}, { 3.5923766059379897`*^9, 3.5923766601948853`*^9}, {3.593700675035776*^9, 3.5937007054628215`*^9}, {3.5960950676401634`*^9, 3.5960950684401646`*^9}, { 3.596095142447279*^9, 3.59609517020732*^9}, {3.5960952227244015`*^9, 3.5960952394224253`*^9}, {3.8067567653335266`*^9, 3.8067568218527255`*^9}, { 3.8222227166970606`*^9, 3.822222717391858*^9}},ExpressionUUID->"8f9af442-96cd-4772-bff6-\ 09a71bc45d4e"], Cell[BoxData[ RowBox[{"(*", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Input", ":", " ", RowBox[{ "set", " ", "of", " ", "K", " ", "partitions", " ", "obtained", " ", "using", " ", "the", " ", RowBox[{"MInc", "[", "]"}], " ", "module", " ", RowBox[{"(", RowBox[{ "each", " ", "partition", " ", "consists", " ", "of", " ", "clusters", " ", "with", " ", "corresponding", " ", "weights", " ", "associated", " ", "to", " ", "their", " ", "elements"}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Output", ":", " ", RowBox[{"MCH", "[", "s", "]"}]}]}], ",", " ", RowBox[{"MDB", "[", "s", "]"}], ",", " ", RowBox[{"MSSWC", "[", "s", "]"}], ",", " ", RowBox[{ RowBox[{"and", " ", RowBox[{"MArea", "[", "s", "]"}], " ", "indexes", " ", "for", " ", "each", "\n", "s"}], "=", "2"}], ",", "...", ",", "K", ",", " ", RowBox[{ RowBox[{"with", " ", "graphs"}], ";", "\[IndentingNewLine]", RowBox[{"Submodules", ":", " ", RowBox[{"dM", "[", "]"}]}]}], ",", " ", RowBox[{"MWF", "[", "]"}], ",", " ", RowBox[{"Hdist", "[", "]"}], ",", " ", RowBox[{"DIRECT", "[", "]"}], ",", " ", RowBox[{"MWKMeans", "[", "]"}]}], "\[IndentingNewLine]", "*)"}]], "Input", CellChangeTimes->{{3.8222227648237658`*^9, 3.8222228783140454`*^9}},ExpressionUUID->"b3a8321d-58d6-4aff-95d5-\ 5ac599347900"] }, Open ]], Cell[CellGroupData[{ Cell["Modules", "Section", CellChangeTimes->{{3.592267492965301*^9, 3.5922674984993086`*^9}},ExpressionUUID->"4ab43163-a9ad-4579-9e58-\ 03db053488e5"], Cell[CellGroupData[{ Cell["Squared Mahalanobis distance", "Subsection", CellChangeTimes->{ 3.5854096427658014`*^9, {3.585636902492603*^9, 3.5856369150848255`*^9}},ExpressionUUID->"dafb582f-fc59-4d24-b27f-\ fba8dcac4aec"], Cell[BoxData[{ RowBox[{ RowBox[{"n", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dM", "[", RowBox[{"xpod_", ",", "ypod_", ",", "S_"}], "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Det", "[", "S", "]"}], "^", RowBox[{"(", RowBox[{"1", "/", "n"}], ")"}]}], ")"}], RowBox[{ RowBox[{"(", RowBox[{"xpod", "-", "ypod"}], ")"}], ".", RowBox[{"Inverse", "[", "S", "]"}], ".", RowBox[{"(", RowBox[{"xpod", "-", "ypod"}], ")"}]}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5866502699378433`*^9, 3.586650277895857*^9}, 3.593699610877139*^9, {3.80675759307401*^9, 3.8067575951690593`*^9}}, CellLabel->"In[71]:=",ExpressionUUID->"31c58d2a-fd13-4f58-85c9-49292444865a"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis objective function", "Subsection", CellChangeTimes->{3.5854097085807176`*^9, 3.585637218737363*^9},ExpressionUUID->"d36b648b-ba45-48dc-940b-\ 54c26491d94c"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"MWF", "[", RowBox[{"PI_", ",", "c_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"m", "=", RowBox[{"Length", "[", RowBox[{"Flatten", "[", RowBox[{"PI", ",", "1"}], "]"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "1"}], "]"}], "]"}], "*", RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", RowBox[{"c", "[", RowBox[{"[", "j", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}], " ", ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], " ", "]"}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5886541897699375`*^9, 3.5886541909877396`*^9}, { 3.588654263023868*^9, 3.588654339794004*^9}, {3.5917588683674088`*^9, 3.591758869047409*^9}, {3.5917590328866577`*^9, 3.591759044816675*^9}, { 3.5923602629384975`*^9, 3.592360265798501*^9}, {3.593700553900591*^9, 3.593700561724603*^9}, 3.593700607408673*^9, 3.5937485374911056`*^9, { 3.8067573402719617`*^9, 3.806757350348838*^9}, {3.8067574110379286`*^9, 3.8067574124867525`*^9}, {3.810879858036764*^9, 3.8108798803874083`*^9}}, CellLabel->"In[73]:=",ExpressionUUID->"73251baf-9953-4113-bb50-712ccc292076"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis Calinski - Harabasz index (F je Area)", "Subsection", CellChangeTimes->{ 3.586745880982668*^9, {3.5886570380848265`*^9, 3.588657043467836*^9}, { 3.592369736897122*^9, 3.592369766942775*^9}, {3.806757572271989*^9, 3.8067575753312635`*^9}},ExpressionUUID->"8ff1266c-46b6-4567-91e8-\ 06d92d292be2"], Cell[BoxData[{ RowBox[{"Clear", "[", "MCH", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"MCH", "[", RowBox[{"PI_", ",", "z_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"k", "=", RowBox[{"Length", "[", "z", "]"}]}], ",", "F", ",", "G"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"m", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"F", "=", RowBox[{"MWF", "[", RowBox[{"PI", ",", "z", ",", "Scov"}], "]"}]}], ";", " ", RowBox[{"(*", " ", "or", " ", "*)"}], RowBox[{"F", "=", RowBox[{"n", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Total", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "1"}], "]"}], "]"}], "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"Det", "[", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], ")"}], "^", RowBox[{"(", RowBox[{"1", "/", "n"}], ")"}]}]}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"G", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "*", RowBox[{"dM", "[", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}], ",", "cc", ",", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{"m", "-", "k"}], ")"}], "*", RowBox[{"G", "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"k", "-", "1"}], ")"}], "F"}], ")"}]}]}]}]}], "\[IndentingNewLine]", " ", "]"}]}]}], "Input", CellChangeTimes->{{3.5867458762860594`*^9, 3.5867458855856757`*^9}, { 3.58865692706863*^9, 3.5886569907453427`*^9}, {3.5886571718280644`*^9, 3.588657180066879*^9}, {3.588659074371236*^9, 3.5886590748548365`*^9}, { 3.592277510524433*^9, 3.5922775159074416`*^9}, 3.592277631556615*^9, { 3.5923633359004774`*^9, 3.592363351765705*^9}, {3.5923634577990913`*^9, 3.5923634812927327`*^9}, 3.592363547343249*^9, {3.5923638848186417`*^9, 3.5923638963782616`*^9}, {3.592363928420718*^9, 3.5923639874200215`*^9}, { 3.592364494889913*^9, 3.5923645241711645`*^9}, {3.5923664586379623`*^9, 3.5923664759705925`*^9}, {3.5943018007507334`*^9, 3.594301802370736*^9}, { 3.806757756441392*^9, 3.8067577616430597`*^9}, 3.8108799029465637`*^9, { 3.822221783468341*^9, 3.8222217912172337`*^9}}, CellLabel->"In[74]:=",ExpressionUUID->"a57c4cc1-6972-4d8a-919e-3ccc245a32c1"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis Davies - Bouldin index", "Subsection", CellChangeTimes->{ 3.5923719354165845`*^9, {3.806757604335451*^9, 3.8067576048899994`*^9}, { 3.806757646861176*^9, 3.8067576476071634`*^9}},ExpressionUUID->"a78eb0ac-1a20-48c3-ad54-\ 3de3bb91e9a6"], Cell[BoxData[{ RowBox[{"Clear", "[", "MDB", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"MDB", "[", RowBox[{"PI_", ",", "z_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"k", "=", RowBox[{"Length", "[", "z", "]"}]}], ",", "FF", ",", "VV", ",", "mj", ",", RowBox[{"sum", "=", "0"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"FF", "=", RowBox[{"VV", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"mj", "=", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"FF", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "i", ",", "1"}], "]"}], "]"}], "*", RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "mj"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"VV", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{ RowBox[{"FF", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sum", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "mj"}], "}"}]}], "]"}]}]}]}], "\[IndentingNewLine]", " ", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"max", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], ",", RowBox[{"max", "=", RowBox[{"Append", "[", RowBox[{"max", ",", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"VV", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"VV", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "/", RowBox[{"dM", "[", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}]}], " ", "]"}]}]}], "]"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"sum", "=", RowBox[{"sum", "+", RowBox[{"Max", "[", "max", "]"}]}]}]}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"sum", "/", "k"}]}]}], "\[IndentingNewLine]", " ", "]"}]}]}], "Input", CellChangeTimes->{{3.5923719395193915`*^9, 3.5923719412041945`*^9}, { 3.5923720988112717`*^9, 3.592372102555278*^9}, {3.5923980453197746`*^9, 3.5923980576497936`*^9}, {3.5923980876578383`*^9, 3.5923981052308645`*^9}, { 3.5923981402759175`*^9, 3.5923981646229544`*^9}, {3.5923990668913393`*^9, 3.592399080523359*^9}, {3.5923991283174324`*^9, 3.592399134126441*^9}, { 3.5923991986715393`*^9, 3.592399230828588*^9}, {3.5923993899258285`*^9, 3.592399407617855*^9}, {3.592399502093996*^9, 3.592399516165018*^9}, { 3.592399553822075*^9, 3.5923995766501102`*^9}, {3.8067578646197467`*^9, 3.8067578652497487`*^9}}, CellLabel->"In[76]:=",ExpressionUUID->"9560d7b0-4a41-4f89-8cce-a870e7d86611"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis Silhouette With Criterion", "Subsection", CellChangeTimes->{{3.5834631917641993`*^9, 3.5834632080712285`*^9}, { 3.5845891037200985`*^9, 3.5845891120549145`*^9}, {3.584590326575663*^9, 3.5845903382786837`*^9}, {3.5845903719027433`*^9, 3.5845903724809446`*^9}, 3.5942968175151563`*^9, {3.8067576135858517`*^9, 3.806757614198058*^9}, 3.8067587209144917`*^9},ExpressionUUID->"44e0a73a-6e67-4c8e-b8b1-\ 2bfd3043e44c"], Cell[BoxData[{ RowBox[{"Clear", "[", "MSWC", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"MSWC", "[", RowBox[{"PI_", ",", "z_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"k", "=", RowBox[{"Length", "[", "PI", "]"}]}], ",", "alfe", ",", "bete", ",", "w", ",", "s"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"w", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Total", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "All", ",", "1"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"alfe", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "s", ",", "1"}], "]"}], "]"}], "*", RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "j", ",", "2"}], "]"}], "]"}], ",", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "s", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"w", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bete", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Min", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"q", ",", "s", ",", "1"}], "]"}], "]"}], "*", RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "j", ",", "2"}], "]"}], "]"}], ",", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"q", ",", "s", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "q", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "q", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"w", "[", RowBox[{"[", "q", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"q", ",", RowBox[{"Delete", "[", RowBox[{ RowBox[{"Range", "[", "k", "]"}], ",", "i"}], "]"}]}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"s", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"bete", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "-", RowBox[{"alfe", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"Max", "[", RowBox[{ RowBox[{"alfe", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", RowBox[{"bete", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"N", "[", RowBox[{"Mean", "[", RowBox[{"Flatten", "[", "s", "]"}], "]"}], "]"}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}]}]}], "Input", CellChangeTimes->{{3.583463220558251*^9, 3.583463254744911*^9}, { 3.583463291727977*^9, 3.5834633202820272`*^9}, {3.5834634235576096`*^9, 3.583463428862619*^9}, 3.58346354254642*^9, 3.584588843561839*^9, 3.584589201132472*^9, {3.5845901320879188`*^9, 3.584590134850124*^9}, { 3.5867436569533305`*^9, 3.586743657748932*^9}, {3.5867437429732823`*^9, 3.586743772759735*^9}, {3.586743824764627*^9, 3.586743839558453*^9}, { 3.586745430168068*^9, 3.586745461891924*^9}, {3.592360606845284*^9, 3.5923606409001436`*^9}, {3.592370010708803*^9, 3.592370011800805*^9}, { 3.5923953503266897`*^9, 3.592395496968911*^9}, 3.592395537761973*^9, { 3.806757905270878*^9, 3.806757941965201*^9}, {3.8067587054784565`*^9, 3.8067587261237345`*^9}}, CellLabel->"In[78]:=",ExpressionUUID->"1743380c-9cc9-454c-9903-9faeeca4b61d"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis Simplify Silhouette With Criterion ", "Subsection", CellChangeTimes->{{3.5834631917641993`*^9, 3.5834632080712285`*^9}, { 3.5845891037200985`*^9, 3.5845891120549145`*^9}, {3.584590326575663*^9, 3.5845903382786837`*^9}, {3.5845903719027433`*^9, 3.5845903724809446`*^9}, {3.592357039074562*^9, 3.5923570525565825`*^9}, 3.5923973781397605`*^9, 3.594296812675149*^9, {3.8067576185228186`*^9, 3.806757619072384*^9}, {3.806758778163951*^9, 3.8067587808213825`*^9}},ExpressionUUID->"6ed3c9ba-db42-4003-a28b-\ 38698789e68a"], Cell[BoxData[ RowBox[{ RowBox[{"MSSWC", "[", RowBox[{"PI_", ",", "c_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"k", "=", RowBox[{"Length", "[", "PI", "]"}]}], ",", "alfe", ",", "bete", ",", "s", ",", "W"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"alfe", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "j", ",", "2"}], "]"}], "]"}], ",", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bete", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Min", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"dM", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "j", ",", "2"}], "]"}], "]"}], ",", RowBox[{"c", "[", RowBox[{"[", "q", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "q", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"q", ",", RowBox[{"Delete", "[", RowBox[{ RowBox[{"Range", "[", "k", "]"}], ",", "i"}], "]"}]}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"s", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "j", ",", "1"}], "]"}], "]"}], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"bete", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "-", RowBox[{"alfe", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"Max", "[", RowBox[{ RowBox[{"alfe", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", RowBox[{"bete", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Return", "[", RowBox[{"N", "[", RowBox[{"Mean", "[", RowBox[{"Flatten", "[", "s", "]"}], "]"}], "]"}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"W", "=", RowBox[{"Total", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Total", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "All", ",", "1"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"N", "[", RowBox[{ RowBox[{"Total", "[", RowBox[{"Flatten", "[", "s", "]"}], "]"}], "/", "W"}], "]"}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}]}]], "Input", CellChangeTimes->{{3.583463220558251*^9, 3.583463254744911*^9}, { 3.583463291727977*^9, 3.5834633202820272`*^9}, {3.5834634235576096`*^9, 3.583463428862619*^9}, 3.58346354254642*^9, 3.584588843561839*^9, 3.584589201132472*^9, {3.5845901320879188`*^9, 3.584590134850124*^9}, { 3.5867436569533305`*^9, 3.586743657748932*^9}, {3.5867437429732823`*^9, 3.586743772759735*^9}, {3.586743824764627*^9, 3.586743839558453*^9}, { 3.586745430168068*^9, 3.586745461891924*^9}, {3.592357062248597*^9, 3.5923570637185993`*^9}, {3.5923571024436584`*^9, 3.592357102753659*^9}, { 3.592357135817709*^9, 3.5923571652547536`*^9}, {3.592357353128038*^9, 3.5923573728030677`*^9}, 3.592357582329384*^9, {3.5923576130564303`*^9, 3.5923576597765017`*^9}, {3.592357710182578*^9, 3.5923577540406437`*^9}, { 3.59235778467169*^9, 3.5923578022587166`*^9}, {3.5923578572888002`*^9, 3.5923579155058875`*^9}, {3.59236983235369*^9, 3.592369873038561*^9}, { 3.5923699811467514`*^9, 3.592370058633087*^9}, {3.59237022530378*^9, 3.592370231325391*^9}, {3.5923705650255766`*^9, 3.592370598097635*^9}, { 3.592370663165349*^9, 3.592370686003789*^9}, {3.592370881830933*^9, 3.592370882283334*^9}, {3.5923709757430983`*^9, 3.592370983792712*^9}, { 3.5923710554280376`*^9, 3.592371155673814*^9}, {3.5923711869830694`*^9, 3.59237118738867*^9}, {3.5923714359293065`*^9, 3.592371454586939*^9}, { 3.5923715623519287`*^9, 3.5923715657995343`*^9}, {3.5923716290264454`*^9, 3.592371663642906*^9}, {3.5923722257798944`*^9, 3.5923723378348913`*^9}, { 3.592372367958544*^9, 3.5923724138226247`*^9}, {3.592372798816701*^9, 3.592372806959915*^9}, {3.5923968273069143`*^9, 3.592396905808032*^9}, { 3.5943031426327653`*^9, 3.594303195793847*^9}, {3.594303254945936*^9, 3.5943032802739744`*^9}, {3.6189816916533003`*^9, 3.618981722759755*^9}, { 3.618981766159031*^9, 3.618981766533432*^9}, {3.6189818387771587`*^9, 3.618981841663164*^9}, 3.806758045114294*^9, 3.80675870178259*^9}, CellLabel->"In[80]:=",ExpressionUUID->"a0965a9a-126a-4fb1-adc4-de251ed20952"] }, Open ]], Cell[CellGroupData[{ Cell["Mahalanobis Area index", "Subsection", CellChangeTimes->{{3.5834631917641993`*^9, 3.5834632080712285`*^9}, { 3.5845891037200985`*^9, 3.5845891120549145`*^9}, {3.584590326575663*^9, 3.5845903382786837`*^9}, {3.5845903719027433`*^9, 3.5845903724809446`*^9}, { 3.592364941893098*^9, 3.5923649476495085`*^9}, {3.8067576278521585`*^9, 3.8067576371502104`*^9}},ExpressionUUID->"dbc85fc2-4ec6-49af-8285-\ c2ade7da3dbb"], Cell[BoxData[{ RowBox[{"Clear", "[", "MArea", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"MArea", "[", RowBox[{"PI_", ",", "z_", ",", "Scov_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Det", "[", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "/", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"Det", "[", RowBox[{"Scov", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}], "/", RowBox[{"Total", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "1"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "]"}]}]}], "Input", CellChangeTimes->{{3.583463220558251*^9, 3.583463254744911*^9}, { 3.583463291727977*^9, 3.5834633202820272`*^9}, {3.5834634235576096`*^9, 3.583463428862619*^9}, 3.58346354254642*^9, 3.584588843561839*^9, 3.584589201132472*^9, {3.5845901320879188`*^9, 3.584590134850124*^9}, { 3.5867436569533305`*^9, 3.586743657748932*^9}, {3.5867437429732823`*^9, 3.586743772759735*^9}, {3.586743824764627*^9, 3.586743839558453*^9}, { 3.586745430168068*^9, 3.586745461891924*^9}, {3.592360606845284*^9, 3.5923606409001436`*^9}, {3.592365663863366*^9, 3.592365776885565*^9}, { 3.5943041512933016`*^9, 3.594304167456326*^9}, {3.594304204796383*^9, 3.594304207877388*^9}, {3.7230994370884895`*^9, 3.723099443201758*^9}, { 3.7231001985974555`*^9, 3.7231002036799326`*^9}, {3.723100321867401*^9, 3.723100330280794*^9}, {3.7231004011513214`*^9, 3.723100425611371*^9}, { 3.7231016860908623`*^9, 3.7231016966699696`*^9}, {3.8067580828086715`*^9, 3.806758084046102*^9}, {3.830083561401351*^9, 3.830083572156329*^9}}, CellLabel->"In[81]:=",ExpressionUUID->"23d3b2b2-0106-467a-a16d-300af3c581e6"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Louding of partitions and calculating centers and covariance matrices\ \>", "Section", CellChangeTimes->{{3.5922675537233934`*^9, 3.5922675637164087`*^9}, { 3.5922676918126016`*^9, 3.592267717667641*^9}, {3.5960952868894935`*^9, 3.5960953699236126`*^9}, {3.6042050556918383`*^9, 3.60420506219705*^9}},ExpressionUUID->"213c3965-022c-457b-b1fc-\ 3806c6c56fbe"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"SetDirectory", "[", RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ StyleBox["PI", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["ReadList", FontColor->GrayLevel[0]], "[", "\"\\"", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "\t\t\t\t\t\t\t", RowBox[{"(*", " ", "Centers", " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"n", "=", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1", ",", "2"}], "]"}], "]"}], "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"centri", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "s", ",", "1"}], "]"}], "]"}], "*", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "s", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"Sum", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "s", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "\t\t\t\t\t", RowBox[{"(*", " ", RowBox[{"Covariance", " ", "Matrices", " ", "Scov"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pom", "=", RowBox[{"mj", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{ RowBox[{"Length", "[", "PI", "]"}], "+", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "-", "1"}]}], "}"}]}], "]"}]}]}], ";", RowBox[{"Scov", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"pom", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "s", ",", "1"}], "]"}], "]"}], "]"}], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "All", ",", "2"}], "]"}], "]"}], "[", RowBox[{"[", "s", "]"}], "]"}], "-", RowBox[{"centri", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"mj", "[", RowBox[{"[", "j", "]"}], "]"}], "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j", ",", "s", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}]}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"Scov", "[", RowBox[{"[", RowBox[{"br", ",", "j"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"pom", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "s"}], "]"}], "]"}], ".", RowBox[{"pom", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "r"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"r", ",", "n"}], "}"}]}], "]"}], "/", RowBox[{"mj", "[", RowBox[{"[", "j", "]"}], "]"}]}]}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "\t\t\t\t\t\t\t", RowBox[{"(*", " ", RowBox[{"Center", " ", "of", " ", "the", " ", "set", " ", "A"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Pod", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "1"}], "]"}]}], ";", " ", RowBox[{"m", "=", RowBox[{"Length", "[", "Pod", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cc", "=", RowBox[{ RowBox[{"Total", "[", RowBox[{ RowBox[{"Pod", "[", RowBox[{"[", RowBox[{"All", ",", "1"}], "]"}], "]"}], "*", RowBox[{"Pod", "[", RowBox[{"[", RowBox[{"All", ",", "2"}], "]"}], "]"}]}], "]"}], "/", RowBox[{"Total", "[", RowBox[{"Pod", "[", RowBox[{"[", RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.592267972600028*^9, 3.5922679728900285`*^9}, { 3.5922766565011435`*^9, 3.59227667421317*^9}, {3.592278977610648*^9, 3.592279002458687*^9}, {3.592279613481639*^9, 3.592279632843668*^9}, { 3.592356114911152*^9, 3.592356115151152*^9}, {3.59235620437529*^9, 3.592356228299327*^9}, 3.592356307965447*^9, 3.592366876579296*^9, 3.592368171573773*^9, {3.592368582712495*^9, 3.5923685982813225`*^9}, { 3.592368746575183*^9, 3.592368757604402*^9}, {3.5923689604369583`*^9, 3.5923689680809717`*^9}, {3.5923690024946322`*^9, 3.592369019342662*^9}, { 3.5923691517244945`*^9, 3.592369151974095*^9}, {3.5923695444083843`*^9, 3.592369544923185*^9}, {3.592371245031771*^9, 3.5923712453281717`*^9}, { 3.592372921776117*^9, 3.592372939029747*^9}, {3.5923732637443175`*^9, 3.592373264103118*^9}, {3.592451395888426*^9, 3.592451396348427*^9}, 3.592534673532527*^9, {3.592537227122818*^9, 3.5925372288838215`*^9}, { 3.5929745899005013`*^9, 3.592974590270502*^9}, {3.59306616000082*^9, 3.593066160260821*^9}, {3.593072364916997*^9, 3.593072365376998*^9}, { 3.593352606044753*^9, 3.593352608674757*^9}, {3.5933538741066823`*^9, 3.5933538748566837`*^9}, {3.59335391206874*^9, 3.5933539389727807`*^9}, { 3.593354049998949*^9, 3.593354052318952*^9}, {3.593354298799322*^9, 3.5933542992193227`*^9}, {3.5936996553672085`*^9, 3.593699655897209*^9}, { 3.5937515762997932`*^9, 3.593751578319796*^9}, {3.594219272539912*^9, 3.5942192741899147`*^9}, {3.594221303420991*^9, 3.594221304770993*^9}, { 3.594301377584095*^9, 3.5943013800860987`*^9}, {3.5943016697015343`*^9, 3.594301670591535*^9}, {3.5943018711568413`*^9, 3.594301873606845*^9}, { 3.5943020888441644`*^9, 3.5943020901341667`*^9}, {3.5943022794354515`*^9, 3.594302280675453*^9}, {3.5943042538044577`*^9, 3.594304278069495*^9}, { 3.5943053209550896`*^9, 3.5943053213450904`*^9}, {3.5945342354778986`*^9, 3.5945342366679*^9}, {3.594534397140143*^9, 3.594534398170144*^9}, { 3.595568462330875*^9, 3.5955684633408766`*^9}, {3.595568746437312*^9, 3.595568748647315*^9}, {3.595942500781396*^9, 3.5959425013613977`*^9}, { 3.5960445849841785`*^9, 3.5960445887841835`*^9}, {3.5960953779866242`*^9, 3.596095422277691*^9}, {3.5961044632194796`*^9, 3.5961044651204824`*^9}, { 3.5962080083801923`*^9, 3.596208008600192*^9}, 3.596384486671713*^9, 3.6014437673922243`*^9, 3.601444918597246*^9, 3.602741778010562*^9, { 3.6027418460110817`*^9, 3.6027418484134855`*^9}, {3.602741987986931*^9, 3.602741988361332*^9}, 3.602742091243512*^9, {3.6027422127053256`*^9, 3.602742213438527*^9}, {3.602742413868679*^9, 3.60274241438348*^9}, 3.6027426094150224`*^9, 3.6027433244582787`*^9, {3.6042053289947195`*^9, 3.6042053352035303`*^9}, {3.7817924420866485`*^9, 3.781792463073543*^9}}, CellLabel->"In[83]:=",ExpressionUUID->"52941722-142c-424c-8f8a-5a1b2f3b207b"] }, Closed]], Cell[CellGroupData[{ Cell[TextData[{ "Indexes\n", StyleBox["MSWC takes a very long time, and the results are similar as by \ MSSWC", FontFamily->"Arial", FontSize->14, FontColor->RGBColor[1, 0.5, 0]] }], "Section", CellChangeTimes->{{3.593072390127034*^9, 3.5930724033190536`*^9}, { 3.593072450385124*^9, 3.5930724508051248`*^9}, {3.5960954603747454`*^9, 3.596095523247838*^9}, 3.5960955749039125`*^9, 3.806757943930112*^9, { 3.8067581321221747`*^9, 3.8067581373878746`*^9}, {3.806758811090575*^9, 3.806758814152806*^9}},ExpressionUUID->"995faa6c-957e-4031-a017-\ 9d7faa5f1a0a"], Cell[BoxData[ RowBox[{"Do", "[", RowBox[{ RowBox[{"i", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", RowBox[{"centri", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "\"\<; k= \>\"", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", "\"\<; mj:\>\"", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"i", ",", "s"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", "\"\<; F=\>\"", ",", StyleBox[ RowBox[{"MWF", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5925347100005813`*^9, 3.592534760661658*^9}, { 3.592534827023759*^9, 3.592534832807768*^9}, {3.5925351927703133`*^9, 3.5925352197353544`*^9}, {3.592535403355632*^9, 3.59253543564568*^9}, { 3.5925354817577505`*^9, 3.592535562356873*^9}, 3.5930662121718984`*^9, 3.810879922443553*^9}, CellLabel->"In[92]:=",ExpressionUUID->"e3ab93be-ae86-4b2a-803d-ba17102a1e7b"] }, Closed]], Cell[CellGroupData[{ Cell["Figures", "Section", CellChangeTimes->{{3.592267645221532*^9, 3.592267652588543*^9}, { 3.592397248146564*^9, 3.5923972656395907`*^9}, {3.593072438522106*^9, 3.593072446613118*^9}},ExpressionUUID->"a8f8d477-73cd-48f2-93f1-\ b4e3cd88b649"], Cell[BoxData[{ RowBox[{ RowBox[{"IndF", "=", RowBox[{"IndArea", "=", RowBox[{"IndMCH", "=", RowBox[{"IndMSSWC", "=", RowBox[{"IndSWC", "=", RowBox[{"IndMDB", "=", RowBox[{"{", "}"}]}]}]}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"br0", "=", RowBox[{ RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "-", "1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ StyleBox["IndF", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndF", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"MWF", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";", " ", RowBox[{"length", "=", RowBox[{"Length", "[", "IndF", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndF", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndF", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], RowBox[{ RowBox[{ RowBox[{ StyleBox["IndArea", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndArea", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"MArea", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndArea", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndArea", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ StyleBox["IndMCH", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndMCH", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"MCH", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndMCH", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndMCH", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ StyleBox["IndMDB", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndMDB", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"MDB", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndMDB", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndMDB", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ StyleBox["IndMSSWC", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndMSSWC", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"MSSWC", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndMSSWC", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndMSSWC", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", "\[IndentingNewLine]", RowBox[{ RowBox[{"time", "=", RowBox[{"Timing", "[", RowBox[{"Do", "[", StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], RowBox[{ RowBox[{ RowBox[{ StyleBox["IndSWC", FontColor->GrayLevel[0]], StyleBox["=", FontColor->GrayLevel[0]], RowBox[{ StyleBox["Append", FontColor->GrayLevel[0]], StyleBox["[", FontColor->GrayLevel[0]], RowBox[{ StyleBox["IndSWC", FontColor->GrayLevel[0]], StyleBox[",", FontColor->GrayLevel[0]], StyleBox[ RowBox[{"SWC", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"centri", "[", RowBox[{"[", "br", "]"}], "]"}], ",", RowBox[{"Scov", "[", RowBox[{"[", "br", "]"}], "]"}]}], "]"}], FontColor->RGBColor[0.5, 0, 0.5]]}], StyleBox["]", FontColor->GrayLevel[0]]}]}], StyleBox[";", FontColor->GrayLevel[0]]}], StyleBox[" ", FontColor->GrayLevel[0]], StyleBox["\[IndentingNewLine]", FontColor->GrayLevel[0]], ",", RowBox[{"{", RowBox[{"br", ",", RowBox[{"Length", "[", "PI", "]"}]}], "}"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"TabIndSWC", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "+", "br0"}], ",", RowBox[{"IndSWC", "[", RowBox[{"[", "j", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "length"}], "}"}]}], "]"}]}], ",", "\"\<; Time=\>\"", ",", RowBox[{"time", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\< sec\>\""}], "]"}], ";"}], "\[IndentingNewLine]", "*)"}], "\[IndentingNewLine]", "\t\t\t\t\t\t\t\t\t", RowBox[{"(*", " ", "Slike", " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slF", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndF", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndF", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slArea", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndArea", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndArea", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMCH", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMCH", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMCH", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMDB", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMDB", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMDB", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMSSWC", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMSSWC", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMSSWC", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"Style", "[", RowBox[{ RowBox[{"GraphicsGrid", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"slF", ",", "slMCH", ",", "slMDB"}], "}"}], ",", RowBox[{"{", RowBox[{"slMSSWC", ",", "slArea"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Magnification", "\[Rule]", "1.2"}]}], "]"}], "]"}], "\t\t\t\t"}]}], "Input", CellChangeTimes->{{3.5922767843803368`*^9, 3.5922767847403374`*^9}, { 3.5922768149843826`*^9, 3.5922768353674135`*^9}, {3.592276877102475*^9, 3.592276994524654*^9}, {3.5922770482317343`*^9, 3.5922770486217346`*^9}, { 3.592277081627785*^9, 3.592277089544797*^9}, {3.592277129862859*^9, 3.5922771509058905`*^9}, {3.5922772042709713`*^9, 3.5922772054949727`*^9}, {3.592277249434039*^9, 3.5922772699900703`*^9}, { 3.592277385496244*^9, 3.5922774001142654`*^9}, {3.592277556759503*^9, 3.5922775572125034`*^9}, {3.5922776752556815`*^9, 3.592277745196787*^9}, { 3.5922777798598385`*^9, 3.5922777964828634`*^9}, {3.592278041512235*^9, 3.5922780419022355`*^9}, {3.592278144680389*^9, 3.5922782428335376`*^9}, { 3.5922782977766204`*^9, 3.5922783733667345`*^9}, {3.5922784508798523`*^9, 3.5922784543908577`*^9}, {3.592279391765274*^9, 3.5922794078452983`*^9}, { 3.592279671248725*^9, 3.592279729317813*^9}, {3.592279823022954*^9, 3.5922798770120344`*^9}, {3.5922799850151987`*^9, 3.5922799921662087`*^9}, {3.592280272662634*^9, 3.5922803157536993`*^9}, { 3.592359437108234*^9, 3.592359468235281*^9}, {3.5923595219113626`*^9, 3.5923595519644084`*^9}, {3.5923596065294914`*^9, 3.5923596839786096`*^9}, 3.5923597814197607`*^9, {3.5923600046391*^9, 3.592360007670105*^9}, { 3.592360679307411*^9, 3.592360679791012*^9}, {3.5923610006835756`*^9, 3.592361006502386*^9}, {3.592362031783987*^9, 3.592362142606581*^9}, 3.5923621847110553`*^9, {3.592363048953573*^9, 3.592363070980812*^9}, { 3.5923635751736975`*^9, 3.5923636997243166`*^9}, {3.5923657961983986`*^9, 3.5923659355690436`*^9}, 3.592365979997922*^9, {3.5923660127579794`*^9, 3.5923660330224147`*^9}, {3.5923661766362667`*^9, 3.5923662523900003`*^9}, {3.592366294962475*^9, 3.5923662958048763`*^9}, { 3.592366932364994*^9, 3.592366943269414*^9}, {3.5923692634362907`*^9, 3.5923694404030013`*^9}, {3.5923728677376223`*^9, 3.592372893914468*^9}, { 3.5923731240772724`*^9, 3.592373228035855*^9}, {3.592373698330881*^9, 3.592373786845436*^9}, {3.5923962188050027`*^9, 3.5923963133711452`*^9}, { 3.592396349962201*^9, 3.592396351212203*^9}, {3.5923964129122963`*^9, 3.592396440324338*^9}, {3.592396481531399*^9, 3.592396589980562*^9}, { 3.592396630297623*^9, 3.5923966594956665`*^9}, {3.5923969717381315`*^9, 3.5923971315743723`*^9}, 3.5923972042904825`*^9, {3.59239733135069*^9, 3.592397340208703*^9}, {3.592398192557997*^9, 3.592398249349083*^9}, { 3.5923982815751333`*^9, 3.5923983295982056`*^9}, 3.593067159438328*^9, { 3.59306734684861*^9, 3.593067347588611*^9}, {3.5930674130527105`*^9, 3.593067529171885*^9}, {3.593067598300989*^9, 3.5930678115703125`*^9}, 3.5930678564623785`*^9, {3.5930725259172387`*^9, 3.5930725338102503`*^9}, { 3.723099651109883*^9, 3.7230996575991077`*^9}, {3.7230998406952267`*^9, 3.7230998824463377`*^9}, {3.72310014793165*^9, 3.7231001859186773`*^9}, 3.723100386558505*^9, {3.723100445378959*^9, 3.723100461086749*^9}, { 3.723101729840216*^9, 3.723101753424996*^9}, {3.7231017856898346`*^9, 3.723101788425127*^9}, {3.806477441462988*^9, 3.8064774420061016`*^9}, { 3.8064774787817545`*^9, 3.8064775135048876`*^9}, {3.806757791200276*^9, 3.806757812495595*^9}, {3.806757867818373*^9, 3.8067578772666516`*^9}, { 3.806757946039962*^9, 3.8067579521565943`*^9}, {3.8067580859722486`*^9, 3.806758095780899*^9}, {3.806758571871671*^9, 3.80675862533539*^9}, { 3.8067588558827767`*^9, 3.8067589087243614`*^9}, 3.8108799327423573`*^9}, CellLabel->"In[93]:=",ExpressionUUID->"e3826702-8d1f-4a49-b7a8-11590686e1f5"], Cell[CellGroupData[{ Cell["Figures", "Subsubsection", CellChangeTimes->{{3.8064775434683375`*^9, 3.806477551379848*^9}},ExpressionUUID->"52b0fdd6-8474-4c80-90fb-\ e93469cb79d9"], Cell[BoxData[{ RowBox[{ RowBox[{"slArea", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndArea", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndArea", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMCH", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMCH", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMCH", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMDB", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMDB", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMDB", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slMSSWC", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMSSWC", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", RowBox[{".75", RowBox[{"Min", "[", "IndMSSWC", "]"}]}]}], "}"}]}], ",", RowBox[{"PlotLabel", "->", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"Style", "[", RowBox[{ RowBox[{"GraphicsGrid", "[", RowBox[{"{", RowBox[{"{", RowBox[{"slMCH", ",", "slMDB", ",", "slMSSWC", ",", "slArea"}], "}"}], "}"}], "]"}], ",", RowBox[{"Magnification", "\[Rule]", "1.2"}]}], "]"}], "]"}], "\t\t"}]}], "Input", CellChangeTimes->{{3.594276001775733*^9, 3.594276023376069*^9}, { 3.594276175824193*^9, 3.594276177210215*^9}, {3.6027428073013697`*^9, 3.602742871620283*^9}, {3.6027435203322225`*^9, 3.6027435208938236`*^9}, { 3.604205534119479*^9, 3.6042055903575783`*^9}, {3.6042056304184484`*^9, 3.6042056460864763`*^9}, {3.781792547702299*^9, 3.7817926210512414`*^9}, { 3.781792652699639*^9, 3.7817927598232574`*^9}, {3.7817927997475233`*^9, 3.781792853313323*^9}, {3.781792884138916*^9, 3.781792916995084*^9}, { 3.7817929534915085`*^9, 3.7817934499702363`*^9}, {3.7817935631916523`*^9, 3.7817935650028076`*^9}, {3.7817936345413723`*^9, 3.781793698923254*^9}, { 3.8064775856033993`*^9, 3.8064776191115847`*^9}, {3.806757827342328*^9, 3.806757831429141*^9}, {3.8067578860613813`*^9, 3.806757888579191*^9}, { 3.8067580072057085`*^9, 3.806758019998727*^9}, {3.8067589309667616`*^9, 3.80675895835153*^9}}, CellLabel-> "In[120]:=",ExpressionUUID->"5a1b8754-92b2-4fdc-a4a3-3b8abeadd878"] }, Open ]], Cell[CellGroupData[{ Cell["Figures corrected", "Subsubsection", CellChangeTimes->{{3.8064775434683375`*^9, 3.806477551379848*^9}, { 3.830085560995288*^9, 3.830085565317997*^9}},ExpressionUUID->"eec0e165-a66b-41ee-b8c1-\ 5caeceef0350"], Cell[BoxData[{ RowBox[{ RowBox[{"slArea", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndArea", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}], ",", RowBox[{"Ticks", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{".01", ",", ".02", ",", ".03"}], "}"}]}], "}"}]}], ",", RowBox[{"GridLines", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{"0", ",", ".01", ",", ".02", ",", ".03", ",", ".04"}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", ".03"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slVCH", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMCH", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", "1000"}], "}"}]}], ",", RowBox[{"Ticks", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{"1000", ",", "2000", ",", "3000", ",", "4000", ",", "5000"}], "}"}]}], "}"}]}], ",", RowBox[{"GridLines", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{"1000", ",", "2000", ",", "3000", ",", "4000", ",", "5000"}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"1000", ",", "5000"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slVDB", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMDB", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}], ",", RowBox[{"Ticks", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{".1", ",", ".2", ",", ".3", ",", ".4"}], "}"}]}], "}"}]}], ",", RowBox[{"GridLines", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{".1", ",", ".2", ",", ".3", ",", ".4"}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", ".4"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"slSWC", "=", RowBox[{"ListPlot", "[", RowBox[{"TabIndMSSWC", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"1", ",", ".8"}], "}"}]}], ",", RowBox[{"GridLines", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{".825", ",", ".85", ",", ".875", ",", ".90"}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{".8", ",", ".9"}], "}"}]}], ",", RowBox[{"Ticks", "\[Rule]", RowBox[{"{", RowBox[{"Automatic", ",", RowBox[{"{", RowBox[{".75", ",", ".8", ",", ".85", ",", ".90"}], "}"}]}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"GraphicsGrid", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"slVCH", ",", "slVDB", ",", "slSWC", ",", "slArea"}], "}"}], "}"}], ",", RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], "]"}], "\t"}]}], "Input", CellChangeTimes->{{3.813119843882594*^9, 3.8131198680792007`*^9}, { 3.8131199354990053`*^9, 3.8131199549856014`*^9}, {3.8131199869629803`*^9, 3.813120150682181*^9}, 3.8131216801464806`*^9, {3.813121714084735*^9, 3.813121795116723*^9}, 3.8300842941639147`*^9, {3.830084337856185*^9, 3.830084348776643*^9}, {3.8300844374489784`*^9, 3.830084450251973*^9}, { 3.8300845070572557`*^9, 3.830084532334301*^9}, {3.8300846037874784`*^9, 3.83008460636095*^9}}, CellLabel-> "In[163]:=",ExpressionUUID->"43131d3f-2dd2-43c0-819f-ba933517ec9a"] }, Open ]] }, Closed]] }, WindowSize->{793, 717}, WindowMargins->{{-837, Automatic}, {-195, Automatic}}, ShowSelection->True, FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 742, 15, 87, "Section",ExpressionUUID->"8f9af442-96cd-4772-bff6-09a71bc45d4e"], Cell[1324, 39, 1437, 33, 162, "Input",ExpressionUUID->"b3a8321d-58d6-4aff-95d5-5ac599347900"] }, Open ]], Cell[CellGroupData[{ Cell[2798, 77, 154, 3, 67, "Section",ExpressionUUID->"4ab43163-a9ad-4579-9e58-03db053488e5"], Cell[CellGroupData[{ Cell[2977, 84, 206, 4, 53, "Subsection",ExpressionUUID->"dafb582f-fc59-4d24-b27f-fba8dcac4aec"], Cell[3186, 90, 781, 21, 48, "Input",ExpressionUUID->"31c58d2a-fd13-4f58-85c9-49292444865a"] }, Open ]], Cell[CellGroupData[{ Cell[4004, 116, 178, 3, 53, "Subsection",ExpressionUUID->"d36b648b-ba45-48dc-940b-54c26491d94c"], Cell[4185, 121, 1885, 45, 86, "Input",ExpressionUUID->"73251baf-9953-4113-bb50-712ccc292076"] }, Open ]], Cell[CellGroupData[{ Cell[6107, 171, 323, 5, 53, "Subsection",ExpressionUUID->"8ff1266c-46b6-4567-91e8-06d92d292be2"], Cell[6433, 178, 3540, 90, 162, "Input",ExpressionUUID->"a57c4cc1-6972-4d8a-919e-3ccc245a32c1"] }, Open ]], Cell[CellGroupData[{ Cell[10010, 273, 264, 5, 53, "Subsection",ExpressionUUID->"a78eb0ac-1a20-48c3-ad54-3de3bb91e9a6"], Cell[10277, 280, 4603, 114, 333, "Input",ExpressionUUID->"9560d7b0-4a41-4f89-8cce-a870e7d86611"] }, Open ]], Cell[CellGroupData[{ Cell[14917, 399, 446, 6, 53, "Subsection",ExpressionUUID->"44e0a73a-6e67-4c8e-b8b1-2bfd3043e44c"], Cell[15366, 407, 5704, 147, 295, "Input",ExpressionUUID->"1743380c-9cc9-454c-9903-9faeeca4b61d"] }, Open ]], Cell[CellGroupData[{ Cell[21107, 559, 558, 8, 53, "Subsection",ExpressionUUID->"6ed3c9ba-db42-4003-a28b-38698789e68a"], Cell[21668, 569, 6134, 144, 200, "Input",ExpressionUUID->"a0965a9a-126a-4fb1-adc4-de251ed20952"] }, Open ]], Cell[CellGroupData[{ Cell[27839, 718, 429, 6, 53, "Subsection",ExpressionUUID->"dbc85fc2-4ec6-49af-8285-c2ade7da3dbb"], Cell[28271, 726, 2356, 51, 105, "Input",ExpressionUUID->"23d3b2b2-0106-467a-a16d-300af3c581e6"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[30676, 783, 377, 7, 91, "Section",ExpressionUUID->"213c3965-022c-457b-b1fc-3806c6c56fbe"], Cell[31056, 792, 10516, 265, 542, "Input",ExpressionUUID->"52941722-142c-424c-8f8a-5a1b2f3b207b"] }, Closed]], Cell[CellGroupData[{ Cell[41609, 1062, 580, 13, 71, "Section",ExpressionUUID->"995faa6c-957e-4031-a017-9d7faa5f1a0a"], Cell[42192, 1077, 1688, 42, 105, "Input",ExpressionUUID->"e3ab93be-ae86-4b2a-803d-ba17102a1e7b"] }, Closed]], Cell[CellGroupData[{ Cell[43917, 1124, 249, 4, 53, "Section",ExpressionUUID->"a8f8d477-73cd-48f2-93f1-b4e3cd88b649"], Cell[44169, 1130, 18824, 510, 903, "Input",ExpressionUUID->"e3826702-8d1f-4a49-b7a8-11590686e1f5"], Cell[CellGroupData[{ Cell[63018, 1644, 160, 3, 44, "Subsubsection",ExpressionUUID->"52b0fdd6-8474-4c80-90fb-e93469cb79d9"], Cell[63181, 1649, 3238, 78, 181, "Input",ExpressionUUID->"5a1b8754-92b2-4fdc-a4a3-3b8abeadd878"] }, Open ]], Cell[CellGroupData[{ Cell[66456, 1732, 219, 4, 44, "Subsubsection",ExpressionUUID->"eec0e165-a66b-41ee-b8c1-5caeceef0350"], Cell[66678, 1738, 4135, 113, 257, "Input",ExpressionUUID->"43131d3f-2dd2-43c0-819f-ba933517ec9a"] }, Open ]] }, Closed]] } ] *)