Helper function for SQUAREM for extracting log likelihood
LL_SQUAREM(theta_vec, ...)
theta_vec
Vector of current parameter values
Negative log-likelihood given current values of parameters