This function is similar to ts.diag in stats package, but with
one more diagnostic plot for the normality of residuals. Also, the default sequence of lags
for a Ljung-Box test is set to be seq(4,24,by = 4) if sample size \(n > 24\),
otherwise seq(1,n,4). This function has been automatically implemented in
estimate function.
Diagnostics are plotted, including the ACF plot, PACF plot, p.value of
white noise checking plot, and Q-Q plot for residuals.