dma_gamma_known: Density function of gamma dfe

In dwinter/dfe: Distributions of Fitness Effects from new mutations

Description

Density function of gamma dfe

Usage

dma_gamma_known(obs, shape, rate, Ve, k, p_neutral, log = TRUE)

Arguments

obs,	numeric, observed fitnesses
shape,	numeric, shape of gamma distribution
rate,	numeric, rate of gamma distribution
Ve,	numeric, experimental variance
k,	numeric, vector of knonw mutation counts
p_neutral	numeric,proportion of all mutations that have no effect
log	logical return log-liklihood (defaults to true)

Value

numeric (log-) liklihood of the specfified model and data

Examples

set.seed(123)
mu <- rpois(50, 10)
w <- rma_known_gamma(shape=1, rate=25, Ve=0.01, k=mu,p_neutral=0.7) 
dma_gamma_known(w, shape=1, rate=25, Ve=0.01, k=mu, p_neutral=0.75)
dma_gamma_known(w, shape=1, rate=25, Ve=0.01, k=mu, p_neutral=0.65)

dwinter/dfe documentation built on May 15, 2019, 6:21 p.m.