hstats: Calculate various statistics for a heteroplasmy set h

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

Calculate various statistics for a heteroplasmy set h

Description

Can enforce an initial h0 or leave as a free parameter. Can can use population or sample statistics. analyticVar offers a simplified version of this function to compute the standard error of the variance.

Usage

hstats(h, h0 = F, usepopn = F)

Arguments

h	A vector containig heteroplasmy measurements. Every observation should be in [0,1].
h0	Logical parameter. A particular h0 value Default is to treat h0 as a fit parameter
usepopn	Logical parameter. Use of population or sample statistics (T and F, respectively)

Value

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

Examples

 X.1 = rnorm(50,0.5,0.1)
hstats(X.1)

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