Get a vector of mutation rates from a gamma distribution by specifying the mean and std. dev. of the distribution
getGammaMutRates(n, gmean = 1e-04, gstd = 1e-05)
Number of loci to select rates for
Mean of the gamma distribution
Std Dev of the gamma distribution
a vector of mutation rates (type numeric)
# NOT RUN { rates = getGammaMutRates(1000,gmean=0.0001, gstd=0.0001) hist(rates) rates = getGammaMutRates(1000,gmean=0.0001, gstd=0.00001) hist(rates) # }
Run the code above in your browser using DataCamp Workspace