kimura_lrt: Likelihood ratio between two heteroplasmy samples

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

Likelihood ratio between two heteroplasmy samples

Description

A function to perform likelihood ratio test exploring difference in bottleneck size between two heteroplasmy samples.

Usage

kimura_lrt(h1, h2, use.h0s = F, h1.h0set = 0, h2.h0set = 0)

Arguments

h1	The first heteroplasmy sample.
h2	The two heteroplasmy sample.
use.h0s	Logical parameter. TBD
h1.h0set	Logical parameter.TBD
h2.h0set	Logical parameter.TBD

Value

The maximum likelihood for the input data according to the Kimura distribution (using bootstrapping)

Examples

 X.1 = rnorm(50,0.5,0.1)
 X.2 = rnorm(50,0.5,0.1)
kimura_lrt(X.1,X.2)

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