CCmed_gene: Estimate cross-condition mediation probability at the...

In kjgleason/CCmed: Cross-condition mediation analysis

Description

For each cis-trans gene pair, estimate the probability of mediation.

Usage

CCmed_gene(cis_res, condcorr_res, cisID, condcorrID, ncond)

Arguments

cis_res	list object returned by running Primo on gene-level cis-association statistics.
condcorr_res	list object returned by running Primo on conditional correlation statistics.
cisID	vector of gene names corresponding that match the identifiers of of cis_res (e.g. rows in list element post_prob).
condcorrID	data.frame with two columns ("cis", "trans") that match the identifiers of condcorr_res (i.e. rows in list element post_prob).
ncond	integer of the number of conditions required for mediation.

Value

A data.frame with the following columns:

cis	cis-gene identifier
trans	trans-gene identifier
ncond	number of conditions required for mediation
prob_cis_all	posterior probability of cis-association in all conditions
prob_condcorr_ncond	posterior probability of conditional correlation in at least ncond conditions
prob_med	probability of mediation in at least ncond conditions

The main element of interest for inference is the column prob_med.

kjgleason/CCmed documentation built on Nov. 14, 2019, 12:14 a.m.