MSMC: very high Ne close to present
Entering edit mode
3.0 years ago
evoecogen ▴ 20


I have a weird problem with MSMC2. I use it to study demography of a non-model species, with populations diverged over the last 100,000-1M years (the times are deeper than for the usual human work).

  1. In almost all cases, MSMC based on SHAPEIT-phased haplotypes recovers ridiculous increases in Ne close to the present, like an order of magnitude boom in the last 10Kyr (last 4-5 slice of time, corresponding to the first part of the SMC command). This seems very unlikely and does not correlate with differences between the populations (plus the inferred Ne values are way to high for this animal). The middle of the curves looks reasonable, and consistent with PSMC on unphased data. The only solution I have had for this is to tie the five most recent slices together, so the Ne stays flat.: -p 1 4+25 1+1 2+1 5 [markodwn: there are normally asterisks instead of the spaces.]

  2. Similar weird thing happens in the deep past - the inferred Ne for the first time slice is ridiculously high. I don't even understand how that happens with the coalescent.

What could cause these behaviours? I don't care about the very recent or the very deep past, but why is this happening?...


PSMC demography popgen MSMC SMC • 975 views
Entering edit mode
2.7 years ago
clinnaeus ▴ 30

Hi evo,

Have you had any thoughts on this since you posted? I'm experiencing very similar problematic patterns with a PSMC analysis; huge upshoots in Ne that surpass 1e06. I'm struggling to find help on the matter, nor can I interpret what I've generated as I'm worried it is a technical artefact of some kind.

Cheers :)

Entering edit mode


This may be too late for you, but last year I found that this is an expected behaviour of the SMC methods. PSMC is not to be trusted when it comes to very recent past. It is common to just ignore the first and the last few time slices (this is what I did for my MSMC)


Login before adding your answer.

Traffic: 1629 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6