8_RABBIT_run_jointModel_OPD.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[      6950,        175]
NotebookOptionsPosition[      6657,        162]
NotebookOutlinePosition[      6996,        177]
CellTagsIndexPosition[      6953,        174]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "User", " ", "needs", " ", "to", " ", "set", " ", "directory", " ", 
     "where", " ", "RABBIT", " ", "packages", " ", "are", " ", "downloaded"}],
     ",", " ", 
    RowBox[{"pull", " ", "required", " ", "package"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
   "SetDirectory", "[", "\"\</directory/RABBIT/RABBIT_Packages\>\"", "]"}], 
   "\n", 
   StyleBox[
    RowBox[{"Needs", "[", "\"\<MagicReconstruct`\>\"", "]"}], "Input"], 
   StyleBox["\[IndentingNewLine]", "Input"], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "This", " ", "line", " ", "will", " ", "change", " ", "directory", " ", 
      "to", " ", "the", " ", "directory", " ", "where", " ", "this", " ", 
      "notebook", " ", "is", " ", "saved"}], ",", " ", 
     RowBox[{
     "needs", " ", "to", " ", "have", " ", "input", " ", "file", " ", "from", 
      " ", "4", "_MACH", "_to", "_RABBIT", 
      RowBox[{"_format", ".", "R"}]}]}], "*)"}], "\[IndentingNewLine]", 
   StyleBox[
    RowBox[{"SetDirectory", "[", 
     RowBox[{"NotebookDirectory", "[", "]"}], "]"}], "Input"], 
   StyleBox[" ", "Input"], 
   StyleBox["\[IndentingNewLine]", "Input"], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"model", " ", "=", " ", "\"\<jointModel\>\""}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"designate", " ", "model"}], ",", " ", 
     RowBox[{"here", " ", "using", " ", "joint", " ", "model"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"estfun", " ", "=", " ", "\"\<origPosteriorDecoding\>\""}], ";", 
    " ", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{"designate", " ", "algorithm"}], ",", " ", 
      RowBox[{"can", " ", "be", " ", "origViterbiDecoding"}], ",", " ", 
      "origPosteriorDecoding", ",", " ", 
      RowBox[{"or", " ", "origPathSampling"}]}], "*)"}], 
    "\[IndentingNewLine]", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{
      "run", " ", "magicReconstruct", " ", "for", " ", "each", " ", 
       "chromosome"}], ",", " ", 
      RowBox[{
      "write", " ", "raw", " ", "output", " ", "and", " ", "creates", " ", 
       "summary", " ", "file", " ", "that", " ", "will", " ", "be", " ", 
       "used", " ", "in", " ", "SPEARS", " ", "pipeline"}]}], "*)"}], 
    "\[IndentingNewLine]", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{"change", " ", "j"}], " ", "\[LessEqual]", " ", 
      RowBox[{
      "10", " ", "to", " ", "number", " ", "corresponding", " ", "to", " ", 
       "the", " ", "number", " ", "of", " ", "chromosomes", " ", "to", " ", 
       "run"}]}], "*)"}], "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", " ", "=", "1"}], ",", " ", 
      RowBox[{"j", " ", "\[LessEqual]", " ", "10"}], ",", " ", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"magicsnp", " ", "=", " ", 
        RowBox[{"Import", "[", 
         RowBox[{
          RowBox[{"\"\<filtered_chrom_\>\"", " ", "<>", " ", 
           RowBox[{"ToString", "[", "j", "]"}], " ", "<>", " ", 
           "\"\<_RABBIT_input.csv\>\""}], ",", "\"\<CSV\>\""}], "]"}]}], ";", 
       " ", "\[IndentingNewLine]", 
       RowBox[{"popdesign", " ", "=", " ", 
        RowBox[{"{", 
         RowBox[{"{", 
          RowBox[{
          "0.156656", ",", "0.000930313", ",", "0.582508", ",", "2.61506", 
           ",", "0.582507", ",", "2.61505"}], "}"}], "}"}]}], ";", " ", 
       RowBox[{"(*", 
        RowBox[{
         RowBox[{"popdesign", " ", "=", " ", 
          RowBox[{
          "prior", " ", "calculated", " ", "based", " ", "on", " ", 
           "specific", " ", "population", " ", "pedigree", " ", 
           RowBox[{"see", ":", " ", 
            RowBox[{"Zheng", " ", "et", " ", 
             RowBox[{"al", ".", " ", "2015"}]}]}]}]}], ",", " ", 
         RowBox[{
         "can", " ", "input", " ", "priors", " ", "or", " ", "designate", " ",
           "pedigree", " ", "here"}], ",", "\[IndentingNewLine]", 
         RowBox[{
         "the", " ", "priors", " ", "designated", " ", "above", " ", "are", 
          " ", "for", " ", "the", " ", "example", " ", "population", " ", 
          "provided"}]}], "*)"}], "\[IndentingNewLine]", 
       RowBox[{"magicReconstruct", "[", 
        RowBox[{"magicsnp", ",", "model", ",", "popdesign", ",", " ", 
         RowBox[{"reconstructAlgorithm", "\[Rule]", "estfun"}], ",", 
         RowBox[{"(*", 
          RowBox[{
          "Define", " ", "parameters", " ", "for", " ", "magicReconstruct"}], 
          "*)"}], "\[IndentingNewLine]", 
         RowBox[{"offspringAllelicError", " ", "\[Rule]", " ", "0.006"}], ",",
          " ", 
         RowBox[{"(*", 
          RowBox[{
           RowBox[{
           "define", " ", "error", " ", "rate", " ", "in", " ", "progeny"}], 
           ",", " ", 
           RowBox[{
           "match", " ", "what", " ", "was", " ", "used", " ", "for", " ", 
            "step", " ", "2"}]}], "*)"}], 
         RowBox[{"outputFileID", "\[Rule]", 
          RowBox[{"\"\<reconstructed_chrom_\>\"", " ", "<>", " ", 
           RowBox[{"ToString", "[", "j", "]"}], " ", "<>", " ", 
           "\"\<_jointModel_OVD\>\""}]}], ",", 
         RowBox[{"isPrintTimeElapsed", "\[Rule]", "True"}]}], "]"}], " ", ";",
        " ", 
       RowBox[{"(*", 
        RowBox[{"name", " ", "of", " ", "output"}], "*)"}], 
       "\[IndentingNewLine]", " ", 
       RowBox[{"saveAsSummaryMR", "[", 
        RowBox[{
         RowBox[{"\"\<reconstructed_chrom_\>\"", " ", "<>", " ", 
          RowBox[{"ToString", "[", "j", "]"}], " ", "<>", 
          "\"\<_jointModel_OPD_magicReconstruct.txt\>\""}], ",", 
         RowBox[{"\"\<reconstructed_chrom_\>\"", " ", "<>", " ", 
          RowBox[{"ToString", "[", "j", "]"}], " ", "<>", 
          "\"\<_jointModel_OPD_summary.csv\>\""}]}], "]"}]}]}], 
     RowBox[{"(*", 
      RowBox[{"output", " ", "summary"}], "*)"}], "\[IndentingNewLine]", 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.8028952741279373`*^9, 
  3.8028952868736486`*^9}},ExpressionUUID->"75aed170-5b34-48b2-8e87-\
03a2f4da5f0c"]
},
WindowSize->{808, 911},
WindowMargins->{{762, Automatic}, {Automatic, 202}},
FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 6095, 140, 661, "Input",ExpressionUUID->"75aed170-5b34-48b2-8e87-03a2f4da5f0c"]
}
]
*)