par(mfrow=c(2,1))
fit2 = arima(sspot,order=c(5,0,0))
plot(fit2$resid,ty="l")



qqnorm(fit2$resid)

abline(a=0,b=1)