MATS | R Documentation |
This function implements Multi-Ancestry TWAS (MATS), based on the framework outlined in Knutson & Pan, 2022
MATS( Y, xhat, groups, C, P, ev, K, categorical.vars = NULL, trait_type = "binary", np = 5 )
Y |
a vector of either continuous or binary trait values for n subjects. |
xhat |
a numeric vector of imputed gene expression for n subjects. |
groups |
a vector giving ancestral group labels for the n subjects. Must include at least 2 groups. |
C |
an optional matrix of covariates for the n subjects (n rows). If any of the covariates are categorical, specify their columnn names in the categorical.vars parameter. Dp not include a column for the intercept. |
P |
a matrix of genetic principal component on the n subjects (n rows). Must have at least nP PCs (at least nP columns). Columns should be orders PC1, PC2, ... and should not include any other columns aside from PCs (i.e., no subject ID column). |
ev |
a numeric vector of ordered eigenvalues for the genetic principal component. Should be the same length as the number of columns in P. |
K |
an n x n genetic covariance matrix. |
categorical.vars |
a character vector of categorical variable column names in C (i.e., c("variable1", "variable2")) |
trait_type |
Either "Continuous" or "Binary", for continous or binary disease traits |
nP |
the number of principal components to include in all models to adjust for population stratification, and the number used in the full PC-interaction model. |
n |
The sample size used to estimate the disease GWAS summary data |
n_case |
The number of cases used in the disease GWAS sample. Required for trait_type = "Binary", default for Continuous trait is NULL |
n_control |
The number of controls used in the disease GWAS sample. Required for trait_type = "Binary", default for Continuous trait is NULL |
MATS()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.