The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.
The lit package implements a kernel-based multivariate
testing procedure, called Latent Interaction Testing (LIT), to test for
latent genetic interactions in genome-wide association studies. See our
manuscript for additional details:
Bass AJ, Bian S, Wingo AP, Wingo TS, Culter DJ, Epstein MP. Identifying latent genetic interactions in genome-wide association studies using multiple traits. Submitted; 2023.
install.packages("devtools")
library("devtools")
# install package
devtools::install_github("ajbass/lit")The vignette can be viewed by typing:
browseVignettes(package = "lit")We provide two ways to use the lit package. For small
GWAS datasets where the genotypes can be loaded in R, the
lit() function can be used:
library(lit)
# set seed
set.seed(123)
# generate SNPs and traits
X <- matrix(rbinom(10 * 10, size = 2, prob = 0.25), ncol = 10)
Y <- matrix(rnorm(10 * 4), ncol = 4)
# test for latent genetic interactions
out <- lit(Y, X)
head(out)
#> wlit ulit alit
#> 1 0.2681410 0.3504852 0.3056363
#> 2 0.7773637 0.3504852 0.6044655
#> 3 0.4034423 0.3504852 0.3760632
#> 4 0.7874949 0.3504852 0.6157108
#> 5 0.8701189 0.3504852 0.7337565
#> 6 0.2352616 0.3504852 0.2847600The output is a data frame of p-values where the rows are SNPs and
the columns are different implementations of LIT to test for latent
genetic interactions: the first column (wlit) uses a linear
kernel, the second column (ulit) uses a projection kernel,
and the third column (alit) maximizes the number of
discoveries by combining the p-values of the linear and projection
kernels.
For large GWAS datasets (e.g., biobank-sized), the lit()
function is not computationally feasible. Instead, the
lit_plink() function can be applied directly to plink
files. To demonstrate how to use the function, we use the example plink
files from the genio package:
# load genio package
library(genio)
# path to plink files
file <- system.file("extdata", 'sample.bed', package = "genio", mustWork = TRUE)
# generate trait expression
Y <- matrix(rnorm(10 * 4), ncol = 4)
# apply lit to plink file
out <- lit_plink(Y, file = file, verbose = FALSE)
head(out)
#> chr id pos alt ref maf wlit ulit alit
#> 1 1 rs3094315 752566 G A 0.3888889 0.7908763 0.3422960 0.6150572
#> 2 1 rs7419119 842013 T G 0.3888889 0.1552580 0.3422960 0.2194972
#> 3 1 rs13302957 891021 G A 0.2500000 0.4088937 0.3325939 0.3687589
#> 4 1 rs6696609 903426 C T 0.3125000 0.5857829 0.3325939 0.4519475
#> 5 1 rs8997 949654 A G 0.4375000 0.6628300 0.3325939 0.4969663
#> 6 1 rs9442372 1018704 A G 0.2500000 0.3192430 0.3325939 0.3258332See ?lit and ?lit_plink for additional
details and input arguments.
Note that a marginal testing procedure for latent genetic
interactions based on the squared residuals and cross products (Marginal
(SQ/CP)) can also be implemented using the marginal and
marginal_plink functions:
# apply Marginal (SQ/CP) to loaded genotypes
out <- marginal(Y, X)
# apply Marginal (SQ/CP) to plink file
out <- marginal_plink(Y, file = file, verbose = FALSE)These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.