Generalized linear models with Gaussian priors in Matlab

I have posted a set of functions to fit GLMs with Gaussian priors on Matlab Central. GLMs are a class of models in which an internal variable is contructed linearly from regressors, ie, \mathbf{r} = \mathbf{Xw}, and the observed variable is generated from the internal variable through a non-Gaussian noise process. GLMs include standard linear regression, Poisson regression, and logistic regression models as special cases. These three models are very useful in the estimation of neuronal receptive fields. The variable \mathbf{w} can be fitted by maximizing the likelihood of the data with respect to the model parameters; the built-in Matlab function glmfit does this . This approach is prone to overfitting, however.

Using a prior on \mathbf{w} and finding the maximum a posteriori estimate of \mathbf{w} alleviates overfitting. Gaussian priors of the form N(\mathbf{w}|0;(\lambda \mathbf{A})^{-1}) allow the expression of prior knowledge that the weights are of small magnitude (A = I) or that the weights are smooth (A = D, a second order finite difference operator), for example. The set of functions allow the determination of the hyperparameter \lambda through k-fold cross-validation.

3 thoughts on “Generalized linear models with Gaussian priors in Matlab

Leave a comment

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s