来自ARIMA（）函数的诊断图

Question

我的时间序列为returns。我发现了ARIMA（p，d，q）的最佳c（p，d，q）（我称之为return.order），所以我终于得到了：

returns.arima <- prima(returns, order=returns.order)

这不是lm文件，因此在计算plot(returns.arima)后我无法获得以下诊断图：

但我得到的只是反单元根：

所以我的问题如下：

1. 有没有办法可以将我的格式转换为lm格式，以便我可以轻松生成此输出？
2. 如果没有，请说我设法通过以下方式获取QQ-Plot的残差：

qqnorm(residuals.arima)

然而，我没有像我发布的第一张图片那样找到报告异常值的方法。 有人知道怎么做吗？
1. 剩余与杠杆，比例位置和库克的距离图怎么样？