First of all, congrats to the R-INLA team for the amazing work.
I have a question regarding the specification of the parameter alpha in the inla.spde2.matern() function.
Alpha is said to be a fractional operator equal to the smoothness of the Matern kernel plus half the number of dimensions. Since the upper bound for alpha is 2, the largest smoothness that one can set in a 2-dimensional problem is apparently 1.
For a comparison study, I would really like to rely on a model based on a Matern 1.5 covariance function; does anybody know how to get the job done in this case?
Thank you very much for your help!
Filippo