powered by
simGammaHaps(pop.size, num.haps, shape, scale, return.freq = TRUE, plot = TRUE)
dgamma
FALSE
haps <- simGammaHaps(1000, 15, 1, 2.5) print(haps)
Run the code above in your browser using DataLab