(* 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[ 34296, 866] NotebookOptionsPosition[ 33001, 822] NotebookOutlinePosition[ 33416, 840] CellTagsIndexPosition[ 33373, 837] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{"(*", RowBox[{ "This", " ", "code", " ", "is", " ", "designed", " ", "to", " ", "numerically", " ", "test", " ", "the", " ", "derivation", " ", "of", " ", "the", " ", "dynamical", " ", "theory", " ", "of", " ", "long", " ", "term", " ", "asexual", " ", "adaptation", " ", "that", " ", "incorporates", " ", "clonal", " ", "interference", " ", "and", " ", "diminishing", " ", "returns", " ", RowBox[{"epistasis", ".", " ", "This"}], " ", "confirms", " ", "the", " ", "accuracy", " ", "of", " ", "the", " ", "derivation", " ", "which", " ", "includes", " ", "approximations", " ", "for", " ", "small", " ", StyleBox["s", FontSlant->"Italic"], " ", "and", " ", "long", " ", RowBox[{"time", ".", " ", "See"}], " ", "simulations", " ", "in", " ", RowBox[{"Fig", ".", " ", "S6"}], " ", "for", " ", "confirmation", " ", "that", " ", "this", " ", "framework", " ", "is", " ", "an", " ", "accurate", " ", "representation", " ", "of", " ", "the", " ", "fitness", " ", RowBox[{"trajectory", "."}]}], "*)"}], "Subsection"]], "Input", CellChangeTimes->{{3.592676533293436*^9, 3.592676534354267*^9}, { 3.592676834141769*^9, 3.592676836389093*^9}, {3.5926777555439663`*^9, 3.592677770601284*^9}, {3.5926797941894617`*^9, 3.592679797844224*^9}, { 3.592680683003478*^9, 3.592680695112646*^9}, {3.5926809949345007`*^9, 3.5926810842925577`*^9}, {3.592681236540517*^9, 3.5926812505439053`*^9}, { 3.592681486467559*^9, 3.5926814933618107`*^9}}], Cell[BoxData[ StyleBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "The", " ", "best", " ", "fit", " ", "of", " ", "the", " ", "power", " ", "law", " ", "model", " ", "to", " ", "the", " ", "six", " ", "populations", " ", "that", " ", "kept", " ", "the", " ", "ancestral", " ", "mutation", " ", "rate", " ", "gives", " ", "a", " ", "value", " ", "of", " ", "the", " ", "diminishing", " ", "returns", " ", "parameter", " ", StyleBox["g", FontSlant->"Italic"], " ", "of", " ", "6.0217", " ", "and", " ", "is", " ", "consistent", " ", "with", " ", "pairs", " ", "of", " ", "values", " ", RowBox[{"(", RowBox[{"\[Alpha]0", " ", "and", " ", "\[Mu]"}], ")"}], " ", "as", " ", "summarized", " ", "in", " ", RowBox[{"Fig", ".", " ", "3"}], RowBox[{"A", ".", " ", "These"}], " ", "pairs", " ", RowBox[{"include", ":"}]}], "*)"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.5926758524728117`*^9, 3.5926758895873137`*^9}, { 3.592676006178081*^9, 3.5926761262776537`*^9}, {3.5926761647400513`*^9, 3.592676209962331*^9}, 3.5926762417011423`*^9, {3.5926765228673077`*^9, 3.592676647003566*^9}, {3.59267676974697*^9, 3.592676794023363*^9}, { 3.5926768401106167`*^9, 3.592676977387388*^9}, {3.592680704880715*^9, 3.592680708335227*^9}, 3.592681281101273*^9, {3.592681508827318*^9, 3.592681511284816*^9}, 3.592682125360487*^9, 3.592682218933366*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "37.3392"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "8"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.592677074505121*^9, 3.592677083415193*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "58.015"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "7"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5926770910664454`*^9, 3.5926770936489887`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "79.643"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "6"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.592677099770459*^9, 3.592677103737167*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Alpha]0", "=", "84.5697"}], ";", " ", RowBox[{"(*", RowBox[{ "this", " ", "pair", " ", "is", " ", "the", " ", "best", " ", "fit", " ", "to", " ", "the", " ", "hypermutator", " ", "trajectories"}], "*)"}], "\[IndentingNewLine]", RowBox[{"\[Mu]", "=", RowBox[{"1.7", "*", RowBox[{"10", "^", RowBox[{"-", "6"}]}]}]}], ";"}]], "Input", CellChangeTimes->{{3.592677161697652*^9, 3.59267718742659*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "101.638"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "5"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.592677110124364*^9, 3.59267711340331*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "123.808"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "4"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5926771189090357`*^9, 3.592677124819723*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "146.074"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"10", "^", RowBox[{"-", "3"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.592679828209852*^9, 3.592679841691132*^9}}], Cell[BoxData[ StyleBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "The", " ", "following", " ", "functions", " ", "represent", " ", "the", " ", "theory", " ", "of", " ", "clonal", " ", "interfence", " ", RowBox[{"(", RowBox[{"Gerrish", " ", "and", " ", "Lenski", " ", "1998"}], ")"}], " ", "with", " ", "no", " ", RowBox[{"approximations", ":"}]}], "*)"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.573391593542989*^9, 3.573391597219928*^9}, { 3.573391696567686*^9, 3.5733916973675833`*^9}, {3.5734006236036787`*^9, 3.573400623902433*^9}, {3.573400681909287*^9, 3.573400691795145*^9}, { 3.573405443864771*^9, 3.573405444544845*^9}, {3.574150234149631*^9, 3.574150243409227*^9}, 3.574180262852131*^9, {3.574470123221428*^9, 3.5744701235032597`*^9}, {3.592669092659823*^9, 3.5926690945011044`*^9}, 3.592671155145281*^9, 3.592671395767634*^9, {3.592672128247942*^9, 3.59267212838291*^9}, 3.592672688716034*^9, {3.5926759212904377`*^9, 3.5926759214727592`*^9}, {3.5926772065336857`*^9, 3.5926772228209133`*^9}, {3.5926772767798157`*^9, 3.592677316724634*^9}, { 3.5926774574258432`*^9, 3.592677479470162*^9}, {3.592677852088139*^9, 3.592677889962743*^9}, {3.59268073978549*^9, 3.592680742681608*^9}, 3.592681277762719*^9, 3.592682123059395*^9, 3.592682213532769*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"s_", ",", "\[Alpha]_"}], "]"}], ":=", RowBox[{"\[Alpha]", " ", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "\[Alpha]"}], " ", "s"}], "]"}]}]}], RowBox[{"(*", RowBox[{ "probability", " ", "density", " ", "of", " ", "available", " ", "beneficial", " ", "mutations", " ", RowBox[{"(", RowBox[{"mean", " ", RowBox[{"1", "/", "\[Alpha]"}]}], ")"}]}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Nu]", "[", "s_", "]"}], ":=", FractionBox[ RowBox[{"4", "s"}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "s"}], ")"}], "2"]]}], RowBox[{"(*", RowBox[{ "probability", " ", "that", " ", "a", " ", "mutation", " ", "with", " ", "effect", " ", "s", " ", "survives", " ", "drift"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Lambda]", "[", RowBox[{"s_", ",", "\[Alpha]_"}], "]"}], ":=", RowBox[{ FractionBox["\[Mu]", RowBox[{"Log", "[", RowBox[{"1", "+", "s"}], "]"}]], "Ne", " ", RowBox[{"Log", "[", "Ne", "]"}], " ", RowBox[{"Integrate", "[", RowBox[{ RowBox[{ RowBox[{"\[Nu]", "[", "sp", "]"}], " ", RowBox[{"f", "[", RowBox[{"sp", ",", "\[Alpha]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"sp", ",", "s", ",", "\[Infinity]"}], "}"}]}], "]"}]}]}], " ", RowBox[{"(*", RowBox[{ "expected", " ", "number", " ", "of", " ", "competing", " ", "mutations"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"p", "[", RowBox[{"s_", ",", "\[Alpha]_"}], "]"}], ":=", FractionBox[ RowBox[{ RowBox[{"\[Nu]", "[", "s", "]"}], " ", RowBox[{"Exp", "[", RowBox[{"-", RowBox[{"\[Lambda]", "[", RowBox[{"s", ",", "\[Alpha]"}], "]"}]}], "]"}], " ", RowBox[{"f", "[", RowBox[{"s", ",", "\[Alpha]"}], "]"}]}], RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"\[Nu]", "[", "sp", "]"}], " ", RowBox[{"Exp", "[", RowBox[{"-", RowBox[{"\[Lambda]", "[", RowBox[{"sp", ",", "\[Alpha]"}], "]"}]}], "]"}], " ", RowBox[{"f", "[", RowBox[{"sp", ",", "\[Alpha]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"sp", ",", "0", ",", "\[Infinity]"}], "}"}]}], "]"}]]}], RowBox[{"(*", RowBox[{ "probability", " ", "that", " ", "a", " ", "mutation", " ", "of", " ", "effect", " ", "s", " ", "is", " ", "fixed"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"sfix", "[", "\[Alpha]_", "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"s", " ", RowBox[{"p", "[", RowBox[{"s", ",", "\[Alpha]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", "0", ",", "\[Infinity]"}], "}"}], ",", RowBox[{"MaxRecursion", "\[Rule]", "20"}]}], "]"}]}], RowBox[{"(*", RowBox[{ "expected", " ", "value", " ", "of", " ", "fixed", " ", "beneficial", " ", "effect"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tfix", "[", "\[Alpha]_", "]"}], ":=", FractionBox["1", RowBox[{ RowBox[{"Log", "[", "2", "]"}], "\[Mu]", " ", "Ne", " ", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"\[Nu]", "[", "sp", "]"}], RowBox[{"Exp", "[", RowBox[{"-", RowBox[{"\[Lambda]", "[", RowBox[{"sp", ",", "\[Alpha]"}], "]"}]}], "]"}], " ", RowBox[{"f", "[", RowBox[{"sp", ",", "\[Alpha]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"sp", ",", "0", ",", "\[Infinity]"}], "}"}]}], "]"}]}]], RowBox[{"(*", RowBox[{"expected", " ", "time", " ", "to", " ", "fixation"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5733894409203377`*^9, 3.573389464871747*^9}, { 3.5733895746605*^9, 3.5733895790681*^9}, {3.573391750929729*^9, 3.573391762457711*^9}, {3.5733918653467913`*^9, 3.5733918672350492`*^9}, { 3.573396763665783*^9, 3.573396787966365*^9}, {3.57340073914016*^9, 3.573400748906588*^9}, {3.5734007807736*^9, 3.573400782057126*^9}, { 3.573401160115759*^9, 3.573401160713336*^9}, {3.573401212299294*^9, 3.573401212819322*^9}, {3.573401698202005*^9, 3.573401699489191*^9}, { 3.5734022222161617`*^9, 3.573402229736616*^9}, {3.573406474687737*^9, 3.5734065197397137`*^9}, 3.573406726459332*^9, {3.573406758665411*^9, 3.573406787023782*^9}, {3.573406880870039*^9, 3.573406993404174*^9}, { 3.573407402707111*^9, 3.573407428770537*^9}, {3.573407622364995*^9, 3.5734076366905823`*^9}, {3.592669832131611*^9, 3.592669833504074*^9}, { 3.592670027580001*^9, 3.592670032082766*^9}, {3.5926713623515453`*^9, 3.592671399257194*^9}, {3.592671501320867*^9, 3.5926715273712482`*^9}, { 3.5926715645511847`*^9, 3.5926715914334307`*^9}, {3.592672167133747*^9, 3.592672169412601*^9}, {3.592672693672345*^9, 3.592672700515389*^9}, { 3.59267275625279*^9, 3.592672759182848*^9}, {3.5926727926262197`*^9, 3.592672950668634*^9}, {3.592673061270187*^9, 3.592673095366528*^9}, { 3.592674138863965*^9, 3.592674143320156*^9}, {3.59267732025176*^9, 3.5926774518744593`*^9}, 3.592677488073292*^9, {3.592677520189649*^9, 3.592677563206328*^9}, {3.592677604152523*^9, 3.592677696828356*^9}, 3.592677902184828*^9, {3.592678109466927*^9, 3.592678110105795*^9}, { 3.5926783691589603`*^9, 3.5926784120921917`*^9}, 3.592678656662311*^9, { 3.592678897254118*^9, 3.592678939880129*^9}, {3.592679150922063*^9, 3.592679186623945*^9}, {3.592679245738606*^9, 3.592679249657916*^9}, { 3.59267931909529*^9, 3.592679319170245*^9}, {3.592679403338517*^9, 3.5926794037414494`*^9}, {3.59267971906982*^9, 3.592679751491281*^9}, { 3.592679870207447*^9, 3.592680009166417*^9}, {3.592681561865315*^9, 3.592681564958968*^9}}], Cell[BoxData[ StyleBox[ RowBox[{"(*", RowBox[{"Set", " ", "parameter", " ", RowBox[{"values", ":"}]}], "*)"}], "Subsection"]], "Input", CellChangeTimes->{{3.592680414050387*^9, 3.592680420126582*^9}, { 3.59268074797613*^9, 3.59268075031293*^9}, 3.592681283946159*^9, 3.5926821596393557`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[Alpha]0", "=", "84.5697"}], ";"}], RowBox[{"(*", RowBox[{ "initial", " ", "distribution", " ", "of", " ", "available", " ", "beneficial", " ", "effects"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"1.7", "*", RowBox[{"10", "^", RowBox[{"-", "6"}]}]}]}], ";"}], RowBox[{"(*", RowBox[{"beneficial", " ", "mutation", " ", "rate"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g", "=", "6.0217"}], ";"}], RowBox[{"(*", RowBox[{"diminishing", " ", "returns", " ", "parameter"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Ne", "=", RowBox[{"3.3", "*", RowBox[{"10", "^", "7"}]}]}], ";"}], RowBox[{"(*", RowBox[{"effective", " ", "population", " ", "size"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Nmutmax", "=", "24"}], ";"}], RowBox[{"(*", RowBox[{ "number", " ", "of", " ", "fixations", " ", "to", " ", "calculate"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.592680020172299*^9, 3.592680029045566*^9}, { 3.59268026836852*^9, 3.592680408326871*^9}, {3.592680442129601*^9, 3.5926804645561323`*^9}}], Cell[BoxData[ StyleBox[ RowBox[{"(*", RowBox[{ "The", " ", "following", " ", "calculates", " ", "the", " ", "series", " ", "of", " ", "\[Alpha]", " ", "values", " ", "and", " ", "fixed", " ", "beneficial", " ", "effects", " ", "and", " ", "fixation", " ", "times", " ", "using", " ", "the", " ", "above", " ", RowBox[{"framework", ":"}]}], "*)"}], "Subsection"]], "Input", CellChangeTimes->{{3.59268009189785*^9, 3.592680139839583*^9}, { 3.592680236220682*^9, 3.592680248938785*^9}, {3.592680422391943*^9, 3.5926804226080523`*^9}, {3.59268075444485*^9, 3.5926807571151543`*^9}, 3.592681285634671*^9, {3.592681586242264*^9, 3.5926815894248257`*^9}, 3.5926821658140373`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[Alpha]new", "[", RowBox[{"{", RowBox[{"\[Alpha]_", ",", "g_"}], "}"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"\[Alpha]", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"g", " ", RowBox[{"sfix", "[", "\[Alpha]", "]"}]}]}], ")"}]}], " ", ",", "g"}], "}"}]}], RowBox[{"(*", RowBox[{"diminishing", " ", "returns", " ", "parameterization"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Alpha]vals", "=", RowBox[{"First", "[", RowBox[{"Transpose", "[", RowBox[{"NestList", "[", RowBox[{"\[Alpha]new", ",", RowBox[{"{", RowBox[{"\[Alpha]0", ",", "g"}], "}"}], ",", RowBox[{"Nmutmax", "-", "1"}]}], "]"}], "]"}], "]"}]}], ";"}], RowBox[{"(*", RowBox[{"series", " ", "of", " ", "\[Alpha]", " ", "values"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"sfixvals", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"sfix", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "\[Alpha]vals"}], "}"}]}], "]"}]}], ";"}], RowBox[{"(*", RowBox[{ "series", " ", "of", " ", "fixed", " ", "beneficial", " ", "effects"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tfixintervals", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"tfix", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "\[Alpha]vals"}], "}"}]}], "]"}]}], ";", RowBox[{"(*", RowBox[{"time", " ", "intervals", " ", "between", " ", "fixations"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tfixvalscalc", "[", RowBox[{"{", RowBox[{"runningtotal_", ",", "index_"}], "}"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"runningtotal", "+", RowBox[{"tfixintervals", "[", RowBox[{"[", "index", "]"}], "]"}]}], ",", RowBox[{"index", "+", "1"}]}], "}"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tfixvals", "=", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"NestList", "[", RowBox[{"tfixvalscalc", ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"Nmutmax", "-", "1"}]}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], RowBox[{"(*", RowBox[{"timepoints", " ", "of", " ", "fixations"}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PlotSettings", "=", RowBox[{"{", RowBox[{ RowBox[{"PlotRange", "\[Rule]", "Full"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Automatic", ",", " ", "None"}], "}"}], ",", RowBox[{"{", RowBox[{"Automatic", ",", "None"}], "}"}]}], "}"}]}], ",", RowBox[{"LabelStyle", "\[Rule]", "13"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"1", "/", "\[Alpha]vals"}], ",", "PlotSettings", ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", " ", ",", " ", "\"\<1/\[Alpha]\>\""}], "}"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{"sfixvals", ",", "PlotSettings", ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ "\"\\"", " ", ",", " ", "\"\\""}], "}"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{"tfixintervals", ",", "PlotSettings", ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ "\"\\"", " ", ",", " ", "\"\