Re: [msprime-users] sampling same population from different time

[Jerome Kelleher]

Hi,

This looks fine to me. There's definitely no problem with having samples
at different times in the same population like this. Drawing the tree
might help a bit for debugging, as you'll see the leaf nodes at
different times.

Cheers,
Jerome

On 09/06/2020 00:38, molkoju wrote:
> Hi,
>
> I want to ask question about if it is possible to simulate samples from
> same population at different time period (i.e. 10 000, 8 000, 6 000, 4
> 000, 2 000 before present). I tried to do that by using samples function
> like below.
>
> |
> samples = (
>     [msprime.Sample(population=0, time=0) for _ in range(Ref)] +
>     [msprime.Sample(population=1, time=0) for _ in range(PopA)] +
>     [msprime.Sample(population=2, time=10000/gen) for _ in
> range(PopB//6)] +
>     [msprime.Sample(population=2, time=8000/gen) for _ in range(PopB//6)] +
>     [msprime.Sample(population=2, time=6000/gen) for _ in range(PopB//6)] +
>     [msprime.Sample(population=2, time=4000/gen) for _ in range(PopB//6)] +
>     [msprime.Sample(population=2, time=2000/gen) for _ in range(PopB//6)] +
>     [msprime.Sample(population=2, time=0) for _ in range(PopB//6)] +
>     [msprime.Sample(population=3, time=0) for _ in range(PopC)])
> |
>
> Unfortunately I could not check whether I am doing right or not neither
> by using debugger or drawing tree.
>
> I am really new with msprime (and all simulations). Can anyone help me
> with this?
>
> Thanks,
>
> M
>
> --
> You received this message because you are subscribed to the Google
> Groups "msprime-users" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to msprime-user...@googlegroups.com
> <mailto:msprime-user...@googlegroups.com>.
> To view this discussion on the web, visit
> https://groups.google.com/d/msgid/msprime-users/bcc02a91-32e3-4322-b989-dd9a8dde289do%40googlegroups.com
> <https://groups.google.com/d/msgid/msprime-users/bcc02a91-32e3-4322-b989-dd9a8dde289do%40googlegroups.com?utm_medium=email&utm_source=footer>.