RDownScriptalphabeta@tmax.nb

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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    185742,       3614]
NotebookOptionsPosition[    178492,       3490]
NotebookOutlinePosition[    178893,       3507]
CellTagsIndexPosition[    178850,       3504]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["RDown alphas & betas", "Title",
 CellChangeTimes->{{3.7171355710697517`*^9, 3.7171356028895473`*^9}, {
  3.8029916730196133`*^9, 
  3.8029916809310827`*^9}},ExpressionUUID->"b53761d4-1b24-4b86-b93e-\
9c44d4b29d65"],

Cell[BoxData[
 RowBox[{"Quit", "[", "]"}]], "Input",
 CellChangeTimes->{{3.772628159141247*^9, 3.772628161249263*^9}, {
  3.785052589392028*^9, 3.785052592823535*^9}},
 CellLabel->
  "In[258]:=",ExpressionUUID->"980e79f6-9c07-4072-bd1e-300ac946e8c3"],

Cell[CellGroupData[{

Cell["Setup", "Section",
 CellChangeTimes->{{3.717145376675954*^9, 3.7171453778419027`*^9}, {
  3.725360092733674*^9, 
  3.725360093603505*^9}},ExpressionUUID->"fe9e23be-b684-46cc-9b54-\
957e00ada2e3"],

Cell[CellGroupData[{

Cell["Load package and set root paths", "Subsection",
 CellChangeTimes->{{3.799119833730782*^9, 3.7991198361011353`*^9}, {
  3.801906007839068*^9, 
  3.8019060131649723`*^9}},ExpressionUUID->"e952a5b9-bbfe-43a9-9668-\
8567ca76bf42"],

Cell[BoxData[
 RowBox[{
  RowBox[{"rootpath", "=", "\"\<~/git/rdstackingproject/\>\""}], 
  ";"}]], "Code",
 CellChangeTimes->{{3.801893836989571*^9, 3.8018938554497147`*^9}, {
  3.8018991673772917`*^9, 3.801899171272743*^9}, {3.801899220625408*^9, 
  3.801899220979004*^9}, {3.801899255251523*^9, 3.8018992567883883`*^9}, {
  3.801900861015358*^9, 3.801900867609252*^9}, {3.8019010062879868`*^9, 
  3.801901006499617*^9}, {3.801905973467177*^9, 3.8019060065780773`*^9}},
 CellLabel->
  "In[203]:=",ExpressionUUID->"ded0ca57-00ba-4008-906d-ee5041d06f91"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"rootpath", "<>", "\"\<code\>\""}], "]"}], ";"}], "\n", 
 RowBox[{"<<", "RDown.m"}], "\n", 
 RowBox[{"<<", "mcmc.m"}]}], "Code",
 CellChangeTimes->{{3.799119881021441*^9, 3.799119911381034*^9}, {