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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.