Re: [simuPOP-list] Bug? in the allele frequency

[Bo P. <be...@gm...>]

Hi, Hugo,

I have a quick look at your script and find that your mating scheme uses a
PyOperator as ops. This operator set offspring info but does not pass
parental genotypes. Please add a MedelianGenoTransmitter or another
genotype transmitter before the PyOperator and let me know if this solves
your problem.

Cheers,
Bo

On Sun, Aug 28, 2016 at 11:59 AM, hmh <hu...@gm...> wrote:

> Hi,
>
>
> I have a bug in the allele frequency, I don't understand what is wrong in
> my script, and the bug persist even with the simplest situation '1 locus, 2
> alleles, no mutation':
> *Every two generations the allele frequency is set to 0*.
>
>
> For every generation, the attached script print the allele frequency of
> the population, and the genotype of the first 5 individuals as well as of
> their 2 parents.
>
> (Genotypes are 0, 0.5 and 1 where 0.5 means heterozygous)
>
>
> Probably because of the same bug, their often inconsistency such as:
>
>
> at the generation N
>
> Off:1.0 Dad:0.0 Mom:0.0
>
>
> at the generation N+1
>
> Off:0.0 Dad:1.0 Mom:1.0
>
>
> I am using this simupop version :
>
> # simuPOP Version 1.1.7.1 : Copyright (c) 2004-2016 Bo Peng# Revision 5003 (May 22 2016) for Python 3.5.1+ (64bit, 1thread)# Random Number Generator is set to mt19937 with random seed 0x91f07829f57cb62.# This is the standard short allele version with 256 maximum allelic states.# For more information, please visit http://simupop.sourceforge.net,# or email sim...@li... (subscription required).
>
>
> Thanks for your help!
>
> hugo
>
>
> --
>
> Hugo Mathé-Hubert
>
> BU-G19
>
> postdoc
>
> eawag (Swiss Federal Institute of Aquatic Science and Technology)
> Evolutionary Ecology
> <http://www.eawag.ch/en/department/eco/main-focus/evolutionary-ecology/>
> - About me
> <http://www.eawag.ch/en/aboutus/portrait/organisation/staff/profile/hugo-mathe-hubert/>
>
>
>
> Überlandstrasse 133
> P.O.Box 611
> 8600 Dübendorf, Switzerland
>
> - - - - - - - - - - - - - - - - - -
>
> Thoughts appear from doubts and die in convictions. Therefore, doubts are
> an indication of strength and convictions an indication of weakness. Yet,
> most people *believe* the opposite.
>
> - - - - - - - - - - - - - - - - - -
>
> Les réflexions naissent dans les doutes et meurent dans les certitudes.
> Les doutes sont donc un signe de force et les certitudes un signe de
> faiblesse. La plupart des gens sont pourtant certains du contraire.
>
>