Nothing
lxy <- function(row, data, pop1, pop2, allele.column, ref.pop=NA)
{
ref.pop <- .subset2(ref.pop,1)
n.ref.pop <- names(ref.pop)
re <- .subset2(data,1)[row]
rat <- .subset2(data,2)[row]
a <- allele.column*2+1
ai <- .subset2(pop1,a)[re]
aj <- .subset2(pop1,a+1)[re]
bi <- .subset2(pop2,a)[rat]
bj <- .subset2(pop2,a+1)[rat]
# Calculations after Lynch and Ritland 1999 (5a)
# Share.RE for colnames/pop2
pxma <- .subset2(ref.pop,which(n.ref.pop==ai))
pymb <- .subset2(ref.pop,which(n.ref.pop==aj))
het <- ai==aj
share.RE <-
(
((ai==bi)+(ai==bj))*pymb
+
((aj==bi)+(aj==bj))*pxma-4*pxma*pymb
)/((1+het)*(pxma+pymb)-4*pxma*pymb)
# weighted for Lynch 1999 equ 7a
share.RE <- share.RE*((1+het)*(pxma+pymb)-4*pxma*pymb)/(2*pxma*pymb)
# Share.RAT for pop2
# het == 1 if Alleles identical otherwise het == 0
pxma <- .subset2(ref.pop,which(n.ref.pop==bi))
pymb <- .subset2(ref.pop,which(n.ref.pop==bj))
het <- bi==bj
share.RAT <-
(
((bi==ai)+(bi==aj))*pymb
+
((bj==ai)+(bj==aj))*pxma-4*pxma*pymb
)/((1+het)*(pxma+pymb)-4*pxma*pymb)
# weighted for Lynch 1999 equ 7a
share.RAT <- share.RAT*((1+het)*(pxma+pymb)-4*pxma*pymb)/(2*pxma*pymb)
return(c(share.RE,share.RAT))
if (length(c(share.RE,share.RAT))>0){return(c(share.RE,share.RAT))}else{return(NA)}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.