cherry pick only user-relevant functions for API

Author: mmkim1210

README.md

@@ -66,8 +66,8 @@ diag(model.Σcov) # m * (binomial(d, 2) + d) parameters
 model.logl
 ```
 ## References
-- H. Zhou, L. Hu, J. Zhou, and K. Lange: **MM algorithms for variance components models** (2019) ([link](https://doi.org/10.1080/10618600.2018.1529601))
-- M. Kim: **Gene regulation in the human brain and the biological mechanisms underlying psychiatric disorders** (2022) ([link](https://escholarship.org/uc/item/9v08q5f7))
+- <u>H. Zhou, L. Hu, J. Zhou, and K. Lange: **MM algorithms for variance components models** (2019) ([link](https://doi.org/10.1080/10618600.2018.1529601))</u>
+- <u>M. Kim: **Gene regulation in the human brain and the biological mechanisms underlying psychiatric disorders** (2022) ([link](https://escholarship.org/uc/item/9v08q5f7))</u>
 
 ## See also
 - J. Kim, J. Shen, A. Wang, D.V. Mehrotra, S. Ko, J.J. Zhou, and H. Zhou: **VCSEL: Prioritizing SNP-set by penalized variance component selection** (2021) ([link](http://doi.org/10.1214/21-aoas1491))

docs/src/api.md

@@ -3,6 +3,9 @@
 ```@index
 ```
 
-```@autodocs
-Modules = [MultiResponseVarianceComponentModels]
+```@docs
+MultiResponseVarianceComponentModels
+MRVCModel
+MRTVCModel
+fit!
 ```

src/parse.jl

@@ -1,5 +1,5 @@
 """
-    lrt(model1::MRVCModel, model0::MRVCModel)
+    lrt(model1::VCModel, model0::VCModel)
 
 Perform a variation of the likelihood ratio test for univariate variance components models as in 
 Molenberghs and Verbeke 2007 with model1 and model0 being the full and nested models, respectively.