Skip to contents

An association test in a tissue-specific manner for a phenotype (trait).

Usage

twas.single.trait(stats, x, list.ts.eQTL)

Arguments

stats

summary statistics for a phenotype. It can be a filename to be read or a data.frame. Suppose it has 5 columns with colnames rsid, a1, a2, chr, z.

x

twas data list

list.ts.eQTL

a list of ts.eQTL, each element corresponds to a tissue

Value

z-statistics and p-values. Both are of size n_gene x n_tissue.