The outlier distribution we use. Right now it is just a beta binomial with mean 1/2 and od 1/3 (so underlying beta is just a uniform from 0 to 1).
doutdist(x, n, logp)
The number of reference counts.
The total number of read-counts.
Return the log density TRUE
or not FALSE
?
A double. The outlier density value.