Distribution.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     68044,       1496]
NotebookOptionsPosition[     63898,       1442]
NotebookOutlinePosition[     64294,       1458]
CellTagsIndexPosition[     64251,       1455]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.7855648473582897`*^9, 3.785564848357767*^9}},
 CellLabel->"In[2]:=",ExpressionUUID->"fe31bd78-5fbf-4a06-86f0-b80159ec6f68"],

Cell[BoxData[
 RowBox[{"NotebookDirectory", "[", "]"}]], "Input",
 CellLabel->"In[3]:=",ExpressionUUID->"2c8e21e4-f654-4aac-8fa2-d6a56104d0b2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{
   "Import", "[", 
    "\"\</home/dlowl/work/stochastic_work_folder/datasets/Gaublomme_GSE75109_\
Trbv.csv\>\"", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7855647734947567`*^9, 3.785564807518098*^9}, {
   3.787911813047371*^9, 3.787911818500882*^9}, {3.7879171517734346`*^9, 
   3.7879171603055706`*^9}, 3.788072188466508*^9, {3.789275707971333*^9, 
   3.789275723700508*^9}, {3.7941322992389307`*^9, 3.794132358229192*^9}, {
   3.7942104488578*^9, 3.79421045792579*^9}, {3.795420158595336*^9, 
   3.795420175517316*^9}, {3.798297007329571*^9, 3.798297021244158*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"5e4db062-34b8-499a-b4d1-fd1bbb88c8f8"],

Cell[BoxData[
 RowBox[{"data", "[", 
  RowBox[{"[", "2", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.7941323733582277`*^9, 3.794132374392075*^9}, {
  3.7941325765177794`*^9, 3.794132617301661*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"64dbf530-b3fc-4d33-b3d8-da822f5c311b"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7941326052309237`*^9, 
  3.794132606524393*^9}},ExpressionUUID->"b41e38b3-a239-42de-9ac0-\
71a31d034b50"],

Cell[BoxData[
 RowBox[{
  RowBox[{"filtered", " ", "=", " ", 
   RowBox[{"Import", "[", "\"\<~/dopaminergic_filtered.csv\>\"", "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.7882516600298777`*^9, 3.788251673812788*^9}, {
  3.789275729346911*^9, 3.78927574108042*^9}},
 CellLabel->"In[2]:=",ExpressionUUID->"7d375fe9-4815-417c-81d9-186f56f11dfe"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GetMixture", "[", "pop_", "]"}], " ", ":=", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", "estimate", "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Print", "[", 
      RowBox[{"pop", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Dbh\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Hsf2bp\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Chaf1a\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Upk1a\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Grin2c\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Wfdc1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Glp1r\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cpt1a\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Ecrg4\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Ifi44\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cd38\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Ppef2\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cpa5\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Wnt5b\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Pon2\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Pus7l\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Pla2g3\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Setmar\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Neurod1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Epgn\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Pmch\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<P2ry14\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Aspg\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Smim3\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Col16a1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Iapp\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
           "\"\<AAdacl4fm3\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<E2f8\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Plaur\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gimap6\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Rpp38\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Shisa3\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Acot4\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Lgals7\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Emcn\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gpr174\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Ttll9\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Thbd\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<C1s2\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Sorbs2os\>\""}], 
          "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm15634\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm12204\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Kif16bos\>\""}], 
          "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cep112os1\>\""}],
           "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
           "\"\<A430108G06Rik\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm12199\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm14664\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm15562\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
           "\"\<4933411K16Rik\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm17669\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
           "\"\<1110020A21Rik\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
           "\"\<1700047M11Rik\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm34583\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm3488\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm43462\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm44148\>\""}]}],
          ",", 
         RowBox[{"Return", "[", 
          RowBox[{"-", "1"}], "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", 
       RowBox[{"dopaminergic", " ", "TPM"}], "*)"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Ang4\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Cdx1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Clec14a\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Dio2\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Esm1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm11627\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm15997\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm3985\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm7173\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Naaladl1\>\""}], 
          "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Olfr1442\>\""}], 
          "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Tcfl5\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Zpld1\>\""}]}], 
         ",", 
         RowBox[{"Return", "[", 
          RowBox[{"-", "1"}], "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", 
       RowBox[{"somatosensory", " ", "tpm"}], "*)"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<5730403I07Rik\>\""}], "  ", "||", "  ", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Clcn1\>\""}], " ", "||", " ", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Cnga4\>\""}], " ", "||", " ", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Gm1123\>\""}], " ", "||", " ", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Gm514\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Hrct1\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Nek5\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Olfr568\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Tekt2\>\""}], "||", 
          RowBox[{
           RowBox[{"pop", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
           "\"\<Uox\>\""}]}], ",", " ", 
         RowBox[{"Return", "[", 
          RowBox[{"-", "1"}], "]"}]}], "]"}], ";"}], " ", 
      RowBox[{"(*", "somatosensory", "*)"}], "*)"}], "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
          "\"\<1110002E22Rik\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
          "\"\<1700034H15Rik\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
          "\"\<2010010A06Rik\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", 
          "\"\<A930002H24Rik\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Abcc4\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Abcg5\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Arhgap6\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cabp4\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Celrr\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cep72\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Colec11\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Creb3l4\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Ctso\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Cwh43\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Dio3\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Eif4ebp1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Esyt1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Fam241a\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Fev\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm11201\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm11978\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm12279\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm13838\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm26510\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm28085\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm28610\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm3830\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm43328\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm44214\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gm7457\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Gnat1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<H2-Aa\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Kctd21\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Leap2\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Lockd\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Mgst1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Mlc1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Mrc2\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Mtbp\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Myl7\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Nsl1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Plekhg2\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Pou3f1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Scgb1c1\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Serpina9\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Stard8\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Tm4sf4\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Tpx2\>\""}], "||", 
         RowBox[{
          RowBox[{"pop", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "\"\<Xcl1\>\""}]}], 
        ",", " ", 
        RowBox[{"Return", "[", 
         RowBox[{"-", "1"}], "]"}]}], "]"}], " ", 
      RowBox[{"(*", "dopaminergic", "*)"}], "*)"}], "\[IndentingNewLine]", 
     RowBox[{"FindDistribution", "[", 
      RowBox[{"pop", "[", 
       RowBox[{"[", 
        RowBox[{"2", ";;"}], "]"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.787918570724206*^9, 3.7879186263243027`*^9}, {
   3.787918663187888*^9, 3.787918693019414*^9}, {3.787918757511053*^9, 
   3.7879188566318398`*^9}, {3.787919672068653*^9, 3.7879196722092943`*^9}, {
   3.787920316202681*^9, 3.7879203237037406`*^9}, 3.7879207268276896`*^9, {
   3.7879207571123943`*^9, 3.7879207991015663`*^9}, {3.787921638313754*^9, 
   3.787921642720581*^9}, {3.787922092601185*^9, 3.78792209786759*^9}, {
   3.7880720864928703`*^9, 3.7880720996506863`*^9}, {3.7880847054336205`*^9, 
   3.7880847229512215`*^9}, {3.788084873196802*^9, 3.7880848899643455`*^9}, {
   3.7880853643778667`*^9, 3.7880853650653877`*^9}, 3.788085425166053*^9, {
   3.788088807730646*^9, 3.7880888188258257`*^9}, {3.788090677529727*^9, 
   3.788090704407821*^9}, 3.7880942903784676`*^9, {3.788094440098759*^9, 
   3.7880944425050507`*^9}, {3.7881529022491283`*^9, 
   3.7881529564641886`*^9}, {3.7881562923616686`*^9, 
   3.7881563094107227`*^9}, {3.7881586678563776`*^9, 3.788158678279437*^9}, {
   3.788162258390897*^9, 3.7881622672984295`*^9}, {3.7881662541733584`*^9, 
   3.7881662614710226`*^9}, {3.7881677717120852`*^9, 3.78816778526029*^9}, {
   3.788168156862338*^9, 3.7881681645505486`*^9}, {3.7881712057327733`*^9, 
   3.7881712436587095`*^9}, {3.788173395136899*^9, 3.788173404559784*^9}, {
   3.789275788266651*^9, 3.789275793733*^9}, {3.789276054871378*^9, 
   3.7892760655773067`*^9}, {3.789276239606235*^9, 3.78927628020903*^9}, {
   3.789276900599433*^9, 3.7892769195662193`*^9}, {3.7892774316289577`*^9, 
   3.7892774646916933`*^9}, {3.789278961093919*^9, 3.789278974762322*^9}, {
   3.789279194049995*^9, 3.789279201127388*^9}, {3.7892795585107307`*^9, 
   3.7892795642253857`*^9}, {3.78928064551243*^9, 3.789280651935339*^9}, {
   3.7892815251387587`*^9, 3.789281532090947*^9}, {3.789282510967246*^9, 
   3.7892825170757837`*^9}, {3.78928271759953*^9, 3.789282723290739*^9}, {
   3.7892833726175413`*^9, 3.789283378962125*^9}, {3.789283647076762*^9, 
   3.789283651520021*^9}, {3.789283935564562*^9, 3.789283939805272*^9}, {
   3.78928403249575*^9, 3.7892840380302687`*^9}, {3.789284519785202*^9, 
   3.7892845311246023`*^9}, {3.789289378829529*^9, 3.789289385184067*^9}, {
   3.7892897931489487`*^9, 3.789289799763404*^9}, {3.789290107073674*^9, 
   3.7892901144756317`*^9}, {3.789290406960224*^9, 3.7892904138212442`*^9}, {
   3.7892922929002533`*^9, 3.789292299153592*^9}, {3.7892931126649847`*^9, 
   3.789293119785273*^9}, {3.7892951629696903`*^9, 3.7892951686911592`*^9}, {
   3.789296043304289*^9, 3.7892960460599003`*^9}, {3.789296079909264*^9, 
   3.789296081978405*^9}, {3.789296918684518*^9, 3.789296923666947*^9}, {
   3.78929768298759*^9, 3.789297689005872*^9}, {3.789298533997549*^9, 
   3.789298539127035*^9}, {3.789298890071073*^9, 3.7892988963932*^9}, {
   3.789299073935729*^9, 3.78929908321598*^9}, {3.789299950644627*^9, 
   3.789299957529258*^9}, {3.789300384941728*^9, 3.789300390791234*^9}, {
   3.789300637106381*^9, 3.789300646297488*^9}, {3.789301041004562*^9, 
   3.789301046854792*^9}, {3.7893612058229923`*^9, 3.789361210906672*^9}, {
   3.789361750002543*^9, 3.78936176017196*^9}, {3.789363012345319*^9, 
   3.789363017253358*^9}, {3.7893635035169888`*^9, 3.789363529370018*^9}, {
   3.78936386053679*^9, 3.789363864721506*^9}, {3.789365682479076*^9, 
   3.789365690657185*^9}, {3.789365932056561*^9, 3.789365937894856*^9}, {
   3.7893661597524652`*^9, 3.789366168020493*^9}, {3.7893664410774593`*^9, 
   3.7893664469496326`*^9}, {3.789366865868002*^9, 3.789366873967388*^9}, {
   3.789367595227816*^9, 3.789367599727377*^9}, {3.7893692710480957`*^9, 
   3.789369276931369*^9}, {3.789369911417919*^9, 3.789369916941128*^9}, {
   3.789374205099449*^9, 3.789374210330028*^9}, {3.789374476142598*^9, 
   3.7893744816544447`*^9}, {3.789376514716964*^9, 3.789376519393956*^9}, {
   3.7893772000099277`*^9, 3.789377204079545*^9}, {3.7893778266183367`*^9, 
   3.7893778326140423`*^9}, {3.78937879189863*^9, 3.789378800335326*^9}, {
   3.794132720933292*^9, 3.7941327678879137`*^9}, {3.794134222243268*^9, 
   3.794134261790318*^9}, {3.794135738466804*^9, 3.794135744650042*^9}, {
   3.7941365060425262`*^9, 3.794136512311205*^9}, {3.794137191037779*^9, 
   3.794137204067133*^9}, {3.7941376436793633`*^9, 3.794137648099907*^9}, {
   3.794138326473496*^9, 3.794138330352372*^9}, {3.794138449806006*^9, 
   3.794138453738491*^9}, {3.794138559423852*^9, 3.794138563661417*^9}, {
   3.7941388139997673`*^9, 3.7941388182622538`*^9}, {3.794140718587081*^9, 
   3.794140724918885*^9}, {3.794141156572933*^9, 3.794141161486472*^9}, {
   3.794143628700659*^9, 3.79414363513277*^9}, {3.794144991276031*^9, 
   3.794144997370873*^9}, 3.7942104961332617`*^9, {3.7942114997199717`*^9, 
   3.794211530215673*^9}, {3.794211865978715*^9, 3.794211923222467*^9}, {
   3.7942122322339563`*^9, 3.794212237140997*^9}, {3.794213661985642*^9, 
   3.794213685426867*^9}, {3.794216133776197*^9, 3.794216139631188*^9}, {
   3.794218038873033*^9, 3.7942180433000183`*^9}, {3.794218532574658*^9, 
   3.79421853856763*^9}, {3.7942192685541973`*^9, 3.794219274486249*^9}, {
   3.794220171463725*^9, 3.794220177549725*^9}, {3.794221835949416*^9, 
   3.7942218407201242`*^9}, {3.794222852566451*^9, 3.794222857422144*^9}, {
   3.794223181869269*^9, 3.79422318751303*^9}, {3.794223464912013*^9, 
   3.79422346964543*^9}, {3.79422367289913*^9, 3.7942236773898773`*^9}, {
   3.7942259781170177`*^9, 3.794225983813163*^9}, {3.7942262353840637`*^9, 
   3.794226257761891*^9}, {3.7942267814801197`*^9, 3.794226787945284*^9}, {
   3.794227189917529*^9, 3.794227194706987*^9}, {3.794227250778002*^9, 
   3.794227256073056*^9}, {3.794227486150888*^9, 3.7942274898805437`*^9}, {
   3.7942276932264423`*^9, 3.7942276975562553`*^9}, {3.7942284968405237`*^9, 
   3.794228501571949*^9}, {3.794228947530232*^9, 3.7942289517410707`*^9}, {
   3.794229181009058*^9, 3.794229188198122*^9}, {3.794230048255272*^9, 
   3.794230076033687*^9}, {3.7942303676547956`*^9, 3.794230374709668*^9}, {
   3.79423043255062*^9, 3.7942304392586823`*^9}, {3.7942313726254663`*^9, 
   3.794231377379468*^9}, {3.7942314181287727`*^9, 3.794231428672526*^9}, {
   3.794231808661038*^9, 3.7942318144699707`*^9}, {3.794232603915566*^9, 
   3.7942326091203833`*^9}, {3.794234658049645*^9, 3.7942346627196417`*^9}, {
   3.7942350686885147`*^9, 3.794235073717525*^9}, {3.79423570481857*^9, 
   3.794235708967889*^9}, {3.7942359672098227`*^9, 3.794235974896982*^9}, {
   3.794239054901988*^9, 3.7942390607500772`*^9}, {3.794239488753891*^9, 
   3.794239495401754*^9}, {3.7942396928714113`*^9, 3.79423971399691*^9}, {
   3.794240317997746*^9, 3.7942403283954268`*^9}, {3.794241591659547*^9, 
   3.794241599365221*^9}, {3.7942422211243258`*^9, 3.794242232252483*^9}, {
   3.79424239416936*^9, 3.794242426342835*^9}, {3.794242628967066*^9, 
   3.794242634865918*^9}, {3.7942427094632683`*^9, 3.794242716409892*^9}, {
   3.7942430491152887`*^9, 3.7942430558000107`*^9}, {3.794243332637001*^9, 
   3.7942433390467443`*^9}, {3.794243697966158*^9, 3.7942437048657303`*^9}, {
   3.7942439310691357`*^9, 3.7942439408392897`*^9}, {3.794244352022147*^9, 
   3.7942443569585657`*^9}, {3.7942950531299677`*^9, 3.794295070264379*^9}, {
   3.7942954789432697`*^9, 3.7942954835742903`*^9}, {3.794304904741172*^9, 
   3.7943049103745117`*^9}, {3.794304996487574*^9, 3.794305000423792*^9}, {
   3.794305686480569*^9, 3.794305691347926*^9}, {3.794305899933074*^9, 
   3.79430590330938*^9}, {3.7954202033568583`*^9, 3.7954202092563868`*^9}},
 NumberMarks->False,
 CellLabel->"In[15]:=",ExpressionUUID->"13bab909-a43e-4633-9293-76f19acad6d1"],

Cell[BoxData[
 RowBox[{"l", " ", "=", " ", 
  RowBox[{"Length", "[", "data", "]"}]}]], "Input",
 CellChangeTimes->{{3.78809452749936*^9, 3.7880945451573577`*^9}, {
   3.788094599585477*^9, 3.7880946153997498`*^9}, {3.7880946971902227`*^9, 
   3.7880948965254707`*^9}, {3.7880949445786114`*^9, 
   3.7880949696752143`*^9}, {3.7880950047106576`*^9, 
   3.7880950068669577`*^9}, {3.7880950618683343`*^9, 
   3.7880951113117094`*^9}, {3.788095160801737*^9, 3.7880951710060263`*^9}, {
   3.7881527944400597`*^9, 3.7881528410388613`*^9}, 3.788152895639106*^9, 
   3.7881529733877316`*^9, 3.7881757098330264`*^9, {3.788261501930628*^9, 
   3.78826150252286*^9}, 3.789279393279681*^9},
 CellLabel->"In[8]:=",ExpressionUUID->"6beb66a9-547a-48e7-8249-6fa038f64fa7"],

Cell[BoxData["23294"], "Input",
 CellChangeTimes->{{3.794288882580977*^9, 
  3.794288884148926*^9}},ExpressionUUID->"84312b97-f8a8-4748-8a3b-\
34a27903b521"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"mixtures", " ", "=", " ", 
    RowBox[{"Array", "[", 
     RowBox[{"Null", ",", " ", "l"}], "]"}]}], ";"}], " "}]], "Input",
 CellChangeTimes->{{3.78809452749936*^9, 3.7880945451573577`*^9}, {
   3.788094599585477*^9, 3.7880946153997498`*^9}, {3.7880946971902227`*^9, 
   3.7880948965254707`*^9}, {3.7880949445786114`*^9, 
   3.7880949696752143`*^9}, {3.7880950047106576`*^9, 
   3.7880950068669577`*^9}, {3.7880950618683343`*^9, 
   3.7880951113117094`*^9}, {3.788095160801737*^9, 3.7880951710060263`*^9}, {
   3.7881527944400597`*^9, 3.7881528410388613`*^9}, 3.788152895639106*^9, 
   3.7881529733877316`*^9, 3.7881757098330264`*^9},
 CellLabel->"In[9]:=",ExpressionUUID->"81fa31a9-cf0c-4d20-8605-45a46c137fd9"],

Cell[BoxData[
 RowBox[{"mixtures", "[", 
  RowBox[{"[", "4001", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.789281775793973*^9, 3.789281794983157*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"e7d6aac5-7f18-45d9-b391-364f1e417212"],

Cell[BoxData[
 RowBox[{"FindDistribution", "[", 
  RowBox[{
   RowBox[{"data", "[", 
    RowBox[{"[", "2", "]"}], "]"}], "[", 
   RowBox[{"[", 
    RowBox[{"2", ";;"}], "]"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.798297197309661*^9, 3.798297213630362*^9}, {
  3.798297281612336*^9, 3.798297307578796*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"66c43638-3b3a-4012-8360-263557794c13"],

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", " ", "=", " ", "2"}], ",", " ", 
    RowBox[{"i", " ", "<=", " ", "l"}], ",", " ", 
    RowBox[{"i", "++"}], ",", " ", 
    RowBox[{"(", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Head", " ", "@", " ", 
          RowBox[{"mixtures", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], " ", "\[Equal]", "  ", "Null"}], 
        ",", "\[IndentingNewLine]", 
        RowBox[{"(", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"Catch", "[", 
           RowBox[{
            RowBox[{"mixtures", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", 
            RowBox[{"GetMixture", "[", 
             RowBox[{"data", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "]"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"Mod", "[", 
              RowBox[{"i", ",", "100"}], "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{"Export", "[", 
             RowBox[{
              RowBox[{"\"\<~/mixtures.\>\"", "<>", 
               RowBox[{"ToString", "[", "i", "]"}], "<>", "\"\<.mx\>\""}], 
              ",", "mixtures", ",", "\"\<MX\>\""}], "]"}]}], "]"}], ";"}], 
         "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]", "]"}], ";"}], 
     "\[IndentingNewLine]", ")"}]}], "]"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.78809452749936*^9, 3.7880945451573577`*^9}, {
   3.788094599585477*^9, 3.7880946153997498`*^9}, {3.7880946971902227`*^9, 
   3.7880948965254707`*^9}, {3.7880949445786114`*^9, 
   3.7880949696752143`*^9}, {3.7880950047106576`*^9, 
   3.7880950068669577`*^9}, {3.7880950618683343`*^9, 
   3.7880951113117094`*^9}, {3.788095160801737*^9, 3.7880951710060263`*^9}, {
   3.7881527944400597`*^9, 3.7881528410388613`*^9}, 3.788152895639106*^9, {
   3.7881529733877316`*^9, 3.788152995109227*^9}, {3.789276554828431*^9, 
   3.78927658803645*^9}, {3.789276991442972*^9, 3.7892769923891687`*^9}, {
   3.7892818603477182`*^9, 3.789281900870123*^9}, 3.794221857506715*^9},
 CellLabel->"In[16]:=",ExpressionUUID->"9b24cbe7-e5a8-4a81-a7ef-0865d25bd28c"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.794244657616744*^9, 3.794244693727565*^9}},
 CellLabel->
  "In[121]:=",ExpressionUUID->"8c9530cd-b15c-4e3f-bec6-11aa140cd24d"],

Cell[BoxData["mixtures"], "Input",
 CellChangeTimes->{{3.7941341896339293`*^9, 3.794134190967795*^9}},
 CellLabel->
  "In[115]:=",ExpressionUUID->"3c22d126-2422-483a-a362-c6e95464878a"],

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{
  "\"\<mixtures_alveolar_tpm.mx\>\"", ",", "mixtures", ",", "\"\<MX\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.7880949776295104`*^9, 3.7880949788015213`*^9}, 
   3.788152236385627*^9, {3.7881743439301453`*^9, 3.788174362885543*^9}, {
   3.7881774335209293`*^9, 3.7881774367087946`*^9}, {3.788177497640091*^9, 
   3.7881775012186236`*^9}, {3.789275827226618*^9, 3.7892758348651953`*^9}, {
   3.789290547199996*^9, 3.789290574841579*^9}, 3.789379495892942*^9, {
   3.794145077994763*^9, 3.7941450875831003`*^9}, {3.794244699857452*^9, 
   3.7942447068635273`*^9}, {3.7954241022737103`*^9, 3.7954241056787977`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"f29668c6-7fbf-4432-81bc-002cca10b38b"],

Cell[BoxData[
 RowBox[{"mixtures", " ", "=", " ", 
  RowBox[{
  "Import", "[", "\"\<mixtures_dopaminergic_tpm.mx\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.7881789889345293`*^9, 3.7881790080148277`*^9}, {
   3.789384729681984*^9, 3.789384731751544*^9}, 3.794145234413423*^9, {
   3.794244716604972*^9, 3.794244720416848*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"a1c4c55a-0aab-48a6-aa63-7f725767072a"],

Cell[BoxData[
 RowBox[{"df", " ", "=", " ", 
  RowBox[{"Transpose", "[", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"Transpose", "[", "data", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}], ",", " ", "mixtures"}], "}"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.7879144745799074`*^9, 3.787914540353039*^9}, {
   3.7879150141942067`*^9, 3.787915021835552*^9}, {3.787915855274321*^9, 
   3.787915878527091*^9}, 3.7879159914773664`*^9, {3.7879172971022863`*^9, 
   3.7879173222147636`*^9}, {3.787917379283657*^9, 3.7879173816586576`*^9}, 
   3.788174530497143*^9, {3.788174563163684*^9, 3.7881745635699425`*^9}, {
   3.788174707148154*^9, 3.7881747285410805`*^9}, {3.788174764199482*^9, 
   3.788174803125371*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"bde6c8c1-95a5-4f83-974d-35c57931d7ae"],

Cell[BoxData[
 RowBox[{"filteredNames", " ", "=", " ", 
  RowBox[{
   RowBox[{"Transpose", "[", "filtered", "]"}], "[", 
   RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7882516917209024`*^9, 3.7882516928460693`*^9}, {
  3.7882517308820777`*^9, 3.7882517391796827`*^9}, {3.788251805967534*^9, 
  3.7882518311424017`*^9}, {3.788252235052998*^9, 3.788252246319801*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"e4fbad93-f975-4f15-a229-c4dc378a9b18"],

Cell[BoxData[
 RowBox[{"df1", " ", "=", " ", "df", 
  RowBox[{"(*", 
   RowBox[{"Select", "[", 
    RowBox[{"df", ",", " ", 
     RowBox[{
      RowBox[{"MemberQ", "[", 
       RowBox[{"filteredNames", ",", " ", 
        RowBox[{"#", "[", 
         RowBox[{"[", "1", "]"}], "]"}]}], "]"}], " ", "&"}]}], "]"}], 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.788251848535366*^9, 3.7882518895400887`*^9}, {
  3.788251970567692*^9, 3.788252017353927*^9}, {3.7882521072591763`*^9, 
  3.788252110962721*^9}, {3.788252206924852*^9, 3.788252209847063*^9}, {
  3.7882522515392423`*^9, 3.788252256633566*^9}, {3.7882523793502035`*^9, 
  3.788252386179367*^9}, {3.7882525656840787`*^9, 3.7882525775132647`*^9}, {
  3.794145451264677*^9, 3.7941454531504498`*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"ab98ed2c-bae2-4e59-a231-1315d4250761"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GetAntimode", "[", 
   RowBox[{"name_", ",", "dist_", ",", " ", "orig_"}], "]"}], " ", ":=", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "i", ",", " ", "pdf", ",", " ", "dPdf", ",", " ", "ddPdf", ",", " ", 
      "zeros"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"dist", " ", "\[Equal]", " ", 
          RowBox[{"-", "1"}]}], ")"}], " ", "||", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Length", "[", 
           RowBox[{"dist", "[", 
            RowBox[{"[", 
             RowBox[{"-", "1"}], "]"}], "]"}], "]"}], " ", "\[Equal]", " ", 
          "0"}], ")"}]}], ",", 
       RowBox[{"Return", "[", "\"\<NotMixture\>\"", "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"name", "\[Equal]", "\"\<Diras1\>\""}], ")"}], ",", 
         RowBox[{"Return", "[", "\"\<HandRemoved\>\"", "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", 
       RowBox[{"Somatosensory", " ", "TPM"}], "*)"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"name", "\[Equal]", "\"\<Dennd1c\>\""}], "||", 
          RowBox[{"name", "\[Equal]", "\"\<Arid5a\>\""}], "||", 
          RowBox[{"name", "\[Equal]", "\"\<Gm11574\>\""}]}], ",", 
         RowBox[{"Return", "[", "\"\<HandRemoved\>\"", "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", 
       RowBox[{"Dopaminergic", " ", "TPM"}], "*)"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"name", " ", "\[Equal]", " ", "\"\<2610204G22Rik\>\""}], 
           " ", "||", " ", 
           RowBox[{"name", "==", "\"\<4931414P19Rik\>\""}]}], ")"}], ",", 
         RowBox[{"Return", "[", "\"\<HandRemoved\>\"", "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", "somatosensory", "*)"}], "*)"}], "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"name", " ", "\[Equal]", " ", "\"\<1700109H08Rik\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<4731419I09Rik\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Adgre1\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Gm10644\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Gm15601\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Gm29585\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Gm5946\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Atp6v1fnb\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Gsx1\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Islr\>\""}], "||", 
           RowBox[{"name", "\[Equal]", "\"\<Rarres2\>\""}]}], ")"}], ",", 
         RowBox[{"Return", "[", "\"\<HandRemoved\>\"", "]"}]}], "]"}], ";"}], 
      RowBox[{"(*", "dopaminergic", "*)"}], "*)"}], "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{"Print", "[", "name", "]"}], ";"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "UniformDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsUniform\>\"", "]"}]}], "]"}], 
         ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "DiscreteUniformDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsDiscreteUniform\>\"", "]"}]}], 
          "]"}], ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "LaplaceDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsLaplace\>\"", "]"}]}], "]"}], 
         ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "ZipfDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsZipf\>\"", "]"}]}], "]"}], 
         ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "LogSeriesDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsLogSeries\>\"", "]"}]}], "]"}],
          ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", " ", 
       RowBox[{"i", " ", "\[LessEqual]", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"dist", "[", 
          RowBox[{"[", 
           RowBox[{"-", "1"}], "]"}], "]"}], "]"}]}], ",", " ", 
       RowBox[{"i", "++"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Head", "  ", "@", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"dist", "[", 
                RowBox[{"[", 
                 RowBox[{"-", "1"}], "]"}], "]"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", " ", 
            "BenfordDistribution"}], ",", " ", 
           RowBox[{"Return", "[", "\"\<ContainsBenford\>\"", "]"}]}], "]"}], 
         ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"pdf", " ", "=", " ", 
      RowBox[{"PDF", "[", 
       RowBox[{"dist", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"dPdf", " ", "=", " ", 
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"PiecewiseExpand", "[", 
         RowBox[{"pdf", ",", " ", 
          RowBox[{"x", ">", "0"}]}], "]"}], ",", "x"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"ddPdf", " ", "=", " ", 
      RowBox[{"D", "[", 
       RowBox[{"dPdf", ",", " ", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{"Print", "[", 
       RowBox[{"Max", "[", 
        RowBox[{"orig", "[", 
         RowBox[{"[", 
          RowBox[{"2", ";;"}], "]"}], "]"}], "]"}], "]"}], "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"zeros", " ", "=", " ", 
      RowBox[{"NSolve", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"dPdf", " ", "\[Equal]", " ", "0"}], ",", " ", 
          RowBox[{"0", "<", "x", "<", 
           RowBox[{"Max", "[", 
            RowBox[{"orig", "[", 
             RowBox[{"[", 
              RowBox[{"2", ";;"}], "]"}], "]"}], "]"}]}]}], "}"}], ",", "x", 
        ",", "Reals"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"evalDdPdf", "[", "x_", "]"}], " ", ":=", " ", 
      RowBox[{"(", "\[IndentingNewLine]", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"ddPdf", "/.", "x"}], ")"}], " ", ">", "0"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"x", ",", " ", "\"\<Antimode\>\""}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"x", ",", " ", "\"\<Mode\>\""}], "}"}]}], "]"}], 
       "\[IndentingNewLine]", ")"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Map", "[", 
      RowBox[{"evalDdPdf", ",", " ", "zeros"}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.7882601590202227`*^9, 3.788260172639021*^9}, {
   3.78826026120148*^9, 3.7882603074627943`*^9}, {3.788260495738029*^9, 
   3.788260640570622*^9}, {3.788260705304728*^9, 3.78826076371565*^9}, {
   3.7882611270505943`*^9, 3.788261127218686*^9}, {3.788261195743717*^9, 
   3.788261330497423*^9}, {3.788261397740982*^9, 3.788261401430027*^9}, {
   3.788261699530663*^9, 3.78826170220295*^9}, {3.788261769618223*^9, 
   3.788261786725203*^9}, {3.788263308668869*^9, 3.78826337778834*^9}, {
   3.788264139244886*^9, 3.788264143113503*^9}, {3.788264438251318*^9, 
   3.788264450962438*^9}, {3.788511209284741*^9, 3.78851122202367*^9}, {
   3.789379607665881*^9, 3.7893796330640097`*^9}, {3.7893810846225433`*^9, 
   3.789381124185883*^9}, 3.789381314039143*^9, 3.7893813824414167`*^9, 
   3.789384655122005*^9, {3.789387790505979*^9, 3.789387816763117*^9}, {
   3.789387849615203*^9, 3.789387856083931*^9}, {3.789389152968321*^9, 
   3.789389167028751*^9}, {3.7893906554740953`*^9, 3.789390666643969*^9}, {
   3.789448921099037*^9, 3.789448926948511*^9}, {3.789456143364398*^9, 
   3.789456147032933*^9}, {3.78945646716094*^9, 3.789456552976347*^9}, {
   3.789456868011952*^9, 3.7894568897690077`*^9}, {3.78945772887554*^9, 
   3.789457739134397*^9}, 3.789457897699584*^9, {3.7894581299605837`*^9, 
   3.7894582230686827`*^9}, {3.789458254749028*^9, 3.789458290558093*^9}, {
   3.789458394774646*^9, 3.789458398955556*^9}, {3.789458538617429*^9, 
   3.789458551073876*^9}, {3.7894631304573107`*^9, 3.789463133865384*^9}, {
   3.789463205046129*^9, 3.789463220457184*^9}, {3.789463952226719*^9, 
   3.789463955567733*^9}, {3.789464081930838*^9, 3.789464086239688*^9}, {
   3.78946476751584*^9, 3.789464772995665*^9}, {3.789465323722101*^9, 
   3.7894653275813017`*^9}, {3.78946754439487*^9, 3.789467548589781*^9}, 
   3.794145346125449*^9, {3.794146029812702*^9, 3.79414604064688*^9}, {
   3.794150002740993*^9, 3.794150025921131*^9}, 3.794244796472526*^9, {
   3.7942455815713453`*^9, 3.794245615223127*^9}, {3.794253886132657*^9, 
   3.794253894659917*^9}, {3.794289150179729*^9, 3.79428918364882*^9}, {
   3.7954241646867123`*^9, 3.795424180250239*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"1ac8b741-f57e-4071-83f1-860cd500bc85"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"antimodes", " ", "=", " ", 
    RowBox[{"Array", "[", 
     RowBox[{"Null", ",", " ", "l"}], "]"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"antimodes", "[", 
   RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", 
  "\"\<antimode\>\""}]}], "Input",
 CellChangeTimes->{{3.788261521362976*^9, 3.7882615251696787`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"dcf64c38-29f3-4915-8365-4546509ce467"],

Cell[BoxData["antimodes"], "Input",
 CellChangeTimes->{{3.794245629675312*^9, 3.794245632310367*^9}},
 CellLabel->
  "In[134]:=",ExpressionUUID->"5969ca59-54f4-414c-9d00-98f5cf6dd2e6"],

Cell[BoxData[
 RowBox[{"antimodes", "=", 
  RowBox[{
  "Import", "[", "\"\<~/antimodes.somatosensory_tpm.20500.mx\>\"", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.788497027149952*^9, 3.788497028684948*^9}, {
  3.788511276553288*^9, 3.788511294833666*^9}, {3.789384663722941*^9, 
  3.789384686021525*^9}, {3.7893879024127483`*^9, 3.789387902770458*^9}, {
  3.789463965511181*^9, 3.789463966071065*^9}, {3.7894641263454227`*^9, 
  3.789464126433343*^9}, {3.789469672979167*^9, 3.789469673934841*^9}, {
  3.7941461127538643`*^9, 3.794146138823436*^9}, {3.794253744770978*^9, 
  3.794253745618363*^9}, {3.794288853656247*^9, 3.7942888544587803`*^9}},
 CellLabel->"In[16]:=",ExpressionUUID->"e50a75cd-5264-4c3a-978b-25e4291eb9f6"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"l", "=", 
   RowBox[{"Length", "[", "df1", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"to", " ", "=", " ", "l"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ProgressIndicator", "[", 
  RowBox[{
   RowBox[{"Dynamic", "[", "i", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "to"}], "}"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", " ", "=", " ", "2"}], ",", " ", 
   RowBox[{"i", " ", "<=", " ", "to"}], ",", " ", 
   RowBox[{"i", "++"}], ",", " ", 
   RowBox[{"(", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Head", " ", "@", " ", 
         RowBox[{"(", 
          RowBox[{"antimodes", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ")"}]}], " ", "\[Equal]", "  ", 
        "Null"}], ",", 
       RowBox[{"(", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Print", "[", 
          RowBox[{"df1", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "]"}], ";", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"antimodes", "[", 
           RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", 
          RowBox[{"GetAntimode", "[", 
           RowBox[{
            RowBox[{"df1", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", " ", 
            RowBox[{"df1", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", " ", 
            RowBox[{"data", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mod", "[", 
             RowBox[{"i", ",", "100"}], "]"}], "\[Equal]", "0"}], ",", 
           RowBox[{"Export", "[", 
            RowBox[{
             RowBox[{"\"\<~/antimodes.somatosensory_tpm.\>\"", "<>", 
              RowBox[{"ToString", "[", "i", "]"}], "<>", "\"\<.mx\>\""}], ",",
              "antimodes", ",", "\"\<MX\>\""}], "]"}]}], "]"}], ";"}], 
        "\[IndentingNewLine]", ")"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
    ")"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.788261531520117*^9, 3.788261617130538*^9}, {
   3.788261713113068*^9, 3.788261714532525*^9}, {3.788262754486804*^9, 
   3.788262790270927*^9}, {3.788267874246996*^9, 3.7882678855997467`*^9}, {
   3.788268002208952*^9, 3.788268036070548*^9}, {3.788268173141788*^9, 
   3.7882681903837843`*^9}, 3.788268409238207*^9, 3.788268475828014*^9, {
   3.788268517352783*^9, 3.788268558201249*^9}, {3.788497192962968*^9, 
   3.788497196142601*^9}, {3.78851162090241*^9, 3.78851163218147*^9}, {
   3.788516059015703*^9, 3.7885160669930143`*^9}, {3.788517612850062*^9, 
   3.788517635565118*^9}, {3.789379790189185*^9, 3.789379794786984*^9}, {
   3.78938113463583*^9, 3.7893811461343412`*^9}, {3.789381276515148*^9, 
   3.789381284749422*^9}, {3.789381323422073*^9, 3.789381376282324*^9}, {
   3.789448998937911*^9, 3.789449010491034*^9}, {3.789449060928458*^9, 
   3.789449066731633*^9}, 3.7894564056102343`*^9, 3.789469641998769*^9, {
   3.794145393752687*^9, 3.794145400735819*^9}, {3.79414543689883*^9, 
   3.79414543884157*^9}, {3.7941455129045258`*^9, 3.794145513241947*^9}, {
   3.794145716075446*^9, 3.794145728845853*^9}, 3.794245639044297*^9},
 CellLabel->"In[16]:=",ExpressionUUID->"a5db6e22-da78-4761-aef5-2b5a4b2d82de"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.794145737645195*^9, 3.794145738934997*^9}, 
   3.794145780376904*^9, {3.7942538262069798`*^9, 
   3.794253827856695*^9}},ExpressionUUID->"235fecd7-4651-4bb3-8a1f-\
8d23103e1674"],

Cell[BoxData["antimodes"], "Input",
 CellChangeTimes->{
  3.794149703043024*^9, {3.794151285163306*^9, 3.79415128877805*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"060dc19b-c311-471b-84bb-db561e8c359f"],

Cell[BoxData["Ftdc"], "Input",ExpressionUUID->"077f37bd-f141-4a7a-b0f7-3793413c3e23"],

Cell[BoxData[
 RowBox[{"antimodes", " ", "=", " ", 
  RowBox[{"antimodes", "[", 
   RowBox[{"[", 
    RowBox[{";;", "l"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.789449385392745*^9, 3.7894493866895638`*^9}, {
  3.789469727072476*^9, 3.789469768775161*^9}, {3.794291348181028*^9, 
  3.7942913484337606`*^9}},
 CellLabel->"In[24]:=",ExpressionUUID->"6abcc5b8-de3e-4e96-a779-8a7ba72bd814"],

Cell[BoxData["l"], "Input",
 CellChangeTimes->{{3.789469755012299*^9, 3.789469757446239*^9}},
 CellLabel->"In[16]:=",ExpressionUUID->"2d046b52-6363-45b9-8b51-755cf5627672"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"PDF", "[", 
    RowBox[{
     RowBox[{"MixtureDistribution", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0.9105960267896854`", ",", "0.08940397321031447`"}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"GeometricDistribution", "[", "0.9892086298770227`", "]"}], 
         ",", 
         RowBox[{"PascalDistribution", "[", 
          RowBox[{"2", ",", "0.033707865074614964`"}], "]"}]}], "}"}]}], 
      "]"}], ",", "x"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0.0", ",", "10.0"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.789380981110139*^9, 3.789381050246477*^9}, {
   3.7893890132270927`*^9, 3.789389061411521*^9}, {3.789389129892315*^9, 
   3.789389130508493*^9}, {3.7893902287479486`*^9, 3.7893902801472*^9}, 
   3.789456430911807*^9, {3.7894646632134953`*^9, 3.789464668756515*^9}},
 CellLabel->"In[15]:=",ExpressionUUID->"c6680cea-3b8f-4165-955b-663da255470b"],

Cell[BoxData["\"\<1700109H08Rik\>\""], "Input",ExpressionUUID->"a8e5dc2d-7ca8-4eef-ba05-7935e3fb3ed0"],

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"\"\<~/antimodes_.mx\>\"", ",", "antimodes", ",", "\"\<MX\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.788261649287261*^9, 3.788261650477108*^9}, {
  3.7884970434688797`*^9, 3.788497056608378*^9}, {3.7885160989106903`*^9, 
  3.7885161210163927`*^9}, {3.789384874271319*^9, 3.789384876779723*^9}, {
  3.795427443358508*^9, 
  3.7954274457840147`*^9}},ExpressionUUID->"59b7e0a2-299f-4f16-b46c-\
f97f85b2a289"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"thresholds", " ", "=", " ", 
    RowBox[{"Array", "[", 
     RowBox[{"Null", ",", " ", "l"}], "]"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"thresholds", "[", 
   RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", 
  "\"\<thresholds\>\""}]}], "Input",
 CellChangeTimes->{{3.788516322765834*^9, 3.7885163292237253`*^9}},
 CellLabel->"In[21]:=",ExpressionUUID->"eac08b2e-77d5-41b4-830f-61dc9827df07"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GetThreshold", "[", "modes_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"i", ",", "j"}], "}"}], ",", 
    RowBox[{"(", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"FoundMode", "=", "False"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "1"}], ",", 
        RowBox[{"i", "\[LessEqual]", 
         RowBox[{"Length", "[", "modes", "]"}]}], ",", 
        RowBox[{"i", "++"}], ",", 
        RowBox[{"(", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"modes", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "\[Equal]", 
             "\"\<Antimode\>\""}], ",", " ", 
            RowBox[{"Return", "[", 
             RowBox[{"Values", "[", 
              RowBox[{"modes", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1", ",", " ", "1"}], "]"}], "]"}], "]"}], 
             "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", ")"}]}], "]"}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", "\"\<NotFound\>\"", "]"}]}], ")"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.788516388572351*^9, 3.788516633483386*^9}, {
  3.788516720235239*^9, 3.788516745973778*^9}, {3.788516809223216*^9, 
  3.788516810549872*^9}, {3.788516843064136*^9, 3.788516847765389*^9}},
 CellLabel->"In[23]:=",ExpressionUUID->"ea6e3fb0-8ea5-49fd-b258-a13c97b84c3b"],

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "0"}], ",", 
   RowBox[{"i", "\[LessEqual]", "to"}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"(", "\[IndentingNewLine]", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Head", "[", 
        RowBox[{"antimodes", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", "List"}], ",", 
      RowBox[{"(", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"thresholds", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "=", 
         RowBox[{"GetThreshold", "[", 
          RowBox[{"antimodes", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";"}], 
       "\[IndentingNewLine]", ")"}]}], "]"}], "\[IndentingNewLine]", ")"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.788516153339023*^9, 3.7885162180263853`*^9}, {
  3.788516277427536*^9, 3.78851629753091*^9}, {3.788516336547123*^9, 
  3.78851637155704*^9}, {3.788516639614751*^9, 3.788516662700115*^9}},
 NumberMarks->False,
 CellLabel->"In[24]:=",ExpressionUUID->"3f24a9db-3250-4060-b476-de1733acf562"],

Cell[BoxData["thresholds"], "Input",
 CellChangeTimes->{{3.788516673483222*^9, 3.788516674200678*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"8eade957-96e8-41aa-9261-bfb58d46b7af"],

Cell[BoxData[{
 RowBox[{"results", "=", 
  RowBox[{"Transpose", "[", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"Transpose", "[", "df1", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}], ",", 
     RowBox[{
      RowBox[{"Transpose", "[", "df1", "]"}], "[", 
      RowBox[{"[", "2", "]"}], "]"}], ",", "antimodes", ",", "thresholds"}], 
    "}"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"results", "[", 
   RowBox[{"[", 
    RowBox[{"1", ",", "2"}], "]"}], "]"}], "=", 
  "\"\<distribution\>\""}]}], "Input",
 CellChangeTimes->{
  3.7882532813289185`*^9, {3.7882545452208614`*^9, 3.788254588553856*^9}, {
   3.7882546575459795`*^9, 3.7882546632809973`*^9}, {3.7882546977536244`*^9, 
   3.7882547048947773`*^9}, {3.7882548430077257`*^9, 
   3.7882548678386393`*^9}, {3.788256541669841*^9, 3.788256545154336*^9}, {
   3.7882589265813932`*^9, 3.788258927422812*^9}, {3.788258970759508*^9, 
   3.788258994743943*^9}, {3.788259152180346*^9, 3.788259152280137*^9}, {
   3.788516774089819*^9, 3.788516834020885*^9}, 3.788516880526482*^9, {
   3.7885169831386337`*^9, 3.788517101928709*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"b58a7c3b-8823-4afb-a15e-b4f8a47d190f"],

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"\"\<mathematica_results_alveolar_tpm.csv\>\"", ",", "results"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.788517108200674*^9, 3.788517126195321*^9}, {
  3.789469839731739*^9, 3.789469845896657*^9}, {3.7941513225659246`*^9, 
  3.794151331889326*^9}, {3.794291414503668*^9, 3.794291418602789*^9}, {
  3.795427475689786*^9, 3.7954274827422867`*^9}},
 CellLabel->"In[28]:=",ExpressionUUID->"c54e0542-f7d5-4198-b657-49e710380d85"],

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{
  "\"\<mathematica_results_alveolar_tpm.mx\>\"", ",", "results", ",", 
   "\"\<MX\>\""}], "]"}]], "Input",
 CellChangeTimes->{{3.788517258221339*^9, 3.788517275091715*^9}, {
  3.78946985227606*^9, 3.7894698551443987`*^9}, {3.7941513339217863`*^9, 
  3.794151340870624*^9}, {3.7942914242368517`*^9, 3.794291426230382*^9}, {
  3.795427488634961*^9, 3.795427492962851*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"4a076f87-1999-4cb5-9483-a543122b6ed8"],

Cell[BoxData["results"], "Input",
 CellChangeTimes->{{3.7942914460458517`*^9, 3.794291447101643*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"44bdf3c2-1c33-4b4d-b4ed-b0aa0f171f4d"]
},
WindowSize->{1437., 777.75},
WindowMargins->{{1.5, Automatic}, {Automatic, 15}},
FrontEndVersion->"12.1 for Linux x86 (64-bit) (June 19, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"49296880-46d6-4d2e-adbe-0f9e9fdc3ab0"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 345, 6, 51, "Input",ExpressionUUID->"fe31bd78-5fbf-4a06-86f0-b80159ec6f68"],
Cell[906, 28, 143, 2, 29, "Input",ExpressionUUID->"2c8e21e4-f654-4aac-8fa2-d6a56104d0b2"],
Cell[1052, 32, 711, 13, 29, "Input",ExpressionUUID->"5e4db062-34b8-499a-b4d1-fd1bbb88c8f8"],
Cell[1766, 47, 277, 5, 29, "Input",ExpressionUUID->"64dbf530-b3fc-4d33-b3d8-da822f5c311b"],
Cell[2046, 54, 154, 3, 29, "Input",ExpressionUUID->"b41e38b3-a239-42de-9ac0-71a31d034b50"],
Cell[2203, 59, 349, 7, 29, "Input",ExpressionUUID->"7d375fe9-4815-417c-81d9-186f56f11dfe"],
Cell[2555, 68, 26769, 566, 484, "Input",ExpressionUUID->"13bab909-a43e-4633-9293-76f19acad6d1"],
Cell[29327, 636, 755, 12, 29, "Input",ExpressionUUID->"6beb66a9-547a-48e7-8249-6fa038f64fa7"],
Cell[30085, 650, 157, 3, 29, "Input",ExpressionUUID->"84312b97-f8a8-4748-8a3b-34a27903b521"],
Cell[30245, 655, 760, 14, 31, "Input",ExpressionUUID->"81fa31a9-cf0c-4d20-8605-45a46c137fd9"],
Cell[31008, 671, 231, 4, 31, "Input",ExpressionUUID->"e7d6aac5-7f18-45d9-b391-364f1e417212"],
Cell[31242, 677, 393, 9, 31, "Input",ExpressionUUID->"66c43638-3b3a-4012-8360-263557794c13"],
Cell[31638, 688, 2253, 48, 216, "Input",ExpressionUUID->"9b24cbe7-e5a8-4a81-a7ef-0865d25bd28c"],
Cell[33894, 738, 175, 3, 31, "Input",ExpressionUUID->"8c9530cd-b15c-4e3f-bec6-11aa140cd24d"],
Cell[34072, 743, 185, 3, 31, "Input",ExpressionUUID->"3c22d126-2422-483a-a362-c6e95464878a"],
Cell[34260, 748, 748, 12, 31, "Input",ExpressionUUID->"f29668c6-7fbf-4432-81bc-002cca10b38b"],
Cell[35011, 762, 408, 7, 31, "Input",ExpressionUUID->"a1c4c55a-0aab-48a6-aa63-7f725767072a"],
Cell[35422, 771, 813, 16, 31, "Input",ExpressionUUID->"bde6c8c1-95a5-4f83-974d-35c57931d7ae"],
Cell[36238, 789, 465, 8, 31, "Input",ExpressionUUID->"e4fbad93-f975-4f15-a229-c4dc378a9b18"],
Cell[36706, 799, 830, 17, 31, "Input",ExpressionUUID->"ab98ed2c-bae2-4e59-a231-1315d4250761"],
Cell[37539, 818, 13252, 310, 906, "Input",ExpressionUUID->"1ac8b741-f57e-4071-83f1-860cd500bc85"],
Cell[50794, 1130, 462, 12, 55, "Input",ExpressionUUID->"dcf64c38-29f3-4915-8365-4546509ce467"],
Cell[51259, 1144, 184, 3, 31, "Input",ExpressionUUID->"5969ca59-54f4-414c-9d00-98f5cf6dd2e6"],
Cell[51446, 1149, 727, 12, 31, "Input",ExpressionUUID->"e50a75cd-5264-4c3a-978b-25e4291eb9f6"],
Cell[52176, 1163, 3407, 73, 239, "Input",ExpressionUUID->"a5db6e22-da78-4761-aef5-2b5a4b2d82de"],
Cell[55586, 1238, 227, 4, 31, "Input",ExpressionUUID->"235fecd7-4651-4bb3-8a1f-8d23103e1674"],
Cell[55816, 1244, 204, 3, 31, "Input",ExpressionUUID->"060dc19b-c311-471b-84bb-db561e8c359f"],
Cell[56023, 1249, 85, 0, 31, "Input",ExpressionUUID->"077f37bd-f141-4a7a-b0f7-3793413c3e23"],
Cell[56111, 1251, 395, 8, 31, "Input",ExpressionUUID->"6abcc5b8-de3e-4e96-a779-8a7ba72bd814"],
Cell[56509, 1261, 172, 2, 31, "Input",ExpressionUUID->"2d046b52-6363-45b9-8b51-755cf5627672"],
Cell[56684, 1265, 995, 23, 34, "Input",ExpressionUUID->"c6680cea-3b8f-4165-955b-663da255470b"],
Cell[57682, 1290, 102, 0, 31, "Input",ExpressionUUID->"a8e5dc2d-7ca8-4eef-ba05-7935e3fb3ed0"],
Cell[57787, 1292, 461, 9, 31, "Input",ExpressionUUID->"59b7e0a2-299f-4f16-b46c-f97f85b2a289"],
Cell[58251, 1303, 466, 12, 55, "Input",ExpressionUUID->"eac08b2e-77d5-41b4-830f-61dc9827df07"],
Cell[58720, 1317, 1533, 37, 147, "Input",ExpressionUUID->"ea6e3fb0-8ea5-49fd-b258-a13c97b84c3b"],
Cell[60256, 1356, 1094, 27, 124, "Input",ExpressionUUID->"3f24a9db-3250-4060-b476-de1733acf562"],
Cell[61353, 1385, 181, 2, 31, "Input",ExpressionUUID->"8eade957-96e8-41aa-9261-bfb58d46b7af"],
Cell[61537, 1389, 1193, 26, 55, "Input",ExpressionUUID->"b58a7c3b-8823-4afb-a15e-b4f8a47d190f"],
Cell[62733, 1417, 476, 8, 31, "Input",ExpressionUUID->"c54e0542-f7d5-4198-b657-49e710380d85"],
Cell[63212, 1427, 499, 9, 31, "Input",ExpressionUUID->"4a076f87-1999-4cb5-9483-a543122b6ed8"],
Cell[63714, 1438, 180, 2, 31, "Input",ExpressionUUID->"44bdf3c2-1c33-4b4d-b4ed-b0aa0f171f4d"]
}
]
*)