[R-sig-genetics] R Package to compute observed heterozygosity within a population

Thibaut Jombart thibautjombart at gmail.com
Mon Jan 16 13:03:14 CET 2017


Hello,

Hs(x) in adegenet will do this for the expected heterozygosity. You can use
'summary' combined with seppop to get the average Hobs per population;
example with sim2pop:

> lapply(seppop(sim2pop), function(e) mean(summary(e)$Hobs, na.rm = TRUE))
$P01
[1] 0.755

$P02
[1] 0.77


Cheers
Thibaut


--
Dr Thibaut Jombart
Lecturer, Department of Infectious Disease Epidemiology, Imperial College
London
Head of RECON: repidemicsconsortium.org
sites.google.com/site/thibautjombart/
github.com/thibautjombart
Twitter: @TeebzR <http://twitter.com/TeebzR>
+44(0)20 7594 3658

On 12 January 2017 at 13:19, Bhuller, Ravneet <
ravneet.bhuller13 at imperial.ac.uk> wrote:

> Dear Members,
>
> Is there any suitable package to compute observed heterozygosity within a
> population?
>
> I am interested to compute observed heterozygosity for a population using
> genind object, rather than at a locus.
>
> Many thanks for all your time.
>
> Kind regards,
>
> Rav
>
> _______________________________________________
> R-sig-genetics mailing list
> R-sig-genetics at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-genetics
>

	[[alternative HTML version deleted]]



More information about the R-sig-genetics mailing list