I'm working on creating a model that examines the effect of ocean characteristics on fishing outcomes. I have spatial data on a 0.5 degree grid and I created the following model:
gam(inverse hyperbolic sine(yvar) ~ s(lat, lon, bs="sos") + s(xvar1) + s(xvar2) + s(xvar3), data = dat, method = "REML")
The QQ plot and histogram of residuals look okay. However, gam.check() produces an odd pattern in the residuals plot. I know that the points should be scattered around 0, but I have a very odd pattern in the residuals. Can anyone provide some insight on the interpretation of this pattern?
A similar question on Cross Validated had a solution suggesting the person asking the question try a Poisson model. I tried this and the fit of the model became much worse.
inverse hyperbolic sinein words inside code? $\endgroup$