admixture.txt

admixtureGeneration = gen(15,2); # this is the generation of the admixture event
admixtureGenPlus1 = gen(admixtureGeneration,1);

nodeEuToAA = node(admixtureGeneration);
nodeAfToAA = node(admixtureGeneration);
nodeAA = node(admixtureGeneration);

sizeEuAtAdmixture = size(100000,1);
sizeAfAtAdmixture = size(70000,1);
ancestralSizeAA = size(10000,2);

EuropeansPostAdmixture = edge(node(G0), nodeEuToAA, size(1000000), size(98000,1));
AfricansPostAdmixture = edge(node(G0), nodeAfToAA, size(700000,1), size(62000,1));
AfricanAmericans = edge(node(G0), nodeAA, size(1000000,2), ancestralSizeAA);
migrationBranchEU = edge(nodeAA, nodeEuToAA, size(2000,2));
migrationBranchAf = edge(nodeAA, nodeAfToAA, size(8000,2));


nodeSplitAfEU = node(gen(400,1));
ancientBranchEU = edge(nodeEuToAA, nodeSplitAfEU, sizeEuAtAdmixture, size(50000,1));
ancientBranchAf = edge(nodeAfToAA, nodeSplitAfEU, sizeAfAtAdmixture, size(50000,1));

nodeAncestralAll = node(gen(600));
ancestralSize = size(10000);
edge(nodeSplitAfEU, nodeAncestralAll, size(100000,1), ancestralSize);
edge(nodeAncestralAll, node(Ginf), ancestralSize);



rateAf-AA=rate(0.001,1);

migration(EuropeansPostAdmixture, AfricansPostAdmixture, rateAf-AA );