I have a model comparison I'm doing using AIC, only the residuals of my models are not normally distributed. I know that in this case, the simple formula for AIC does not work.
I'm thinking I may need to use bootstrapping to calculate the likelihood of various models. However, I see several pitfalls in this approach, the biggest being that the bootstrapped likelihood will, in essence, be coming up with a different distribution of residuals for each model fitted, totally invalidating the comparison. How do I do this bootstrapping properly?