get_llh: Get the log-likelihood of data while incorporating sampling...

In benmarwick/signatselect: signatselect: Identifying signatures of selection

Description

get_llh uses Laplace approximation to compute the log-likelihood of the gaussian approximation assuming the data is drawn from a series of binomial choices.

Usage

get_llh(prior, lgth, tvec, nvec, bvec, prec, maxiter, verbose, s, alpha, f0)

Arguments

prior	a function that returns the prior distribution (default: function(x)1)
lgth	number of data points in the time series
tvec	time points for the time series (start with 0)
nvec	total number of samples for each time point
bvec	number of successful samples at each time point
prec	minimum precision for optimisation
maxiter	maximum number of iterations for optimisation
verbose	print intermediate output (TRUE or FALSE)
s	selection coefficient
alpha	population size, bounded [0,Inf)
f0	initial frequency for the underlying logistic change

Value

A log-likelihood estimate.

benmarwick/signatselect documentation built on Dec. 4, 2019, 4:22 p.m.