Bowling Green maxes, mins, and their standard deviations

Andy Long
Spring, 2020 -- year of the Covid-19 virus.

Let’s read in the data, which was created with BGdata.lsp:

(* 
iris=Import["/Users/longa/Desktop/NKU/classes/000 2020Spring/mat375/Fletcher/analysis/SummaryStatsBG.csv","CSV"];
*)
iris=Import["http://ceadserv1.nku.edu/longa//classes/mat375/data/Fletcher/SummaryStatsBG.csv","CSV"];
header=iris[[1]]
BG=iris[[2;;]];
{date,minmeans,minsds,maxmeans,maxsds}=Transpose[BG];
xs=date/366;

GraphicsGrid[
{{ListPlot[Transpose[{date,minmeans}],PlotLabel->"Mean Minima"],ListPlot[Transpose[{date,maxmeans}],PlotLabel->"Mean Maxima"]},{ListPlot[Transpose[{date,minsds}],PlotLabel->"StdDev Minima"],ListPlot[Transpose[{date,maxsds}],PlotLabel->"StdDev Maxima"]}}
]

data=Transpose[{xs,maxsds}];
dtrplot=ListPlot[data,PlotLabel->"Standard Deviation of Maximum Temperature"]
regMaxSD=LinearModelFit[data, {
Cos[2 Pi x], Sin[2Pi x],
Sin[4 Pi x],Cos[4 Pi x],
Sin[6 Pi x],Cos[6 Pi x],
Sin[8 Pi x],Cos[8 Pi x],
Sin[10 Pi x],Cos[10 Pi x],
Sin[12 Pi x],Cos[12 Pi x]
}
,x]
regMaxSDfits=Table[regMaxSD[x],{x,xs}];
regMaxSD["ParameterTable"]
regMaxSD["RSquared"]
regMaxSD["AdjustedRSquared"]
regMaxSD["ANOVATable"]
regMaxSD["ParameterConfidenceIntervals"]
maxsdsWmodel=Show[dtrplot,Plot[regMaxSD[x],{x,0,xs[[Length[xs]]]}]]

data=Transpose[{xs,minsds}];
dtrplot=ListPlot[data,PlotLabel->"Standard Deviation of Minimum Temperature"]
regminSD=LinearModelFit[data, {
Cos[2 Pi x], Sin[2Pi x],
Sin[4 Pi x],Cos[4 Pi x],
Sin[6 Pi x],Cos[6 Pi x],
Sin[8 Pi x],Cos[8 Pi x],
Sin[10 Pi x],Cos[10 Pi x],
Sin[12 Pi x],Cos[12 Pi x]
}
,x]
regminSDfits=Table[regminSD[x],{x,xs}];
regminSD["ParameterTable"]
regminSD["RSquared"]
regminSD["AdjustedRSquared"]
regminSD["ANOVATable"]
regminSD["ParameterConfidenceIntervals"]
minsdsWmodel=Show[dtrplot,Plot[regminSD[x],{x,0,xs[[Length[xs]]]}]]

data=Transpose[{xs,minmeans}];
dtrplot=ListPlot[data,PlotLabel->"Mean of Minimum Temperature"]
regminmeans=LinearModelFit[data, {
Cos[2 Pi x], Sin[2Pi x],
Sin[4 Pi x],Cos[4 Pi x],
Sin[6 Pi x],Cos[6 Pi x],
Sin[8 Pi x],Cos[8 Pi x],
Sin[10 Pi x],Cos[10 Pi x],
Sin[12 Pi x],Cos[12 Pi x]
}
,x]
regminmeansfits=Table[regminSD[x],{x,xs}];
regminmeans["ParameterTable"]
regminmeans["RSquared"]
regminmeans["AdjustedRSquared"]
regminmeans["ANOVATable"]
regminmeans["ParameterConfidenceIntervals"]
minmeansWmodel=Show[dtrplot,Plot[regminmeans[x],{x,0,xs[[Length[xs]]]}]]

data=Transpose[{xs,maxmeans}];
dtrplot=ListPlot[data,PlotLabel->"Mean of maximum Temperature"]
regmaxmeans=LinearModelFit[data, {
Cos[2 Pi x], Sin[2Pi x],
Sin[4 Pi x],Cos[4 Pi x],
Sin[6 Pi x],Cos[6 Pi x],
Sin[8 Pi x],Cos[8 Pi x],
Sin[10 Pi x],Cos[10 Pi x],
Sin[12 Pi x],Cos[12 Pi x]
}
,x]
regmaxmeans["ParameterTable"]
regmaxmeans["RSquared"]
regmaxmeans["AdjustedRSquared"]
regmaxmeans["ANOVATable"]
regmaxmeans["ParameterConfidenceIntervals"]
maxmeansWmodel=Show[dtrplot,Plot[regmaxmeans[x],{x,0,xs[[Length[xs]]]}]]

GraphicsGrid[{{minmeansWmodel,maxmeansWmodel},{minsdsWmodel,maxsdsWmodel}}]