Fast inference over mean and covariance parameters for Generalised Linear Mixed Models.
It implements the mathematical tricks of FaST-LMM for the special case of Linear Mixed Models with a linear covariance matrix and provides an interface to perform inference over millions of covariates in seconds. The Generalised Linear Mixed Model inference is implemented via Expectation Propagation and also makes use of several mathematical tricks to handle large data sets with thousands of samples and millions of covariates.
Via pip:
pip install glimix-coreHere it is a very simple example to get you started:
>>> from numpy import array, ones
>>> from numpy_sugar.linalg import economic_qs_linear
>>> from glimix_core.lmm import LMM
>>>
>>> X = array([[1, 2], [3, -1], [1.1, 0.5], [0.5, -0.4]], float)
>>> QS = economic_qs_linear(X, False)
>>> X = ones((4, 1))
>>> y = array([-1, 2, 0.3, 0.5])
>>> lmm = LMM(y, X, QS)
>>> lmm.fit(verbose=False)
>>> lmm.lml()
-2.2726234086180557We also provide an extensive documentation about the library.
This project is licensed under the MIT License.