Strengths: A well-motivated project trying a range of techniques. The conclusions noted previous projects addressing similar questions with different data, and discussed what new things were learned from their analysis.
Points for consideration:
“data after first difference show more stationary than before” does not formally have meaning - a model is stationary or not, but we have no definition for “more stationary”.
The peaks in the estimated spectrum are reported wrongly. For example, 0.044 is reported as 4.044. In fact, the spectral analysis is consistent with the wavelets.
The conclusion “This leads us to conclude that there is no meaningful seasonality in our Detroit PM2.5 data” is incorrect. You find seasonal frequencies in the spectrum and wavelet analysis. ARIMA(2,1,3) can model seasonal fluctuations.
Typo: \(\mu=E[Y_n]\) not \(\mu=E[X_n]\).
Inconsistent values in the AIC table (not mathematically possible under accurate likelihood maximization and evaluation) should be noted and discussed. The discrepancy is large for the SARIMA table.
Looking at the SARIMA roots might help to understand numerical instability.
The degrees of freedom for the likelihood ratio test is incorrectly take to be 3 rather than 2.