Skip to content

Commit

Permalink
MM-89: incorporate CRAN feedback (#95)
Browse files Browse the repository at this point in the history
  • Loading branch information
jdstamp authored Dec 2, 2022
1 parent c6e8133 commit c8794f2
Show file tree
Hide file tree
Showing 2 changed files with 26 additions and 18 deletions.
38 changes: 23 additions & 15 deletions DESCRIPTION
Original file line number Diff line number Diff line change
Expand Up @@ -8,21 +8,29 @@ Authors@R: c(
role = c("cre", "aut"), comment = c(ORCID = "0000-0003-3014-6249")),
person("Lorin", "Crawford", email = "[email protected]",
role = "aut", comment = c(ORCID = "0000-0003-0178-8242")))
Description: Epistasis, commonly defined as the interaction between genetic loci, is known to play an important role in
the phenotypic variation of complex traits. As a result, many statistical methods have been developed to identify
genetic variants that are involved in epistasis, and nearly all of these approaches carry out this task by focusing on
analyzing one trait at a time. Previous studies have shown that jointly modeling multiple phenotypes can often
dramatically increase statistical power for association mapping. In this study, we present the
'multivariate MArginal ePIstasis Test' (mvMAPIT) – a multi-outcome generalization of a recently proposed epistatic
detection method which seeks to detect marginal epistasis or the combined pairwise interaction effects
between a given variant and all other variants. By searching for marginal epistatic effects, one can identify genetic
variants that are involved in epistasis without the need to identify the exact partners with which the variants
interact – thus, potentially alleviating much of the statistical and computational burden associated with
conventional explicit search based methods. Our proposed mvMAPIT builds upon this strategy by taking advantage of
correlation structure between traits to improve the identification of variants involved in epistasis. We formulate mvMAPIT as a
multivariate linear mixed model and develop a multi-trait variance component estimation algorithm for efficient
parameter inference and P-value computation. Together with reasonable model approximations, our proposed approach is
scalable to moderately sized GWA studies.
Description: Epistasis, commonly defined as the interaction between genetic
loci, is known to play an important role in the phenotypic variation of
complex traits. As a result, many statistical methods have been developed to
identify genetic variants that are involved in epistasis, and nearly all of
these approaches carry out this task by focusing on analyzing one trait at a
time. Previous studies have shown that jointly modeling multiple phenotypes
can often dramatically increase statistical power for association mapping. In
this package, we present the 'multivariate MArginal ePIstasis Test'
('mvMAPIT') – a multi-outcome generalization of a recently proposed epistatic
detection method which seeks to detect marginal epistasis or the combined
pairwise interaction effects between a given variant and all other variants.
By searching for marginal epistatic effects, one can identify genetic variants
that are involved in epistasis without the need to identify the exact
partners with which the variants interact – thus, potentially alleviating
much of the statistical and computational burden associated with conventional
explicit search based methods. Our proposed 'mvMAPIT' builds upon this
strategy by taking advantage of correlation structure between traits to
improve the identification of variants involved in epistasis.
We formulate 'mvMAPIT' as a multivariate linear mixed model and develop a
multi-trait variance component estimation algorithm for efficient parameter
inference and P-value computation. Together with reasonable model
approximations, our proposed approach is scalable to moderately sized
genome-wide association studies.
Crawford et al. (2017) <doi:10.1371/journal.pgen.1006869>.
Stamp et al. (2022) <doi:10.1101/2022.11.30.518547>.
License: GPL (>= 3)
Expand Down
6 changes: 3 additions & 3 deletions vignettes/tutorial-docker-mvmapit.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -30,8 +30,6 @@ docker build -t mvmapit .

This will produce an image named `mvmapit` that contains `Rstudio`, `mvMAPIT`, and all dependencies.

Alternatively, the image will be made available through Dockerhub.

## Run the mvMAPIT Image

With a local copy of the docker image `mvmapit` available, run the following code.
Expand All @@ -49,5 +47,7 @@ container, `mvMAPIT` is already installed and can be imported and run in the R c

```{r run_mvmapit}
library(mvMAPIT)
mvmapit(t(simulated_data$genotype[1:100,1:10]), t(simulated_data$trait[1:100,]), cores = 4, logLevel = "DEBUG")
mvmapit(t(simulated_data$genotype[1:100,1:10]),
t(simulated_data$trait[1:100,]),
cores = 2, logLevel = "DEBUG")
```

0 comments on commit c8794f2

Please sign in to comment.