inst/doc/test/distance-functions.R

In stranda/kernelPop2: Spatially explicit population genetic simulations

pollination.dist <- function(l)
  {
    sqrt((l$ind[,6]-l$ind[,8])^2+(l$ind[,7]-l$ind[,9])^2)
  }
seed.dist <- function(l)
  {
    sv <- l$ind[,6]|l$ind[,7]|l$ind[,8]|l$ind[,9]
    sqrt((l$ind[sv,6]-l$ind[sv,4])^2+(l$ind[sv,7]-l$ind[sv,5])^2)
  }
realized.selfing <- function(l)
  {
    pd <- pollination.dist(l)
    sum(pd==0)/length(pd)
  }