R/interceptToM.R

In KamilSJaron/ConjunctionStats: The parser of simulations of hybrid zones performed by simulator Conjunction (http://github.com/KamilSJaron/Conjunction)

#' @title interceptToM
#'
#' @description
#' \code{interceptToM} converts intercept of line deifning loglog distribution
#' of blocks to number of migrants needed for distribution in 0D
#'
#' @param intercept
#'
#' @param selection
#'
#' @param theta
#'
#' @author Kamil Jaron \email{kamiljaron at gmail.com}
#'
#' @export

interceptToM <- function(intercept, selection, theta){
  # intercept = log((M * theta * 2) / (S * (1 + theta)^2))
  # exp(intercept) = (M * theta * 2) / (S * (1 + theta)^2)
  # (exp(intercept) * (S * (1 + theta)^2)) / (theta * 2) = M
  return( (exp(intercept) * (selection * (1 + theta)^2)) / (theta * 2) )
}