fix fit_reml! to return 6 values as before

juhyun-kim-ucla · Hua-Zhou · commit 1cbf9b0ab8d4 · 2019-01-21T20:27:54.000-08:00
diff --git a/src/two_variance_component.jl b/src/two_variance_component.jl
@@ -1379,8 +1379,7 @@ function fit_reml!(
   copyto!(Bse, sqrt.(diag(Bcov)))
 
   ## output
-  # logpdf(vcmodel, vcdatarot), vcmodel, Σse, Σcov, Bse, Bcov
-  logpdf(vcmodel, vcdatarot), Σcov, Bse, Bcov
+  logpdf(vcmodel, vcdatarot), vcmodel, (Σse[1], Σse[2]), Σcov, Bse, Bcov
 end
 
 #---------------------------------------------------------------------------#
