11.07.2015 Views

The genetics Package

The genetics Package

The genetics Package

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

12 diseqr is the correlation coefficient between two alleles, and can be computed byr =−D√(pi ∗ (1 − p i )p(j)(1 − p j ))where- p i defined as the observed probability of allele ’i’,- p j defined as the observed probability of allele ’j’, and- p ij defined as the observed probability of the allele pair ’ij’.When there are more than two alleles, the summary values for these statistics are obtained by computinga weighted average of the absolute value of each allele pair, where the weight is determinedby the expected frequency. For example:D overall = ∑ i≠j|D ij | ∗ p ijBootstrapping is used to generate confidence interval in order to avoid reliance on parametric assumptions,which will not hold for alleles with low frequencies (e.g. D ′ following a a Chi-squaredistribution).See the function HWE.test for testing Hardy-Weinberg Equilibrium, D = 0.Valuediseq returns an object of class diseq with componentscalldataD.hatTODOfunction call used to create this object2-way table of allele pair countsmatrix giving the observed count, expected count, observed - expected difference,and estimate of disequilibrium for each pair of alleles as well as an overalldisequilibrium value.more slots to be documentednormal-bracket100bracket-normaldiseq.ci returns an object of class bootciAuthor(s)Gregory R. Warnes 〈warnes@bst.rochester.edu 〉See Alsogenotype, HWE.test, boot, bootci

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!