function logL = log_like_norm(gam,y,X)
%log likelihood for NLRM
%gam = [beta' sigma^2]', y [n x 1],  X [n x k]

k = length(y);
b = gam(1:end-1);
s2 = gam(end);

pom = -k/2*log(2*pi)-k/2*log(s2)+(-1/2*1/s2*(y-X*b)'*(y-X*b)); 
logL = -pom;

