estimate_parameters_ks: Max likelihood estimation of Kimura parameters that minimize...

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

Max likelihood estimation of Kimura parameters that minimize the KS statistic

Description

Using maximum likelihood to estimate the parameters of a fitted Kimura distribution to the input sample values that minimizes the KS statistic. It is used to showcase that using the KS statistic to support selection needs extra caution.

Usage

estimate_parameters_ks(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 that minimze the KS statistic of a KS test.

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_ks(mouseData1)

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