Random generation for the \(\pi_0\) distribution.

Random generation for the \(\pi_0\) distribution.

rpi0(n, n_targets, n_decoys)

Arguments

n	number of observations.
n_targets	number of observed target PSMs.
n_decoys	number of observed decoy PSMs.

Value

vector of random deviates. The length equals `n'.

Examples

## visualize the pi0 distribution when observing 10 targets and 3 decoys
x = rpi0(100000, 10 ,3 )
hist(x, breaks = 50 , xlab = 'pi0', ylab = 'counts')