MATS: MATS: Performing Multi-Ancestry TwaS

View source: R/MATSfunction.R

MATSR Documentation

MATS: Performing Multi-Ancestry TwaS

Description

This function implements Multi-Ancestry TWAS (MATS), based on the framework outlined in Knutson & Pan, 2022

Usage

MATS(
  Y,
  xhat,
  groups,
  C,
  P,
  ev,
  K,
  categorical.vars = NULL,
  trait_type = "binary",
  np = 5
)

Arguments

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

Examples

MATS()

kathalexknuts/MATS documentation built on July 5, 2022, 10:36 p.m.