ComputePValues: Compute p-values for affinity scores.

In kinsigne/atSNP_modified: Affinity test for identifying regulatory SNPs.

Description

This function computes the p-values for allele-specific affinity scores and between-allele affinity score changes using the importance sampling technique.

Usage

ComputePValues(motif.lib, snp.info, motif.scores, ncores = 1, figdir = NULL)

Arguments

motif.lib	A list object with the output format of function LoadMotifLibrary.
snp.info	A list object with the output format of function LoadSNPData.
motif.scores	A data.table object containing at least the following columns: motif The name of the motif. log_lik_ref The log-likelihood score for the reference allele. log_lik_snp The log-likelihood score for the SNP allele.
ncores	An integer for the number of parallel process. Default: 1.
figdir	A string for the path to print p-value plots for monitoring results. Default: NULL (no figure).

Value

A data.table extending motif.scores by the following additional columns:

pval_ref	P-values for scores on the reference allele.
pval_snp	P-values for scores on the SNP allele.
pval_cond_ref	Conditional p-values for scores on the reference allele.
pval_cond_snp	Conditional p-values for scores on the SNP allele.
pval_diff	P-values for the difference in scores between the reference and the SNP alleles.
pval_rank	P-values for the log rank ratio between the reference and the SNP alleles.

Author(s)

Chandler [e-mail]

Examples

data(example)
ComputePValues(motif_library, snpInfo, motif_scores$motif.scores, ncores = 4)

kinsigne/atSNP_modified documentation built on May 23, 2019, 4:24 p.m.