Calculate influential genes for a given trait and cell type using DFBETAS.

Calculate influential genes for a given trait and cell type using DFBETAS.

Usage

find_inf_genes(
  ct,
  sscore,
  magma,
  magma_gene_col = "GENE",
  magma_z_col = "ZSTAT"
)

Arguments

ct

A character string containing a valid cell type name in sscore.

sscore

A dgeMatrix of seismic specificity scores where each column is a cell type and row names are gene identifiers. (Note: the identifiers used should match those used in the MAGMA input)

magma

A data.frame or file path to MAGMA output for a particular GWAS with at least 2 columns: gene identifiers and z-scores.

magma_gene_col

A character string corresponding to the column name in the MAGMA data containing gene identifiers. Defaults to 'GENE' from the MAGMA output.

magma_z_col

A character string corresponding to the column name in the MAGMA data containing z-scores. Defaults to 'ZSTAT' from the MAGMA output.

Value

A data.frame containing genes, seismic specificity scores, magma z-stats, the dfbetas influential gene score, and a Boolean value indicating if the gene is influential.