This function simulates 95% confidence level Z-score based detection power of allele biases for a given number of samples and a range of depths
power.bias(
Dlist = c(2, 4, 8, 16),
sam = 100,
intensity = 0.005,
nsims = 1000,
p = 0.5,
plot = TRUE
)Returns a list of detection probability values for the given range of samples and depth
numerical. vector of depths values to be tested
numerical. number of samples
numerical. frequency of bias
numerical. number of simulations to be done for each sample
numerical. expected allele ratio (0.5 for data with known sequencing biases)
logical. plot the output
Pascal Milesi, Piyal Karunarathne