My data is about the number of earthquakes with magnitude greater than or equal to 7 that have occurred from 1912-2013. The Jarque-Bera test shows that the residuals of arima models are not equal. What should I do? I am using E-views to analyze my data. Below is the descriptive stat of ARI(4,1)
From my perspective you have already picked data that tends to be not normal as you are focussing on somewhat extreme cases. How do the residuals look if you include the number of earthquakes of all magnitudes? Additionally you could think about control variables that capture some of the structure in the residuals.
Another thing I would recommend you to do would be a closer analysis in which way the resdiuals deviate from normal. Are the skew? Do they have heavy tails? You can make a time-series plot of your residuals in order to see, whether there is a general structure in your data that your model does not cover. You can create a qqplot in order to look for heavy tails and you can plot a kernel density estimate of your residuals against a normal distribution with $\mu=0$ and $\sigma^2$ equal to the variance of your residuals in order to get some general impression how your residuals deviate.
I you can't get rid of the non-normality of the residuals you still can interpret the parameters but you should be very carefull with making inference and interpreting the p-values.