R/forwardBackward.r

In genostats/FEstim: Genetic Data Analyses in Presence of Inbreeding

##################################################################################
#This function allows the computation of HBD probabilities                       #
#                                                                                #
#!!! logEmission : a dataframe with log.emission computed                        #                                       
#!!! delta.Distance : a vector with delta.dist computed                          #
#!!! a : the value of a for the submaps                                          #
#!!! f : the value of f for the submaps                                          #
#                                                                                #
#*** return a dataframe with HBD probabilities computed                          #
##################################################################################

forwardBackward <- function(logEmission, delta.Distance, a, f) {
  if( a < 0 ) stop("a should be >= 0")
  if( f < 0 | f > 1) stop("f should be between 0 and 1")
  if( f == 0) 
    return( matrix( c(1,0), nrow = 2, ncol = ncol(logEmission) ) )
  if( f == 1) 
    return( matrix( c(0,1), nrow = 2, ncol = ncol(logEmission) ) )
  .Call('_Fantasio_forward_backward', PACKAGE = "Fantasio", logEmission, delta.Distance, a, f)
}