Length[Q]
1164
rpita=TimeSeries[Q,{{1926,1,1},{2022,12,1},"Month"}];DateListPlot[rpita,PlotStyleBlack]
linear=Predict[Thread[tQ],Method"LinearRegression"]rf=Predict[Thread[tQ],Method"RandomForest"]nn=Predict[Thread[tQ],Method"NearestNeighbors"]sarima=TimeSeriesModelFit[Data,{"SARIMA",{5,0,4,2,0,1,12}}]
PredictorFunction
PredictorFunction
PredictorFunction
TimeSeriesModel
erL=Table[Q[[i]]-linear[i],{i,949,1164}];rmseL=Sqrt[Mean[erL^2]]maeL=Mean[Abs[erL]]
902.733
724.992
SST=Total[(Q2-Mean[Q2])^2]Qpredl=Table[linear[i],{i,949,1164}];SSE=Total[(Q2-Qpredl)^2]R2l=1-(SSE/SST)
1.76528×
8
10
1.76024×
8
10
0.00285025
erN=Table[Q2[[i]]-nn[i],{i,949,1164}];rmseN=Sqrt[Mean[erN^2]]maeN=Mean[Abs[erN]]
399.033
305.846
SST=Total[(Q2-Mean[Q2])^2]QpredN=Table[nn[i],{i,949,1164}];SSE=Total[(Q2-QpredN)^2]R2N=1-(SSE/SST)
1.76528×
8
10
3.43931×
7
10
0.805169
erR=Table[Q2[[i]]-rf[i],{i,949,1164}];rmseR=Sqrt[Mean[erR^2]]maeR=Mean[Abs[erR]]
438.637
322.643
SST=Total[(Q2-Mean[Q2])^2]QpredR=Table[rf[i],{i,949,1164}];SSE=Total[(Q2-QpredR)^2]R2R=1-(SSE/SST)
1.76528×
8
10
4.1559×
7
10
0.764575
erS=Table[Q[[i]]-sarima[i],{i,949,1164}];rmseS=Sqrt[Mean[erS^2]]maeS=Mean[Abs[erS]]
0.
0.
SST=Total[(Q2-Mean[Q2])^2]QpredR=Table[rf[i],{i,949,1164}];SSE=Total[(Q2-QpredR)^2]R2R=1-(SSE/SST)
Show[ListPlot[testingData],Plot[{sarima[x],linear[x],rf[x],nn[x]},{x,949,1164}]]
Plot[{rf[x],nn[x],linear[x]},{x,0,110}]