Download All

Code covered by the BSD License  

Highlights from
Fit GLM with quadratic penalty

image thumbnail
from Fit GLM with quadratic penalty by Patrick Mineault
Fits GLM with a quadratic penalty, determines hyperparams through cross-validation or evidence

qfsmooth1D(numx) 
function [qf] = qfsmooth1D(numx)
    %[qf] = qfsmooth1D(numx)
    %Create a quadratic form for smoothness regularization based on
    %second-order derivative operator, for a one-dimensional signal
    D = zeros(numx+1,numx);
    for ii = 1:numx-1
        D(ii,ii:ii+1) = [-1 1];
    end
    
    qf = D'*D;
end

Contact us at files@mathworks.com