R-based Log-Likelihood Implementation
loglik_R(theta, W, X, Y, weights, nReg, ...)
loglik_cpp, opsr.fit
loglik_cpp
opsr.fit