bugNumerical Library for Haploid Genetics - Bugs: bug #34879, ld_sub_geno () should calculate...

 
 

bug #34879: ld_sub_geno () should calculate recursively

Submitted by:  Joel J. Adamson <trashbird1240>
Submitted on:  Tue 22 Nov 2011 02:57:50 PM UTC  
 
Category: LibrarySeverity: 3 - Normal
Item Group: NoneStatus: None
Privacy: PublicAssigned to: Joel J. Adamson <trashbird1240>
Open/Closed: Open

Add a New Comment(Rich Markup)
   

You are not logged in

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

 

Tue 22 Nov 2011 03:18:41 PM UTC, comment #1:

Check to make sure this reduces to Bennett's Equation.

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.
Tue 22 Nov 2011 02:57:50 PM UTC, original submission:

ld_sub_geno () does not calculate the full linkage disequilibrium among a set of loci, it calculates a subset of it. We can use the following recursive algorithm to calculate the right value (where U is the set of loci we're interested in, and |U| is its size):

  1. If |U| = 0, or |U| = 1 then return 0
  2. If |U| = 2 then use the standard linkage disequilibrium
  3. If |U| > 2 then split off a locus and use the following formula

E[X_i \Pi_{j \in U/i} (X_j - p_j) ] - p_i D_U/i

with recursion on the last term. U/i is the set difference of U and i, X_i is an indicator variable giving the status to the locus (on or off) and E represents expectation (genotypic average).

For each such calculation we will need the genotypic frequency (for averaging) and the allele frequencies of the constituent alleles (can be calculated once).

Joel J. Adamson <trashbird1240>
Project AdministratorIn charge of this item.

 

(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 (Submitted the item)
  •  

    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):

     

     

    No Changes Have Been Made to This Item

    Back to the top


    Powered by Savane 3.1-cleanup1