Re: [R-sig-phylo] Mismatch distribution (pegas): expected versus empirical and other questions

2013-12-11 Thread Emmanuel Paradis
david.win...@gmail.com Sender: r-sig-phylo-boun...@r-project.org Date: Tue, 10 Dec 2013 12:47:39 To: Alastair Pottspott...@gmail.com Cc: mailman, r-sig-phylor-sig-phylo@r-project.org Subject: Re: [R-sig-phylo] Mismatch distribution (pegas): expected versus empirical and other questions Hi Alistair, I

[R-sig-phylo] Mismatch distribution (pegas): expected versus empirical and other questions

2013-12-10 Thread Alastair Potts
Hi all, The mismatch distribution function in pegas (MMD) produces the familiar histogram and then also a line that I always took as the theoretical expected distribution under a scenario of constant population size. I had a quick look at the MMD function and noticed that this line is not

Re: [R-sig-phylo] Mismatch distribution (pegas): expected versus empirical and other questions

2013-12-10 Thread David Winter
Hi Alistair, I wrote a function to make the standard expectation under constant population size quite a while ago. It's up here https://github.com/dwinter/Pegas/blob/new_functions/R/MMD.stable.R It's been a little while since I've had to think about mismatch distributions, so, though I'm sure