Hi community,
I have a dataset of individuals of a species (several populations) genotyped at 9 microsatellites.
I would like to ask if (and how) I could simulate that STR data forward in time, so to implement different scenarious (e.g. cut-off of one population, varying migration rates between populations, etc)?
Ideally, the simulated data (forwarded for X years/generations) would be in allele frequencies or even STR data format as well, so that I can compare the initial and the forward-simulated data using standard procedures of population genetics. It seems to me that some programs are able to simulate STR-data forward in time (e.g. R package skeleSim using Rmetasim), yet as initial data they also use simulated data. So I think it should be able to use real STR data instead of the simulated one? But please correct me if I am wrong here.
Your help is appreciated,
Florian