joint_neg_log_lik: Joint negative log-likelihood function for several...

In kostasgian21/heteroplasmy: Calculation of the standard error of the variance for heteroplasmy data

Joint negative log-likelihood function for several heteroplasmy measurements

Description

Joint negative log-likelihood function for several families' heteroplasmy measurements. theta = [ b, h0.1, h0.2, ... ] (use h values if use.h0s=F, otherwise initial heteroplasmies are enforced via h0s).

Usage

joint_neg_log_lik(theta, hlist, use.h0s = F, h0s = -1)

Arguments

theta	Kimura parameters p (or h0 here) and b.
hlist	TBDD list of different sets of heteroplasmy measurements
use.h0s	Logical parameter. TBD
h0s	TBDD

Value

The negative log likelihood for the list of inputs.

Examples

 X.1 = rnorm(50,0.5,0.1)
joint_neg_log_lik(c(0.5,0.91),X.1)

kostasgian21/heteroplasmy documentation built on Jan. 30, 2024, 12:30 a.m.