estimate_parameters_ml: Max likelihood estimation for Kimura distribution parameters

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

Max likelihood estimation for Kimura distribution parameters

Description

Using maximum likelihood to estimate the parameters of a fitted Kimura distribution to the input sample values.

Usage

estimate_parameters_ml(h)

Arguments

h	A vector containig heteroplasmy measurements. Every observation should be in [0,1].

Value

The maximum likelihood estimates for a fitted Kimura distribution parameters.

Author(s)

Kostas and Iain, [e-mail]

References

Site or paper

See Also

readHeteroplasmyData

Examples

# size of the sample
n=50
#generate a random sample of size n from a normal distribution
data_ex=rnorm(n,0.5,0.1)
estimate_parameters_ml(data_ex)

mouseData=readHeteroplasmyData("LE")
mouseData1 = mouseData[which(!is.na(mouseData[,1])),1]
estimate_parameters_ml(mouseData1)

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