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

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

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     28879,        648]
NotebookOptionsPosition[     23133,        558]
NotebookOutlinePosition[     23563,        575]
CellTagsIndexPosition[     23520,        572]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[{
 StyleBox["Covid-19 Case Fitting Curves", "Title"],
 "\n",
 StyleBox["Esmerindo Bernardes", "Author"],
 "\n",
 StyleBox["sousa@ifsc.usp.br", "Affiliation"],
 "\n04/2020"
}], "Text",
 CellChangeTimes->{{3.796547237632004*^9, 3.7965472587694883`*^9}, 
   3.7965473940725307`*^9, {3.7965474294289565`*^9, 3.7965474671091633`*^9}, {
   3.7965475246711864`*^9, 3.7965475427308784`*^9}},
 TextAlignment->Center,
 TextJustification->1.,ExpressionUUID->"79c1882b-12d4-4336-94a2-cac77e8e48ae"],

Cell[TextData[StyleBox["Based on data from Covid-19 cases from China, a \
universal fitting curve is presented. Only five parameters are required to be \
determined by a given Covid-19 case data over a number of days. Under the \
special condition of reliable mass testing, this curve can indicate the \
flattening of a given case data after its inflection point is reached.", \
"Abstract"]], "Text",
 CellChangeTimes->{{3.796547237632004*^9, 3.7965472587694883`*^9}, 
   3.7965473940725307`*^9, {3.7965474294289565`*^9, 3.7965474671091633`*^9}, {
   3.7965475246711864`*^9, 3.796547545985223*^9}},
 TextJustification->1.,ExpressionUUID->"f4dd83a3-7b45-45b9-bd69-59e6ea175a3d"],

Cell["Model function: ", "Text",
 CellChangeTimes->{{3.7959822802307677`*^9, 3.7959822841782045`*^9}, {
   3.795982411522565*^9, 3.7959824302035985`*^9}, {3.795982527250079*^9, 
   3.7959825338354397`*^9}, 
   3.7959826124252586`*^9},ExpressionUUID->"500a9d89-2205-4e4e-80fe-\
c200b71e7f15"],

Cell[BoxData[
 RowBox[{"model", "=", 
  RowBox[{
   RowBox[{"b", "*", 
    RowBox[{"Tanh", "[", 
     RowBox[{
      RowBox[{"a", "*", "z"}], "-", "c"}], "]"}]}], "+", "d"}]}]], "Input",
 CellChangeTimes->{{3.7959823370727167`*^9, 3.7959823705870595`*^9}, {
  3.796547053509556*^9, 3.796547055302727*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"0e3220ee-ded7-4261-88a5-6254bb13ed07"],

Cell[CellGroupData[{

Cell["Generic Case Data:", "Text",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.796564078585355*^9, 3.7965641094856515`*^9}, {
  3.7965650354323792`*^9, 
  3.796565037592555*^9}},ExpressionUUID->"0009df83-a2c5-4963-92e4-\
bf892bc1e3bd"],

Cell[BoxData[
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "covid", "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.7965635533523207`*^9, 3.796563556549786*^9}, 
   3.7965641094856515`*^9},ExpressionUUID->"19882364-e372-47ba-afe8-\
04328fbfe1a3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"covid", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.796563558271174*^9, 3.796563578219798*^9}, 
   3.796564109486689*^9},ExpressionUUID->"a25ebe61-0f17-4682-ac71-\
18337cd507bc"],

Cell[BoxData[
 RowBox[{"lsf", "=", 
  RowBox[{"NonlinearModelFit", "[", 
   RowBox[{"data", ",", "model", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"a", ",", "0.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"c", ",", "4.0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"b", ",", 
        RowBox[{"10", "^"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"d", ",", 
        RowBox[{"10", "^"}]}], "}"}]}], "}"}], ",", 
    RowBox[{"{", "z", "}"}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.7965635823786397`*^9, 3.796563589584407*^9}, 
   3.796564109486689*^9},ExpressionUUID->"75516e81-37e7-44c4-bcc9-\
a53123c726cb"],

Cell[BoxData[
 RowBox[{"lsf", "[", "\"\<BestFit\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.796563592852621*^9, 3.796563593116957*^9}, 
   3.796564109486689*^9},ExpressionUUID->"f6a8a8ec-23ec-4953-9714-\
ec84f60fdaaa"],

Cell[BoxData[
 RowBox[{"lsf", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.79656359620271*^9, 3.796563596501899*^9}, 
   3.796564109486689*^9},ExpressionUUID->"9c92a5df-0056-4ef6-9d8f-\
fbeb9a77cb90"],

Cell[BoxData[
 RowBox[{"lsf", "[", "\"\<ANOVATable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.796563598223304*^9, 3.7965635984885473`*^9}, 
   3.796564109486689*^9},ExpressionUUID->"bf4fc05d-4ab8-4ed7-b623-\
b8ae5e8c5c80"],

Cell[BoxData[
 RowBox[{"lsf", "[", "\"\<AdjustedRSquared\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.7965638698766265`*^9, 3.7965638701100206`*^9}, 
   3.796564109487692*^9},ExpressionUUID->"5094aa3d-d56f-47ab-a5e8-\
25fa9b265eba"],

Cell[BoxData[
 RowBox[{"gr", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"lsf", "[", "i", "]"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "90"}], "}"}], ",", 
    RowBox[{"Epilog", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Medium", "]"}], ",", 
       RowBox[{"Point", "[", "data", "]"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", "Red", "}"}], "}"}]}], ",", 
    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
    RowBox[{"GridLines", "->", "Automatic"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<N (days)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(Z\), \(N\)]\)\>\""}],
       "}"}]}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.7965636001840415`*^9, 3.7965636197686625`*^9}, {
   3.7965640418595915`*^9, 3.796564042026146*^9}, 
   3.796564109487692*^9},ExpressionUUID->"fc26a5db-5b8b-48b8-82a0-\
d02e8525b3ac"],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"lsf", "[", "\"\<FitResiduals\>\"", "]"}], ",", 
   RowBox[{"Filling", "\[Rule]", "Axis"}]}], "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.796563613785625*^9, 3.796563614044936*^9}, 
   3.796564109487692*^9},ExpressionUUID->"67c63e37-377a-40fb-b0e1-\
80c5e5b5df98"]
}, Closed]],

Cell[CellGroupData[{

Cell["US Case Data:", "Text",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.7965477061952963`*^9, 3.7965477119748354`*^9}, 
   3.796547937474633*^9, {3.796564277909169*^9, 
   3.796564278440709*^9}},ExpressionUUID->"0a927a1e-9ff5-4cd4-884d-\
c8b250a073ce"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"covidus", "=", 
   RowBox[{"{", 
    RowBox[{
    "5", ",", "53", ",", "57", ",", "60", ",", "63", ",", "68", ",", "75", 
     ",", "100", ",", "124", ",", "158", ",", "221", ",", "319", ",", "435", 
     ",", "541", ",", "704", ",", "994", ",", "1301", ",", "1630", ",", 
     "2183", ",", "2770", ",", "3613", ",", "4596", ",", "6344", ",", "9197", 
     ",", "13779", ",", "19367", ",", "24192", ",", "33592", ",", "43781", 
     ",", "54856", ",", "68211", ",", "85435", ",", "104126", ",", "123578", 
     ",", "143491", ",", "163788", ",", "188530", ",", "215003", ",", 
     "244877", ",", "277161", ",", "311357", ",", "336673", ",", "367004", 
     ",", "400335", ",", "434927", ",", "468566", ",", "502876", ",", 
     "532879", ",", "560300", ",", "586941", ",", "617661", ",", "648003", 
     ",", "677570", ",", "709735", ",", "738792", ",", "764636", ",", 
     "792759", ",", "818744"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"nus", "=", 
  RowBox[{"Length", "[", "covidus", "]"}]}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.796547740668082*^9, 3.7965477819457135`*^9}, {
   3.796547812449071*^9, 3.7965478230746775`*^9}, {3.7965478686248426`*^9, 
   3.7965478700440083`*^9}, 3.796547937474633*^9},
 CellLabel->"In[1]:=",ExpressionUUID->"117929e0-92c5-47fc-8fb2-4abf5f6df02e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"dataus", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"n", ",", 
       RowBox[{"covidus", "[", 
        RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "nus"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.7965478480099864`*^9, 3.7965478616863637`*^9}, {
  3.7965479060856004`*^9, 3.796547955019704*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"ce842777-c97d-47f9-b426-f246e49f0209"],

Cell["Nonlinear fitting:", "Text",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.7959829937613907`*^9, 3.7959830116724834`*^9}, 
   3.79656418037705*^9},ExpressionUUID->"cd235f9e-3847-4a67-97a0-\
4ab15bb18d54"],

Cell[BoxData[
 RowBox[{"lsfus", "=", 
  RowBox[{"NonlinearModelFit", "[", 
   RowBox[{"dataus", ",", "model", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"a", ",", "0.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"c", ",", "4.0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"b", ",", 
        RowBox[{"10", "^", "5"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"d", ",", 
        RowBox[{"10", "^", "5"}]}], "}"}]}], "}"}], ",", 
    RowBox[{"{", "z", "}"}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.795983903500692*^9, 3.795983906745016*^9}, {
   3.796548040044265*^9, 3.7965480426702886`*^9}, {3.796548080338478*^9, 
   3.7965480970068955`*^9}, {3.7965481666496396`*^9, 
   3.7965481696176643`*^9}, {3.7965482779548655`*^9, 
   3.7965483227231116`*^9}, {3.796548683130039*^9, 3.796548683199852*^9}, {
   3.7965615835628633`*^9, 3.79656164461755*^9}, 3.79656418037705*^9},
 CellLabel->"In[30]:=",ExpressionUUID->"539c17f3-4103-4069-ac52-c1b059f74754"],

Cell[BoxData[
 RowBox[{"lsfus", "[", "\"\<BestFit\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.796548352916387*^9, 3.79654836165697*^9}, {
   3.7965485105706277`*^9, 3.7965485173654957`*^9}, {3.796548690339754*^9, 
   3.7965486904295163`*^9}, 3.7965641938130774`*^9},
 CellLabel->"In[31]:=",ExpressionUUID->"1b52d6f7-db8a-4d9a-b955-428a99adcf73"],

Cell[BoxData[
 RowBox[{"lsfus", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.796548141238573*^9, 3.7965481420883455`*^9}, {
   3.796548675358828*^9, 3.796548698249633*^9}, 3.7965641938130774`*^9},
 CellLabel->"In[19]:=",ExpressionUUID->"f6292780-e6a1-43e4-b540-e29a86963b0b"],

Cell[BoxData[
 RowBox[{"lsfus", "[", "\"\<ANOVATable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.796548746977254*^9, 3.7965487470640182`*^9}, 
   3.79656419381411*^9},
 CellLabel->"In[20]:=",ExpressionUUID->"c3e40960-9850-4995-851f-0e97d8017c33"],

Cell[BoxData[
 RowBox[{"lsfus", "[", "\"\<AdjustedRSquared\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.7965626159692554`*^9, 3.7965626456628065`*^9}, {
   3.796562898826573*^9, 3.796562917855715*^9}, 3.796564207749848*^9},
 CellLabel->"In[34]:=",ExpressionUUID->"d84eace6-5bd0-4a1b-a37c-5cda5b381595"],

Cell[BoxData[
 RowBox[{"grus", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"lsfus", "[", "n", "]"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "90"}], "}"}], ",", 
    RowBox[{"Epilog", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Medium", "]"}], ",", 
       RowBox[{"Point", "[", "dataus", "]"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", "Red", "}"}], "}"}]}], ",", 
    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
    RowBox[{"GridLines", "->", "Automatic"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<N (days)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(Z\), \(N\)]\)\>\""}],
       "}"}]}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.7965487933721905`*^9, 3.79654882306572*^9}, 
   3.796548872637116*^9, {3.796548979761569*^9, 3.796549049098093*^9}, 
   3.796563681761838*^9, 3.7965642077508373`*^9},
 CellLabel->"In[40]:=",ExpressionUUID->"a0c0f0cb-4840-4e57-8e4b-a7405d36b0a6"],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"lsfus", "[", "\"\<FitResiduals\>\"", "]"}], ",", 
   RowBox[{"Filling", "\[Rule]", "Axis"}]}], "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.796549082453868*^9, 3.7965490840276556`*^9}, 
   3.7965642077508373`*^9},
 CellLabel->"In[25]:=",ExpressionUUID->"83411e6c-41ac-4e30-a174-560c4d70fca5"]
}, Closed]],

Cell[CellGroupData[{

Cell["Iran Case Data:", "Text",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.79656305413314*^9, 3.796563061185311*^9}, 
   3.7965642416959944`*^9, {3.796564273455077*^9, 
   3.7965642825038834`*^9}},ExpressionUUID->"a489925e-ff4d-4688-9d5d-\
69a6e4a834aa"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"covidir", "=", 
   RowBox[{"{", 
    RowBox[{
    "2", ",", "5", ",", "18", ",", "29", ",", "43", ",", "61", ",", "95", ",",
      "139", ",", "245", ",", "388", ",", "593", ",", "978", ",", "1501", ",",
      "2336", ",", "2922", ",", "3513", ",", "4747", ",", "5823", ",", "6566",
      ",", "7161", ",", "8042", ",", "9000", ",", "10075", ",", "11364", ",", 
     "12729", ",", "13938", ",", "14991", ",", "16169", ",", "17361", ",", 
     "18407", ",", "19644", ",", "20610", ",", "21638", ",", "23049", ",", 
     "24811", ",", "27017", ",", "29406", ",", "32332", ",", "35408", ",", 
     "38309", ",", "41495", ",", "44605", ",", "47593", ",", "50468", ",", 
     "53183", ",", "55743", ",", "58226", ",", "60500", ",", "62589", ",", 
     "64586", ",", "66220", ",", "68192", ",", "70029", ",", "71686", ",", 
     "73303", ",", "74877", ",", "76389", ",", "77995", ",", "79494", ",", 
     "80868", ",", "82211", ",", "83505", ",", "84802"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"nir", "=", 
  RowBox[{"Length", "[", "covidir", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"datair", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"n", ",", 
       RowBox[{"covidir", "[", 
        RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "nir"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.796563070957132*^9, 3.796563086759915*^9}, {
   3.7965631238676414`*^9, 3.7965631838342247`*^9}, 3.7965642416959944`*^9},
 CellLabel->"In[35]:=",ExpressionUUID->"efa0cc09-ec0e-4a38-b5a2-4b4fc5bc5b62"],

Cell[BoxData[
 RowBox[{"lsfir", "=", 
  RowBox[{"NonlinearModelFit", "[", 
   RowBox[{"datair", ",", "model", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"a", ",", "0.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"c", ",", "4.0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"b", ",", 
        RowBox[{"10", "^", "4"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"d", ",", 
        RowBox[{"10", "^", "4"}]}], "}"}]}], "}"}], ",", 
    RowBox[{"{", "z", "}"}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.7965632302426224`*^9, 3.7965632471364403`*^9}, 
   3.79656424169699*^9},
 CellLabel->"In[38]:=",ExpressionUUID->"af46e476-a026-45bc-b664-eee02cfceb32"],

Cell[BoxData[
 RowBox[{"lsfir", "[", "\"\<BestFit\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.7965632961363597`*^9, 3.79656329704694*^9}, 
   3.79656424169699*^9},
 CellLabel->"In[39]:=",ExpressionUUID->"dc602c4a-64a0-418f-bcd2-ed1112aa7955"],

Cell[BoxData[
 RowBox[{"lsfir", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.79656359620271*^9, 3.796563596501899*^9}, {
   3.796563778115055*^9, 3.7965637788920155`*^9}, 3.7965642600608683`*^9},
 CellLabel->"In[41]:=",ExpressionUUID->"4237c66b-92a6-460e-abe3-42b27f53dce1"],

Cell[BoxData[
 RowBox[{"lsfir", "[", "\"\<ANOVATable\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.796563598223304*^9, 3.7965635984885473`*^9}, {
   3.7965637814172707`*^9, 3.796563782262002*^9}, 3.7965642600608683`*^9},
 CellLabel->"In[42]:=",ExpressionUUID->"42c87934-9e81-4502-9252-af5cefcedafc"],

Cell[BoxData[
 RowBox[{"lsfir", "[", "\"\<AdjustedRSquared\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.796563882457983*^9, 3.796563883384471*^9}, 
   3.796564260061897*^9},
 CellLabel->"In[43]:=",ExpressionUUID->"f0361e05-0b25-476b-acd0-a566ce50e706"],

Cell[BoxData[
 RowBox[{"grir", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"lsfir", "[", "i", "]"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "90"}], "}"}], ",", 
    RowBox[{"Epilog", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Medium", "]"}], ",", 
       RowBox[{"Point", "[", "datair", "]"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", "Red", "}"}], "}"}]}], ",", 
    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
    RowBox[{"GridLines", "->", "Automatic"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<N (days)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(Z\), \(N\)]\)\>\""}],
       "}"}]}]}], "]"}]}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.7965636001840415`*^9, 3.7965636197686625`*^9}, {
   3.796563785227079*^9, 3.796563789715025*^9}, {3.7965639080624814`*^9, 
   3.796563909246317*^9}, 3.796564260061897*^9},
 CellLabel->"In[45]:=",ExpressionUUID->"c2e42463-18f6-45ac-ac2f-7ce63896226f"],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"lsfir", "[", "\"\<FitResiduals\>\"", "]"}], ",", 
   RowBox[{"Filling", "\[Rule]", "Axis"}]}], "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.796563613785625*^9, 3.796563614044936*^9}, {
   3.796563797948036*^9, 3.796563798935358*^9}, 3.796564260061897*^9, 
   3.7965643786399665`*^9, {3.7965644524924173`*^9, 3.796564452855442*^9}},
 CellLabel->"In[48]:=",ExpressionUUID->"62cb5cd5-8f0e-441e-a23d-6a88c5388116"],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"lsfir", "[", "\"\<StandardizedResiduals\>\"", "]"}], ",", 
   RowBox[{"Filling", "\[Rule]", "Axis"}]}], "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{3.7965644255973134`*^9, 3.796564733577778*^9},
 CellLabel->"In[47]:=",ExpressionUUID->"8a918b27-f84c-49a7-bd27-884bf38d3de4"]
}, Closed]],

Cell["", "Text",ExpressionUUID->"510f8ec8-8475-43e2-81ec-9a249d9e834f"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"{", 
   RowBox[{"grus", ",", "grir"}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.7959704254067717`*^9, 3.795970476877093*^9}, 
   3.795972289291597*^9, {3.7959724105971127`*^9, 3.7959724470894966`*^9}, {
   3.7959819848821697`*^9, 3.795982015345639*^9}, {3.7959820698708305`*^9, 
   3.79598211446255*^9}, {3.7959837623057365`*^9, 3.7959837649726343`*^9}, {
   3.7965648551984515`*^9, 3.796564886806932*^9}},
 CellLabel->"In[49]:=",ExpressionUUID->"686e02b2-ecdc-45ca-b245-e9640233da53"],

Cell[BoxData[
 RowBox[{"grw", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"lsf", "[", "n", "]"}], ",", 
      RowBox[{"lsf2", "[", 
       RowBox[{"n", "+", "22"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "90"}], "}"}], ",", 
    RowBox[{"Epilog", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Medium", "]"}], ",", 
       RowBox[{"Point", "[", "data", "]"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", "Red", "}"}], ",", "Blue"}], "}"}]}], ",", 
    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
    RowBox[{"GridLines", "->", "Automatic"}], ",", 
    RowBox[{"PlotLegends", "\[Rule]", 
     RowBox[{"Placed", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
        "\"\<\!\(\*SubscriptBox[\(Z\), \(83\)]\)\>\"", ",", 
         "\"\<\!\(\*SubscriptBox[\(Z\), \(61\)]\)\>\""}], "}"}], ",", 
       "Below"}], "]"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<N (days)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(Z\), \(N\)]\)\>\""}],
       "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7959701651128826`*^9, 3.795970231738864*^9}, 
   3.79597028129426*^9, {3.795970512062011*^9, 3.7959705141812973`*^9}, {
   3.7959725626394095`*^9, 3.795972568614464*^9}, {3.7959732460683804`*^9, 
   3.795973253297045*^9}, {3.7959732837834496`*^9, 3.795973310922859*^9}, {
   3.7959733691869984`*^9, 3.7959734162201915`*^9}, {3.7959734711342955`*^9, 
   3.795973486957021*^9}, {3.7959735430618925`*^9, 3.795973616650051*^9}, {
   3.795973760203054*^9, 3.79597380585697*^9}, {3.7959814809159117`*^9, 
   3.7959814837922173`*^9}, {3.7959816483720737`*^9, 3.79598164848573*^9}, {
   3.795983112999463*^9, 3.7959831142710757`*^9}, {3.795983805468328*^9, 
   3.7959838588191743`*^9}, {3.7959846258331966`*^9, 
   3.7959847253300486`*^9}, {3.79598491219718*^9, 3.7959849205867333`*^9}, {
   3.795985046275531*^9, 3.795985111005371*^9}, {3.7959851537779346`*^9, 
   3.795985202862626*^9}, {3.7959858321850224`*^9, 3.795985848065854*^9}, 
   3.795985894759939*^9, {3.7959859795212536`*^9, 3.795985983346019*^9}, {
   3.7959861302380896`*^9, 3.7959861400777254`*^9}, {3.7959861787283697`*^9, 
   3.795986196874827*^9}, {3.795986232802726*^9, 3.795986232806716*^9}, {
   3.7959863421033354`*^9, 3.795986448019844*^9}, {3.795986510460026*^9, 
   3.7959865146398516`*^9}, {3.795986886312641*^9, 3.7959868929109993`*^9}, {
   3.79598707606621*^9, 3.7959871787955227`*^9}, {3.795987707270095*^9, 
   3.7959877124462605`*^9}, {3.795987934959015*^9, 3.795987962630041*^9}},
 CellLabel->
  "In[127]:=",ExpressionUUID->"c4dca211-4f86-4e0f-bed2-1d9f2b140cba"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.795986898990738*^9, 3.79598691277786*^9}, {
   3.795987490162591*^9, 3.7959874995166035`*^9}, 
   3.796564870996189*^9},ExpressionUUID->"62e6b5e2-5b87-48f4-bd86-\
105f4c15ff34"]
},
WindowSize->{1217, 711},
WindowMargins->{{Automatic, 66}, {Automatic, 0}},
Magnification:>1.25 Inherited,
FrontEndVersion->"12.1 for Microsoft Windows (64-bit) (March 18, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"60c325b8-744c-4a6c-a8c6-2307c6c037e9"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 500, 12, 224, "Text",ExpressionUUID->"79c1882b-12d4-4336-94a2-cac77e8e48ae"],
Cell[1061, 34, 677, 9, 121, "Text",ExpressionUUID->"f4dd83a3-7b45-45b9-bd69-59e6ea175a3d"],
Cell[1741, 45, 291, 5, 53, "Text",ExpressionUUID->"500a9d89-2205-4e4e-80fe-c200b71e7f15"],
Cell[2035, 52, 383, 9, 44, "Input",ExpressionUUID->"0e3220ee-ded7-4261-88a5-6254bb13ed07"],
Cell[CellGroupData[{
Cell[2443, 65, 268, 5, 53, "Text",ExpressionUUID->"0009df83-a2c5-4963-92e4-bf892bc1e3bd",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[2714, 72, 294, 6, 44, "Input",ExpressionUUID->"19882364-e372-47ba-afe8-04328fbfe1a3",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[3011, 80, 514, 14, 44, "Input",ExpressionUUID->"a25ebe61-0f17-4682-ac71-18337cd507bc",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[3528, 96, 722, 20, 44, "Input",ExpressionUUID->"75516e81-37e7-44c4-bcc9-a53123c726cb",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[4253, 118, 274, 5, 44, "Input",ExpressionUUID->"f6a8a8ec-23ec-4953-9714-ec84f60fdaaa",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[4530, 125, 280, 5, 44, "Input",ExpressionUUID->"9c92a5df-0056-4ef6-9d8f-fbeb9a77cb90",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[4813, 132, 279, 5, 44, "Input",ExpressionUUID->"bf4fc05d-4ab8-4ed7-b623-b8ae5e8c5c80",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[5095, 139, 287, 5, 44, "Input",ExpressionUUID->"5094aa3d-d56f-47ab-a5e8-25fa9b265eba",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[5385, 146, 1047, 27, 74, "Input",ExpressionUUID->"fc26a5db-5b8b-48b8-82a0-d02e8525b3ac",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[6435, 175, 377, 8, 44, "Input",ExpressionUUID->"67c63e37-377a-40fb-b0e1-80c5e5b5df98",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}]
}, Closed]],
Cell[CellGroupData[{
Cell[6849, 188, 287, 5, 44, "Text",ExpressionUUID->"0a927a1e-9ff5-4cd4-884d-c8b250a073ce",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[7139, 195, 1385, 23, 190, "Input",ExpressionUUID->"117929e0-92c5-47fc-8fb2-4abf5f6df02e",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[8527, 220, 570, 14, 44, "Input",ExpressionUUID->"ce842777-c97d-47f9-b426-f246e49f0209",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[9100, 236, 241, 4, 53, "Text",ExpressionUUID->"cd235f9e-3847-4a67-97a0-4ab15bb18d54",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[9344, 242, 1060, 24, 44, "Input",ExpressionUUID->"539c17f3-4103-4069-ac52-c1b059f74754",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[10407, 268, 401, 6, 44, "Input",ExpressionUUID->"1b52d6f7-db8a-4d9a-b955-428a99adcf73",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[10811, 276, 355, 5, 44, "Input",ExpressionUUID->"f6292780-e6a1-43e4-b540-e29a86963b0b",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[11169, 283, 302, 5, 44, "Input",ExpressionUUID->"c3e40960-9850-4995-851f-0e97d8017c33",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[11474, 290, 357, 5, 44, "Input",ExpressionUUID->"d84eace6-5bd0-4a1b-a37c-5cda5b381595",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[11834, 297, 1116, 27, 74, "Input",ExpressionUUID->"a0c0f0cb-4840-4e57-8e4b-a7405d36b0a6",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[12953, 326, 405, 8, 44, "Input",ExpressionUUID->"83411e6c-41ac-4e30-a174-560c4d70fca5",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[13395, 339, 288, 5, 44, "Text",ExpressionUUID->"a489925e-ff4d-4688-9d5d-69a6e4a834aa",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[13686, 346, 1703, 33, 219, "Input",ExpressionUUID->"efa0cc09-ec0e-4a38-b5a2-4b4fc5bc5b62",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[15392, 381, 759, 20, 44, "Input",ExpressionUUID->"af46e476-a026-45bc-b664-eee02cfceb32",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[16154, 403, 298, 5, 44, "Input",ExpressionUUID->"dc602c4a-64a0-418f-bcd2-ed1112aa7955",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[16455, 410, 354, 5, 44, "Input",ExpressionUUID->"4237c66b-92a6-460e-abe3-42b27f53dce1",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[16812, 417, 353, 5, 44, "Input",ExpressionUUID->"42c87934-9e81-4502-9252-af5cefcedafc",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[17168, 424, 307, 5, 44, "Input",ExpressionUUID->"f0361e05-0b25-476b-acd0-a566ce50e706",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[17478, 431, 1121, 27, 74, "Input",ExpressionUUID->"c2e42463-18f6-45ac-ac2f-7ce63896226f",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[18602, 460, 523, 9, 44, "Input",ExpressionUUID->"62cb5cd5-8f0e-441e-a23d-6a88c5388116",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[19128, 471, 384, 7, 44, "Input",ExpressionUUID->"8a918b27-f84c-49a7-bd27-884bf38d3de4",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}]
}, Closed]],
Cell[19527, 481, 71, 0, 44, "Text",ExpressionUUID->"510f8ec8-8475-43e2-81ec-9a249d9e834f"],
Cell[19601, 483, 538, 9, 44, "Input",ExpressionUUID->"686e02b2-ecdc-45ca-b245-e9640233da53"],
Cell[20142, 494, 2758, 56, 103, "Input",ExpressionUUID->"c4dca211-4f86-4e0f-bed2-1d9f2b140cba"],
Cell[22903, 552, 226, 4, 44, "Input",ExpressionUUID->"62e6b5e2-5b87-48f4-bd86-105f4c15ff34"]
}
]
*)