maxlik: Compute maximum likelihood parameters and confidence...

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

Compute maximum likelihood parameters and confidence intervals for heteroplasmy data

Description

Compute maximum likelihood parameters and confidence intervals for a given heteroplasmy set. We can do this while imposing a specific h0 as an argument or allowing a search over h0 values.

Usage

maxlik(h, conf.level = 0.95, h0 = F)

Arguments

h	A vector containig heteroplasmy measurements. Every observation should be in [0,1].
conf.level	The preferred confidence interval calculation, Default value is 0,95 (95%).
h0	Logical parameter. A particular h0 value Default is to treat h0 as a fit parameter

Value

The maximum likelihood for the input data according to the Kimura distribution

Examples

 X.1 = rnorm(50,0.5,0.1)
maxlik(X.1,conf.level=0.95)

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