integrand: Function generates the required integrand for the test.

Description Usage Arguments Value

View source: R/ANEVA_DOT.R

Description

Function generates the required integrand for the test.

Usage

1
integrand(dE, eh1, eh2, Eg_std, r0, p0, log_BinCoeff)

Arguments

dE

Variable representing the true change in gene dosage (variable of integration)

eh1

Integer value for count of expression of haplotype1

eh2

Integer value for count of expression of haplotype2

Eg_std

Vector containing genetic standard deviation, the square root of genetic variation in regulation in general population in natural log scale (numeric value).

r0

The ratio of the eh1 allele (i.e., eh1/(eh1+eh2)) in the absence of any regulatory difference (reference bias due to alignment). The simplest way to get such an estimate would be to get the median ratio between eh1 and eh2 across the entire library.

p0

An average noise rate p(R->A) or p(A->R), i.e., the probability of seeing an allele due to noise when it is essentially not there. (For v7: LAMP = 0.0003).

log_BinCoeff

a vector of natural log (ln) transformed binomial coefficients 0:N=eh1+eh2.

Value

The integrand over which we wish to integrate to get desired p-values


PejLab/ASEtestTools documentation built on Oct. 11, 2021, 12:10 p.m.