I have a dataset where I suspect the residuals are approximately Laplace-distributed.
There are three continuous predictors. When I split up the data into many bins, based on the values of these predictors, and plot the density for each subset, it seems to me that the data is always centered at 0, but lower values of the predictors lead to wide laplace densities and higher values lead to narrower laplace densities.
My intuition is to use some kind of linear regression where the dispersion, rather than the mean, is modeled.
- Can a GLM handle the situation where the residuals are laplace-distributed?
- Is there a mature package in R that can do this? Seems like stats::glm cannot.
- I am familiar with the DGLM (double generalized linear model), so I suspect if I could overcome 1 and 2 I could do something to model the dispersion.
gamlss
can do that. You can specifically model the scale and location separately. The accompanying packagegamlss.dist
includes the Laplace distribution viaPE
. $\endgroup$