tsdiag for ar?

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...