Linear mixed model to study multivariate genotype-environment interactions
Project description
Struct-LMM
Structured Linear Mixed Model (StructLMM) is a computationally efficient method to test for and characterize loci that interact with multiple environments [1].
This a standalone module that implements the basic functionalities of StructLMM. However, we recommend using StructLMM via LIMIX2 as this additionally implements:
- Multiple methods for GWAS;
- Methods to characterize GxE at specific variants;
- Command line interface.
Install
From terminal, it can be installed using pip:
pip install struct-lmm
Usage
>>> from numpy import ones, concatenate
>>> from numpy.random import RandomState
>>>
>>> from struct_lmm import StructLMM
>>>
>>> random = RandomState(1)
>>> n = 20
>>> k = 4
>>> y = random.randn(n, 1)
>>> E = random.randn(n, k)
>>> M = ones((n, 1))
>>> x = 1.0 * (random.rand(n, 1) < 0.2)
>>>
>>> lmm = StructLMM(y, M, E)
>>> lmm.fit(verbose=False)
>>> # Association test
>>> pv = lmm.score_2dof_assoc(x)
>>> print(pv)
0.8470017313426488
>>> # Association test
>>> pv, rho = lmm.score_2dof_assoc(x, return_rho=True)
>>> print(pv)
0.8470017313426488
>>> print(rho)
0
>>> M = concatenate([M, x], axis=1)
>>> lmm = StructLMM(y, M, E)
>>> lmm.fit(verbose=False)
>>> # Interaction test
>>> pv = lmm.score_2dof_inter(x)
>>> print(pv)
0.6781100453132024
Problems
If you encounter any problem, please, consider submitting a new issue.
Authors
License
This project is licensed under the MIT License.
[1] Moore, R., Casale, F. P., Bonder, M. J., Horta, D., Franke, L., Barroso, I., & Stegle, O. (2018). A linear mixed-model approach to study multivariate gene–environment interactions (p. 1). Nature Publishing Group.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for struct_lmm-0.3.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 3a2d8560e716ce4c6be0f6f2d9adf3015f9f04e456f97e01efd7ee368dd8f047 |
|
MD5 | 5d35737678c3c50b7f9a258388dd7a96 |
|
BLAKE2b-256 | 5bebc7076b118dd99778d2619b0d3408e1dbe96616b765bcecf66222539194e7 |