dist_cavalli: Estimation of Bray-Curtis distance

In dyerlab/gstudio: Tools Related to the Spatial Analysis of Genetic Marker Data

Estimation of Bray-Curtis distance

Description

This function returns a measure of genetic distance based upon the Cavalli-Sforza & Edwards (1967) distance metric. Assuming drift is the only source or differences observed among strata.

Usage

dist_cavalli(x, stratum = "Population")

Arguments

x	The genetic data, either as a single locus or multilocus (data.frame) object.
stratum	The groups among which you are going to estimate genetic distances.

Value

A matrix of Cavalli-Sforza Genetic distance estimates.

Author(s)

Rodney J. Dyer [e-mail]

Examples

  AA <- locus( c("A","A") )
  AB <- locus( c("A","B") )
  BB <- locus( c("B","B") )
  AC <- locus( c("A","C") )
  AD <- locus( c("A","D") )
  BC <- locus( c("B","C") )
  BD <- locus( c("B","D") )
  CC <- locus( c("C","C") )
  CD <- locus( c("C","D") )
  DD <- locus( c("D","D") )
  loc1 <- c(AA,AB,AC,AD,BB,BC,BD,CC,CD,DD)
  loc2 <- c(AA,AA,AC,AA,CC,CC,AC,CC,AA,AC)
  df <- data.frame( Population=c(rep("Pop-A",5),rep("Pop-B",5)), TPI=loc1, PGM=loc2 )

dyerlab/gstudio documentation built on Feb. 2, 2024, 8:24 p.m.