Is there a way to have the ar function work with tsdiag for on-the-fly visualization of ar fits? I have to fit a great many models of varying order and would like to save the diagnostic graphs. For instance, tsdiag(ar(lh)) tsdiag(arima(lh, order = c(3,0,0))) Thanks...