function [newlik] = likelihood(lambda, lambda_t, psi, inv_psi, lambda_psi_inv, s)

newlik = -trace(lambda_psi_inv*s);

newlik = newlik - (sum(log(diag(psi)))+log(det(inv_psi*lambda*lambda'+eye(560))));