R/morans.fin.R

In Demerelate: Functions to Calculate Relatedness on Diploid Genetic Data

morans.fin <- function(row, data, pop1, pop2, allele.column, ref.pop=NA)
{
  # Morans I from Spagedi manual
  re <- .subset2(data,1)[row]
  rat <- .subset2(data,2)[row]
  a <- allele.column*2+1

  n.ref.pop <- names(ref.pop)
  
  ai <- .subset2(pop1,a)[re]
  aj <- .subset2(pop1,a+1)[re]
  bi <- .subset2(pop2,a)[rat]
  bj <- .subset2(pop2,a+1)[rat]
  
  ai<-n.ref.pop==ai
  aj<-n.ref.pop==aj
  bi<-n.ref.pop==bi
  bj<-n.ref.pop==bj
  
  r.return<-sum((ai/2+aj/2-ref.pop)*(bi/2+bj/2-ref.pop))

  
  return(r.return)
  
}