[adegenet-forum] proportion of shared alleles

jombart at biomserv.univ-lyon1.fr jombart at biomserv.univ-lyon1.fr
Fri Apr 25 14:05:50 CEST 2008


Quoting Grant Gillis <grant.j.gillis at gmail.com>:

> List and Thibaut,
>
> I posted this question on the R-sig-phylo list and Thibaut suggested that
> this is a more appropriate list.  Great looking package by the way.   I am
> looking for a way to calculate the proportions of shared alleles between
> individuals.  I have been looking through the documentation but am having a
> hard time getting started (a result of my lack of familiarity with R).  Any
> ideas or nudges in the right direction are appreciated.
>
> Grant
>

Dear Grant,

here is a function that should do what you want (called propShared), and that
will eventually become a part of adegenet.

I am out for vacations until tuesday, and all I have here is an old mac with
R 2.2 installed... So I could not test the function
as much as I would have liked to, but this is a start. The input is a genind
object, i.e. the representation of a set of genotypes in adegenet. It computes
the proportions of shared alleles for each locus, (function f1) and then
averages theses proportions over loci to get the overall measure. The handling
of missing data as it is here inderestimates the proportion of shared alleles.
I'll fix this in the final implementation of the method in the package, which
will likely be in C anyway.

Best,

Thibaut.



--------------------------------------------------------------------------
Ce message a été envoyé depuis le webmail IMP (Internet Messaging Program)
-------------- next part --------------
A non-text attachment was scrubbed...
Name: propShared.R
Type: application/octet-stream
Size: 1226 bytes
Desc: not available
Url : http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20080425/13082470/attachment.obj 


More information about the adegenet-forum mailing list