More Drossel fine print regarding fixation time: MatLab demo

Submitted by David Liao on Tue, 2010/02/02 - 22:50

You can try out FixationTimeDemor001.tgz using these commands.

>> [rme,tfix,s] = WrightEnsemble(100,20,100,0.1);
>> WrightReplot(rme,17,s);

In FixationTimeDemoR002.tgz, I added a command to look at metapopulations with subpopulation interbreeding:

>> WrightMetapopulation(10,20,100,0,0.003,0.2,0)

Here is some example output:

Figure: Ensemble of populations from WrightEnsemble

Figure: Zoom-in on a particular population using WrightReplot

Figure: Metapopulation using WrightMetapopulation

Attachment	Size
2010 June 26: FixationTimeDemoR003.tgz	3.15 KB
"Alpha" draft 2010 Feb 04 01h46: FixationTimeDemoR002.tgz	3.15 KB
"Alpha" draft 2010 Feb 02 22h51: FixationTimeDemor001.tgz	1.89 KB

lookatphysics.com

	Copyright 2008-2011 by David Liao
	You can view some content by visiting this website without logging in. You are welcome to distribute and reuse these materials using a Creative Commons Attributions-Share Alike 3.0 US License (http://creativecommons.org/licenses/by-sa/3.0/us/). However, content requiring a login to be viewed on this site includes unpublished data and notes from David Liao and other authors. These communications are private. To login or logout, click login or logout. The main index is here.