function sfran_loglike(mu,sigma2_mu,sigma2,y) d, ni = size(y) Sigmai = sigma2*diagm(ones(ni)) .+ sigma2_mu*ones(ni,ni) l = -(ni*d)/2*log(2*pi) - d/2*log(det(Sigmai)) for i=1:d yi = y[i,:]; l = l - .5*(yi .- mu)'*(Sigmai\(yi .- mu)) end return l end