powered by
wrightFst(Ne, dispersal, gen.time, ploidy)
dispersal <- seq(0.05, 0.8, by = 0.05) fst <- wrightFst(100, dispersal, 20, 2) plot(dispersal, fst, type = "l")
Run the code above in your browser using DataLab