mvtBifact_nllk function

negative log-likelihood for the bi-factor Gaussian/t model