gammaIntegral:
gammaIntegral: get the conditional expectation of a convolved gamma RV
Description
Conditional expectation of one (out of two) gamma RVs given their sum & params
Usage
gammaIntegral(total, params, offset=50, minx=1)
Arguments
total
a vector of sums (convolutions)
params
parameters for the signal and background distributions
offset
an offset to be added to the estimate
minx
a minimum value below which we will replace x
Value
The conditional expectation(s) of the signal given the total.
Details
See Triche et al., Nucleic Acids Research (submitted).