bugNumerical Library for Haploid Genetics - Bugs: bug #30720, Allele frequencies changing...

 
 

bug #30720: Allele frequencies changing without selection or mutation

Submitted by:  None
Submitted on:  Tue Aug 10 00:02:57 2010  
 
Category: LibrarySeverity: 3 - Normal
Item Group: NoneStatus: In Progress
Privacy: PublicAssigned to: Joel J. Adamson <trashbird1240>
Open/Closed: Open

Add a New Comment (Rich MarkupRich Markup):
   

You are not logged in

Please log in, so followups can be emailed to you.

 

(Jump to the original submission Jump to the original submission)

Mon Dec 6 16:38:34 2010, comment #17:

-- mode: compilation; default-directory: "~/Public/src/haploid/" --
Compilation started at Thu Dec 2 13:28:16

make -k check
make sim_stop pop_ck sparse_test diseq rec_test
make[1]: Entering directory `/home/joel/Public/src/haploid'
make[1]: `pop_ck' is up to date.
make[1]: `sparse_test' is up to date.
make[1]: `diseq' is up to date.
make[1]: `rec_test' is up to date.
make[1]: Leaving directory `/home/joel/Public/src/haploid'
make check-TESTS
make[1]: Entering directory `/home/joel/Public/src/haploid'
PASS: sim_stop
PASS: pop_ck
PASS: sparse_test
Change in allele frequency detected:
Starting frequency: p[0] = 0.50000000000000000000000000000000000000000000000000000
New frequency: p[0] = 0.50000000000001088018564132653409615159034729003906250
[...]
/bin/sh: line 5: 11068 Aborted (core dumped) ${dir}$tst
FAIL: rec_test
Value by hand: 0.170442
Value by ld_from_geno: 0.170442
PASS: diseq
=========================================
1 of 5 tests failed
Please report to -unavailable-
=========================================
make[1]: *** [check-TESTS] Error 1
make[1]: Leaving directory `/home/joel/Public/src/haploid'
make: *** [check-am] Error 2
make: Target `check' not remade because of errors.

Compilation exited abnormally with code 2 at Fri Dec 3 04:31:32

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Thu Dec 2 18:35:59 2010, comment #16:

With 9 loci there is a difference at the ninth decimal place. At this point it looks like a floating point error and we will have to conduct further tests to determine if it is significant.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Thu Sep 2 12:33:22 2010, comment #15:

With 9 loci (= 512 genotypes) the test failed, printing

We can conclude that this is a floating-point error.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Wed Sep 1 19:51:05 2010, comment #14:

The test as it is goes pretty close the limits of published recombination algorithms (with r = 0.5), so we may be waiting a while for the test to finish.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Wed Sep 1 19:46:20 2010, comment #13:

Changes in revision 1112 fix this. Revision 1113 includes a thorough test. We can close this as soon as the test finishes on my machine.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Mon Aug 30 00:15:45 2010, comment #12:

We need to track the parental genome we are starting with. Use the branching array design to keep track of genuine branch points (true "maybes").

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 24 00:32:18 2010, comment #11:

Test fails with three loci and r = 0.3. Why?

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Mon Aug 23 21:35:39 2010, comment #10:

Four-locus simulation chokes; my suspicion is that we are ending up with more than two expressions to sum, thereby not subtracting enough. We need to find a way to add the expressions on the fly so we get as many as we need.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Mon Aug 23 18:55:04 2010, comment #9:

Now producing correct tables for two and three loci; we should not close this bug until we have tested larger numbers of loci.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 20:27:45 2010, comment #8:

Revising the mating table calculations for the random mating examples produces allele frequencies that do not change. See mating.c in revision 1086.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 19:59:04 2010, comment #7:

Revision 1084 fixes the problem of not assigning to crosses that can produce certain offspring.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 18:18:06 2010, comment #6:

Rec_gen_table is producing incorrect recombination tables!

For example, the table for offspring 3 in a tlta with r = 0.5:

but instead it is

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 18:17:05 2010, comment #5:

Rec_gen_table is producing incorrect recombination tables!

For example, the table for offspring 3 in a tlta with r = 0.5:

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 18:15:22 2010, comment #4:

Rec_gen_table is producing incorrect recombination tables!

For example, the table for offspring 3 in a tlta with r = 0.5:

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 17:11:11 2010, comment #3:

Removing the denominator from rec_mating is logical (the denominator should be taken care of in the mating table calculation). However, small errors lead to a large deviation from frequencies that sum to 1.0:

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 15:06:33 2010, comment #2:
  • Possible Programming Errors
  1. Uncaught errors
  2. floating point errors
  3. adding over the wrong indices
Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 15:04:46 2010, comment #1:

== Possible Programming Errors ==

  1. Uncaught errors
  1. floating point errors
  1. adding over the wrong indices
Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue Aug 10 00:02:57 2010, original submission:

See the following output of nrm:

According to Emacs calc, the frequency of the focal allele at locus 1 changes from 0.95874309 to 0.99005369.

Anonymous

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by trashbird1240 (Posted a comment)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only logged-in users can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 8 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Thu Dec 2 18:35:59 2010trashbird1240Severity5 - Blocker=>3 - Normal
      StatusReady For Test=>In Progress
      Open/ClosedClosed=>Open
    Thu Sep 2 15:22:55 2010trashbird1240CategoryNone=>Library
      Open/ClosedOpen=>Closed
    Wed Sep 1 19:46:20 2010trashbird1240StatusNone=>Ready For Test
    Tue Aug 10 18:15:22 2010trashbird1240Severity3 - Normal=>5 - Blocker
      Assigned toNone=>trashbird1240

    Back to the top


    Powered by Savane 3.1-cleanup1